Создайте вектор из 100 равноотстоящих чисел от $$0$$ до $$\mathrm{}\mathrm{3\delta }$$. Генерируйте квадратную волну с периодом $$\mathrm{}\mathrm{2\delta }$$.

t = linspace(0,3*pi)';
x = square(t);

Постройте график квадратной волны и наложите синус. Нормализуйте ось X на $$\delta$$. Генерируемая квадратная волна имеет значение $$1$$ для интервалов [$\mathit{}\mathrm{n\lambda },\mathit{(}n\mathrm{}+$1) λ) с $\mathit{}$чётным n и значение $$\mathrm{}$$-1 для $\mathit{}\mathrm{интервалов}[\mathit{}n\lambda \mathrm{,}($n + 1) λ) $\mathit{}$с нечётным n. Волна никогда не имеет $$\mathrm{}$$значения 0.

plot(t/pi,x,'.-',t/pi,sin(t))
xlabel('t / \pi')
grid on

Повторите расчет, но теперь оцените square(2*t) при 121 равноотстоящих числах между $$-\lambda$$ и $$\mathrm{}\mathrm{2\lambda }$$. Измените амплитуду на $$\mathrm{}\mathrm{}\mathrm{1,15}$$. Постройте график волны и наложите синус с теми же параметрами. Эта новая волна отрицательна при $$t\mathrm{=}$$ 0 и положительна в конечных точках$$,$$-δ $$\mathrm{и}$$2δ.

t = linspace(-pi,2*pi,121);
x = 1.15*square(2*t);

plot(t/pi,x,'.-',t/pi,1.15*sin(2*t))
xlabel('t / \pi')
grid on

Создайте квадратную волну 30 Гц, отобранную на частоте 1 кГц для 70 мс. Укажите рабочий цикл 37%. Добавьте белый гауссов шум с дисперсией 1/100.

t = 0:1/1e3:0.07;
y = square(2*pi*30*t,37)+randn(size(t))/10;

Вычислите рабочий цикл волны. Постройте график формы сигнала и аннотируйте рабочий цикл.

dutycycle(y,t)

ans = 0.3639