久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

標(biāo)題: 2ASK 2PSK 2FSK在不同信噪比下三種調(diào)制方式誤碼性能MATLAB仿真 [打印本頁]
作者: BinaryIvan    時間: 2018-3-18 21:56
標(biāo)題: 2ASK 2PSK 2FSK在不同信噪比下三種調(diào)制方式誤碼性能MATLAB仿真
信號檢測與估計中，在不同信噪比下對比2ASK 2PSK 2FSK三種不同調(diào)制方式的誤碼性能對比
              采用相干解調(diào)時，在誤碼率相同的情況下，所需要的信噪比2ASK比2FSK高3dB，2FSK比2PSK高3dB,2ASK比2PSK高6dB。反過來，若信噪比r相同，2PSK系統(tǒng)的誤碼率比2FSK的小，2FSK系統(tǒng)的誤碼率比2ASK的小，相干解調(diào)的2PSK系統(tǒng)的誤碼率最小。由此看來，在抗加性高斯白噪聲方面，相干2PSK性能最好，2FSK次之，2ASK最差。

源程序：
  1. clear,clc;
  2. L=10000;
  3. t=0:0.01*pi:2*pi;
  4. SNR=-4:0.5:10;
  5. snr=SNR/10;
  6. for ii=1:length(snr)
  7.     Sigma(ii)=1/sqrt(2*(10^snr(ii)));%標(biāo)準(zhǔn)差
  8.     psk_error1=0; fsk_error1=0; ask_error1=0;   %誤碼數(shù)初始值為0
  9.     psk_error2=0; fsk_error2=0; ask_error2=0;
  10.     S=randi(L);
  11.     for k=1:L
  12.         n=Sigma(ii)*randn(1,length(t)); %一定信噪比下的高斯白噪聲
  13.        for i=1:length(t)
  14.             psk1(i)=sin(t(i)); psk2(i)=-sin(t(i));    %2PSK的S0(t)和S1(t)
  15.             fsk1(i)=sin(3*t(i)); fsk2(i)=sin(5*t(i));   %2FSK的S0(t)和S1(t)
  16.             ask1(i)=sin(t(i)); ask2(i)=0;          %2ASK的S0(t)和S1(t)      
  17.          end
  18.         %------------------2PSK
  19.         overlay11=sum(psk1.^2); overlay12=sum(psk2.^2);
  20.         psk1=psk1/sqrt(overlay11); psk2=psk2/sqrt(overlay12);%進(jìn)行歸一化
  21.         xp1=psk1+n; %接收信號（含有噪聲的）
  22.         xp2=psk2+n ;
  23.         %-----------------2FSK
  24.         overlay21=sum(fsk1.^2); overlay22=sum(fsk2.^2);
  25.         fsk1=fsk1/sqrt(overlay21); fsk2=fsk2/sqrt(overlay22);
  26.         xf1=fsk1+n;
  27.         xf2=fsk2+n;
  28.         %-----------------2ASK
  29.         overlay31=sum(ask1.^2);
  30.         ask1=ask1/sqrt(overlay31);
  31.         xa1=ask1+n;
  32.         xa2=ask2+n;                                    
  33.         if  S[k]==0  % 假設(shè)為H0的情況下
  34.            p1=0; f1=0; a1=0;
  35.            men1=0; men2=0; men3=0;
  36.             for i=1:length(t)
  37.                 p1=p1+xp1(i)*(psk2(i)-psk1(i));
  38.                 f1=f1+xf1(i)*(fsk2(i)-fsk1(i));
  39.                 a1=a1+xa1(i)*(ask2(i)-ask1(i));
  40.                 men1=men1+(psk2(i)^2-psk1(i)^2)/2;
  41.                 men2=men2+(fsk2(i)^2-fsk1(i)^2)/2;
  42.                 men3=men3+(ask2(i)^2-ask1(i)^2)/2;
  43.             end
  44.             if p1>men1
  45.                psk_error1=psk_error1+1;
  46.             end
  47.             if f1>men2
  48.                 fsk_error1=fsk_error1+1;
  49.             end
  50.             if a1>men3
  51.                 ask_error1=ask_error1+1;
  52.             end
  53.         else  %假設(shè)為 H1的情況下
  54.            p2=0; f2=0; a2=0; men1=0; men2=0; men3=0;
  55.             for i=1:length(t)
  56.                 p2=p2+xp2(i)*(psk2(i)-psk1(i));
  57.                 f2=f2+xf2(i)*(fsk2(i)-fsk1(i));
  58.                 a2=a2+xa2(i)*(ask2(i)-ask1(i));
  59.                 men1=men1+(psk2(i)^2-psk1(i)^2)/2;
  60.                 men2=men2+(fsk2(i)^2-fsk1(i)^2)/2;
  61.                 men3=men3+(ask2(i)^2-ask1(i)^2)/2;
  62.             end
  63.            if p2<men1
  64.                psk_error2=psk_error2+1;
  65.             end
  66.             if f2<men2
  67.                 fsk_error2=fsk_error2+1;
  68.             end
  69.             if a2<men3
  70.                 ask_error2=ask_error2+1;
  71.             end
  72.         end
  73.     end
  74.     %誤碼率
  75.     psk_pe(ii)=(psk_error1+psk_error2)/L;
  76.     fsk_pe(ii)=(fsk_error1+fsk_error2)/L;
  77.     ask_pe(ii)=(ask_error1+ask_error2)/L;
  78. end
  79. semilogy(SNR,psk_pe,'*'); hold on;
  80. semilogy(SNR,fsk_pe,'+'); hold on;
  81. semilogy(SNR,ask_pe,'o'); hold on;
  82. %理論值
  83. snr1=10.^(SNR/10);
  84. psk_pe_ideal=0.5*erfc(sqrt(snr1));
  85. fsk_pe_ideal=0.5*erfc(sqrt(0.5*snr1));
  86. ask_pe_ideal=0.5*erfc(sqrt(0.25*snr1));
  87. semilogy(SNR,psk_pe_ideal,'-b'); hold on;
  88. semilogy(SNR,fsk_pe_ideal,'-m'); hold on;
  89. semilogy(SNR,ask_pe_ideal,'-k'); hold on;
  90. axis([-4 10 10^(-5) 1]);
  91. legend('2PSK','2FSK','2ASK','2PSK-ideal','2FSK-ideal','2ASK-ideal');
  92. grid off;
  93. xlabel('E/n_0 dB');ylabel('P_e 誤碼率');
  94. title('不同信噪比下三種調(diào)制方式誤碼性能仿真');
