正弦信号谱分析和提取.docVIP

西南科技大学 课 程 设 计 报 告 课程名称： 数字通信课程设计 设计名称： 正弦信号的谱分析及提取 姓 名： 学 号: 班 级： 指导教师： 起止日期： 2010.6.28-2010.7.9 西南科技大学信息工程学院制课 程 设 计 任 务 书 学生班级： 学生姓名： 学号： 设计名称： 正弦信号的谱分析及提取 起止日期：2010.6.28-2010.7.9 指导教师： 设计要求： 1、采用matlab或者其他软件工具产生不同频率，不同幅度的两种正弦波信号与，并将这两个信号叠加为一个信号，观察信号的波形。 2、对叠加后的信号用FFT作谱分析。 3、选择合适的指标，设计FIR数字滤波器，从信号中提取信号与。 课 程 设 计 学 生 日 志 时间 设计内容 2010.6.29 查阅相关资料,看《matlab》、《现代通信原理》、《数字信号处理》等书籍 2010.7.3 根据设计要求，确定实现原理，实现方案 2010.7.5 根据前阶段所做的工作开始书写程序，并调试 2010.7.8 完成设计报告的书写，准备答辩 2010.7.9 去老师处答辩 课 程 设 计 考 勤 表 周 星期一 星期二 星期三 星期四 星期五 课 程 设 计 评 语 表 指导教师评语： 成绩： 指导教师： 年 月 日 一、正弦信号的谱分析及提取 设计目的和意义 1、熟练掌握MATLAB在数字信号处理中的应用； 2、加深对FFT变换的理解，并掌握FFT变换； 3、掌握信号的合成和提取； 4、熟练运用各种数字滤波器恢复原始信号； 5、理解MATLAB工具箱的应用。 设计原理 如题，首先要产生两个正弦信号，而且按扩展要求，信号的幅度和频率还要用户可以自定义，这里可以直接用命令产生需要的信号。然后将这两个正弦信号进行叠加得到目标信号，就可以直接用MATLAB自带的fft(x,n)做FFT变换进行频谱分析。最后按要求，还要进行滤波恢复出原始信号，显然地，这里需要两个不同的滤波器来进行滤波，分别得到两个原始信号x1和x2。 整个实验原理的原理图如下： 图1 同样MATLAB自带有fftfilt(b,x)可以进行滤波仿真，设计FIR数字滤波器时用窗函数法，关键的是要选取合适的窗函数，并对其参数进行合理设置，按要求要用户输入数字滤波器的技术指标，这样就可以仿真恢复出原始信号。 详细设计步骤 1、产生两个正弦波信号，信号幅度和频率要求用户可以自定义。 通过前期做数字信号处理实验时就已经知道，可以利用输入命令input( )来输入信号指标，同时MATLAB自带有可以直接产生正弦信号波形的函数sin(x)。所以第一步用于产生所需幅度和频率的信号比较简单。 2、产生FFT变换所需波形，即x1(n)和x2(n)MATLAB自带有可以进行快速傅里叶变换所需的函数y=fft(x)，或y=fft(x，n)。x为所要进行FFT变换的信号。这里n的长度若为2的整数次幂，则MATLAB将自动进行基2的FFT算法；否则，将采用较慢的分列基算法。所以这里需要注意进行FFT变换的点数需要为2的整数次幂。对x3(n)进行FFT变换后，求其幅度做出其频谱图。在这里我将采样频率简化为2倍信号的最高的频率，满足奈奎斯特抽样定律。 4、设计满足要求的FIR数字滤波器恢复原始信号。 用窗函数法设计FIR数字滤波器。通过查阅相关资料后知道可以用函数fftfilt模拟信号经过相关滤波器，即信号的恢复。还有用ceil函数用以求滤波器的阶数。这里选用汉