12通信实验指导书.doc

实验 连续时间信号与系统的频域分析 一、 实验目的 1、学会用MATLAB实现连续时间信号傅里叶变换 2、学会用MATLAB分析LTI系统的频域特性 3、学会用MATLAB分析LTI系统的输出响应 二、实验原理 1．傅里叶变换的MATLAB求解 MATLAB的symbolic Math Toolbox 提供了直接求解傅里叶变换及逆变换的函数fourier()及ifourier()两者的调用格式如下。 Fourier 变换的调用格式 F=fourier(f)：它是符号函数f的fourier变换默认返回是关于w的函数。 F=fourier(f，v):它返回函数F是关于符号对象v的函数，而不是默认的w，即 Fourier逆变换的调用格式 f=ifourier(F):它是符号函数F的fourier逆变换，默认的独立变量为w，默认返回是关于x的函数。 f=ifourier(f,u):它的返回函数f是u的函数，而不是默认的x. 注意：在调用函数fourier()及ifourier()之前，要用syms命令对所用到的变量（如t,u,v,w）进行说明,即将这些变量说明成符号变量。 例-1 求的傅立叶变换 解: 可用MATLAB解决上述问题： syms t Fw=fourier(exp(-2*abs(t))) 例-2 求的逆变换f(t) 解： 可用MATLAB解决上述问题 syms t w ft=ifourier(1/(1+w^2),t) 2．连续时间信号的频谱图 例-3 求调制信号的频谱，式中 解：MATLAB程序如下所示 ft=sym(4*cos(2*pi*6*t)*(Heaviside(t+1/4)-Heaviside(t-1/4))); Fw=simplify(fourier(ft)) subplot(121) ezplot(ft,[-0.5 0.5]),grid on subplot(122) ezplot(abs(Fw),[-24*pi 24*pi]),grid fourier( )和ifourier( )函数的一个局限性是，如果返回函数中有诸如单位冲激函数等项，则用ezplot()函数无法作图。对某些信号求变换时，其返回函数可能包含一些不能直接用符号表达的式子，因此不能对返回函数作图。此外，在很多实际情况中，尽管信号是连续的，但经过抽样所获得的信号则是多组离散的数值量，因此无法表示成符号表达式，此时不能应用fourier()函数对f(n)进行处理，而只能用数值计算方法来近似求解。 从傅里叶变换定义出发有， 当足够小时，上式的近似情况可以满足实际需要。对于时限信号，或者在所研究的时间范围内让衰减到足够小，从而近似地看成时限信号，则对于上式可以考虑有限n的取值。假设是因果信号，则有 傅里叶变换后在域用MATLAB进行求解，对上式的角频率进行离散化。假设离散化后得到N个样值，即 －1，是频率取样点 时间信号取样间隔应小于奈奎斯特取样时间间隔，若不是带限信号可根据计算精度要求确定一个频率 W0为信号的带宽。 因此有 。采用行向量，用矩阵表示为 。其要点是要正确生成的M个样本向量与向量。当足够小时，上式的内积运算（即相乘求和运算）结果即为所求的连续时间信号傅里叶变换的数值解。 例-4 用数值计算法求信号的傅里叶变换 解信号频谱是，第一个过零点是，一般将此频率视为信号的带宽，若将精度提高到该值的50倍，既W0=50，据此确定取样间隔， 3．用MATLAB分析LTI系统的频率特性 当系统的频率响应H（jw）是jw的有理多项式时，有 MATLAB信号处理工具箱提供的freqs函数可直接计算系统的频率响应的数值解。其调用格式如下 H=freqs(b,a,w) 其中，a和b分别是H(jw)的分母和分子多项式的系数向量，w为形如w1:p:w2的向量，定义系统频率响应的频率范围，w1为频率起始值，w2为频率终止值,p为频率取样间隔。H返回w所定义的频率点上，系统频率响应的样值。 例如，运行如下命令，计算0~2pi频率范围内以间隔0.5取样的系统频率响应的样值 a=[1 2 1]; b=[0 1]; h=freqs(b,a,0:0.5:2*pi) 例 -5 三阶归一化的butterworth 低通滤波器的频率响应为 试画出该系统的幅度响应和相位响应。 解 其MATLAB程序及响应的波形如下 w=0:0.025:5; b=[1];a=[1,2,2,1]; H=freqs(b,a,w); subplot(2,1,1); plot(w,abs(H));grid; xlabel(\omega(rad/s)); ylabel(|H(j\omega)|); title(H(jw)的幅频特性); subplot(2,1,2);