迫零均衡器--代码

1、程序 1clear allN=10000;x=randn(1,N); %产生均值 0, 方差为 1的随机序列for i=1:Nif x(i)>b(i)=1;elseb(i)=-1;endend%调制 xx=pskmod(x,M);y= 0 ;%衰减信道A1=zeros(5,5);for m=1:5for n=1:5if n=mA1(m,n)=;else if n=m+1A1(m,n)=;else if n=m-1A1(m,n)=;else if n=m+2A1(m,n)=;else if n=m-2A1(m,n)=0;elseA1(m,n)=0;endendendendendendend

2、B1=0;0;1;0;0;C1=pinv(A1)*B1;A2=zeros(17,17);for i=1:17for j=1:17if j=iA2(i,j)=;else if j=i-1A2(i,j)=;else if j=i-2A2(i,j)=0;elseifj=i-3A2(i,j)elseifj=i-4A2(i,j)elseifj=i-5A2(i,j)elseifj=i+1A2(i,j)elseifj=i+2A2(i,j)elseifj=i+3A2(i,j)elseifj=i+4A2(i,j)elseifj=i+5A2(i,j)elseA2(i,j)=0;endendendendendend

3、endendendendendendendB2=0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;C2=pinv(A2)*B2;A=zeros(31,31);for i=1:31for j=1:31if j=iA(i,j)=;else if j=i-1A(i,j)=;else if j=i-2 A(i,j)=0;else if j=i-3 A(i,j)=;else if j=i-4A(i,j)=;else if j=i-5 A(i,j)=;else if j=i+1 A(i,j)=;else if j=i+2 A(i,j)=;else if j=i+3 A(i,j)=;els

4、e if j=i+4 A(i,j)=;else if j=i+5 A(i,j)=;elseA(i,j)=0;end end end end end end end end end end end end endB=0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0; C=pinv(A)*B;SNR=1:2:21;for m=1:length(SNR)for i=1:N noise=randn; x1noise(i)=b(i)*sqrt(2*10(SNR(m)/10)+noise;if x1noise(i)<0;t1(i)

5、=-1;elset1(i)=1;endend%解调： X1=pskdemod(x1noise,M);numbers,pm1(m)=symerr(b,t1); xxy=conv(b,y);for i=1:N+6noise=randn; x2noise(i)=xxy(i)*sqrt(2*10(SNR(m)/10)+noise;if x2noise(i)<0; t2(i)=-1;else t2(i)=1;endend %解调： X2=pskdemod(x2noise,M); numbers,pm2(m)=symerr(b,t2(6:length(x)+5); x2C1=conv(x2noise

6、,C1);for i=1:N+8 if x2C1(i)<0;t3(i)=-1;else t3(i)=1;endend%解调 X2C1=pskdemod(x2C1,M)numbers,pm3(m)=symerr(b,t3(8:length(x)+7); x2C2=conv(x2noise,C2);for i=1:N+14 if x2C2(i)<0;t4(i)=-1;else t4(i)=1;endend %解调 X2C2=pskdemod(x2C2,M); numbers,pm4(m)=symerr(b,t4(14:length(x)+13); x2C=conv(x2noise,C)

7、;for i=1:N+21 if x2C(i)<0;t5(i)=-1;else t5(i)=1;endend %X2C=pskdemod(x2C,M) numbers,pm5(m)=symerr(b,t5(21:length(x)+20);end figure(1) semilogy(SNR,pm1);grid on; title( ' 无串扰的误码率 ' ); axis(1 16 10-5 1);xlabel( 'SNR(dB)' );ylabel( '误码率 ' ); figure(2)semilogy(SNR,pm2);grid on

8、; title( ' 经过衰减信道的误码率 ' );xlabel( 'SNR(dB)' );ylabel( '误码率 ' ); figure(3)semilogy(SNR,pm2, 'r-' ,SNR,pm3, 'ms:' );grid on; title( '5 抽头迫零均衡器 ' );axis（1 16 10-2 1）;xlabel（ 'SNR（dB）' ）;ylabel（ '误码率 ' ）;legend（ ' 多径串扰信道 ','5 抽头

9、'） figure（4）semilogy（SNR,pm2, 'r-' ,SNR,pm4, 'ms:' ）;grid on;title（ '17 抽头迫零均衡器 ' ）;axis（1 16 10-5 1）;xlabel（ 'SNR（dB）' ）;ylabel（ '误码率 ' ）;legend（ '多径串扰信道 ','17 抽头'） figure（5）semilogy（SNR,pm2, 'r-' ,SNR,pm5, 'ms:' ）;grid on;

10、title（ '31 抽头迫零均衡器 ' ）;axis（1 16 10-5 1）;xlabel（ 'SNR（dB）' ）;ylabel（ '误码率 ' ）;legend（ '多径串扰信道 ','31 抽头'） figure（6）on;semilogy（SNR,pm1, 'r-' ,SNR,pm2, 'ms:' ,SNR,pm3, 'cd-' ,SNR,pm4, 'm*-' ,SNR,pm5, 'ko-' ）;grid title（ '