基于MATLAB的dpsk系统仿真代码.docVIP

下载本文档

28
0
约3.07千字
约 5页
2017-08-22 发布于重庆
举报
版权申诉

基于MATLAB的dpsk系统仿真代码.doc

1、本文档共5页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

基于MATLAB的dpsk系统仿真代码.doc

差分编码 function Y = difference(X) %UNTITLED Summary of this function goes here % Detailed explanation goes here a=length(X); Y=zeros(1,a); b=0.01:0.01:a; for(i=1:a) is((i-1)*100+1:i*100)=X(i); end figure plot(b,is); axis([0,a+1,-0.1,1.1]); title(原始序列) grid if (X(1)==1) Y(1)=1; else Y(1)=0; end for i=2:a if(X(i)-X(i-1)~=0) Y(i)=1; else Y(i)=0; end end Y; t=0.01:0.01:a; for(i=1:a) st((i-1)*100+1:i*100)=Y(i); end figure plot(t,st); axis([0,a+1,-0.1,1.1]); title(差分编码) grid end PSK调制 function z=dpsk(X) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here fs=200; w=2; n=length(X); t=1:100; nn=1:100; n1=0.01:0.01:n; %os=zeros(1,100); for(i=1:n) is((i-1)*100+1:i*100)=X(i); st((i-1)*100+1:i*100)=0; end for(nn=1:100) oss(nn)=sin(2*pi*w*(t(nn)/100)); osc(nn)=sin(2*pi*w*(t(nn)/100)+pi); end %figure %subplot(2,1,1) %plot(t/100,oss); %axis([0,n+1,-1.1,1.1]); %title(信号波形1) %grid %subplot(2,1,2) %plot(t/100,osc); %axis([0,n+1,-1.1,1.1]); %title(信号波形2) %grid for(i=1:n) if(X(i)==1) for(nn=1:100) st(100*(i-1)+nn)=oss(nn); end else for(nn=1:100) st(100*(i-1)+nn)=osc(nn); end end end st; figure subplot(2,1,1) plot(n1,is); axis([0,n+1,-0.1,1.1]); title(原始信号) grid subplot(2,1,2) plot(n1,st); axis([0,n+1,-1.1,1.1]); title(调制信号) grid %%%%%%%%%%%%%功率谱 %f1=50; %figure %psd(spectrum.welch,st,fs,f1),title(调制信号功率谱密度); %grid tt=0.01:0.01:n; %%%%%%%%%%%%%%%%%%%%%高斯白噪声信道 stt=awgn(st,5); figure plot(tt,stt); axis([0,n+1,-1.5,1.5]); title(通过高斯信道后的信号); grid dt=stt.*sin(2*pi*w*tt); figure subplot(2,1,1) plot(tt,dt); axis([0,n+1,-1.5,1.5]); title(相乘后的波形); grid %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%低通滤波器 [N,Wn]=buttord(2*pi*25,2*pi*50,3,25,s);%临界频率采用角频率表示? [b,a]=butter(N,Wn,s); [bz,az]=impinvar(b,a,fs);%映射为数字的?? dt=filter(bz,az,dt); subplot(2,1,2) %15? plot(tt,dt); axis([0,n+1,-1.5,1.5]); title(低通滤波后的波形); grid %%%%%%%%%%%%%%%%%%%