通信原理课程设计1.docVIP

下载本文档

23
0
约8.83万字
约 10页
2017-02-11 发布于河南
举报
版权申诉

通信原理课程设计1.doc

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

通信原理课程设计1

附录Ⅰ：数字基带信号的仿真代码与结果 function[f,sf]=T2F(t,st) dt=t(2)-t(1); T=t(end); df=1/T; N=length(st); f=-N/2*df:df:N/2*df-df; sf=fft(st); sf=T/N*fftshift(sf); end %在把序列d插入到序列M中,得到一个新序列 function [out]=sigexpand(d,M) N=length(d); out=zeros(M,N); out(1,:)=d; out=reshape(out,1,M*N); %数字基带信号的功率谱密度 digit_baseband.m clear all; close all; Ts=1;%每个码元的长度 N_sample=8;%每个码元的抽样点数 dt=Ts/N_sample;%抽样时间间隔 N=1000;%码元数 t=0:dt:(N*N_sample-1)*dt;%1x8000 T=t(end); gt1=ones(1,N_sample);%1x8 gt2=ones(1,N_sample/2); % 1x4?1x8? gt2=[gt2 zeros(1,N_sample/2)];%1x8 mt3=sinc((t-5)/Ts);% 1x8000 sin(pi*t/Ts)/(pi*t/Ts)波形 gt3=mt3(1:10*N_sample);%截段取10个码元 d=(sign(randn(1,N))+1)/2;%1x1000 d=?? data=sigexpand(d,N_sample);%插入N_sample-1个0 gt4=ones(1,N_sample); gt5=ones(1,N_sample/2); gt5=[gt5 zeros(1,N_sample/2)]; st1=conv(data,gt1);%调用Matlab的卷积函数conv st2=conv(data,gt2); d=2*d-1; % 变成双极性序列 data=sigexpand(d,N_sample); st3=conv(data,gt3); st4=conv(data,gt4); st5=conv(data,gt5); [f,st1f]=T2F(t,[st1(1:length(t))]); [f,st2f]=T2F(t,[st2(1:length(t))]); [f,st3f]=T2F(t,[st3(1:length(t))]); [f,st4f]=T2F(t,[st4(1:length(t))]); [f,st5f]=T2F(t,[st5(1:length(t))]); figure(1) subplot(321) plot(t,[st1(1:length(t))]); axis([0 20 -1.5 1.5]); xlabel(单极性NRZ波形);grid subplot(322); plot(f,10*log10(abs(st1f).^2/T));grid axis([-5 5 -40 10]); xlabel(单极性NRZ功率谱密度(dB/Hz)); subplot(323) plot(t,[st2(1:length(t))]);grid axis([0 20 -1.5 1.5]); xlabel(单极性RZ波形); subplot(324); plot(f,10*log10(abs(st2f).^2/T));grid axis([-5 5 -40 10]); xlabel(单极性RZ功率谱密度(dB/Hz)); subplot(325) plot(t-5,[st3(1:length(t))]);grid axis([0 20 -2 2]); ylabel(双极性sinc波形);xlabel(t/Ts); subplot(326); plot(f,10*log10(abs(st3f).^2/T));grid axis([-5 5 -40 10]); ylabel(sinc波形功率谱密度(dB/Hz));xlabel(f*Ts); figure(2) subplot(221) plot(t,[st4(1:length(t))]); axis([0 20 -1.5 1.5]); xlabel(双极性NRZ波形);grid subplot(222); plot(f,10*log10(abs(st4f).^2/T));grid axis([-5 5 -40 10]); xlabel(双极性NRZ功率谱密度(dB/Hz)); subplot(223) plot(t,[st5(1:length(t))]);grid axis([0 20 -1.5 1.5]); xlabel(双极性RZ波