exa2_60.m

来自「有关matlab r2007的实用教程」· M 代码 · 共 14 行

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
figure
plot(f,p*180/pi),ylabel('phase');
grid on