compare.m

关于信号传输,信号处理和滤波器生成的matlab源代码.(压缩包内附详细的readme和程序注释) 找不到关于"信号处理"的分类,所以我归在了"matlab例程"之下,望谅解. 感谢您们关注我

%预加重(预滤波) 信号在噪声信道内传播 及 解码 模拟算法 与 不加重信号传输的结果比较
%程序设计:刘辉
[e-mail] (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)');