资源描述：
《matlab语言编程实现回波信号的产生和消除》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、山东建筑大学信电学院课程设计说明书目录1·摘要…………………………………………………………………………22·正文…….......................................................................................................3(1)设计目的和要求………………………………………………….…....3(2)设计原理………………………………………………………….…….3(3)设计内容………………………………………………………………..33·总结与致谢……………………………………
2、…………………………....114·参考文献……………………………………………………………………12-12-山东建筑大学信电学院课程设计说明书摘要回波是现代通信中影响通信质量的噪声干扰，也有测距估算等的使用价值。MATLAB语言具备高效、可及推理能力强等特点,在数值计算方法、图形功能、用户界面设计编程手段和工具等方面有着重要的应用，随着MATLAB应用领域不断扩大，使其成为一种影响大、流行广的科学计算语言。本文基于MATLAB以回波这一最简单的线性模型进行分析，实现了回波的产生、消除、测量反射物的距离。本设计方案用MATLAB自带的函数功能读取语音信号，通过延迟叠加产
3、生回波信号来模拟回声，通过计算分析输出信号的自相关函数利用基于自相关函数的方法来估计衰减系数和延迟，并利用快速傅里叶变换将回波信号滤除回声，得到原声音信号，同时根据延迟估计反射物的距离。关键字：回波自相关函数MATLAB延迟距离-12-山东建筑大学信电学院课程设计说明书正文（1）设计目的和要求本课程设计要求用MATLAB语言编程实现回波信号的产生和消除。具体要求：（1）利用声音信号x产生带有回声的声音文件y。（2）从带有回声的文件y中消除回声。（3）从y中估计反射物的距离。（2）设计原理1·根据设计要求分析系统功能，掌握设计中所需理论（采样频率、采样位数的概念，采样定
4、理：回声信号，相关性分析；数字滤波器设计原理和方法）。采样频率：采样频率，也称为采样速度或者采样率，定义了每秒从连续信号中提取并成离散信号的采样个数。采样位数：即采样值或取样值，用来衡量声音波动变化的参数，是指声卡在采集和播放声音文件时所使用数字声音信号的二进制位数。采样定理：在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中最高频率fmax的2倍时(fs.max>=2fmax)，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。相关函数：X(t)为随机过程，a(t)=EX(
5、t)为期望，则相关函数定义为：B(s,t)=E(X(s)-a(t))(X(t)-a(t))。若X(t)=Y(t)+i*Z(t)，Y,Z为实过程，则称X(t)为复随机过程，相关函数定义为：B(s,t)=E(X(s)-a(t))(X(t)-a(t))（后一个括号取共轭）相关函数两个本质特性：1）共轭对称：B(s,t)=B(t,s)的共轭2）非负定：对任意的n>=1，t1……tn属于T,n个复数Z1,Z2……，Zn，有。-12-山东建筑大学信电学院课程设计说明书2·熟悉MATLAB的编程语言。（3）设计内容1·信号的采集录制一段音频（如歌曲，说话声等），绘制其时域波形，对此
6、音频信号用FFT作频谱分析。我截取了一段自己的手机铃声，转换为WAV格式，重命名为gwm，输入代码并得到图谱如下：[x,fs]=wavread('gwm');figure(1);subplot(3,1,1);plot(x(1:4000));title('原始信号');y=fft(x);subplot(3,1,2);plot(abs(y));title('幅值');-12-山东建筑大学信电学院课程设计说明书2·回声信号的产生在录制信号基础上产生带有回声的声音信号并绘制其时域波形，对其进行FFT频谱分析，绘制频谱图。本段输入代码及得到图谱如下：[x,fs]=wavread
7、('gwm');x1=x(1:4000);x2=x(1:4000);x3=x(1:4000);x1=[x1,zeros(1,1800)];x2=[zeros(1,600),0.7*x2,zeros(1,1200)];x3=[zeros(1,1200),0.4*x3,zeros(1,600)];y=x1+x2+x3;figure(2);subplot(3,1,1);-12-山东建筑大学信电学院课程设计说明书plot(y(1:5800));title('回声');y1=fft(y);subplot(3,1,2);plot(abs(y1));title('
