实验：连续系统的频域分析.pdfVIP

实验4 ：连续系统的频域分析 一、实验目的 （1 ）掌握连续时间信号的傅里叶变换和傅里叶逆变换的实现方法。 （2 ）掌握傅里叶变换的数值计算方法和绘制信号频谱的方法。 二、实验原理 f (t) 根据傅里叶级数的原理，任何周期信号都可以分解为三角级数的组合——称为 的傅里叶级数。在误差确定的前提下，可以由一组三角函数的有限项叠加而得到。 例如一个方波信号可以分解为： 9% 合成波形所包含的谐波分量越多，除间断点附近外，它越接近于原波形，在间断点附近，即使合成的波形所含谐波次数足够多，也任存在约 的偏差，这就是吉布 斯现象（Gibbs ）。 由傅里叶变换的公式： f (t) n N 当 为时限信号时，上式中的 取值可以认为是有限项 ，则有： N 2      j n   k F (k )  f (n )e k ,0 k N ， 其 中 k N n0 Y () H ()  连 续 LTI 系统的频率特性称为频率响应特性 ，是指在正弦信号激励作用下稳态响应随激励信号频率的变化而变化的情况 ，表示为 X () 三、实验内容与方法 50Hz 【例 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( ‘信号叠加后’); 产生的波形如图所示 ： 求傅里叶变换，可以调用 fourier 函数，调用格式为 F=fourier(f,u,v) ，是关于 u 的函数 f 的傅里叶变换，返回函数 F 是关于 v 的函数。 求傅里叶逆变换，可以调用 ifourier 函数，调用格式为 f=ifourier(F,v,u) ，是关于 v 的函数 F 的傅里叶逆变换，返回函数 f 是关于 u 的函数。 【例 2】已知连续信号 f (t) e2 t ，通过程序完成其傅里叶变换。 MATLAB 程序如下： syms