2PSK调制解调_2psk调制_2psk调制解调_wonderfuljl6_双极性码_调制解调2PSK_

clear close all i=10;%基带信号码元数 j=5000; t=linspace(0,5,j);%0-10之间产生5000个点行矢量，即将[0,5]分成5000份 fc=4;%载波频率 fm=i/5;%码元速率 B=2*fm;%信号带宽 %产生基带信号 a=round(rand(1,i)); %a=inputdlg('请输入十位二级制基带信号','输入信号'); %a=[1,1,1,,0,0,0,1,0,1,0];%随机序列，基带信号 figure(3);stem(a); st1=t; for n=1:10 if a(n)<1; for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end figure(1); subplot(311); plot(t,st1); title('基带信号st1'); axis([0,5,-1,2]); %基带信号求反 %由于PSK中的是双极性信号，因此对上面所求单极性信号取反来与之一起构成双极性码 st2=t; for k=1:j; if st1(k)>=1; st2(k)=0; else st2(k)=1; end end subplot(312); plot(t,st2); title('基带信号反码st2'); axis([0,5,-1,2]); %基带信号变成双极性 st3=st1-st2; subplot(313); plot(t,st3); title('双极性基带信号st3'); axis([0,5,-2,2]); %载波信号 s1=sin(2*pi*fc*t); figure(2); subplot(311); plot(s1); title('载波信号s1'); %调制 e_psk=st3.*s1; subplot(312); plot(t,e_psk); title('调制后波形e-2psk'); %加噪 noise=rand(1,j); psk=e_psk+noise; subplot(313) plot(t,psk) title('加噪后波形') %相干解调 psk=psk.*s1;%与载波相乘 figure(3) subplot(311) plot(t,psk); title('与载波相乘后波形'); [f,af]=T2F(t,psk);%傅里叶变换 %subplot(614); %plot(t,psk); %title('傅里叶变换后波形'); [t,psk] = lpf(f,af,B);%通过低通滤波器 subplot(312); plot(t,psk); title('低通滤波后波形'); %抽样判决 for m=0:i-1 if psk(1,m*500+250)<0; for j=m*500+1:(m+1)*500; psk(1,j)=0; end else for j=m*500+1:(m+1)*500; psk(1,j)=1; end end end subplot(313); plot(t,psk); axis([0,5,-1,2]); title('抽样判决后波形');