信号与系统实验二.doc

实验二、连续时间系统的频域分析 一、实验目的 (1) 理解周期信号的傅里叶分解，掌握傅里叶系数的计算方法； 深刻理解和掌握非周期信号的傅里叶变换及其计算方法； (3) 熟悉傅里叶变换的性质，并能应用其性质实现信号的幅度调制； (4) 理解连续时间系统的频域分析原理和方法，掌握连续系统的频率响应求解方法，并画出相应的幅频、相频响应曲线。 二、实验原理 (1) 周期信号的傅里叶分解 设有连续时间周期信号，它的周期为T，角频率，且满足狄里赫利条件，则该周期信号可以展开成傅里叶级数，即可表示为一系列不同频率的正弦或复指数信号之和。傅里叶级数有三角形式和指数形式两种。 三角形式的傅里叶级数： 式中系数，称为傅里叶系数，可由下式求得： 指数形式的傅里叶级数： 式中系数称为傅里叶复系数，可由下式求得： 周期信号的傅里叶分解用Matlab进行计算时，本质上是对信号进行数值积分运算。Matlab中进行数值积分运算的函数有quad函数和int函数。其中int函数主要用于符号运算，而quad函数（包括quad8，quadl）可以直接对信号进行积分运算。因此利用Matlab进行周期信号的傅里叶分解可以直接对信号进行运算，也可以采用符号运算方法。quadl函数(quad系)的调用形式为：y＝quadl(‘func’,a,b)或y＝quadl(@myfun,a,b)。其中func是一个字符串，表示被积函数的.m文件名（函数名）；a、b分别表示定积分的下限和上限。第二种调用方式中”@”符号表示取函数的句柄，myfun表示所有限定义的函数的文件名。 示例1：给定一个周期为4、脉冲宽度为2的矩形信号，用Matlab计算其傅里叶系数并绘图。 程序代码： T=4;tao=2;w=2*pi/T; singrect=inline((abs(t)=1).*1); a0=quadl(singrect,-2,2)/T;%计算a0 N=10;an=zeros(1,N);bn=zeros(1,N); for k=1:N an(k)=quadl(@rectcos,-2,2,[],[],k,w)*2/T; %计算an; end for k=1:N bn(k)=quadl(@rectsin,-2,2,[],[],k,w)*2/T; %计算bn; end; n=1:1:N; figure(1); subplot(1,2,1);plot(n,an);grid on; subplot(1,2,2);plot(n,bn);grid on; % 有限项级数逼近 t=-6:0.01:6; A0=a0; AN=zeros(1,N);fiN=zeros(1,N); AN=sqrt(an.^2+bn.^2); fiN=-atan(bn./an); figure(2);subplot(6,2,1); x=pulstran(t,-8:4:8,rectpuls,2); plot(t,x);axis([-8,8,-1,2]);grid on; subplot(6,2,2);plot(t,A0/2);grid on; wave=a0/2; for k=1:12 wave=wave+an(k)*cos(k*w*t+fiN(k)); subplot(6,2,k+2);plot(t,wave);grid on; end 程序中“rectcos”和“rectsin”分别为所预先定义的函数文件。其中rectcos.m文件为： function y=rectcos(t,n,w); y=(abs(t)=1).*1.*cos(n*w*t); 定义了函数与余弦函数之间的乘积。 rectsin.m文件为： function y=rectsin(t,n,w); y=(abs(t)=1).*1.*sin(n*w*t); 定义了函数与正弦弦函数之间的乘积。 运行结果：傅里叶系数图形： 有限项数逼近图形，采用前10次谐波： (2)周期信号的频谱 周期信号经过傅里叶分解可表示为一系列正弦或复指数信号之和。为了直观地表示出信号所含各分量的振幅，以频率（或角频率）为横坐标，以各谐波的振幅或虚指数函数的幅度为纵坐标，可画出幅度－频率关系图，称为幅度频谱或幅度谱。类似地，可画出各谐波初相角与频率的关系图，称为相位频谱或相位谱。 在计算出信号的傅里叶分解系数后，就可以直接求出周期信号的频谱并画出其频谱图。 (3)非周期信号的傅里叶变换和性质 非周期信号的傅里叶变换定义为： 称为频谱密度函数，一般需要用幅度谱和相位谱两个图形才能将它完全表示出来。 Matlab中symbolic工具箱提供了直接求解信号的傅里叶变换和逆变换的函数fourier()和ifourier()。这两个函数采用符号运算方法，在调用之前要