通信原理matlab課程設計--2ask、2fsk、2psk、2dpsk調(diào)制解調(diào)matlab仿真-文庫吧
2025-05-15 17:27 本頁面
【正文】 it1=ones(1,N)。 end bita=[bita,bit1]。 c=sin(2*pi*t*fc)。 sl=[sl c]。 end figure(1)。 subplot(K,1,1)。plot(bita,39。LineWidth39。,),title(39。 基 帶 信 號 39。),grid on。axis([0,N*length(a),])。 tz=bita*6.*sl。 subplot(K,1,2)。plot(tz,39。LineWidth39。,)。title(39。ASK 調(diào)制后信號 39。)。grid on。 signal=awgn(tz,80,39。measured39。)。 subplot(K,1,3)。plot(signal,39。LineWidth39。,),grid on。title(39。信號 +噪聲39。) Fs=3e3。 [b,a]=ellip(4,40,[,]*2/Fs)。%設計 IIR 帶通濾波器，階數(shù)為4，通帶紋波 ，阻帶衰減 40DB sf=filter(b,a,signal)。%信號通過該濾波器 figure(2)。 K1=4。 subplot(K1,1,1)。plot(sf,39。LineWidth39。,),grid on。title(39。BPF39。) signal2=abs(sf)。 %乘同頻同相 sin subplot(K1,1,2)。plot(signal2,39。LineWidth39。,),grid on。title(39。全波整流器 39。)。 Fs=3e3。%抽樣頻率 400HZ [b,a]=ellip(4,40,[50]*2/Fs)。%設計 IIR 低通濾波器 sf1=filter(b,a,signal2)。%信號通過該濾波器，輸出信號 sf subplot(K1,1,3)。plot(sf1,39。LineWidth39。,),grid on。title(39。LPF39。)。 sf2=[]。 LL=fc/bitRate*N。 5 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]。 end figure(1)。 subplot(K,1,4)。 plot(bitb,39。LineWidth39。,),grid on。title(39。解調(dià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信號波形 subplot(2,1,2),plot(l,f2),axis([0 pi ]),xlabel(39。t39。),ylabel(39。f239。)。%畫出 f2信 號波形 a=[0 1 0 0 0 1 1 0 1] for i=1:9 %2pSK 編碼 if a(i)==0 for k=1:50 %如果二進制原碼為 0 則輸出 f1 波形 out(k+50*(i1))=f1(k)。 end else 6 for j=1:50 out(j+50*(i1))=f2(j)。 %r 如果二進制原碼為 1則輸出 f2 波形 end end end for 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。 end end b for i=1:9 %畫出解碼后的波形 ,包括原碼和解碼出的碼 ,進行對比 for j=1+50*(i1):50*i f(j)=a(i)。 g(j)=b(i)。 end end figure(2)。 subplot(3,1,1),plot(t,f),axis([0 9*pi ]),xlabel(39。t39。),ylabel(39。數(shù)字基帶調(diào)制原碼 3
畢業(yè)設計相關推薦
鄂ICP備17016276號-1