畢業(yè)設(shè)計----基于matlab的qpsk系統(tǒng)仿真設(shè)計與實現(xiàn)(參考版)
2024-12-06 23:42本頁面
【正文】 % since 2N bits are transmitted 。 end ps1=numofsymbolerror/(N*counter)。 end。 end。 if (decis2~=dsource2(i)), numofbiterror=numofbiterror+1。 symbolerror=1。 % Increment the error counter, if the decision is not correct symbolerror=0。 decis2=1。 decis2=0。 decis2=1。 decis2=0。 % The decision on the ith symbol is made next c_max=max([c00,c01,c10,c11])。 c10=dot(r,s10)。 % The correlation metrics are puted below c00=dot(r,s00)。 else r=s11+n。 elseif ((dsource1(i)==1) amp。 elseif ((dsource1(i)==0) amp。 % 2 normal distributed with 0, variance sgma if ((dsource1(i)==0) amp。 end。 dsource2(i)=1。 dsource2(i)=0。 dsource2(i)=1。 dsource2(i)=0。 while(numofbiterror100) for i=1:N, temp=rand。 numofbiterror=0。 s10=[0 1]。 s01=[0 1]。 % signal to noise ratio sgma=sqrt(E/snr)/2。 E=1。 % since there are totally N symbols pb=numofbiterror/(2*N*counter)。 end counter=counter+1。 if (symbolerror==1), numofsymbolerror=numofsymbolerror+1。 symbolerror=1。 end。 if (decis1~=dsource1(i)), numofbiterror=numofbiterror+1。 end。 else decis1=1。 elseif (c10==c_max), decis1=1。 elseif (c01==c_max), decis1=0。 if (c00==c_max), decis1=0。 c11=dot(r,s11)。 c01=dot(r,s01)。 end。 (dsource2(i)==0)), r=s10*ray+n。 (dsource2(i)==1)), r=ray*s01+n。 (dsource2(i)==0)), r=ray*s00+n。 n=sgma*randn(1,2)。 end。 dsource2(i)=1。 dsource2(i)=0。 dsource2(i)=1。 dsource2(i)=0。 % signal mapping % generation of the data source while(numofbiterror100) for i=1:N, temp=rand。 s11=[1 0]。 % noise variance s00=[1 0]。 snr=10^(snr_in_dB/10)。 numofbiterror=0。 E=1。)。,39。,39。,39。 legend(39。BER39。)。 xlabel(39。*39。hold on。hold on。r*39。hold on。誤碼率曲線 39。r39。)。,39。 legend(39。)。 plot(Ray_ichsum(i),Ray_qchsum(i),39。hold on。QPSK 信號星座圖（ Awgn+Ray） 39。r+39。 subplot(2,2,3)。調(diào)制信號功率譜密度 (Ray)39。fs39。 subplot(2,2,2)。)。 end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %通過高斯信道再通過瑞利衰落信道 figure(3) subplot(2,2,1) plot(Ray_s),title(39。實際值（接收端） 39。理論值（發(fā)射端） 39。hold on。*39。 axis([2 2 2 2])。)。),title(39。 for i=1:nb/2 plot(idata(i),qdata(i),39。 axis([0 2021 5 5])。高斯噪聲曲線 39。 subplot(2,2,3)。調(diào)制信號功率譜密度 (Awgn)39。fs39。 subplot(2,2,2)。)。 plot(s1),title(39。)。,fs),title(39。 psd(h,demodata1,39。 axis([0 20210 2 2])。解調(diào)輸出 39。 subplot(3,2,5)。調(diào)制信號功率譜密度 39。fs39。 subplot(3,2,4)。)。 plot(s),title(39。)。,fs),title(39。 psd(h,data1,39。 axis([0 20210 2 2])?；鶐盘?39。 % 類似于 C 語言的宏定義，方便以下的調(diào)用 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%% 輸出顯示部分 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 第一部分（理想） figure(1) subplot(3,2,1)。 end。 % 高斯噪聲理論誤碼率 theo_err_ray_prb(i)=*(11/sqrt(1+1/SNR))。 for i=1:length(SNRindB2), SNR=exp(SNRindB2(i)*log(10)/10)。 end。 disp([ps1,pb1])。 smld_bit_awgn_err_prb(i)=pb1。 % 高斯信道 echo on。 echo off。 smld_symbol_ray_err_prb(i)=ps。 % 瑞利衰落信道 for i=1:length(SNRindB1), [pb,ps]=cm_sm32(SNRindB1(i))。和 cm_sm33(）函數(shù) %%聲明： 函數(shù)聲明在另外倆個 M 文件中 %%作用： cm_sm32()用于瑞利信道誤碼率的計算 %% cm_sm33()用于高斯信道誤碼率的計算 %% ecoh on/off 作用在于決定是否顯示指令內(nèi)容 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% SNRindB1=0:1:6。 % 創(chuàng)建一個 1*nb/delta_T 的零矩陣 for q=1:nb Ray_demodata1((q1)/delta_T+1:q/delta_T)=Ray_demodata(q)。 % 存放奇數(shù)位 Ray_demodata(2:ml
公司管理相關(guān)推薦
鄂ICP備17016276號-1