发布时间 : 星期六 文章通信原理课程设计 - （基于MATLAB的 - 2PSK - 2DPSK仿真）更新完毕开始阅读
江西农业大学课程设计报告 st2(k)=1; end end;
subplot(412); plot(t,st2);
title('基带信号反码st2'); axis([0,5,-1,2]); st3=st1-st2; subplot(413); plot(t,st3);
title('双极性基带信号st3'); axis([0,5,-2,2]); s1=sin(2*pi*fc*t); subplot(414); plot(s1);
title('载波信号s1'); e_psk=st3.*s1; figure(2); subplot(511); plot(t,e_psk); title('e_2psk'); noise=rand(1,j);
psk=e_psk+noise; %加入噪声 subplot(512); plot(t,psk);
title('加噪后波形');
psk=psk.*s1; %与载波相乘 subplot(513); plot(t,psk);
title('与载波s1相乘后波形');
[f,af] = T2F(t,psk); %通过低通滤波器 [t,psk] = lpf(f,af,B); subplot(514); 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(515); plot(t,psk);
axis([0,5,-1,2]);
title('抽样判决后波形')
2DPSK调制解调程序及注释: clear all close all i=10; j=5000;
fc=4.6; fm=i/5; B=2*fm;
t=linspace(0,5,j); a=round(rand(1,i)); figure(4); 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(321); plot(t,st1); title('绝对码'); axis([0,5,-1,2]);
b=zeros(1,i); b(1)=a(1); for n=2:10
if a(n)>=1;
江西农业大学课程设计报告 %载波频率 %码元速率 %全零矩阵 江西农业大学课程设计报告 if b(n-1)>=1 b(n)=0; else
b(n)=1; end else
b(n)=b(n-1); end end st1=t;
for n=1:10
if b(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
subplot(323); plot(t,st1);
title('相对码st1'); axis([0,5,-1,2]); st2=t; for k=1:j;
if st1(k)>=1; st2(k)=0; else
st2(k)=1; end end;
subplot(324); plot(t,st2);
title('相对码反码st2'); axis([0,5,-1,2]); s1=sin(2*pi*fc*t); subplot(325); plot(s1);
title('载波信号s1'); s2=sin(2*pi*fc*t+pi); subplot(326); plot(s2);
title('载波信号s2'); d1=st1.*s1; d2=st2.*s2; figure(2); subplot(411); plot(t,d1);
title('st1*s1'); subplot(412); plot(t,d2);
title('st2*s2'); e_dpsk=d1+d2; subplot(413); plot(t,e_dpsk);
title('调制后波形'); noise=rand(1,j);
dpsk=e_dpsk+noise; subplot(414); plot(t,dpsk);
title('加噪声后信号');
dpsk=dpsk.*s1; figure(3); subplot(411); plot(t,dpsk);
title('与载波相乘后波形');
[f,af]=T2F(t,dpsk); [t,dpsk]=lpf(f,af,B); subplot(412); plot(t,dpsk);
title('低通滤波后波形');
st=zeros(1,i); for m=0:i-1;
if dpsk(1,m*500+250)<0; st(m+1)=0;
for j=m*500+1:(m+1)*500; dpsk(1,j)=0; end else
for j=m*500+1:(m+1)*500; st(m+1)=1; dpsk(1,j)=1; end end end
subplot(413);
江西农业大学课程设计报告 %加入噪声 %与载波s1相乘 %通过低通滤波器 %全零矩阵 江西农业大学课程设计报告 plot(t,dpsk); axis([0,5,-1,2]);
title('抽样判决后波形')
dt=zeros(1,i); %全零矩阵 dt(1)=st(1); for n=2:10;
if (st(n)-st(n-1))<=0&&(st(n)-st(n-1))>-1; dt(n)=0; else
dt(n)=1; end end st=t;
for n=1:10
if dt(n)<1;
for m=j/i*(n-1)+1:j/i*n st(m)=0; end else
for m=j/i*(n-1)+1:j/i*n st(m)=1; end end end
subplot(414); plot(t,st);
axis([0,5,-1,2]);
title('码反变换后波形');