欢迎来到天天文库
浏览记录
ID:52717318
大小:37.00 KB
页数:4页
时间:2020-03-29
《数字通信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
此文档下载收益归作者所有