square

创建一个由 $$0$$ 到 $$3\pi$$ 之间的 100 个等间距数字组成的向量。生成周期为 $$2\pi$$ 的方波。

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

绘制方波并叠加一个正弦波。用 $$\pi$$ 归一化 x 轴。对于具有偶数 $\mathit{n}$ 的区间 [$\mathit{n}\pi ,\left(\mathit{n}+1\right)\pi$)，生成的方波的值为 $$1$$，而对于具有奇数 $\mathit{n}$ 的区间 [$\mathit{n}\pi ,\left(\mathit{n}+1\right)\pi$)，生成的方波的值为 $$-1$$。该波的值永远不会是 $$0$$。

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

重复计算，但现在，基于 $$-\pi$$ 和 $$2\pi$$ 之间的 121 个等间距数字计算 square(2*t)。将振幅更改为 $$1.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 Hz 方波，其采样率为 1 kHz，持续时间为 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