原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数字通信qam调制及其matlab仿真程序
数字通信QAM调制MATLAB仿真程序
%数字通信QAM仿真
clear
%产生基带符号流
M=16;? ?? ? %进制 16QAM
N=40;? ?? ? %信源比特流长度,为保证对齐,选择是log2(M)的整倍数
bitstream=zeros(1,N);
symbolstream=zeros(1,N/log2(M));
for n1=1:N
? ? bitstream(n1)=binornd(1,0.5,1,1);
end
%格雷码编码
for n1=1:length(symbolstream)
? ? pack=[bitstream(4*(n1-1)+1) bitstream(4*(n1-1)+2) bitstream(4*(n1-1)+3) bitstream(4*(n1-1)+4)];
? ? if??pack==[0 0 0 0]
? ?? ???symbolstream(n1)=1+j;
? ? elseif??pack==[0 0 1 0]
? ?? ???symbolstream(n1)=2+j;
? ? elseif??pack==[0 0 1 1]
? ?? ???symbolstream(n1)=2+2j;
? ? elseif??pack==[0 0 0 1]
? ?? ???symbolstream(n1)=1+2j;
? ? elseif??pack==[0 1 0 1]
? ?? ???symbolstream(n1)=-1+2j;
? ? elseif??pack==[0 1 0 0]
? ?? ???symbolstream(n1)=-1+j;
? ? elseif??pack==[0 1 1 1]
? ?? ???symbolstream(n1)=-2+2j;
? ? elseif??pack==[0 1 1 0]
? ?? ???symbolstream(n1)=-2+j;
? ? elseif??pack==[1 1 1 0]
? ?? ???symbolstream(n1)=-2-j;
? ? elseif??pack==[1 1 1 1]
? ?? ???symbolstream(n1)=-2-2j;
? ? elseif??pack==[1 1 0 1]
? ?? ???symbolstream(n1)=-1-2j;
? ? elseif??pack==[1 1 0 0]
? ?? ???symbolstream(n1)=-1-j;
? ? elseif??pack==[1 0 0 0]
? ?? ???symbolstream(n1)=1-j;
? ? elseif??pack==[1 0 1 0]
? ?? ???symbolstream(n1)=2-j;
? ? elseif??pack==[1 0 0 1]
? ?? ???symbolstream(n1)=1-2j;
? ? else symbolstream(n1)=2-2j;
? ? end
end
figure(1)
hold on
for n1=-2:2
? ? for n2=-2:2
? ?? ???if ~((n1==0)|(n2==0))
? ?? ?? ?? ?scatter(n1,n2);
? ?? ???end
? ? end
end
axis([-3 3 -3 3]);
hold off
title(16QAM星座图);
xlabel(同相支路)
ylabel(正交支路)
%产生载波
fc=10;? ?? ?? ???%载波频率
fs=100;? ?? ?? ? %采样频率
T=1;? ?? ?? ?? ? %每符号持续时间
N_samples=T*fs;??%每符号内的采样点数
t=0:T/N_samples:(T-T/N_samples);
carrier=exp(j*2*pi*fc*t);
%产生信号脉冲g(t)
%gt=ones(1,length(carrier));
gt=1-cos(2*pi*(1/T)*t);
%生成调制信号S(t)
St_complex=zeros(1,length(carrier)*length(symbolstream));
for n1=1:length(symbolstream)
? ? St_complex((N_samples*(n1-1)+1):(N_samples*(n1-1)+N_samples))=(symbolstream(n1)*carrier).*gt;
end
figure(2)
St_real=real(St_complex);
plot(St_real)
title(QAM仿真波形图 载波10Hz g(t)为升余弦脉冲);
xla
文档评论（0）