【文章內(nèi)容簡介】 該依靠高版本的優(yōu)化功能把在低版本不能實(shí)現(xiàn)的功能強(qiáng)制執(zhí)行，所以我還有很多地方需要學(xué)習(xí)，包括對誤碼率進(jìn)行計(jì)算的兩個函數(shù)我也是借鑒的網(wǎng)上的程序。 這次方向設(shè)計(jì)讓我學(xué)會了很多，也認(rèn)識到了自己還有很多方面的欠缺。 八、 參考文獻(xiàn) [1] 李人厚 、 張平安 .精通 MATLAB， 西安交通大學(xué) [2] 肖明波，通信系統(tǒng)仿真原理與無線應(yīng)用 機(jī)械工業(yè)出版社 [3] 萬永革編著 ，通信系統(tǒng)仿真原理與無線應(yīng)用 ， 北京 科學(xué)出版社 [4] 郭文彬 、 桑林 ，通信原理 — 基于 MATLAB 的計(jì)算機(jī)仿真， 北京郵電大學(xué)出版社 ， [5]網(wǎng)上資料 附件： 主程序： T=1。 % 基帶信號寬度，也就是頻率 fc=10/T。 % 載波頻率 ml=2。 % 調(diào)制信號類型的一個標(biāo)志位 nb=100。 % 傳輸?shù)谋忍財(cái)?shù) delta_T=T/200。 % 采樣間隔 fs=1/delta_T。 % 采樣頻率 SNR=0。 % 信噪比 t=0:delta_T:nb*Tdelta_T。 % 限定 t 的取值范圍 N=length(t)。 % 采樣數(shù) % 調(diào)制部分 % 基帶信號的產(chǎn)生 data=randn(1,nb)。 datanrz=data.*21。 data1=zeros(1,nb/delta_T)。 for q=1:nb data1((q1)/delta_T+1:q/delta_T)=datanrz(q)。 end % 串并轉(zhuǎn)換，將奇偶位數(shù)據(jù)分開 idata=datanrz(1:ml:(nb1))。 qdata=datanrz(2:ml:nb)。 % QPSK 信號的調(diào)制 ich=zeros(1,nb/delta_T/2)。 for i=1:nb/2 ich((i1)/delta_T+1:i/delta_T)=idata(i)。 end 15 for ii=1:N/2 a(ii)=sqrt(2/T)*cos(2*pi*fc*t(ii))。 end idata1=ich.*a。 qch=zeros(1,nb/2/delta_T)。 for j1=1:nb/2 qch((j11)/delta_T+1:j1/delta_T)=qdata(j1)。 end for jj=1:N/2 b(jj)=sqrt(2/T)*sin(2*pi*fc*t(jj))。 end qdata1=qch.*b。 s=idata1+qdata1。 ss=abs(fft(s))。 % 瑞利衰落信道和高斯信道 % 瑞利衰落信道 ray_ich=raylrnd(,1,nb/2/delta_T)。 ray_qch=raylrnd(,1,nb/2/delta_T)。 Ray_idata=idata1.*ray_ich。 Ray_qdata=qdata1.*ray_qch。 Ray_s=Ray_idata+Ray_qdata。 % 高斯信道 s1=awgn(s,SNR)。 s11=abs(fft(s1))。 s111=s1s。 Awgn_s=awgn(Ray_s,SNR)。 % 通過高斯信道再通過瑞利衰落信道 % QPSK 解調(diào)部分 % 解調(diào)部分（高斯信道） idata2=s1.*a。 qdata2=s1.*b。 idata3=zeros(1,nb/2)。 qdata3=zeros(1,nb/2)。 % 抽樣判決的過程，與 0 作比較， data=0,則置 1，否則置 0 for n=1:nb/2 if sum(idata2((n1)/delta_T+1:n/delta_T))=0 idata3(n)=1。 else idata3(n)=0。 end 16 if sum(qdata2((n1)/delta_T+1:n/delta_T))=0 qdata3(n)=1。 else qdata3(n)=0。 end end % 為了顯示星座圖 ,將信號進(jìn)行處理 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 % 將判決之后的數(shù)據(jù)存放進(jìn)數(shù)組 demodata=zeros(1,nb)。 demodata(1:ml:(nb1))=idata3。 demodata(2:ml:nb)=qdata3。 %為了顯示，將它變成波形信號（即傳輸一個 1代表單位寬度的高電平） demodata1=zeros(1,nb/delta_T)。 for q=1:nb demodata1((q1)/delta_T+1:q/delta_T)=demodata(q)。 end % 累計(jì)誤碼數(shù) % abs(demodatadata)求接收端和發(fā)射端 % 數(shù)據(jù)差的絕對值，累計(jì)之后就是誤碼個數(shù) Awgn_num_BER=sum(abs(demodatadata)) % 解調(diào)部分（瑞利 +高斯） Ray_idata2=Ray_s.*a。 Ray_qdata2=Ray_s.*b。 % 為了顯示星座圖 ,將信號進(jìn)行處理 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。 17 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 % 將判決之后的數(shù)據(jù)存放進(jìn)數(shù)組 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)。 for q=1:nb Ray_demodata1((q1)/delta_T+1:q/delta_T)=Ray_demodata(q)。 end % 累計(jì)誤碼數(shù) % abs(demodatadata)求接收端和發(fā)射端 % 數(shù)據(jù)差的絕對值，累計(jì)之后就是誤碼個數(shù) Ray_num_BER=sum(abs(Ray_demodatadata)) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 誤碼率計(jì)算 %