傅里叶的分析及应用.docVIP

实验二 傅里叶分析及应用 一、实验目的 （一）掌握使用Matlab进行周期信号傅里叶级数展开和频谱分析 1、学会使用Matlab分析傅里叶级数展开，深入理解傅里叶级数的物理含义 2、学会使用Matlab分析周期信号的频谱特性 （二）掌握使用Matlab求解信号的傅里叶变换并分析傅里叶变换的性质 1、学会运用Matlab求连续时间信号的傅里叶变换 2、学会运用Matlab求连续时间信号的频谱图 3、学会运用Matlab分析连续时间信号的傅里叶变换的性质 掌握使用Matlab完成信号抽样并验证抽样定理 1、学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析 2、学会运用MATLAB改变抽样时间间隔，观察抽样后信号的频谱变化 3、学会运用MATLAB对抽样后的信号进行重建 二、实验条件 安装winXP系统的电脑一台、matlab 7.0软件 三、实验内容 1、已知周期三角信号如下图所示[注：图中时间单位为：毫秒(ms)]： （1）试求出该信号的傅里叶级数[自己求或参见课本P112或P394]，利用Matlab编程实现其各次谐波[如1、3、5、13、49]的叠加，并验证其收敛性； 解： 命令文件： clear all;close all;clc; t=-10:0.01:10; omega=pi; y=abs(sawtooth(pi*0.5*t,0.5)); plot(t,y),grid on; axis([-10,10,0,3]); n_max=[1,3,5,13,49]; N=length(n_max); for k=1:N n=1:2:n_max(k); b=4./((pi*n).^2); x=b*cos(omega*n*t); figure; plot(t,y); hold on; x=x+1/2; plot(t,x); hold off; axis([-10,10,0,3]); title([最大谐波数=,num2str(n_max(k))]); end 图像： （2）用Matlab分析该周期三角信号的频谱[三角形式或指数形式均可]。当周期三角信号的周期（如由2ms(1ms或由2ms(4ms）和宽度（如2ms(1ms）分别变化时，试观察分析其频谱的变化。 解： 周期为2ms： 命令文件： clear all;close all;clc; dt=0.01; t=-4:dt:4; ft=(t=-1t0).*(t+1)+(t0t=1).*(1-t); %subplot(2,1,1) %plot(t,ft);grid on n=2000; k=-n:n; w=pi*k/(n*dt); f=dt*ft*exp(-i*t*w); f=abs(f); %subplot(2,1,2) plot(w,f); axis([-20 20 0 1.1]),grid on; 图像： 当周期为1ms时： 命令文件： clear all;close all;clc; dt=0.01; t=-4:dt:4; ft=(t=-0.5t0).*(t+1)+(t0t=0.5).*(1-t); %subplot(2,1,1) %plot(t,ft);grid on n=2000; k=-n:n; w=pi*k/(n*dt); f=dt*ft*exp(-i*t*w); f=abs(f); %subplot(2,1,2) plot(w,f); axis([-20 20 0 1.1]),grid on; 图像： 宽度与周期均变为1ms，所以只有两个图 2、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT，并画出其频谱图（包括幅度谱和相位谱）[注：图中时间单位为：毫秒(ms)]。 解： 符号运算求解法： 命令文件： %符号运算求解法 clear all;close all;clc; ft=sym((t+2)*(Heaviside(t+2)-Heaviside(t+1))+Heaviside(t+1)-Heaviside(t-1)+(2-t)*(Heaviside(t-1)-Heaviside(t-2))); Fw=fourier(ft); subplot(2,1,1) axis([-8,8,0,3]); ezplot(abs(Fw)),grid on title(幅度谱) phase=atan(imag(Fw)/real(Fw)); subplot(2,1,2) ezplot(phase),grid on title(相位谱) 图像： 数值计算法： 命令文件： clear all;close all;clc; dt=0.01;