下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、程序1clear allN=10000;x=randn(1,N);%产生均值0,方差为1的随机序列for i=1:N if x(i)0.5 b(i)=1; else b(i)=-1; endend %调制xx=pskmod(x,M); y=0.04 -0.05 0.07 -0.21 -0.5 0.72 0.36 0 0.21 0.03 0.07; %衰减信道A1=zeros(5,5);for m=1:5 for n=1:5 if n=m A1(m,n)=0.72; else if n=m+1 A1(m,n)=-0.5; else if n=m-1 A1(m,n)=0.36; else if n
2、=m+2 A1(m,n)=-0.21; else if n=m-2 A1(m,n)=0; else A1(m,n)=0; end end end end end endendB1=0;0;1;0;0;C1=pinv(A1)*B1;A2=zeros(17,17);for i=1:17 for j=1:17 if j=i A2(i,j)=0.72; else if j=i-1 A2(i,j)=0.36; else if j=i-2 A2(i,j)=0; else if j=i-3 A2(i,j)=0.21; else if j=i-4 A2(i,j)=0.03; else if j=i-5 A2(
3、i,j)=0.07; else if j=i+1 A2(i,j)=-0.5; else if j=i+2 A2(i,j)=-0.21; else if j=i+3 A2(i,j)=0.07; else if j=i+4 A2(i,j)=-0.05; else if j=i+5 A2(i,j)=0.04; else A2(i,j)=0; end end endendendendendendendendendendendB2=0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;C2=pinv(A2)*B2;A=zeros(31,31);for i=1:31 for j=1:31 i
4、f j=i A(i,j)=0.72; else if j=i-1 A(i,j)=0.36; else if j=i-2 A(i,j)=0; else if j=i-3 A(i,j)=0.21; else if j=i-4 A(i,j)=0.03; else if j=i-5 A(i,j)=0.07; else if j=i+1 A(i,j)=-0.5; else if j=i+2 A(i,j)=-0.21; else if j=i+3 A(i,j)=0.07; else if j=i+4 A(i,j)=-0.05; else if j=i+5 A(i,j)=0.04; else A(i,j)=
5、0; end end endendendendendendendendendendend B=0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;C=pinv(A)*B;SNR=1:2:21;for m=1:length(SNR) for i=1:N noise=randn; x1noise(i)=b(i)*sqrt(2*10(SNR(m)/10)+noise; if x1noise(i)0; t1(i)=-1; else t1(i)=1; end end %解调:X1=pskdemod(x1noise,M); numbers
6、,pm1(m)=symerr(b,t1); xxy=conv(b,y); for i=1:N+6 noise=randn; x2noise(i)=xxy(i)*sqrt(2*10(SNR(m)/10)+noise; if x2noise(i)0; t2(i)=-1; else t2(i)=1; end end %解调:X2=pskdemod(x2noise,M); numbers,pm2(m)=symerr(b,t2(6:length(x)+5); x2C1=conv(x2noise,C1); for i=1:N+8 if x2C1(i)0; t3(i)=-1; else t3(i)=1; e
7、nd end %解调 X2C1=pskdemod(x2C1,M) numbers,pm3(m)=symerr(b,t3(8:length(x)+7); x2C2=conv(x2noise,C2); for i=1:N+14 if x2C2(i)0; t4(i)=-1; else t4(i)=1; end end %解调X2C2=pskdemod(x2C2,M); numbers,pm4(m)=symerr(b,t4(14:length(x)+13); x2C=conv(x2noise,C); for i=1:N+21 if x2C(i)0; t5(i)=-1; else t5(i)=1; en
8、d end %X2C=pskdemod(x2C,M) numbers,pm5(m)=symerr(b,t5(21:length(x)+20);endfigure(1)semilogy(SNR,pm1);grid on;title(无串扰的误码率);axis(1 16 10-5 1);xlabel(SNR(dB));ylabel(误码率);figure(2)semilogy(SNR,pm2);grid on;title(经过衰减信道的误码率);xlabel(SNR(dB));ylabel(误码率);figure(3)semilogy(SNR,pm2,r-,SNR,pm3,ms:);grid on
9、;title(5抽头迫零均衡器);axis(1 16 10-2 1);xlabel(SNR(dB));ylabel(误码率);legend(多径串扰信道,5抽头)figure(4)semilogy(SNR,pm2,r-,SNR,pm4,ms:);grid on;title(17抽头迫零均衡器);axis(1 16 10-5 1);xlabel(SNR(dB));ylabel(误码率);legend(多径串扰信道,17抽头)figure(5)semilogy(SNR,pm2,r-,SNR,pm5,ms:);grid on;title(31抽头迫零均衡器);axis(1 16 10-5 1);xlabel(SNR(dB));ylabel(误码率);legend(多径串扰信道,31抽头)figure(6)semilogy
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 医疗仪器在妇产科的应用实例
- 医药研发合作协议书
- 城市肠道传染病应急处置策略
- 文案策划劳务合同
- 珠宝行业卫生防疫规定
- 春季流感防治预案
- 农业买卖合同模版
- 商业用租赁合同
- 物业满意度提升方案
- 2024春六年级下语文期中归类整合复习-字词、句子品析、课文积累
- 物理自主招生考试试卷
- 2023年出版专业资格考试基础中级真题及答案
- 少儿美术绘本教案课件-3-6岁 《100层巴士》
- 部编版道德与法治四年级下册第二单元《做聪明的消费者》单元作业设计
- 申报绿色施工示范工程资料
- 带式输送机传动装置中的一级圆柱齿轮减速器设计说明书(毕业论文)
- 二太阳能前期开发工作流程
- 苏科版三年级下册劳动第3课《纸黏土浮雕》课件
- 2021年度全矿性反风技术方案
- 2022年福建省厦门市集美区初中毕业年级模拟考试英语试题及答案
- GB/T 5836.1-1992建筑排水用硬聚氯乙烯管材
评论
0/150
提交评论