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

標(biāo)題: 2ASK 2PSK 2FSK在不同信噪比下三種調(diào)制方式誤碼性能MATLAB仿真 [打印本頁(yè)]
作者: BinaryIvan    時(shí)間: 2018-3-18 21:56
標(biāo)題: 2ASK 2PSK 2FSK在不同信噪比下三種調(diào)制方式誤碼性能MATLAB仿真
信號(hào)檢測(cè)與估計(jì)中，在不同信噪比下對(duì)比2ASK 2PSK 2FSK三種不同調(diào)制方式的誤碼性能對(duì)比
              采用相干解調(diào)時(shí)，在誤碼率相同的情況下，所需要的信噪比2ASK比2FSK高3dB，2FSK比2PSK高3dB,2ASK比2PSK高6dB。反過(guò)來(lái)，若信噪比r相同，2PSK系統(tǒng)的誤碼率比2FSK的小，2FSK系統(tǒng)的誤碼率比2ASK的小，相干解調(diào)的2PSK系統(tǒng)的誤碼率最小。由此看來(lái)，在抗加性高斯白噪聲方面，相干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; %接收信號(hào)（含有噪聲的）
  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    時(shí)間: 2018-3-19 14:52
由蒙特卡洛實(shí)驗(yàn)實(shí)現(xiàn)MATLAB仿真



歡迎光臨 (http://www.zg4o1577.cn/bbs/) Powered by Discuz! X3.1
主站蜘蛛池模板： 欧美精品一区二区在线观看 | av网站免费看 | 国产成人综合一区二区三区 | 午夜精品一区二区三区在线观看 | 亚洲国产精品成人久久久 | 91综合网 | 国产免费一区 | 欧美日韩成人 | 久久久青草婷婷精品综合日韩 | 欧美成人一级视频 | 97成人精品 | 亚洲一区中文字幕 | 新91视频网| 亚洲最大福利网 | 成人精品视频在线观看 | 亚洲高清视频一区二区 | 一区在线免费视频 | 成人毛片视频免费 | 欧美国产精品 | 青草青草久热精品视频在线观看 | 精品日韩一区二区三区av动图 | av在线播放一区二区 | 98久久 | 天天干视频 | a级黄色片视频 | 欧美一区二区免费视频 | 国产欧美一区二区三区在线播放 | 国产91亚洲精品一区二区三区 | 国产精品久久久久久久久免费高清 | 日本又色又爽又黄又高潮 | 日韩精品在线免费 | 午夜国产羞羞视频免费网站 | 91精品久久久久 | 亚洲精品一区国语对白 | 欧美一区成人 | 国产精品高潮呻吟久久 | www.天天操.com | 91免费观看 | 日本三级日产三级国产三级 | 操久久| 国产精品色综合 |