%預加重(預濾波) 信號在噪聲信道內傳播 及 解碼 模擬算法 與 不加重信號傳輸的結果比較
%程序設計:劉輝
%pangelaca@yahoo.com.cn (email/msn)
%由于我在柏林工業大學學習,源程序的一些變量名由德語縮寫成.(但其實閱讀上沒有影響)

%參數說明:
%預加重和輸入反制的系數,與PreAndDeemphae相同
%程序比較了經預加重,和不經預加重直接產生的信號,在輸出時的差異.從聲音上可以明顯聽出,從波形和誤差比較上也可以得出.

%Hpre(z)=(5z-3)/2z
Bpre=[5 -3];
Apre=[2 0];

[x,sampl]=wavread ('ftzm1.wav');
subplot(3,2,1);
plot(x);
title ('原始聲波 x(n) (ftzm1.wav)')
sound(x);

pause(3);

xpre=filter(Bpre,Apre,x);
xvar=var(x);
[m,n]=size(x);
r=xvar*randn(m,1);
Mixsignal=xpre+r;

%Hde(z)=2z/5z-3
Bde=[2 0];
Ade=[5 -3];

y=filter(Bde,Ade,Mixsignal);
subplot(3,2,3);
plot(y);
title ('通過預加重,反制加重所產生的輸出聲波 y(n)')
sound(y);

Fehler=y-x;
subplot(3,2,5);
plot(Fehler);
title ('誤差  e(n)=y(n)-x(n)')

pause(3);


%Ohne Pre_De
subplot(3,2,2);
plot(x);
title ('原始聲波 x(n) (ftzm1.wav)')
sound(x);

pause(3);

y2=x+r;
subplot(3,2,4);
plot(y2);
title ('直接通過噪聲信道,得到的輸出聲波 y(n)')
sound(y2);

subplot (3,2,6);
plot (r);
title ('誤差 e(n)=y(n)-x(n)=r(n)');