汉明编码+matlab,matlab通信系统仿真汉明编码.doc

matlab通信系统仿真汉明编码.doc

此代码实现了汉明编码，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 message 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 ; %===========================信号扩码================================= %信号 msg1 扩码 for mm1=1:175 for nn1=1:200u1((mm1-1)*200+nn1)=DS_chips12(mm1);end end %信号 msg2 扩码 for mm2=1:len*1.75for nn2=1:200u2((mm2-1)*200+nn2)=DS_chips22(mm2);end end %信号 msg3 扩码 for mm3=1:len*1.75for nn3=1:200u3((mm3-1)*200+nn3)=DS_chips32(mm3);end end %===========================调制前信号频谱 ================================= figure(1); subplot(3,1,1); stem(t,abs(fft(u1)), . );xlabel( 单位:Hz );ylabel( 幅度 ); title( 调制前 3 路信号各自的频谱图 ); subplot(3,1,2); stem(t,abs(fft(u2)), . );xlabel( 单位:Hz );ylabel( 幅度 ); subplot(3,1,3); stem(t,abs(fft(u3)), . );xlabel( 单位:Hz );ylabel( 幅度 ); %================================FDMA 调制 ================================= x1=4*u1.*cos(2*pi*f1*t/fs); %频谱搬移 x2=4*u2.*cos(2*pi*f2*t/fs); x3=4*u3.*cos(2*pi*f3*t/fs); s=x1+x2+x3; % 复用信号频谱为各信号频谱的叠加 %===========================信号调制后 3 路信号频谱 ================================= figure(2); subplot(3,1,1); 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 );% Transmit signal through an AWGN channel. %=============================带通滤波器设计 ================================= %通过带通滤波器 Rp1=0.5;Rs1=40; % 用切比雪夫设计带通滤波器 1； Wp1=[4000 8000]/22050; %fs/2=22050 Ws1=[3800 8500]/22050; [n1,Wn1]=cheb2ord(Wp1,Ws1,Rp1,Rs1); [b1,a1]=cheby2(n1,Rs1,Wn1); [h1,w1]=freqz(b1,a1); mag1=abs(h1); % 求实部与虚部的算术平方根 db1=20*log10((mag1+eps)/max(mag1)); %用切比雪夫设计带通滤波器 2； Wp2=[9000 13000]/22050; %用切比雪夫设计带通滤波器 2； 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)); %用切比雪夫设计带通滤波器 3； Wp3=[14500 18500]/22050; %用切比雪夫设计带通滤波器 3； Ws3=[14000 19000]/22050; [n3,Wn3]=ch