- 1、本文档共4页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 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
您可能关注的文档
- 广东省佛山市南海区黄岐高中学业水平测试生物模拟卷三(附答案).doc
- 常州市教育学会学业水平监测七年级数学试题2015年2月.doc
- 工程数值方法上(2015-3-17讲稿).doc
- 建设数字化示范校园的实践思考.doc
- 微分方程数值解上机报告第二章报告.doc
- 微分方程数值解打靶法求边值问题.doc
- 开题报告黄智渊.doc
- 广东省深圳市20152016学年高一物理上学期期中试题.doc
- 幼儿园数学教育目标体系的再认识模式(.docx
- 微机实验五0809AD转换器实验.doc
- 2024年山东德州夏津县融媒体中心招聘4人历年高频考题难、易错点模拟试题(共500题)附带答案详解.docx
- 2024年山东济南历城区卫生健康局所属事业单位招聘50人历年高频考题难、易错点模拟试题(共500题)附带答案详解.docx
- 2024江苏徐州锡沂高新技术产业开发区招聘50人历年高频难、易点(行政职业能力测验共200题含答案解.docx
- 2024年山东德州职业技术学院招聘49人历年高频考题难、易错点模拟试题(共500题)附带答案详解.docx
- 2024年山东德州高级师范学校招聘工作人员12人高频考题难、易错点模拟试题(共500题)附带答案详解.docx
- 2023年12月恩施巴东县“才聚荆楚·智汇恩施”第二次专项事业单位工作人员7名笔试近6年高频考题难、易错点荟萃答案带详解附后.docx
- 2023年12月广东广州市黄埔区机关事务管理局招考聘用政府雇员笔试近6年高频考题难、易错点荟萃答案带详解附后.docx
- 2023年12月广西马山县纪委监委机关面向社会公开招考基层廉洁监督员及防贫监测员笔试近6年高频考题难、易错点荟萃答案带详解附后.docx
- 2023年12月广东广州市增城区招考聘用教师(上海设点)92人笔试近6年高频考题难、易错点荟萃答案带详解附后.docx
- 2023年12月广东省茂名市电白区赴高校现场公开招聘62名教师314笔试近6年高频考题难、易错点荟萃答案带详解附后.docx
文档评论(0)