调幅信号的解调.docxVIP

《信号与系统》课程工程设计问题设计报告题目：_____调幅信号的解调_____分组号： _____________组 长： __________组 员：__________________________________时 间： _______ 指导教师：报告日期：年月日报告撰写说明按照参考模板的内容和格式撰写报告理论模型部分须结合本课程知识分析问题、建立模型程序设计部分应给出设计思路、主要流程图和关键函数的说明；结果分析不能只是简单给出结论，应结合具体问题，对关键参数或算法在不同取值条件下对结果的影响情况进行分析和总结。如果可能，还应进行误差分析组内互评分A、B、C、D四个贡献等级，最终评价应区分出前三个等级在规定时间内，完成叙述并回答问题。目录1 问题描述22 理论模型22.1 原理分析与设计思路22.2 数学模型23 程序设计23.1 编程思路24 组内互评35 总结与体会3参考文献3附录 MATLAB程序主要代码41 问题描述文件project.wav中包含了一段错误解调的音频信息的采样值，原始信号是以fc作为载波频率进行调制的，但在解调是却用了错误的频率fc≠fc，进行相干解调原始信号的带宽为fb=4，project.wav是对错误解调后得到的连续时间信号进行采样得到的。采样频率为fs，且fs远大于fb和|fc-fc|。本文则根据上述信息恢复出正确语音信号。2 理论模型2.1 原理分析与设计思路相干解调是指利用乘法器，输入一路与载频相干（同频同相）的参考信号与载频相乘。比如原始信号 A 与载频 cos(ωt + θ) 调制后得到信号 Acos(ωt + θ)；解调时引入相干（同频同相）的参考信号 cos(ωt + θ)，则得到：Acos(ωt+θ)cos(ωt+θ)，利用积化和差公式可以得到A*1/2*[cos(ωt+θ+ωt+θ)+cos(ωt+θ-ωt-θ)]=A*1/2*[cos(2ωt+2θ)+cos(0)]=A/2*[cos(2ωt+2θ)+1]=A/2+A/2cos(2ωt+2θ)利用低通滤波器将高频信号cos(2ωt+2θ)滤除，即得原始信号 A。而本题采用cos(ω0t)乘以原始信号X1(t),则X2(t)=1/2X1(t)+1/2X1(t)cos(2ω0t)。如果我们在进行解调时使用了错误的频率ω0≠ω0，设=ω0-ω0，则解调得出的错误信号乘以cos(ωSt)进行相干解调再通过低通滤波器则可以得出原信号，如果有混叠现象则需要先通过高通滤波器进行处理，所以本题的关键在于寻找到ωS。2.2 数学模型巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。贝赛尔(Bessel)滤波器是具有最大平坦的群延迟（线性相位响应）的线性过滤器。切比雪夫滤波器是在通带或阻带上频率响应幅度等波纹波动的滤波器，振幅特性在通带内是等波纹。当滤波器具有相同阶数时：　巴特沃斯滤波器通带最平坦，阻带下降慢。　切比雪夫滤波器通带等纹波，阻带下降较快。　贝塞尔滤波器通带等纹波，阻带下降慢。也就是说幅频特性的选频特性最差。但是，贝塞尔滤波器具有最佳的线性相位特性。　　此外，还有椭圆滤波器，椭圆滤波器在通带等纹波（阻带平坦或等纹波），阻带下降最快。 本题采用巴特沃斯滤波器，根据查表求出设计的高通与低通巴特沃斯滤波器均为8阶，其中归一化的巴特沃斯多项式因子Bn（s）为3 程序设计3.1 编程思路信号的频谱对题中所给的声音信号进行傅里叶变换得到频谱图。由图可知零振幅点的平移量3000hz，则在原理分析中所要求的ωS可得出，fd=3000hz。因为平移量为fc和fc差的绝对值，所以不能判断两者的大小，也不会对错误解调的结果产生影响。（2）滤波器设计使用matlab设计滤波器，其中高阶滤波器截止频率为3000hz，低通滤波器截止频率为4000hz。使用函数tf构建上述两个滤波器的线性系统，得到系统传递函数：高通滤波器H_H = s^8 --------------------------------------------------------------------- s^8 + 9.662e04 s^7 + 4.668e09 s^6 + 1.463e14 s^5 + 3.243e18 s^4 + 5.199e22 s^3 + 5.893e26 s^2 + 4.334e30 s + 1.594e34低通滤波器H_L = 5.47e34 --------------------------------