畢業(yè)設(shè)計(jì)基于matlab的qpsk系統(tǒng)仿真-文庫(kù)吧資料
2025-06-25 20:33本頁(yè)面
【正文】 the given value of snr_in_dB, signal to noise ratio in dB.N=100。仿真Rayleigh39。理論Rayleigh39。仿真AWGN39。理論AWGN39。)。ylabel(39。Eb/No39。)。 semilogy(SNRindB1,smld_bit_ray_err_prb,39。 semilogy(SNRindB2,theo_err_ray_prb)。)。 semilogy(SNRindB1,smld_bit_awgn_err_prb,39。)。),title(39。end subplot(2,2,4) semilogy(SNRindB2,theo_err_awgn_prb,39。實(shí)際值（接收端）39。理論值（發(fā)射端）39。hold on。*39。axis([2 2 2 2])。)。),title(39。for i=1:nb/2plot(idata(i),qdata(i),39。)。,fs),title(39。psd(h,Ray_s,39。axis([0 500 5 5])。調(diào)制信號(hào)(Ray+Awgn)39。)。,39。legend(39。)。plot(Awgn_ichsum(i),Awgn_qchsum(i),39。hold on。QPSK信號(hào)星座圖（Awgn）39。r+39。subplot(2,2,4)。)。plot(s111),title(39。)。,fs),title(39。psd(h,s1,39。axis([0 500 5 5])。調(diào)制信號(hào)(Awgn)39。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 通過(guò)高斯信道figure(2)subplot(2,2,1)。解調(diào)輸出功率譜密度39。fs39。subplot(3,2,6)。)。plot(demodata1),title(39。)。,fs),title(39。psd(h,s,39。axis([0 500 3 3])。調(diào)制信號(hào)39。subplot(3,2,3)?；鶐盘?hào)功率譜密度39。fs39。subplot(3,2,2)。)。plot(data0),title(39。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%h = 。 % 高斯噪聲理論誤碼率 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ù)聲明在另外倆個(gè)M文件中%%作用： cm_sm32()用于瑞利信道誤碼率的計(jì)算%% cm_sm33()用于高斯信道誤碼率的計(jì)算%% ecoh on/off 作用在于決定是否顯示指令內(nèi)容%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%SNRindB1=0:1:6。 % 創(chuàng)建一個(gè)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:nb)=Ray_qdata4。 endend % 將判決之后的數(shù)據(jù)存放進(jìn)數(shù)組Ray_demodata=zeros(1,nb)。 if Ray_qchsum(n)=0 Ray_qdata4(n)=1。 else Ray_idata4(n)=0。for n=1:nb/2 Ray_ichsum(n)=sum(idata2((n1)/delta_T+1:n/delta_T))*delta_T。% end% end % 為了顯示星座圖,將信號(hào)進(jìn)行處理Ray_idata4=zeros(1,nb/2)。 % end% if Ray_sum(Ray_qdata2((n1)/delta_T+1:n/delta_T))=0% Ray_qdata3(n)=1。% 抽樣判決的過(guò)程，與0作比較，data=0,則置1，否則置0% for n=1:nb/2% if Ray_sum(Ray_idata2((n1)/delta_T+1:n/delta_T))=0% Ray_idata3(n)=1。 % 對(duì)應(yīng)的信號(hào)與正余弦信號(hào)相乘% Ray_idata3=zeros(1,nb/2)。 % 將極性碼變成對(duì)應(yīng)的波形信號(hào)end % 累計(jì)誤碼數(shù)% abs(demodatadata)求接收端和發(fā)射端% 數(shù)據(jù)差的絕對(duì)值，累計(jì)之后就是誤碼個(gè)數(shù)Awgn_num_BER=sum(abs(demodatadata)) %%%%%%%%%%%%%%%%%%%% 解調(diào)部分（瑞利+高斯）Ray_idata2=Ray_s.*a。 % 存放偶數(shù)位%為了顯示，將它變成波形信號(hào)（即傳輸一個(gè)1代表單位寬度的高電平）demodata1=zeros(1,nb/delta_T)。demodata(1:ml:(nb1))=idata3。 else qdata4(n)=0。 end Awgn_qchsum(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))*delta_T。 if Awgn_ichsum(n)=0 idata4(n)=1。qdata4=zeros(1,nb/2)。 else qdata3(n)=0。 end% A2(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))。 if sum(idata2((n1)/delta_T+1:n/delta_T))=0 idata3(n)=1。 % 建立1*nb數(shù)組，以存放解調(diào)之后的信號(hào)qdata3=zeros(1,nb/2)。 % 這里面其實(shí)隱藏了一個(gè)串并轉(zhuǎn)換的過(guò)程qdata2=s1.*b。 % 高斯噪聲曲線(xiàn)%%%%%%%%%%%%%%%%%%%%%%%%%%Awgn_s=awgn(Ray_s,SNR)。 % 通過(guò)高斯信道之后的信號(hào)s11=abs(fft(s1))。Ray_s=Ray_idata+Ray_qdata。Ray_idata=idata1.*ray_ich。 % 快速傅里葉變換得到頻譜%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 瑞利衰落信道和高斯信道%%%%%%%%%%%%%%%%%%%%%%%
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號(hào)-1