復(fù)制代碼

完整的Word格式文檔51黑下載地址：
MATLAB仿真程序.zip (47.72 KB, 下載次數(shù): 26)



作者: BinaryIvan    時間: 2018-3-19 14:52
由蒙特卡洛實(shí)驗(yàn)實(shí)現(xiàn)MATLAB仿真



歡迎光臨 (http://m.zg4o1577.cn/bbs/) Powered by Discuz! X3.1
主站蜘蛛池模板： 日韩手机看片 | 手机看片1024日韩 | 欧美专区第一页 | 国产成人精品毛片 | 欧美久久一区二区 | 亚洲精品911 | 欧美午夜在线观看 | 免费视频一区二区 | 免费av观看 | 欧美一区二区三区在线观看 | 中文字幕亚洲一区 | 欧美在线性爱视频 | 美女无遮挡网站 | 一区二区三区高清 | 99色在线 | 成人激情视频 | 免费理论片 | 国产成人免费观看 | 夜夜操网站 | 国产伦精品一区二区三区88av | 亚洲成人中文字幕 | 日本成人中文字幕 | 欧美日韩大片 | 日日干av| 日韩一区二区视频在线观看 | 91精品国产日韩91久久久久久 | 日韩久久久久久久 | 99只有精品 | 欧美亚洲一区二区三区 | 欧美极品一区 | 日本免费一级片 | 久草视频免费看 | 日本黄网站 | 亚洲国产成人精品女人久久久 | 亚洲视频在线免费观看 | 精品九九九| 伊人网在线视频 | 国产三级黄色片 | 天堂在线中文资源 | 亚洲专区在线 | 91亚洲视频|