最新实验4：连续系统的频域分析
- 格式：doc
- 大小：147.00 KB
- 文档页数：8
实验4：连续系统的频域分析 精品好文档，推荐学习交流
仅供学习与交流，如有侵权请联系网站删除 谢谢7 实验4：连续系统的频域分析
一、实验目的
（1）掌握连续时间信号的傅里叶变换和傅里叶逆变换的实现方法。
（2）掌握傅里叶变换的数值计算方法和绘制信号频谱的方法。
二、实验原理
1.周期信号的分解
根据傅里叶级数的原理，任何周期信号都可以分解为三角级数的组合——称为()ft的傅里叶级数。在误差确定的前提下，可以由一组三角函数的有限项叠加而得到。
例如一个方波信号可以分解为：
11114111()sinsin3sin5sin7357Efttttt
合成波形所包含的谐波分量越多，除间断点附近外，它越接近于原波形，在间断点附近，即使合成的波形所含谐波次数足够多，也任存在约9%的偏差，这就是吉布斯现象（Gibbs）。
2.连续时间信号傅里叶变换的数值计算
由傅里叶变换的公式：
0()()lim()jtjnnFjftedtfne
当()ft为时限信号时，上式中的n取值可以认为是有限项N，则有：
0()(),0kNjnnFkfnekN，其中2kkN
3.系统的频率特性
连续LTI系统的频率特性称为频率响应特性，是指在正弦信号激励作用下稳态响应随激励信号频率的变化而变化的情况，表示为()()()YHX
三、实验内容与方法
1.周期信号的分解
【例1】用正弦信号的叠加近似合成一个频率为50Hz的方波。
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; 精品好文档，推荐学习交流
仅供学习与交流，如有侵权请联系网站删除 谢谢7 sum=sum+4/pi*1/n*sin(2*pi*n*f0*t);
end
plot(t,sum,’k’);
title(‘信号叠加后’);
产生的波形如图所示：
00.010.020.030.040.050.060.070.080.090.1-2-1012信号叠加前00.010.020.030.040.050.060.070.080.090.1-2-1012信号叠加后
2.傅里叶变换和逆变换的实现
求傅里叶变换，可以调用fourier函数，调用格式为F=fourier(f,u,v)，是关于u的函数f的傅里叶变换，返回函数F是关于v的函数。
求傅里叶逆变换，可以调用ifourier函数，调用格式为f=ifourier(F,v,u)，是关于v的函数F的傅里叶逆变换，返回函数f是关于u的函数。
【例2】已知连续信号2()tfte，通过程序完成其傅里叶变换。
MATLAB程序如下：
syms t;
f=fourier(exp(-2*abs(t)));
ezplot(f) ;
得到的傅里叶变换如图所示： 精品好文档，推荐学习交流
仅供学习与交流，如有侵权请联系网站删除 谢谢7 -6-4-2024600.10.20.30.40.50.60.70.80.91w4/(4+w2)
【例3】已知连续信号21()1Fj，通过程序完成其傅里叶逆变换。
MATLAB程序如下：
syms t w
ifourier(1/(1+w^2),t)
得到的结果为：ans =1/2*exp(-t)*heaviside(t)+1/2*exp(t)*heaviside(-t)
图形如图所示：
-3-2-1012300.050.10.150.20.250.30.350.40.450.5t1/2 exp(-t) heaviside(t)+1/2 exp(t) heaviside(-t)
3.傅里叶变换的性质
举例验证傅里叶变换的时移特性和频移特性。
【例4】分别绘出信号21()()2tftet和(1)ft的频谱，求21()()2tftet的频谱。
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); 精品好文档，推荐学习交流
仅供学习与交流，如有侵权请联系网站删除 谢谢7 plot(w,p1*180/pi);grid;xlabel('w');ylabel('相位(度)');
-5-4-3-2-101234500.5tf(t)f(t)-8-6-4-20246800.20.4wF(jw)-8-6-4-202468-1000100w相位(度)
再求信号(1)ft的频谱，MATLAB程序如下：
%求(1)ft的频谱
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('相位(度)'); 精品好文档，推荐学习交流
仅供学习与交流，如有侵权请联系网站删除 谢谢7 -5-4-3-2-101234500.5tf(t)f(t-1)-8-6-4-20246800.20.4wF(jw)的模-8-6-4-202468-2000200w相位(度)
【例5】傅里叶变换的频移特性：信号()()ftgt为门信号，绘出信号101()()jtftfte和信号102()()jtftfte的频谱，并与原信号的频谱图进行比较。
（1）()()(1)(1)ftgttt，求其频谱可以采用数值就算得方法。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=f*exp(-j*t'*W)*R;%求F（jw）
F=real(F);
W=[-fliplr(W),W(2:501)];%形成负半轴及正半轴的2N+1个频率点W
F=[fliplr(F),F(2:501)];%形成对应于W的F（jw）的值
subplot(2,1,1);plot(t,f);
xlabel('t');ylabel('f(t)');axis([-2,2,-0.5,2]);
title('f(t)=u(t+1)-u(t-1)');subplot(2,1,2);plot(W,F);
xlabel('W');ylabel('F(W)'); title('f(t)的傅里叶变换'); 精品好文档，推荐学习交流
仅供学习与交流，如有侵权请联系网站删除 谢谢7 -2-1.5-1-0.500.511.52-0.500.511.52tf(t)f(t)=u(t+1)-u(t-1)-40-30-20-10010203040-0.500.511.52WF(W)f(t)的傅里叶变换
（2）得到101()()jtftfte，102()()jtftfte的频谱的MATLAB程序如下：
R=0.02;t=-2:R:2;
f=stepfun(t,-1)-stepfun(t,1);
f1=f.*exp(-j*10*t) ;f2=f.exp*(j*10*t)
W1=2*pi*5;
N=500;
k=-N :N ;
W=k*W1/N;
F1=f1*exp(-j*t'*W)*R;
F2=f2*exp(-j*t'*W)*R;
F1=real(F1);F2=real(F2);
subplot(2,1,1);plot(W,F1);
xlabel('W');ylabel('F1(W)');title('频谱F1(jw)');
subplot(2,1,2);plot(W,F2);
xlabel('W');ylabel('F2(W)');title('频谱F2(jw)');
得到的傅里叶变换的频移特性如图所示： 精品好文档，推荐学习交流
仅供学习与交流，如有侵权请联系网站删除 谢谢7 -40-30-20-10010203040-0.500.511.52WF1(W)频谱F1(jW)-40-30-20-10010203040-0.500.511.52WF2(W)频谱F2(jW)
四、程序设计实验
（1）方波的合成实验。用5项谐波合成一个频率为50Hz，幅值为3的方波，写出MATLAB程序，给出实验的结果。
（2）编写程序，画出信号3()()tftet，(4)ft以及信号4()jtfte的频谱图。
五、实验预习要求
（1）预习实验原理。
（2）熟悉实验程序。
（3）思考课程设计实验部分程序的编写。
六、实验报告要求
（1）在MATLAB中输入程序，验证实验结果，并将实验结果存入指定存储区域中。
（2）对于程序设计实验，要求通过对验证性实验的练习，自行编制完整的实验程序，实现对信号的模拟，并得出实验结果。
（3）在实验报告中写出完整的自编程序，并给出实验结果。
七、思考题
（1）傅里叶级数是什么？非周期傅里叶变换的定义是什么？
（2）将信号进行分解成谐波函数，n次谐波时能否得到原波形？如不能会存在多少误差？
（3）常数和阶跃函数是否能够直接利用傅里叶变换定义公式进行变换？为什么不能？
