程控實驗報告cc(共6頁)

1、實驗報告實驗(shyn)內容簡介電話(dinhu)交換系統中的號碼(hom)分析交換機時隙交換回波消除器小組成員及分工實驗記錄回波消除器實驗是我主要負責的內容。我們花了比較長的一段時間來確定究竟實驗要達到的效果是什么。A向B發(fā)送話音信號，需要將A的回聲與B的話音信號疊加嗎？是應該把A的話音信號與回聲疊加，然后消除回聲只留下原話音信號，還是只需要把話音信號延遲與回聲相減？A是否知道回聲的延遲時間(也就是那個RTT)？詢問過老師以后，我們決定采用假設：A不知道RTT，通過將原話音信號與回聲信號進行相關函數分析確定延時，回聲信號不與任何其他信號疊加，僅將原話音信號延遲后減去回聲信號。實驗流程為：采集

2、(cij)的原始(yunsh)話音(huyn)信號與收到的回聲信號：原始(yunsh)話音(huyn)信號與收到的回聲(hushng)信號的互相關函數：將原始(yunsh)話音信號延時RTT時間(shjin)后與收到的回聲(hushng)信號相減，正好能完全抵消：存儲的文件信息：實驗代碼fs = 8000; %采樣頻率N=2400; %0.3秒的延時x = wavrecord(3*fs, fs, double); %wacrecord( 采樣點數， 采樣頻率， 采樣點的數據類型 )wavwrite(x, 原始信號); %存儲音頻subplot(211);x1 = x;zeros(N,1);pl

3、ot(x1);grid on;xlabel(時間(shjin);ylabel(幅值);title(原始信號(xnho)時域波形);y=zeros(N,1);x; %y=x(n-N)wavwrite(y, 回聲(hushng)信號);subplot(212);plot(y);grid on;xlabel(時間);ylabel(幅值);title(回聲信號時域波形);R = xcorr(x1,y);figureplot(R);value,index = max(R);d = abs(index - (3*fs+N) %延遲點數delay = d/fs %延遲時間（RTT）grid on;title

4、(原始信號與回聲信號的互相關函數);x2 = zeros(d,1);x; %x2=x(n-d)z = x2 - y; %z是原信號延遲d個單位后與回聲相減的信號figureplot(z);grid on;xlabel(時間);ylabel(幅值);title(原信號延遲RTT時間后與回聲信號相減的時域波形);實驗心得這次實驗做出結果其實并不難，難的就在對題意的理解是否正確。開始時聽一些同學說就設計兩個數組，按位相減就可以了，但我一直認為不對，因為這樣完全沒有使用價值，接收端并不能知道延遲時間，也就是說從哪一位開始按位相減，它是不知道的。理解到了題目的意思后，開始編寫程序，首先遇到的困難是對ma

5、tlab內置庫函數不熟悉，不知道怎么調用，每次都是按照功能在百度上搜索出函數，再查看這個函數的參數說明，最后調用測試是否完成了我們需要的功能。最大的困擾出現在調用互相關函數xcorr()時，得到的結果總是不是我們認為正確的那個，所以嘗試了很多其他互相關函數比如corrcoef()、更改xcorr()參數選項、還有按步分解成傅里葉變換的乘積的逆傅里葉變換等等方法，做出來結果總是覺得不對，弄了一晚上，發(fā)現在進行互相關函數時是將原函數周期拓展以后再與另一個函數進行相關性比較的，所以直接使用xcorr()得出的結果應該是正確的，所以問題就解決了。使用MATLAB都是用序列、矩陣進行運算，長度（點數）是

6、否一致直接導致程序是不是能跑通，在這個問題上也有過一點阻礙，不過還好都及時發(fā)現解決了。這個實驗讓我更清楚了回聲消除器的原理，讓我更加熟悉了MATLAB的操作，基本完成了回聲消除的功能，不過還是有點瑕疵就是沒有考慮衰減，因為相關函數是不反折的卷積，所以衰減后相關函數與原來未衰減的相關函數之間有什么系數關系弄不清楚，所以沒有做出來。如果在數學方面多進行一些學習應該能做出來，不過對于實驗來說應該測出延時就夠了。內容總結（1）實驗報告實驗內容簡介電話交換系統中的號碼分析交換機時隙交換回波消除器小組成員及分工實驗記錄回波消除器實驗是我主要負責的內容（2）我們花了比較長的一段時間來確定究竟實驗要達到的效果是什么（3）grid o