1、程序1clear allN=10000;x=randn(1,N);%產(chǎn)生均值0,方差為1的隨機序列for i=1:N if x(i)0.5 b(i)=1; else b(i)=-1; endend %調(diào)制xx=pskmod(x,M); y=0.04 -0.05 0.07 -0.21 -0.5 0.72 0.36 0 0.21 0.03 0.07; %衰減信道A1=zeros(5,5);for m=1:5 for n=1:5 if n=m A1(m,n)=0.72; else if n=m+1 A1(m,n)=-0.5; else if n=m-1 A1(m,n)=0.36; else if n

2、=m+2 A1(m,n)=-0.21; else if n=m-2 A1(m,n)=0; else A1(m,n)=0; end end end end end endendB1=0;0;1;0;0;C1=pinv(A1)*B1;A2=zeros(17,17);for i=1:17 for j=1:17 if j=i A2(i,j)=0.72; else if j=i-1 A2(i,j)=0.36; else if j=i-2 A2(i,j)=0; else if j=i-3 A2(i,j)=0.21; else if j=i-4 A2(i,j)=0.03; else if j=i-5 A2(

3、i,j)=0.07; else if j=i+1 A2(i,j)=-0.5; else if j=i+2 A2(i,j)=-0.21; else if j=i+3 A2(i,j)=0.07; else if j=i+4 A2(i,j)=-0.05; else if j=i+5 A2(i,j)=0.04; else A2(i,j)=0; end end endendendendendendendendendendendB2=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:31 for j=1:31 i

4、f j=i A(i,j)=0.72; else if j=i-1 A(i,j)=0.36; else if j=i-2 A(i,j)=0; else if j=i-3 A(i,j)=0.21; else if j=i-4 A(i,j)=0.03; else if j=i-5 A(i,j)=0.07; else if j=i+1 A(i,j)=-0.5; else if j=i+2 A(i,j)=-0.21; else if j=i+3 A(i,j)=0.07; else if j=i+4 A(i,j)=-0.05; else if j=i+5 A(i,j)=0.04; else A(i,j)=

5、0; end end endendendendendendendendendendend B=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)=-1; else t1(i)=1; end end %解調(diào)：X1=pskdemod(x1noise,M); numbers

6、,pm1(m)=symerr(b,t1); xxy=conv(b,y); for i=1:N+6 noise=randn; x2noise(i)=xxy(i)*sqrt(2*10(SNR(m)/10)+noise; if x2noise(i)0; t2(i)=-1; else t2(i)=1; end end %解調(diào)：X2=pskdemod(x2noise,M); numbers,pm2(m)=symerr(b,t2(6:length(x)+5); x2C1=conv(x2noise,C1); for i=1:N+8 if x2C1(i)0; t3(i)=-1; else t3(i)=1; e

7、nd end %解調(diào) 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; end end %解調(diào)X2C2=pskdemod(x2C2,M); numbers,pm4(m)=symerr(b,t4(14:length(x)+13); x2C=conv(x2noise,C); for i=1:N+21 if x2C(i)0; t5(i)=-1; else t5(i)=1; en

8、d end %X2C=pskdemod(x2C,M) numbers,pm5(m)=symerr(b,t5(21:length(x)+20);endfigure(1)semilogy(SNR,pm1);grid on;title(無串擾的誤碼率);axis(1 16 10-5 1);xlabel(SNR（dB）);ylabel(誤碼率);figure(2)semilogy(SNR,pm2);grid on;title(經(jīng)過衰減信道的誤碼率);xlabel(SNR（dB）);ylabel(誤碼率);figure(3)semilogy(SNR,pm2,r-,SNR,pm3,ms:);grid on

9、;title(5抽頭迫零均衡器);axis(1 16 10-2 1);xlabel(SNR（dB）);ylabel(誤碼率);legend(多徑串擾信道,5抽頭)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;title(31抽頭迫零均衡器);axis(1 16 10-5 1);xlabel(SNR（dB）);ylabel(誤碼率);legend(多徑串擾信道,31抽頭)figure(6)semilogy