程控实验报告cc.docxVIP

实验报告 实验内容简介 电话交换系统中的号码分析 交换机时隙交换 回波消除器 小组成员及分工 实验记录 回波消除器实验是我主要负责的内容。我们花了比较长的一段时间来确定究竟实验要达到的效果是什么。A向B发送话音信号，需要将A的回声与B的话音信号叠加吗？是应该把A的话音信号与回声叠加，然后消除回声只留下原话音信号，还是只需要把话音信号延迟与回声相减？A是否知道回声的延迟时间(也就是那个RTT)？ 询问过老师以后，我们决定采用假设：A不知道RTT，通过将原话音信号与回声信号进行相关函数分析确定延时，回声信号不与任何其他信号叠加，仅将原话音信号延迟后减去回声信号。 实验流程为： 采集的原始话音信号与收到的回声信号： 原始话音信号与收到的回声信号的互相关函数： 将原始话音信号延时RTT时间后与收到的回声信号相减，正好能完全抵消： 存储的文件信息： 实验代码 fs = 8000; %采样频率 N=2400; %0.3秒的延时 x = wavrecord(3*fs, fs, double); %wacrecord( 采样点数， 采样频率， 采样点的数据类型 ) wavwrite(x, 原始信号); %存储音频 subplot(211); x1 = [x;zeros(N,1)]; plot(x1); grid on; xlabel(时间);ylabel(幅值);title(原始信号时域波形); y=[zeros(N,1);x]; %y=x(n-N) wavwrite(y, 回声信号); subplot(212); plot(y); grid on; xlabel(时间);ylabel(幅值);title(回声信号时域波形); R = xcorr(x1,y); figure plot(R); [value,index] = max(R); d = abs(index - (3*fs+N)) %延迟点数 delay = d/fs %延迟时间（RTT） grid on; title(原始信号与回声信号的互相关函数); x2 = [zeros(d,1);x]; %x2=x(n-d) z = x2 - y; %z是原信号延迟d个单位后与回声相减的信号 figure plot(z); grid on; xlabel(时间);ylabel(幅值);title(原信号延迟RTT时间后与回声信号相减的时域波形); 实验心得 这次实验做出结果其实并不难，难的就在对题意的理解是否正确。开始时听一些同学说就设计两个数组，按位相减就可以了，但我一直认为不对，因为这样完全没有使用价值，接收端并不能知道延迟时间，也就是说从哪一位开始按位相减，它是不知道的。理解到了题目的意思后，开始编写程序，首先遇到的困难是对matlab内置库函数不熟悉，不知道怎么调用，每次都是按照功能在百度上搜索出函数，再查看这个函数的参数说明，最后调用测试是否完成了我们需要的功能。最大的困扰出现在调用互相关函数xcorr()时，得到的结果总是不是我们认为正确的那个，所以尝试了很多其他互相关函数比如corrcoef()、更改xcorr()参数选项、还有按步分解成傅里叶变换的乘积的逆傅里叶变换等等方法，做出来结果总是觉得不对，弄了一晚上，发现在进行互相关函数时是将原函数周期拓展以后再与另一个函数进行相关性比较的，所以直接使用xcorr()得出的结果应该是正确的，所以问题就解决了。使用MATLAB都是用序列、矩阵进行运算，长度（点数）是否一致直接导致程序是不是能跑通，在这个问题上也有过一点阻碍，不过还好都及时发现解决了。这个实验让我更清楚了回声消除器的原理，让我更加熟悉了MATLAB的操作，基本完成了回声消除的功能，不过还是有点瑕疵就是没有考虑衰减，因为相关函数是不反折的卷积，所以衰减后相关函数与原来未衰减的相关函数之间有什么系数关系弄不清楚，所以没有做出来。如果在数学方面多进行一些学习应该能做出来，不过对于实验来说应该测出延时就够了。