通信原理matlab課程設(shè)計(jì)--2ask、2fsk、2psk、2dpsk調(diào)制解調(diào)matlab仿真-文庫吧
2024-12-29 19:18 本頁面
【正文】 gnal2=abs(sf)。 %乘同頻同相sinsubplot(K1,1,2)。plot(signal2,39。LineWidth39。,),grid on。title(39。全波整流器39。)。 Fs=3e3。%抽樣頻率400HZ[b,a]=ellip(4,40,[50]*2/Fs)。%設(shè)計(jì)IIR低通濾波器sf1=filter(b,a,signal2)。%信號(hào)通過該濾波器，輸出信號(hào)sfsubplot(K1,1,3)。plot(sf1,39。LineWidth39。,),grid on。title(39。LPF39。)。 sf2=[]。LL=fc/bitRate*N。i=LL/2。bitb=[]。while (i=length(sf1)) %判決 sf2=[sf2,sf1(i)=]。 i=i+LL。 end for i=1:length(sf2) if sf2(i)==0 bit1=zeros(1,N)。 else bit1=ones(1,N)。 end bitb=[bitb,bit1]。endfigure(1)。subplot(K,1,4)。 plot(bitb,39。LineWidth39。,),grid on。title(39。解調(diào)后信號(hào)39。)。 axis([0,N*length(sf2),])。 case 2,scolor=39。2PSK39。 l=linspace(0,pi,50)。% 數(shù)據(jù)初始化t=linspace(0,9*pi,450)。b=1:1:9。out=1:1:450。f=1:1:450。g=1:1:450。w1=2 %正弦波f1的頻率,可以根據(jù)自己想要的頻率在此改寫 %正弦波f2的頻率,可以根據(jù)自己想要的頻率在此改寫f1=sin(w1*l)。figure(1)。f2=sin(w1*l+pi)。figure(1)。subplot(2,1,1),plot(l,f1),axis([0 pi ]),xlabel(39。t39。),ylabel(39。f139。)。%畫出f1信號(hào)波形subplot(2,1,2),plot(l,f2),axis([0 pi ]),xlabel(39。t39。),ylabel(39。f239。)。%畫出f2信號(hào)波形a=[0 1 0 0 0 1 1 0 1]for i=1:9 %2pSK編碼 if a(i)==0 for k=1:50 %如果二進(jìn)制原碼為0則輸出f1波形 out(k+50*(i1))=f1(k)。 end else for j=1:50 out(j+50*(i1))=f2(j)。 %r如果二進(jìn)制原碼為1則輸出f2波形 end endendfor i=1:9 %2PSK解碼 n=0。m=0。 for j=1:50 if out(j+50*(i1))f1(j)==0 n=n+1。 else if out(j+50*(i1))f2(j)==0 m=m+1。 end end end if nm b(i)=0。 else b(i)=1。 endend bfor i=1:9 %畫出解碼后的波形 ,包括原碼和解碼出的碼,進(jìn)行對(duì)比 for j=1+50*(i1):50*i f(j)=a(i)。 g(j)=b(i)。 endendfigure(2)。subplot(3,1,1),plot(t,f),axis([0 9*pi ]),xlabel(39。t39。),ylabel(39。數(shù)字基帶調(diào)制原碼39。)。subplot(3,1,2),plot(t,out),axis([0 9*pi ]),xlabel(39。t39。),ylabel(39。調(diào)制好的波形39。)。subplot(3,1,3),plot(t,g),axis([0 9*pi ]),xlabel(39。t39。),ylabel(39。解碼得到的碼39。)。[Num,Rat]=biterr(f,g) case 3,scolor=39。2FSK39。%==生成隨機(jī)碼元、基帶信號(hào)、調(diào)制================%n=8。%隨機(jī)碼元個(gè)數(shù)N=100。%模擬一個(gè)碼元
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1