资源描述:
《迫零均衡器--代码(毕业设计).doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、程序1clearallN=10000;x=randn(1,N);%产生均值0,方差为1的随机序列fori=1:Nifx(i)>0.5b(i)=1;elseb(i)=-1;endend%调制xx=pskmod(x,M);y=[0.04-0.050.07-0.21-0.50.720.3600.210.030.07];%衰减信道A1=zeros(5,5);form=1:5forn=1:5ifn==mA1(m,n)=0.72;elseifn==m+1A1(m,n)=-0.5;elseifn==m-1A1(m,n)=0.36;elseifn==m+2A1(m,n)=-0.21;els
2、eifn==m-2A1(m,n)=0;elseA1(m,n)=0;endendendendendendendB1=[0;0;1;0;0];C1=pinv(A1)*B1;A2=zeros(17,17);fori=1:17forj=1:17ifj==iA2(i,j)=0.72;elseifj==i-1A2(i,j)=0.36;elseifj==i-2A2(i,j)=0;elseifj==i-3A2(i,j)=0.21;elseifj==i-4A2(i,j)=0.03;elseifj==i-5A2(i,j)=0.07;elseifj==i+1A2(i,j)=-0.5;elseif
3、j==i+2A2(i,j)=-0.21;elseifj==i+3A2(i,j)=0.07;elseifj==i+4A2(i,j)=-0.05;elseifj==i+5A2(i,j)=0.04;elseA2(i,j)=0;endendendendendendendendendendendendendB2=[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);fori=1:31forj=1:31ifj==iA(i,j)=0.72;elseifj==i-1A(i,j)=0.36;elseifj==i-2
4、A(i,j)=0;elseifj==i-3A(i,j)=0.21;elseifj==i-4A(i,j)=0.03;elseifj==i-5A(i,j)=0.07;elseifj==i+1A(i,j)=-0.5;elseifj==i+2A(i,j)=-0.21;elseifj==i+3A(i,j)=0.07;elseifj==i+4A(i,j)=-0.05;elseifj==i+5A(i,j)=0.04;elseA(i,j)=0;endendendendendendendendendendendendendB=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;
5、0;0;0;0;0;0;0;0;0;0;0;0;0;0];C=pinv(A)*B;SNR=1:2:21;form=1:length(SNR)fori=1:Nnoise=randn;x1noise(i)=b(i)*sqrt(2*10^(SNR(m)/10))+noise;ifx1noise(i)<0;t1(i)=-1;elset1(i)=1;endend%解调:X1=pskdemod(x1noise,M);[numbers,pm1(m)]=symerr(b,t1);xxy=conv(b,y);fori=1:N+6noise=randn;x2noise(i)=xxy(i)*sq
6、rt(2*10^(SNR(m)/10))+noise;ifx2noise(i)<0;t2(i)=-1;elset2(i)=1;endend%解调:X2=pskdemod(x2noise,M);[numbers,pm2(m)]=symerr(b,t2(6:length(x)+5));x2C1=conv(x2noise,C1);fori=1:N+8ifx2C1(i)<0;t3(i)=-1;elset3(i)=1;endend%%解调X2C1=pskdemod(x2C1,M)[numbers,pm3(m)]=symerr(b,t3(8:length(x)+7));x2C2=con
7、v(x2noise,C2);fori=1:N+14ifx2C2(i)<0;t4(i)=-1;elset4(i)=1;endend%解调X2C2=pskdemod(x2C2,M);[numbers,pm4(m)]=symerr(b,t4(14:length(x)+13));x2C=conv(x2noise,C);fori=1:N+21ifx2C(i)<0;t5(i)=-1;elset5(i)=1;endend%X2C=pskdemod(x2C,M)[numbers,pm5(m)]=symerr(b,t5(21:len