畢業設計----基于matlab的qpsk系統仿真設計與實現(編輯修改稿)

　

【文章內容簡介】 ata3(n)=1。 else idata3(n)=0。 end % A2(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))。 if sum(qdata2((n1)/delta_T+1:n/delta_T))=0 qdata3(n)=1。 else qdata3(n)=0。 end end % 為了顯示星座 圖 ,將信號進行處理 idata4=zeros(1,nb/2)。 qdata4=zeros(1,nb/2)。 for n=1:nb/2 Awgn_ichsum(n)=sum(idata2((n1)/delta_T+1:n/delta_T))*delta_T。 if Awgn_ichsum(n)=0 idata4(n)=1。 else idata4(n)=0。 end Awgn_qchsum(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))*delta_T。 if Awgn_qchsum(n)=0 qdata4(n)=1。 else qdata4(n)=0。 end end % 將判決之后的數據存放進數組 demodata=zeros(1,nb)。 demodata(1:ml:(nb1))=idata3。 % 存放奇數位 demodata(2:ml:nb)=qdata3。 % 存放偶數位 %為了顯示，將它變成波形信號（即傳輸一個 1 代表單位寬度的高電平） demodata1=zeros(1,nb/delta_T)。 % 創建一個 1*nb/delta_T 的零矩陣 for q=1:nb demodata1((q1)/delta_T+1:q/delta_T)=demodata(q)。 % 將極性碼變成對應的波形信號 end % 累計誤碼數 % abs(demodatadata)求接收端和發射端 % 數據差的絕對值，累計之后就是誤碼個數 Awgn_num_BER=sum(abs(demodatadata)) %%%%%%%%%%%%%%%%%%% % 解 調部分（瑞利 +高斯） Ray_idata2=Ray_s.*a。 % 這里面其實隱藏了一個串并轉換的過程 Ray_qdata2=Ray_s.*b。 % 對應的信號與正余弦信號相乘 % Ray_idata3=zeros(1,nb/2)。 % 建立 1*nb 數組，以存放解調之后的信號 % Ray_qdata3=zeros(1,nb/2)。 % 抽樣判決的過程，與 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。 % else Ray_idata3(n)=0。 % end % if Ray_sum(Ray_qdata2((n1)/delta_T+1:n/delta_T))=0 % Ray_qdata3(n)=1。 % else Ray_qdata3(n)=0。 % end % end % 為了顯示星座圖 ,將信號進行處理 Ray_idata4=zeros(1,nb/2)。 Ray_qdata4=zeros(1,nb/2)。 for n=1:nb/2 Ray_ichsum(n)=sum(idata2((n1)/delta_T+1:n/delta_T))*delta_T。 if Ray_ichsum(n)=0 Ray_idata4(n)=1。 else Ray_idata4(n)=0。 end Ray_qchsum(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))*delta_T。 if Ray_qchsum(n)=0 Ray_qdata4(n)=1。 else Ray_qdata4(n)=0。 end end % 將判決之后的數據存放進數組 Ray_demodata=zeros(1,nb)。 Ray_demodata(1:ml:(nb1))=Ray_idata4。 % 存放奇數位 Ray_demodata(2:ml:nb)=Ray_qdata4。 % 存放偶數位 %為了顯示，將它變成波形信號（即傳輸一個 1 代表單位寬度的高電平） Ray_demodata1=zeros(1,nb/delta_T)。 % 創建一個 1*nb/delta_T 的零矩陣 for q=1:nb Ray_demodata1((q1)/delta_T+1:q/delta_T)=Ray_demodata(q)。 % 將極性碼變成對應的波形信號 end % 累計誤碼數 % abs(demodatadata)求接收端和發射端 % 數據差的絕對值，累計之后就是誤碼個數 Ray_num_BER=sum(abs(Ray_demodatadata)) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 誤碼率計算 %% 調用了 cm_sm32()。和 cm_sm33(）函數 %%聲明： 函數聲明在另外倆個 M 文件中 %%作用： cm_sm32()用于瑞利信道誤碼率的計算 %% cm_sm33()用于高斯信道誤碼率的計算 %% ecoh on/off 作用在于決定是否顯示指令內容 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% SNRindB1=0:1:6。 SNRindB2=0::6。 % 瑞利衰落信道 for i=1:length(SNRindB1), [pb,ps]=cm_sm32(SNRindB1(i))。 % 比特誤碼率 smld_bit_ray_err_prb(i)=pb。 smld_symbol_ray_err_prb(i)=ps。 disp([ps,pb])。 echo off。 end。 % 高斯信道 echo on。 for i=1:length(SNRindB1), [pb1,ps1]=cm_sm33(SNRindB1(i))。 smld_bit_awgn_err_prb(i)=pb1。 smld_symbol_awgn_err_prb(i)=ps1。 disp([ps1,pb1])。 echo off。 end。 % 理論曲線 echo on。 for i=1:length(SNRindB2), SNR=exp(SNRindB2(i)*log(10)/10)。 % 信噪比 theo_err_awgn_prb(i)=*erfc(sqrt(SNR))。 % 高斯噪聲理論誤碼率 theo_err_ray_prb(i)=*(11/sqrt(1+1/SNR))。 % 瑞利衰落信道理論誤碼率 echo off。 end。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% h = 。 % 類似于 C 語言的宏定義，方便以下的調用 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%