数字通信 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.for n1=1：N9. bitstream（n1)=binornd（1,0。5，1，1）；10.end11.%格雷码编码12.for n1=1：length（symbolstream）13. pack=［bitstream（4*（n1—1)+1） bitstream(4*（n1—1)+2) bitstream(4*(n1—1）+3) bitstream(4*（n1-1）+4）]；14. if pack==[0 0 0 0］15. symbolstream(n1）=1+j;16. elseif pack==［0 0 1 0]17. symbolstream(n1）=2+j;18. elseif pack==［0 0 1 1]19. symbolstream(n1)=2+2j;20. elseif pack==［0 0 0 1］21. symbolstream（n1）=1+2j;22. elseif pack==[0 1 0 1]23. symbolstream(n1)=-1+2j；24. elseif pack==［0 1 0 0］25. symbolstream(n1）=—1+j；26. elseif pack==[0 1 1 1］27. symbolstream(n1)=-2+2j；28. elseif pack==［0 1 1 0]29. symbolstream(n1）=-2+j;30. elseif pack==［1 1 1 0]31. symbolstream（n1）=—2-j；32. elseif pack==［1 1 1 1］33. symbolstream(n1)=-2-2j；34. elseif pack==[1 1 0 1］35. symbolstream（n1）=-1-2j；36. elseif pack==［1 1 0 0]37. symbolstream(n1)=-1—j；38. elseif pack==[1 0 0 0］39. symbolstream(n1)=1—j;40. elseif pack==[1 0 1 0］41. symbolstream（n1）=2—j;42. elseif pack==[1 0 0 1]43. symbolstream(n1）=1-2j;44. else symbolstream(n1）=2-2j；45. end46.end47.figure（1）48.hold on49.for n1=-2：250. for n2=—2：251. if ~（（n1==0)｜(n2==0）） 52. scatter（n1,n2);53. end54. end55.end56.axis(［—3 3 —3 3］）；57.hold off58.title('16QAM 星座图')；59.xlabel（'同相支路’）60.ylabel（'正交支路')61.%产生载波62.fc=10; %载波频率63.fs=100； %采样频率64.T=1; ％每符号持续时间65.N_samples=T＊fs； ％每符号内的采样点数66.t=0:T/N_samples：(T—T/N_samples）；67.carrier=exp(j*2＊pi*fc*t）；68.%产生信号脉冲 g(t）69.％gt=ones(1,length（carrier）);70.gt=1—cos（2*pi*(1/T)＊t）;71.72.％生成调制信号 S(t）73.St_complex=zeros（1...
