matlab通信系统仿真-汉明编码,FDMA调制

1、此代码实现了汉明编码，FDMA调制，高斯白噪声信道传输的通信系统matlab仿真。%=所需数据 len = 100;%产生信号长度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;%= 随 机 2 进 制 信 号 产 生msg1 = randint(1,len,2); % Random binary mes

2、sage of 2-bit symbols 信源 msg2 = randint(1,len,2);msg3 = randint(1,len,2);%=汉明编码= %信号 1 编码DS_chips11=encode(msg1,7,4,hamming/fmt); DS_chips12=DS_chips11;%信号2编码DS_chips21=encode(msg2,7,4,hamming/fmt); DS_chips22=DS_chips21;%信号 3 编码DS_chips31=encode(msg3,7,4,hamming/fmt);DS_chips32=DS_chips31;%=信号扩码= %

3、信号 msg1 扩码for mm1=1:175for nn1=1:200u1(mm1-1)*200+nn1)=DS_chips12(mm1);endend%信号 msg2 扩码for mm2=1:len*1.75for nn2=1:200 u2(mm2-1)*200+nn2)=DS_chips22(mm2);endend%信号 msg3 扩码for mm3=1:len*1.75for nn3=1:200 u3(mm3-1)*200+nn3)=DS_chips32(mm3);endend%=调制前信号频谱%=调制前信号频谱figure(1);subplot(3,1,1);stem(t,abs(f

4、ft(ul),.);xlabel(单位:Hz);ylabel (幅度); title(调制前3路信号各自的频谱图)； 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调制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);信 号 调 制 后 3 路 信 号 频 谱s

5、=xl+x2+x3;信 号 调 制 后 3 路 信 号 频 谱%= figure(2);subplot(3,l,l);stem(t,abs(fft(x1),.);xlabel(单位:Hz);ylabel (幅度);title(调制后3路信号各自的频谱图)； subplot(3,1,2);stem(t,abs(fft(x2),.);xlabel(单位:Hz);ylabel(幅度); subplot(3,1,3)；高斯白噪声信号stem(t,abs(fft(x3),.);xlabel(单位:Hz);ylabel(幅度);高斯白噪声信号%=snoisy = awgn(s,snr,measured);

6、% Transmit signal through an AWGN channel.%=带通滤波器%=带通滤波器%=带通滤波器设计%通过带通滤波器Rp1=0.5;Rs1=40;%用切比雪夫设计带通滤波器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);%求实部与虚部的算术平方根%用切比雪夫设计带通

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);%用切比雪夫设计带通滤波器2；%用切比雪夫设计带通滤波器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);%用切比雪夫设计带通滤波器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型设计三个带通滤波器)；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);%滤出三路未解调信号%=信号解调%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);%各个已调信号分别乘以各自%的高频载波信号%= 解 调 后3路 信 号 的频谱figure(5);%图五

10、为解调后3路信号各自的频谱图subplot(3,1,1)；stem(t,abs(fft(y01),.);xlabel(单位:Hz);ylabel(幅度); title(解调后3路信号各自的频谱图)； subplot(3,1,2);stem(t,abs(fft(y02),.);xlabel(单位:Hz);ylabel(幅度); subplot(3,1,3);stem(t,abs(fft(yO3),.);xlabel(单位:Hz);ylabel(幅度);%= 低 通 滤 波 器 设 计Rp2=0.5;%低通滤波器参数选择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);%求实部与虚部的算术平方根db4=20*log10(mag4+eps)/max(mag4);figure(6);%图六为低通滤波器的频率响应plot(w4/pi,db4);axis(0 1 -50 20);xlabel(w/pi);ylabel(20lg|H(ejw)|);title(低通滤波器的频率响应)；=低通滤波= yyl=filter(dl,cl,y01

12、);%通过低通滤波器yy2=filter(dl,cl,y02);%通过低通滤波器yy3=filter(d1,c1,yO3);%通过低通滤波器%= 信 号 抽 样 、 判 决%每200个抽一次，抽第l00个，即中间那个 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、明解码%= 汉%信号 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);%信号 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);%信号1 解码 ddd3=reshape(rt3,7,len*1.75/7);abc3=ddd3;ccc3= decode(abc3,7,4,hamming