matlab通信系統(tǒng)仿真-漢明編碼,FDMA調(diào)制

1、此代碼實(shí)現(xiàn)了漢明編碼，F(xiàn)DMA調(diào)制，高斯白噪聲信道傳輸?shù)耐ㄐ畔到y(tǒng)matlab仿真。%=所需數(shù)據(jù) len = 100;%產(chǎn)生信號(hào)長度fs=44100;%抽樣頻率t=0:len*200*1.75-1; %t 值范圍 mm1=1:len;mm2=1:len;mm3=1:len;nn1=1:200;nn2=1:200;nn3=1:200;uu=1:300;snr=10;%高斯信道信噪比fazhi=0.5;%判決閥值f1=5200;%載波頻率f2=10200;f3=16200;%= 隨 機(jī) 2 進(jìn) 制 信 號(hào) 產(chǎn) 生msg1 = randint(1,len,2); % Random binary mes

2、sage of 2-bit symbols 信源 msg2 = randint(1,len,2);msg3 = randint(1,len,2);%=漢明編碼= %信號(hào) 1 編碼DS_chips11=encode(msg1,7,4,hamming/fmt); DS_chips12=DS_chips11;%信號(hào)2編碼DS_chips21=encode(msg2,7,4,hamming/fmt); DS_chips22=DS_chips21;%信號(hào) 3 編碼DS_chips31=encode(msg3,7,4,hamming/fmt);DS_chips32=DS_chips31;%=信號(hào)擴(kuò)碼= %

