m0202.m

来自「MATLAB仿真技术与实例应用教程」· M 代码 · 共 10 行

t=0:0.1:1;
s=square(2*pi*25*t);  %周期为0.04
x=sin(t);
y=cos(t);
f=abs(fft(s,128));    %对方波作FFT变换后取其幅值
a=0:1:127;
x=sin(a/128);
y=cos(a/128);
stem3(x,y,f,'b','fill');  %用菱形杆头标记，用蓝色填色