通信原理数字信号传输系统实验代码.doc

通信原理数字信号传输系统实验代码.doc

  1. 1、本文档共4页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
通信原理数字信号传输系统实验代码

clear %清除工作区所有变量,以防止已经存在的变量对程序执行的影响 clc %清空命令窗口 %数字信号传输系统实验xx.m文件 %xxx创建于xx/xx/xxxx length_of_bit=input('请输入比特流长度:'); while length_of_bit<=0 length_of_bit=input('比特流长度不能为负,请重新输入比特流长度:'); end dr=input('请输入占空比((0.0-1.0之间):'); while ~(dr>=0.0&&dr<1.0) dr=input('占空比取值必须在0.0-1.0之间,请重新输入占空比:'); end snr=input('请输入信噪比:'); bit=randint(1,length_of_bit);%产生一个输入长度的bit流 T=10;%码元宽度 tao=T*dr;%脉冲宽度 t=0:length_of_bit*T-1;%信号持续时长 fs=1e6;%采样频率 % %产生基带信号 sig=zeros(1,length_of_bit*T); A=5;%单极性归零码信号的幅度 for k0=1:length_of_bit if bit(k0)==1 sig(((k0-1)*T+1):(k0*T))=[ones(1,tao)*A,zeros(1,T-tao)]; end end SIG=fft(sig);%产生基带信号频谱 % %基带信号功率谱计算 sig_cxn=xcorr(sig,'unbiased'); nfft=1024; sig_CXk=fft(sig_cxn,nfft); sig_Pxx=abs(sig_CXk); index=0:round(nfft/2-1); k=index*fs/nfft; plot_sig_Pxx=10*log10(sig_Pxx(index+1)); % %画出基带信号波形,频谱及其功率谱 figure('Name','基带信号波形,频谱,功率谱','NumberTitle','off') subplot(3,1,1);plot(t,sig);title('基带信号波形');ylim([0,A+1]);xlabel('Time(sec)');ylabel('Voltage(v)'); subplot(3,1,2);plot(t,SIG);title('基带信号频谱');xlabel('Frequency(Hz)');ylabel('Amplitude(v)'); subplot(3,1,3);plot(k,plot_sig_Pxx);title('基带信号功率谱'); % %添加噪声 sig_noise=awgn(sig,snr); SIG_NOISE=fft(sig_noise); % %添加噪声后功率谱计算 sig_noise_cxn=xcorr(sig_noise,'unbiased'); nfft=1024; sig_noise_CXk=fft(sig_cxn,nfft); sig_noise_Pxx=abs(sig_noise_CXk); index=0:round(nfft/2-1); k=index*fs/nfft; plot_sig_noise_Pxx=10*log10(sig_noise_Pxx(index+1)); % %画出添加噪声后的信号波形 figure('Name','基带信号加噪声后的波形,频谱,功率谱','NumberTitle','off') subplot(3,1,1);plot(t,sig_noise);title('基带信号加噪声后波形');ylim([0,A+1]);xlabel('Time(sec)');ylabel('Voltage(v)'); subplot(3,1,2);plot(t,SIG_NOISE);title('基带信号加噪声后频谱');xlabel('Frequency(Hz)');ylabel('Amplitude(v)'); subplot(3,1,3);plot(k,plot_sig_noise_Pxx);title('基带信号加噪声后功率谱'); % %信号滤波 %低通滤波器设计 Wp=(1/(2*T)); Rp=1; Rs=50; n=1; [b,a]=ellip(n,Rp,Rs,Wp); sig_filter=filter(b,a,sig_noise);%进行滤波 SIG_FILTER=fft(sig_filter);%产生滤波后的频谱 % %滤波后信号功率谱计算 sig_filter_cxn=xcorr(sig_filter,'unbiased'); nfft=1024; sig_f

您可能关注的文档

文档评论(0)

haihang2017 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档