3、信號(hào) msg1 擴(kuò)碼for mm1=1:175for nn1=1:200u1(mm1-1)*200+nn1)=DS_chips12(mm1);endend%信號(hào) msg2 擴(kuò)碼for mm2=1:len*1.75for nn2=1:200 u2(mm2-1)*200+nn2)=DS_chips22(mm2);endend%信號(hào) msg3 擴(kuò)碼for mm3=1:len*1.75for nn3=1:200 u3(mm3-1)*200+nn3)=DS_chips32(mm3);endend%=調(diào)制前信號(hào)頻譜%=調(diào)制前信號(hào)頻譜figure(1);subplot(3,1,1);stem(t,abs(f

4、ft(ul),.);xlabel(單位:Hz);ylabel (幅度); title(調(diào)制前3路信號(hào)各自的頻譜圖)； subplot(3,l,2);stem(t,abs(fft(u2),.);xlabel(單位:Hz);ylabel(幅度); subplot(3,l,3);stem(t,abs(fft(u3),.);xlabel(單位:Hz);ylabel(幅度);%=FDMA調(diào)制xl=4*ul.*cos(2*pi*fl*t/fs);%頻譜搬移x2=4*u2.*cos(2*pi*f2*t/fs);x3=4*u3.*cos(2*pi*f3*t/fs);信 號(hào) 調(diào) 制 后 3 路 信 號(hào) 頻 譜s

5、=xl+x2+x3;信 號(hào) 調(diào) 制 后 3 路 信 號(hào) 頻 譜%= figure(2);subplot(3,l,l);stem(t,abs(fft(x1),.);xlabel(單位:Hz);ylabel (幅度);title(調(diào)制后3路信號(hào)各自的頻譜圖)； subplot(3,1,2);stem(t,abs(fft(x2),.);xlabel(單位:Hz);ylabel(幅度); subplot(3,1,3)；高斯白噪聲信號(hào)stem(t,abs(fft(x3),.);xlabel(單位:Hz);ylabel(幅度);高斯白噪聲信號(hào)%=snoisy = awgn(s,snr,measured);

6、% Transmit signal through an AWGN channel.%=帶通濾波器%=帶通濾波器%=帶通濾波器設(shè)計(jì)%通過帶通濾波器Rp1=0.5;Rs1=40;%用切比雪夫設(shè)計(jì)帶通濾波器1；Wp1=4000 8000/22050;Ws1=3800 8500/22050;n1,Wn1=cheb2ord(Wp1,Ws1,Rp1,Rs1);b1,a1=cheby2(n1,Rs1,Wn1);h1,w1=freqz(b1,a1);%fs/2=22050mag1=abs(h1);db1=20*log10(mag1+eps)/max(mag1);%求實(shí)部與虛部的算術(shù)平方根%用切比雪夫設(shè)計(jì)帶通

7、濾波器2；Wp2=9000 13000/22050;Ws2=8000 14000/22050; n2,Wn2=cheb2ord(Wp2,Ws2,Rp1,Rs1);b2,a2=cheby2(n2,Rs1,Wn2);h2,w2=freqz(b2,a2);mag2=abs(h2);db2=20*log10(mag2+eps)/max(mag2);%用切比雪夫設(shè)計(jì)帶通濾波器2；%用切比雪夫設(shè)計(jì)帶通濾波器3；Wp3=14500 18500/22050;Ws3=14000 19000/22050; n3,Wn3=cheb2ord(Wp3,Ws3,Rp1,Rs1);b3,a3=cheby2(n3,Rs1,W

8、n3); h3,w3=freqz(b3,a3);mag3=abs(h3); db3=20*log10(mag3+eps)/max(mag3);%用切比雪夫設(shè)計(jì)帶通濾波器3；figure(4);subplot(3,1,1);plot(w1/pi,db1);axis(0 1 -50 20);xlabel(w/pi);ylabel(20lg|H(ejw)|); title(用切比雪夫2型設(shè)計(jì)三個(gè)帶通濾波器)；subplot(3,1,2);plot(w2/pi,db2)；axis(0 1 -50 20)；xlabel(w/pi)；ylabel(20lg|H(ejw)|)；subplot(3,1,3)；

9、plot(w3/pi,db3)；axis(0 1 -50 20)；xlabel(w/pi)；ylabel(20lg|H(ejw)|)；%=帶通濾波y1=filter(b1,a1,snoisy);y2=filter(b2,a2,snoisy);y3=filter(b3,a3,snoisy);%濾出三路未解調(diào)信號(hào)%=信號(hào)解調(diào)%fs=44100y01=y1.*cos(2*pi*f1*t/fs);y02=y2.*cos(2*pi*f2*t/fs);y03=y3.*cos(2*pi*f3*t/fs);%各個(gè)已調(diào)信號(hào)分別乘以各自%的高頻載波信號(hào)%= 解 調(diào) 后3路 信 號(hào) 的頻譜figure(5);%圖五

10、為解調(diào)后3路信號(hào)各自的頻譜圖subplot(3,1,1)；stem(t,abs(fft(y01),.);xlabel(單位:Hz);ylabel(幅度); title(解調(diào)后3路信號(hào)各自的頻譜圖)； subplot(3,1,2);stem(t,abs(fft(y02),.);xlabel(單位:Hz);ylabel(幅度); subplot(3,1,3);stem(t,abs(fft(yO3),.);xlabel(單位:Hz);ylabel(幅度);%= 低 通 濾 波 器 設(shè) 計(jì)Rp2=0.5;%低通濾波器參數(shù)選擇Rs2=40;Wp4=3400/(22050);Ws4=4000/(22050

11、);n4,Wn4=cheb2ord(Wp4,Ws4,Rp2,Rs2); %采用切比雪夫 2 型(cheby2)帶通濾波器 d1,c1=cheby2(n4,Rs2,Wn4);h4,w4=freqz(d1,c1);mag4=abs(h4);%求實(shí)部與虛部的算術(shù)平方根db4=20*log10(mag4+eps)/max(mag4);figure(6);%圖六為低通濾波器的頻率響應(yīng)plot(w4/pi,db4);axis(0 1 -50 20);xlabel(w/pi);ylabel(20lg|H(ejw)|);title(低通濾波器的頻率響應(yīng))；=低通濾波= yyl=filter(dl,cl,y01

12、);%通過低通濾波器yy2=filter(dl,cl,y02);%通過低通濾波器yy3=filter(d1,c1,yO3);%通過低通濾波器%= 信 號(hào) 抽 樣 、 判 決%每200個(gè)抽一次，抽第l00個(gè)，即中間那個(gè) for uu=l:len*l.75if yyl(uu-l)*200+l00)=fazhi rtl(uu)=l;else rtl(uu)=0;endif yy2(uu-l)*200+l00)=fazhi rt2(uu)=l;else rt2(uu)=0;endif yy3(uu-l)*200+l00)=fazhi rt3(uu)=l;else rt3(uu)=0;end明解碼end

13、明解碼%= 漢%信號(hào) l 解碼dddl=reshape(rtl,7,len*l.75/7);abcl=dddl;cccl= decode(abcl,7,4,hamming/fmt);cll=cccl;zhl=reshape(cll,l,len);%信號(hào) l 解碼ddd2=reshape(rt2,7,len*1.75/7);abc2=ddd2;ccc2= decode(abc2,7,4,hamming/fmt);cl2=ccc2;zh2=reshape(cl2,1,len);%信號(hào)1 解碼 ddd3=reshape(rt3,7,len*1.75/7);abc3=ddd3;ccc3= decode(abc3,7,4,hamming