代码与截图.docx

下载文档 降价啦

22
0
约7.8千字
约 16页
2017-07-06 发布于湖北
举报
版权申诉
保障服务

代码与截图.docx

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

代码与截图

%主函数Ts=1e-3; %采样间隔t=0:Ts:50*Ts; %仿真时间序列x=sin(2*pi*50*t); %信号figure(1);subplot(2,1,1);plot(t,x);title(原始信号);subplot(2,1,2);stem(t,x);title(抽样信号);y1=PCM(x); figure(2);subplot(2,1,1);stairs(y1);axis([0 30 0 1.2]);title(PCM调制);y2=CMIcoding(y1); subplot(2,1,2);stairs(y2);axis([0 30 0 1.2]);title(CMI编码);y3=Cycliccoding(y2); figure(3);subplot(211);stairs(y2);axis([0 30 0 1.2]);title(循环编码); fc1=2000; %信号为0时的频率 fc1 fc2=5000; %信号为1时的频率 fc2 fb=1000; %波特率 f fs=32000; %采样频率64kHz(中频采样定理) fsflag_k=1;[fsksig, fs, fb, fc1, fc2] = FSK(y3, fs, fb, fc1, fc2,flag_k);%fsk调制disp(调制后码长:);l_fsksig=length(fsksig) %24n*fs/fbRaylei=Rayleighchannel(fsksig);%幅值放大，加噪声figure(5)stairs(Raylei)title(加噪声后信号)axis([1700 2000,-1.5 1.5])[y4, x_p] = FSKdet(fsksig, fs,fb, fc1,fc2,flag_k);disp(解调后码长:);l_newCode=length(y4)%24ny5=Cyclicdecoding(y4);figure(3);subplot(212);stairs(y5);axis([0 30 0 1.2]);title(循环解码);y6=CMIdecoding(y5);figure(7);subplot(2,1,1);stairs(y6);axis([0 60 0 1.2]);title(CMI解码);y7=D_PCM(y6);subplot(2,1,2);plot(t,y7);axis([0 50*Ts,-1 1]);title(PCM解调);da=0;fori=length(t)dc=(x(i)-y7(i))^2/length(t);da=da+dc;enddisp(失真度为:);da%FSK 调制 fc1=2000; %信号为0时的频率 fc1 fc2=5000; %信号为1时的频率 fc2 fb=1000; %波特率 f fs=32000; %采样频率64kHz(中频采样定理) fsflag_k=1;[fsksig, fs, fb, fc1, fc2] = FSK(y3, fs, fb, fc1, fc2,flag_k);%fsk调制disp(调制后码长:);l_fsksig=length(fsksig) %24n*fs/fb%噪声Raylei=Rayleighchannel(fsksig);%幅值放大，加噪声figure(5)stairs(Raylei)title(加噪声后信号)axis([1700 2000,-1.5 1.5])%FSK 解调[y4, x_p] = FSKdet(fsksig, fs,fb, fc1,fc2,flag_k);%fsk解调disp(解调后码长:);l_newCode=length(y4)%24ny5=Cyclicdecoding(y4);figure(7);stairs(y5);axis([0 30 0 1.2]);title(循环解码);y6=CMIdecoding(y5);figure(8);subplot(2,1,1);stairs(y6);axis([0 60 0 1.2]);title(CMI解码);y7=D_PCM(y6);subplot(2,1,2);plot(t,y7);axis([0 50*Ts,-1 1]);title(PCM解调);%计算失真度da=0;fori=length(t)dc=(x(i)-y7(i))^2/length(t);da=da+dc;enddisp(失真度为