ex1510.m

数字信号处理的MATLAB实现

%例15-10 傅立叶分析
t=1:20;
y=sin(0.6*t)+cos(3*t);
plot(t,y)	%绘制y关于t的函数图像，如图15-13
fy=fft(y)	%对y进行快速傅立叶变换
mfy=abs(fy);	%求y经傅立叶变换后的幅值
plot(t,mfy)	%如图15-14
afy=angle(fy);	%求y经傅立叶变换后的相位角
uafy=unwrap(afy);	%调整间距超过pi的相位角
plot(t,afy,'r*',t,afy,'r')	%用红色*标记数据点并连线画图
hold on	%打开图形的叠加绘制模式
plot(t,uafy,'go',t,uafy,'g')	%用绿色o标记数据点并连线画图，如图15-15
hold off	%关闭图形的叠加绘制模式
ifft(fy)-y