实验4：连续系的频域分析.docVIP

实验4：连续系统的频域分析 一、实验目的 （1）掌握连续时间信号的傅里叶变换和傅里叶逆变换的实现方法。 （2）掌握傅里叶变换的数值计算方法和绘制信号频谱的方法。 二、实验原理 1.周期信号的分解 根据傅里叶级数的原理，任何周期信号都可以分解为三角级数的组合——称为的傅里叶级数。在误差确定的前提下，可以由一组三角函数的有限项叠加而得到。 例如一个方波信号可以分解为： 合成波形所包含的谐波分量越多，除间断点附近外，它越接近于原波形，在间断点附近，即使合成的波形所含谐波次数足够多，也任存在约的偏差，这就是吉布斯现象（Gibbs）。 2.连续时间信号傅里叶变换的数值计算 由傅里叶变换的公式： 当为时限信号时，上式中的取值可以认为是有限项，则有： ，其中 3.系统的频率特性 连续LTI系统的频率特性称为频率响应特性，是指在正弦信号激励作用下稳态响应随激励信号频率的变化而变化的情况，表示为 三、实验内容与方法 1.周期信号的分解 【例1】用正弦信号的叠加近似合成一个频率为的方波。 MATLAB程序如下： clear all; fs=10000; t=[0:1/fs:0.1]; f0=50;sum=0; subplot(211) for n=1:2:9 plot(t,4/pi*1/n*sin(2*pi*n*f0*t),’k’); hold on; end title(‘信号叠加前’); subplot(212) for n=1:2:9; sum=sum+4/pi*1/n*sin(2*pi*n*f0*t); end plot(t,sum,’k’); title(‘信号叠加后’); 产生的波形如图所示： 2.傅里叶变换和逆变换的实现 求傅里叶变换，可以调用fourier函数，调用格式为F=fourier(f,u,v)，是关于u的函数f的傅里叶变换，返回函数F是关于v的函数。 求傅里叶逆变换，可以调用ifourier函数，调用格式为f=ifourier(F,v,u)，是关于v的函数F的傅里叶逆变换，返回函数f是关于u的函数。 【例2】已知连续信号，通过程序完成其傅里叶变换。 MATLAB程序如下： syms t; f=fourier(exp(-2*abs(t))); ezplot(f)?; 得到的傅里叶变换如图所示： 【例3】已知连续信号，通过程序完成其傅里叶逆变换。 MATLAB程序如下： syms t w? ifourier(1/(1+w^2),t)? 得到的结果为：ans =1/2*exp(-t)*heaviside(t)+1/2*exp(t)*heaviside(-t) 图形如图所示： 3.傅里叶变换的性质 举例验证傅里叶变换的时移特性和频移特性。 【例4】分别绘出信号和的频谱，求的频谱。 MATLAB程序如下： r=0.02;t=-5:r:5;N=200;W=2*pi;k=-N:N;w=k*W/N; f1=1/2*exp(-2*t).*stepfun(t,0); F=r*f1*exp(-j*t*w); F1=abs(F);p1=angle(F);subplot(3,1,1);plot(t,f1);grid xlabel(t);ylabel(f(t));title(f(t));subplot(3,1,2); plot(w,F1);xlabel(w);grid;ylabel(F(jw));subplot(3,1,3); plot(w,p1*180/pi);grid;xlabel(w);ylabel(相位(度)); 再求信号的频谱，MATLAB程序如下： %求的频谱 r=0.02;t=-5:r:5;N=200;W=2*pi;k=-N:N;w=k*W/N; f1=1/2*exp(-2*(t-1)).*stepfun(t,1); F=r*f1*exp(-j*t*w); F1=abs(F);p1=angle(F);subplot(3,1,1);plot(t,f1);grid xlabel(t);ylabel(f(t));title(f(t-1));subplot(3,1,2); plot(w,F1);xlabel(w);grid;ylabel(F(jw)的模);subplot(3,1,3); plot(w,p1*180/pi);grid;xlabel(w);ylabel(相位(度)); 【例5】傅里叶变换的频移特性：信号为门信号，绘出信号和信号的频谱，并与原信号的频谱图进行比较。 （1），求其频谱可以采用数值就算得方法。MATLAB程序如下： R=0.02;t=-2:R:2; f=stepfun(t,-1)-stepfun(t,1); W1=2*pi*5;%频率宽度 N=500;k=0:N;W=k*W1/N;%采样数为N，W为频率正半轴的采样点 F