通信原理matlab課程設(shè)計(jì)--2ask、2fsk、2psk、2dpsk調(diào)制解調(diào)matlab仿真(編輯修改稿)
2024-07-10 17:27 本頁面
【文章內(nèi)容簡介】 9。)。 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ī)碼元、基帶信號、調(diào)制 ================% n=8。%隨機(jī)碼元個數(shù) N=100。%模擬一個碼元的點(diǎn)數(shù) K=4。%figure1 畫 四個小圖 a=randint(1,n)%碼元生成 bita=[]。%定義空數(shù)組，存放基帶信號 for i=1:length(a) 7 if a(i)==0 bit1=zeros(1,N)。 else bit1=ones(1,N)。 end bita=[bita,bit1]。%基帶信號 end figure(1)。 subplot(K,1,1)。plot(bita,39。LineWidth39。,),title(39。 基 帶 信 號 39。),grid on。axis([0,N*length(a),])。 bitRate=1e3。%每一個碼元中采樣點(diǎn)的間隔寬度 fc=1e3。%載頻 1KHZ t=linspace(0,1/bitRate,N)。 tz=[]。 c1=sin(2*pi*t*fc)。%載波 c2=sin(2*pi*t*fc*2)。%載波 for i=1:length(a) if a(i)==1 tz=[tz,c1]。 else tz=[tz,c2]。 end end subplot(K,1,2)。plot(tz,39。LineWidth39。,)。title(39。2FSK 已調(diào) 信號 39。)。grid on。 signal=awgn(tz,20,39。measured39。)。%加噪 subplot(K,1,3)。plot(signal,39。LineWidth39。,),grid on。title(39。信號 +噪聲39。) %===========解調(diào) ============% Fs=5e3。%采樣頻率 [b1,a1]=ellip(4,40,[,]*2/Fs)。%設(shè)計(jì) IIR 帶通濾波器，階數(shù)為 4，通帶紋波 ，阻帶衰減 40DB [b2,a2]=ellip(4,40,[,]*2/Fs)。 sa=filter(b1,a1,signal)。%信號通過該濾波器 sb=filter(b2,a2,signal)。 figure(2)。 K1=3。%figure2 畫（ 3*2)幅圖表示解調(diào)過程 subplot(K1,2,1)。plot(sa,39。LineWidth39。,),grid on。title(39。BPF39。) subplot(K1,2,2)。plot(sb,39。LineWidth39。,),grid on。title(39。BPF39。) %======相乘器 =============% t=linspace(0,1/bitRate,N)。 c1=sin(2*pi*t*fc)。 c2=sin(2*pi*t*fc*2)。 8 sia=[]。 sib=[]。 for i=1:n sia=[sia,c1]。 sib=[sib,c2]。 end siga=sa.*sia。%乘同頻同相 sia sigb=sb.*sib。%乘同頻同相 sib subplot(K1,2,3)。plot(siga,39。LineWidth39。,),grid on。title(39。相乘器 39。)。 subplot(K1,2,4)。plot(sigb,39。LineWidth39。,),grid on。title(39。相乘器 39。)。 %========LPF================% Fs=5e3。%抽樣頻率 400HZ [b1,a1]=ellip(4,40,[50]*2/Fs)。%設(shè)計(jì) IIR低通濾波器 [b2,a2]=ellip(4,40,[50]*2/Fs)。%設(shè)計(jì) IIR低通濾波器 sfa=filter(b1,a1,siga)。%信號通過該濾波器，輸出信號 sfa sfb=filter(b2,a2,sigb)。%信號通過該濾波器，輸出信號 sfb subplot(K1,2,5)。plot(sfa,39。LineWidth39。,),grid on。title(39。LPF39。)。 subplot(K1,2,6)。plot(sfb,39。LineWidth39。,),grid on。title(39。LPF39。)。 %===========抽樣判決 ==================% s2a=[]。s2b=[]。 LL=fc/bitRate*N。 i1=LL/2。i2=LL。 bitb=[]。 while (i1=length(sfa)) %判決過程 s2a=[s2a,sfa(i1)=0]。 i1=i1+LL。 end while (i2=length(sfb)) s2b=[s2b,sfb(i2)=0]。 i2=i2+LL。 end for i=1:n if s2a(i)s2b(i) bit1=zeros(1,N)。 else
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號-1