例2-60.m

来自「上传最新Matlab7基础与提高,与各位爱好者分享」· M 代码 · 共 13 行

>> x=[1 ;3; 5; 7; 9];
>> y=fft(x)
>> t=0:0.02:1;                                %波形方程
>> x=sin(2*pi*20*t)+sin(2*pi*12*t);
>> y=fft(x);                                  %经过傅立叶变换观察其模数和相位
>> m=abs(y);
>> p=unwrap(angle(y));
>> f=(0:length(y)-1)*50/length(y);                %以频率为x轴画图
>> plot(f,m),ylabel('Abs');
>> grid on
>> plot(f,p*180/pi),ylabel('phase');
>> grid on