【導(dǎo)讀】matlab/simulink軟件仿真出DPSK信號(hào)。畫(huà)出信號(hào)的波形、眼圖，調(diào)制。誤碼率，使得產(chǎn)生的信號(hào)誤碼率盡可能達(dá)到最小，得出信號(hào)的仿真圖。圖，星座圖，誤碼率與信噪比曲線的形式給出。課程設(shè)計(jì)說(shuō)明書(shū)中應(yīng)附仿真結(jié)果圖及。真模型圖，需提交相應(yīng)模塊的參數(shù)設(shè)置情況。[4]劉學(xué)勇.詳解MATLAB/Simulink通信系統(tǒng)建模與仿真[M].電子工業(yè)出版社.2020年。[6]邵玉斌.MATLAB/Simulink通信系統(tǒng)建模與仿真實(shí)例分析[M].清華大學(xué)出版社.2020. 加上各種噪聲源,用誤碼測(cè)試模塊測(cè)量誤碼率;最后根據(jù)運(yùn)行。結(jié)果和波形來(lái)分析該系統(tǒng)性能。

　　

【正文】 。)。 %%%%%%%%%%%%%%% 179。233。209。249。197。208。189。246。 %%%%%%%%%%%%%%%%%%%% st=zeros(1,i)。 for m=0:i1。 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(4,1,3)。 plot(t,dpsk)。 axis([0,5,1,2])。 title(39。179。233。209。249。197。208。189。246。185。243。178。168。 208。206。39。) %%%%%%%%%%%%%% 194。235。180。177。 228。186。186。 %%%%%%%%%%%%%%%%%% dt=zeros(1,i)。 dt(1)=st(1)。 for n=2:10。 if (st(n)st(n1))=0amp。amp。(st(n)st(n1))1。 dt(n)=0。 else dt(n)=1。 35 35 end end st=t。 for n=1:10 if dt(n)1。 for m=j/i*(n1)+1:j/i*n st(m)=0。 end else for m=j/i*(n1)+1:j/i*n st(m)=1。 end end end subplot(4,1,4)。 plot(t,st)。 title(39。194。235。180。177。 228。186。186。185。243。178。168。 208。206。39。)。 axis([0,5,1,2])。 2 178。238。 214。207。224。201。168。 clear all close all i=10。 %194。235。212。170。181。196。246。202。253。 j=5000。 fc=4。 fm=i/5。 B=2*fm。 t=linspace(0,5,j)。 % 188。171。05199。248。187。228。198。188。189。249。 214。206。170。j 221。 %%%%%%%%%%%%% 178。201。186。249。180。248。208。197。185。197。 %%%%%%%%%%%%%%%% a =round(rand(1,i))。 %178。201。10246。203。230。186。194。235。163。172。187。199。206。170。a st1=t。 36 36 for n=1:10 if a(n)1。 for m=j/i*(n1)+1:j/i*n % j/i206。170。195。190。246。194。235。212。170。181。196。 st1(m)=0。 end else for m=j/i*(n1)+1:j/i*n st1(m)=1。 end end end figure(1)。 subplot(411)。 plot(t,st1)。 title(39。189。248。182。212。194。235。39。)。 axis([0,5,1,2])。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%178。238。 214。177。 228。186。186。 b=zeros(1,i)。 b(1)=a(1)。 for n=2:10 if a(n)=1。 if b(n1)=1 b(n)=0。 else b(n)=1。 end else b(n)=b(n1)。 end end st1=t。 for n=1:10 if b(n)1。 37 37 for m=j/i*(n1)+1:j/i*n st1(m)=0。 end else for m=j/i*(n1)+1:j/i*n st1(m)=1。 end end end subplot(412)。 plot(t,st1)。 title(39。207。224。182。212。194。235。39。)。 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(413) plot(t,st2) title(39。207。224。182。212。194。235。181。196。180。 194。235。39。) axis([0 5 1 2]) %%%%%%%%%%%% 212。216。178。168。 208。197。185。197。 %%%%%%%%%%%%%%%%% s1=sin(2*pi*fc*t)。 subplot(414)。 plot(s1)。 title(39。212。216。178。168。 208。197。185。197。39。)。 %%%%%%%%%%%%%% 181。247。214。198。 %%%%%%%%%%%%%%%%%%%%%% d1=st1.*s1。 d2=st2.*(s1)。%207。224。210。198。180 figure(2)。 38 38 subplot(4,1,1)。 plot(t,d1)。 title(39。st1*s139。)。 subplot(4,1,2)。 plot(t,d2)。 title(39。st2*s239。)。 e_dpsk=d1+d2。 subplot(4,1,3)。 plot(t,e_dpsk)。 title(39。181。247。214。198。185。243。178。168。 208。206。39。)。 noise=rand(1,j)。 dpsk=e_dpsk+*noise。%187。211。200。235。212。235。201。249。 subplot(4,1,4)。 plot(t,dpsk)。 title(39。187。211。212。235。201。249。208。197。185。197。39。)。 %%%%%%%%%%%% 209。211。179。217。181。165。212。170。 %%%%%%%%%%%%%%%%%%%% if dpsk(65)0 dpsk_delay(1:j/i)=dpsk(1:j/i)。 else dpsk_delay(1:j/i)=dpsk(1:j/i) end dpsk_delay(j/i+1:j)=dpsk(1:jj/i)。 %%%%%%%%%%%% 211。235。206。180。209。211。179。217。208。197。185。197。207。224。179。203。 %%%%%%%%%%%%%% dpsk=dpsk.*dpsk_delay。 figure subplot(3,1,1) plot(t,dpsk)。 title(39。209。211。179。217。207。224。179。203。185。243。178。168。 208。206。39。)。 %%%%%%%%%%%%% 181。205。205。168。 194。203。178。168。 %%%%%%%%%%%%%%%%%%%%%%%% [f,af] = T2F(t,dpsk)。%205。168。184。 253。181。205。205。168。 194。203。178。168。 198。247。 [t,dpsk] = lpf(f,af,B)。 subplot(3,1,2)。 plot(t,dpsk)。 title(39。205。168。184。 253。181。205。205。168。 194。203。178。168。 198。247。178。168。 208。206。39。)。 39 39 %%%%%%%%%%%%%%% 179。233。209。249。197。208。189。246。 %%%%%%%%%%%%%%%%%%%% st=zeros(1,i)。 for m=0:i1。 if dpsk(1,m*500+250)0。 st(m+1)=0。 for j=m*500+1:(m+1)*500。 dpsk(1,j)=1。 end else for j=m*500+1:(m+1)*500。 st(m+1)=1。 dpsk(1,j)=0。 end end end subplot(3,1,3)。 plot(t,dpsk)。 axis([0,5,1,2])。 title(39。179。233。209。249。197。208。189。246。185。243。178。168。 208。206。39。)