Создайте вектор 100 равномерно распределенных чисел от $$0$$ к $$3\pi$$. Сгенерируйте прямоугольную волну с периодом $$2\pi$$.

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

Постройте прямоугольную волну и наложите синус. Нормируйте ось X $$\pi$$. Сгенерированная прямоугольная волна имеет значение $$1$$ для интервалов [$\mathit{n}\pi ,\left(\mathit{n}+1\right)\pi$) с даже $\mathit{n}$ и значение $$-1$$ для интервалов [$\mathit{n}\pi ,\left(\mathit{n}+1\right)\pi$) с нечетным $\mathit{n}$. Волна никогда не имеет значение $$0$$.

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

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

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