%---------------设计和使用均衡器改善信道----------------
%-------改变k的值分别完成对k=1和k=3时的均衡器设计-------
clc;
clear all;
f=[0.407 0.815 0.407];
w=-pi:2*(1/512)*pi:pi;
e=[exp(-j*w*1);exp(-j*w*2);exp(-j*w*3)];
F0=f*e;%傅立叶变换
F1=20*log10(abs(F0));%dB转换，信道频谱
D0=(1/(f(2)))*((abs(f(1))+abs(f(3))));%均衡器输入峰值失真
%-------设计均衡器抽头系数-------((m-n)>=2)
k=1;%抽头系数为2k+1
for m=1:(2*k+1)
for n=1:(2*k+1)
if (((m-n)<-1)|((m-n)>=2))
F2(m,n)=0;
else F2(m,n)=f(m-n+2);
end
end
end
F3=F2^(-1);%求逆矩阵或用inv（）函数
q=zeros(2*k+1,1);
q(k+1)=1;
C=F3*q;
Q=conv(C,f);%等效信道模型
total=0;
for kk=1:length(Q)
total=total+(abs(Q(kk)));
end
D1=(1/(abs(Q(k+2))))*(total-(abs(Q(k+2))));%均衡器输出峰值失真
F4=0;
for n1=1:(2*k+1)
F4=F4+(exp(-j*w*n1))*C(n1);%傅立叶变换
end
F5=20*log10(abs(F4));%均衡器频谱
F6=20*log10(abs(F0.*F4));%等效信道频谱
plot(w,F1,':b',w,F5,'-.r',w,F6,'g-');
xlabel('频率(w)');ylabel('幅度谱/dB');
title('抽头系数为3（K=1）均衡器的幅度谱特性曲线');
legend('信道','均衡器','等效信道')
评论0