BPSK系统仿真—MATLAB

1、程序代码：主程序：%基带信号s(t)clearallcloseallclcnum=10;%码元个数tnum=200;%码元长度N=num*tnum;a=randint(1,num,2);%产生1行num列的矩阵，矩阵内0和1随机出现fc=0.1;%载波频率t=0:0.05:9.99;s=;c=;fori=1:numif(a(i)=0)A=zeros(1,tnum);elseA=ones(1,tnum);ends=sA;cs=sin(2*pi*fc*t);c=ccs;%产生载波信号endfigure(1);subplot(4,1,1);plot(s);gridon;xlabel(基带信号s(t)

2、;axis(0N-22);%BPSK已调信号。采用模拟调制方法s_NRZ=;fori=1:numif(a(i)=0)A=ones(1,tnum);elseA=-1*ones(1,tnum);ends_NRZ=s_NRZ,A;ende=s_NRZ.*c;figure(1);subplot(4,1,2);plot(e);gridon;xlabel(BPSK调制信号);axis(0N-22);%加高斯噪声am=0.7;%输入信号经信道后振幅由1衰减为0.7SNR=5;snr=10人(SNR/10);N0=(am*am)/2/snr;%计算噪声功率N0_db=10*log10(N0);%将噪声功率转换

3、为dBWni=wgn(1,N,N0_db);%产生1行N列的高斯噪声yi=e+ni;%BSK已调信号中加入白噪声，输入信噪比为SNRfigure(1);subplot(4,1,3);plot(yi);gridon;xlabel(加入高斯白噪声的已调信号yi(t);%带通滤波器b1,a1=BUTTER(3,2*pi*0.0001,2*pi*0.01);%计算带通滤波器的H(z)系数y=filter(b1,a1,yi);%对信号yi进行滤波，得到信号yfigure(1);subplot(4,1,4);plot(y);gridon;xlabel(经带通滤波器后信号)；%与恢复载波相乘x1=2*c.*

4、y;figure(2);subplot(4,1,1);plot(x1);gridonxlabel(与恢复载波相乘后的信号x1(t);%低通滤波器b2,a2=butter(2,0.005);%计算H系数，频率为(1/200)x=filter(b2,a2,x1);%对信号x1滤波，得到信号xfigure(2);subplot(4,1,2);plot(x);gridonaxis(0N-22);xlabel(经低通滤波器后信号波形)x=fun_panjue(x);%调用函数，进行抽样判决figure(2);subplot(4,1,3);plot(x);gridon;xlabel(加噪后解调信号x(t)

5、;axis(0N-22);x=fun_yanc(x);%调用函数，进行消除延迟figure(2);subplot(4,1,4);plot(x);gridon;xlabel(加噪后去掉延迟的解调信号x(t);axis(0N-22);%误码率计算Err1=length(find(x=s)%计算解调信号中错误码元个数Pe_test1=Err1/N%计算实际误码率Pe1=(1/2)*erfc(sqrt(snr)%计算系统理论误码率%理论误码率曲线Pe=;forSNR=1:10am=0.7;%输入信号经信道后振幅由1衰减为0.7E=am*am/2;snr=10人(SNR/10);N0=(am*am)/(

6、2*snr);no=N0/(2*200);%计算噪声功率N0_db=10*log10(N0);%将噪声功率转换为dBWni=wgn(1,N,N0_db);%产生1行N列的高斯噪声yi=e+ni;%BSK已调信号中加入白噪声，输入信噪比为SNRy=filter(b1,a1,yi);%对yi进行滤波(带通滤波器)，得到信号yx1=2*c.*y;%与恢复载波相乘xx=filter(b2,a2,x1);%经低通滤波器滤波xx=fun_panjue(xx);%抽样判决xx=fun_yanc(xx);%消除延迟snr=10A(SNR/10);Pe=Pe,(1/2)*erfc(sqrt(snr);%计算理论

7、误码率endPe;figure;SNR=1:10;semilogy(SNR,Pe,b-);holdon%以log10(Pe)为纵坐标画图gridon子程序一：抽样判决(由于存在“倒pi”现象，故以0为基准进行判决，大于0,判为0;小于0,判为1)functionw=fun_panjue(w)N=length(w);ifw(100)0w(1:100)=0;elsew(1:100)=1;endfori=101:Nifw(i)0;w(i)=0;elsew(i)=1;endend子函数二：消除延迟(由于抽样判决后会产生延迟，故认为对解调信号进行时移，以减小甚至消除延迟，便于与基带信号进行比较，计算实际

8、误码率)functionm=fun_yanc(m)N=length(m);leng=0;ifm(1)=0fori=1:Nifm(i)=1leng=i;break;endendelsefori=1:Nifm(i)=0leng=i;break;endendendleng1=leng-(floor(leng/200)*200;fori=1:(N-leng1)m(i)=m(i+leng1);endfori=(N-leng1):Nm(i)=m(N-200+10);end运行结果：输入信噪比为5dB时，计算的错误码元，实际误码率，理论误码率Err1=15Pe_test1=0.0075Pe1=0.0060各波形图如下所示FileEditViewInsertToolsDesktopWindow旦即口Qa|fl|B|s5厂|ii|it0_giiiiiiiiiTOC o 1-5 h z0200400600800100012001400160018002000与恢复载減相乘后的信号id(t)_2iiiiiiiii0200400600800100012001400150018002000经低通滤波器后信号波形_2iiiiiiii