通信原理matlab課程設(shè)計--2ask、2fsk、2psk、2dpsk調(diào)制解調(diào)matlab仿真-在線瀏覽
2024-07-31 17:27本頁面
【正文】 s([0,N*length(sf2),])。2PSK39。% 數(shù)據(jù)初始化 t=linspace(0,9*pi,450)。 out=1:1:450。 g=1:1:450。 figure(1)。 figure(1)。t39。f139。%畫出 f1信號波形 subplot(2,1,2),plot(l,f2),axis([0 pi ]),xlabel(39。),ylabel(39。)。 end else 6 for j=1:50 out(j+50*(i1))=f2(j)。m=0。 else if out(j+50*(i1))f2(j)==0 m=m+1。 else b(i)=1。 g(j)=b(i)。 subplot(3,1,1),plot(t,f),axis([0 9*pi ]),xlabel(39。),ylabel(39。)。t39。調(diào)制好的波形 39。 subplot(3,1,3),plot(t,g),axis([0 9*pi ]),xlabel(39。),ylabel(39。)。2FSK39。%隨機(jī)碼元個數(shù) N=100。%figure1 畫 四個小圖 a=randint(1,n)%碼元生成 bita=[]。 else bit1=ones(1,N)。%基帶信號 end figure(1)。plot(bita,39。,),title(39。),grid on。 bitRate=1e3。%載頻 1KHZ t=linspace(0,1/bitRate,N)。 c1=sin(2*pi*t*fc)。%載波 for i=1:length(a) if a(i)==1 tz=[tz,c1]。 end end subplot(K,1,2)。LineWidth39。title(39。)。 signal=awgn(tz,20,39。)。plot(signal,39。,),grid on。信號 +噪聲39。%采樣頻率 [b1,a1]=ellip(4,40,[,]*2/Fs)。 sa=filter(b1,a1,signal)。 figure(2)。%figure2 畫（ 3*2)幅圖表示解調(diào)過程 subplot(K1,2,1)。LineWidth39。title(39。) subplot(K1,2,2)。LineWidth39。title(39。) %======相乘器 =============% t=linspace(0,1/bitRate,N)。 c2=sin(2*pi*t*fc*2)。 sib=[]。 sib=[sib,c2]。%乘同頻同相 sia sigb=sb.*sib。plot(siga,39。,),grid on。相乘器 39。 subplot(K1,2,4)。LineWidth39。title(39。)。%抽樣頻率 400HZ [b1,a1]=ellip(4,40,[50]*2/Fs)。%設(shè)計 IIR低通濾波器 sfa=filter(b1,a1,siga)。%信號通過該濾波器，輸出信號 sfb subplot(K1,2,5)。LineWidth39。title(39。)。plot(sfb,39。,),grid on。LPF39。 %===========抽樣判決 ==================% s2a=[]。 LL=fc/bitRate*N。i2=LL。 while (i1=length(sfa)) %判決過程 s2a=[s2a,sfa(i1)=0]。 end while (i2=length(sfb)) s2b=[s2b,sfb(i2)=0]。 end for i=1:n if s2a(i)s2b(i) bit1=zeros(1,N)。 end bitb=[bitb,bit1]。 subplot(K,1,4)。LineWidth39。title(39。)。 9 case 4,scolor=39。 % 2DPSK 調(diào)制與解調(diào) % %Initial_Part % fs = 30000。 Num_Unit = fs * Time_Hold_On。 Low_Level = zeros ( 1, Num_Unit )。 A = 1。 st = zeros ( 1, Num_Unit * Lenth_Of_Sign )。 sign_result = zeros ( 1, Num_Unit * Lenth_Of_Sign )。 % %產(chǎn)生基帶信號 % for I = 1 : Lenth_Of_Sign if Sign_Set(I) == 1 sign_orign( (I1)*Num_Unit + 1 : I*Num_Unit) = High_Level。 end end % %調(diào)制部分 % for I = 1 : Lenth_Of_Sign if Sign_Set(I) == 1 st( (I1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w * t( (I1)*Num_Unit + 1 : I*Num_Unit ) + ( pi / 2 ) )。 end end figure subplot ( 2,1,1 ) plot(t, sign_orign)。 title ( 39。 )。 plot
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1