信号与系统实验 连续信号的频域特性分析.doc

信号与系统 课程实验报告实验三实验名称连续信号的频域特性分析系 别教师姓名实验地点5309实验日期2011-06-20学生姓名学号一、实验内容1．求门函数的傅立叶变换，并画出其频率特性曲线图2．已知频率特性函数为：求其幅频特性和相频特性二、实验目的1．通过实验掌握傅立叶变换的Matlab实现，加深对傅立叶变换基本定义的理解2．掌握用Matlab软件绘制信号频谱的方法三、涉及实验的相关情况介绍（包含使用软件或实验设备等情况）计算机一台(安装MATLAB6.5版本或以上版本)四、实验试做记录（含程序、数据记录及分析）在Matlab中求解时域信号的傅立叶变换常用函数QUADL实现：（1）调用格式：Q = QUADL(‘F’,A,B)（2）调用格式：Q = QUADL (‘F’, A,B,[ ],[ ],P)函数QUADL可以实现F函数的积分运算，F为一个字符串，它表示被积函数的文件名A,B分别表示定积分的下限和上限，[ ]为空表示积分递归过程中的绝对误差选用默认值（默认值为），P表示被积函数中的一个参数，对每个P值计算一次积分1．求门函数的傅立叶变换，并画出其频率特性曲线图程序：%fly.mfunction y=fly(t,w)y=(heaviside(t+2)-heaviside(t-2)).*exp(-j*w*t);clearw=linspace(-5*pi,5*pi,512);N=length(w);F=zeros(1,N);for k=1:N F(k)=quadl('fly',-2,2,[],[],w(k));endplot(w,real(F));xlabel('w');ylabel('F(w)');title('f(t)的傅里叶变换F(w)');结果：2．已知频率特性函数为：求其幅频特性和相频特性。

程序：clearw=linspace(0,5,250);b=[0 2 0 1 4];a=[1 3 2 5 2];H=freqs(b,a,w);subplot(211);plot(w, abs(H));set(gca, 'xtick', [0 1 2 3 4 5]);grid;xlabel('\omega');ylabel('|H(j\omega)|');subplot(212);plot(w, angle(H));set(gca, 'xtick', [0 1 2 3 4 5]);grid;xlabel('\omega');ylabel('phi(j\omega)');结果：五、实验总结通过本实验，我了解了傅里叶变换以及信号的幅频曲线和相频曲线的由来，特别是通过调用quadl函数，从定义上理解了傅里叶变换的过程实验表明，熟练使用Matlab便于我们对“信号与系统”这门课程加深理解，同时也有利于高效的处理信号变换。