数字通信QAM调制及其MATLAB仿真程序.doc

《数字通信QAM调制及其MATLAB仿真程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数字通信QAM调制及其MATLAB仿真程序1.%数字通信QAM仿真2.clear3.%产生基带符号流4.M=16;    %进制16QAM5.N=40;    %信源比特流长度,为保证对齐,选择是log2(M)的整倍数6.bitstream=zeros(1,N);7.symbolstream=zeros(1,N/log2(M));8.forn1=1:N9.  bitstream(n1)=binornd(1,0.5,1,1);10.end11.%格雷码编码12.forn1=1:length(symbolstream)13.  pack=[bitstream(4*(n1-

2、1)+1)bitstream(4*(n1-1)+2)bitstream(4*(n1-1)+3)bitstream(4*(n1-1)+4)];14.  if  pack==[0000]15.      symbolstream(n1)=1+j;16.  elseif  pack==[0010]17.      symbolstream(n1)=2+j;18.  elseif  pack==[0011]19.      symbolstream(n1)=2+2j;20.  elseif  pack==[0001]21.      symbolstream(n1)=1+2j

3、;22.  elseif  pack==[0101]23.      symbolstream(n1)=-1+2j;24.  elseif  pack==[0100]25.      symbolstream(n1)=-1+j;26.  elseif  pack==[0111]27.      symbolstream(n1)=-2+2j;28.  elseif  pack==[0110]29.      symbolstream(n1)=-2+j;30.  elseif  pack==[1110]1.      symbolstream(n1)=-2-j;2.  e

4、lseif  pack==[1111]3.      symbolstream(n1)=-2-2j;4.  elseif  pack==[1101]5.      symbolstream(n1)=-1-2j;6.  elseif  pack==[1100]7.      symbolstream(n1)=-1-j;8.  elseif  pack==[1000]9.      symbolstream(n1)=1-j;10.  elseif  pack==[1010]11.      symbolstream(n1)=2-j;12.  elseif  pack==[

5、1001]13.      symbolstream(n1)=1-2j;14.  elsesymbolstream(n1)=2-2j;15.  end16.end17.figure(1)18.holdon19.forn1=-2:220.  forn2=-2:221.      if~((n1==0)

6、(n2==0))22.        scatter(n1,n2);23.      end24.  end25.end26.axis([-33-33]);27.holdoff28.title('16QAM星座图');29.xlabel('同相支路')30.ylabel(

7、'正交支路')31.%产生载波32.fc=10;        %载波频率1.fs=100;      %采样频率2.T=1;        %每符号持续时间3.N_samples=T*fs;  %每符号内的采样点数4.t=0:T/N_samples:(T-T/N_samples);5.carrier=exp(j*2*pi*fc*t);6.%产生信号脉冲g(t)7.%gt=ones(1,length(carrier));8.gt=1-cos(2*pi*(1/T)*t);9.10.%生成调制信号S(t)11.St_complex=zeros(1,length(carri

8、er)*length(symbolstream));12.forn1=1:length(symbolstream)13.  St_complex((N_samples*(n1-1)+1):(N_samples*(n1-1)+N_samples))=(symbolstream(n1)*carrier).*gt;14.end15.figure(2)16.St_real=real(St_complex);17.plot(St_real)18.title('QAM仿真波形图载波10Hzg(t)为升余弦脉冲');19.xlabel('采样点')20.ylabe