数字信号处理实验三.docVIP

数字信号处理实验三 实验三 离散时间信号的频域分析 实验室名称： 信息学院2204 实验时间：2015年10月15日 姓 名：蒋逸恒 学号：20131120038 专业： 通信工程 指导教师：陶大鹏 成绩 教师签名： 年 月 日 的实数部分，XEF的虚部应该为零。你能验证它们吗？你怎样解释仿真结果？ 三、实验器材及软件 1. 微型计算机1台 MATLAB 7.0软件 实验原理 3.1；3.2；3.3；3.4 离散时间傅立叶变换的结果是关于w的连续函数，对于系统函数的离散时间傅立叶变换的求法是，其中，B是f序列傅立叶变换的系数，A是y序列傅立叶变换的系数。离散时间傅立叶变换的结果是w的周期函数，在（2k+1）π附近为高频，在2kπ附近为低频（k=0，+1，-1，+2，-2。。。。） 3.6 离散时间傅立叶变换的时移特性： 3.10 离散时间傅立叶变换的频移特性： 3.14；3.15 离散时间傅立叶变换的卷积性质： 3.17 离散时间傅立叶变换的调制特性： 3.20 离散时间傅立叶变换的反转特性： 3.23 在matlab中，fft（）函数可以快速的计算有限长序列的离散傅立叶变换，ifft（）函数可以快速的计算离散傅立叶逆变换，对于计算中的不同序列长度N，若把时间当作1s，则N相当于采样率Fs，L是傅立叶变换后的序列的长度。此时，采样点的频率可表示为Fn=(n-1)*Fs/L，当N与L越接近，Fs/L越小，Fn的变化速度越慢，此时相位谱也就相应的变化减慢，因为相位是频率f的一次函数。 3.26；3.27；3.28；3.29 圆周移位函数和圆周卷积函数都是在“圆周”上循环的，该圆周的长度就是序列的长度。 3.30；3.31；3.32；3.33 圆周时移实际上是把一个序列的后面的点按顺序搬到前面来， 这里与反转和线性时移有着完全的区别。圆周时移实际上的移动范围不会超过序列长 度值。圆周时移性质：若，则，其中，。 3.36；3.38 由实验我们可以知道一个圆周卷积性质：线性卷积可通过圆周卷积得到。 3.41；3.42 由教材可知：，即序列的偶部分的傅立叶变换是序列的傅立叶变换的实部。 实验步骤 进行本实验，首先必须熟悉matlab的运用，所以第一步是学会使用matlab。 学习相关基础知识，根据《数字信号处理》课程的学习理解实验内容和目的。 在充分熟悉基础知识的情况下进行实验，利用matlab完成各种简单的波形产生和观察，理解各种波形产生的原理和方法。 从产生的图形中学习新的知识，掌握实验的目的，充分学习数字信号处理的运用。 最后需要思考各种波形的联系和建立完整的知识体系，如整理噪声和原波形之间的叠加关系等。 六、实验记录（数据、图表、波形、程序等） 3.2 w=-4*pi:8*pi/511:4*pi; num=[2 1];den=[1 -0.6]; h=freqz(num,den,w); subplot(2,1,1); plot(w/pi,real(h));grid; title(H(e^{j\omega})的实部); xlabel(\omega/\pi); ylabel(振幅); subplot(2,1,2); plot(w/pi,imag(h));grid; title(H(e^{j\omega})的虚部); xlabel(\omega/\pi); ylabel(振幅); pause; subplot(2,1,1); plot(w/pi,abs(h));grid; title(|H(e^{j\omega}|幅度谱); xlabel(\omega/\pi); ylabel(振幅);subplot(2,1,2); plot(w/pi,angle(h));grid; title(相位谱arg[H(e^{j\omega})]); xlabel(\omega/\pi); ylabel(以弧度为单位的相位); 3.3 clf; w=0:8*pi/511:pi; num=[0.7 -0.5 0.3 1];den=[1 0.3 -0.5 0.7]; h=freqz(num,den,w); subplot(2,1,1); plot(w/pi,real(h));grid; title(H(e^{j\omega})的实部); xlabel(\omega/\pi); ylabel(振幅); subplot(2,1,2); plot(w/pi,imag(h));grid; title(H(e^{j\omega})的虚部); xlabel(\omega/\pi); ylabel(振幅); pause; subp