基于matlab的移動衰落信道仿真畢業設計說明書-資料下載頁

　　

【正文】 復確定性過程的自相關函數完全趨近于參考模型的自相關函數。當越大，范圍一定時，值越大。我們知道，從發射端到接收端之一過程中，信道的路徑數越少，信道越發雜，路徑數越多，信道越好，所以為了改變這一情況，我們可以通過改變諧波函數個數來解決。當諧波函數個數值越大，參考模型波形越接近仿真模型的波形，兩者波形非常吻合，信道模型越好。最后，本文還有一些不足，例如只在仿真過程中，只考慮了JAKES法的使用等，在以后的學習過程中有待進一步的研究與提高。參考文獻[1][M]，2010.[2]立寧，樂光新，[M].人民郵電出版社，1999.[3]劉保柱，蘇彥華，[M].人民郵電出版社，2010.[4] 王欣, 酆廣增. 多徑非相關瑞利信道生成的改進[J].通信學報，2007，28（5）:123125.[5][M]..[6]汪安春，周世東，姚彥. 無線衰落信道中信道預測的誤差分析[J].清華大學學報，2004，44（10）:14311436.[7]周富相，[J].河北北方學院學報，2010，26（3）:2124.[8] 代光發，[J].系統仿真學報，2005，17（1）:215237.[9][M]..[10][D]..[11]黃舒，嚴云保，[J].保山師專學報，2009，28（2）:1314.[12]王繪宇. 移動信道建模與仿真及其在點到多點通信系統中的應用研究[D].重慶大學，2007.[13]趙海濤，[J].電力系統通信，2008，29（9）：3543.[14][J].中國傳媒大學學報，2007，14（4）：5984.[15][J].黑龍江科技信息，2011，（12）：94125.致謝大學生活一晃而過，回首走過的歲月，心中倍感充實，本論文的順利完成，離不開各位老師、同學和朋友的關心和幫助。先誠摯的感謝我的論文指導老師賈子彥老師。他在忙碌的教學工作中擠出時間來審查、修改我的論文。還有教過我給過我幫助的所有老師們，你們嚴謹細致、一絲不茍的作風一直是我工作、學習中的榜樣；他們循循善誘的教導和不拘一格的思路給予我無盡的啟迪。感謝四年中陪伴在我身邊的同學和朋友，感謝他們為我提出的有益的建議和意見，有了他們的支持、鼓勵和幫助，我才能充實的度過了四年的學習生活。附錄本次設計運用的MATLAB仿真程序如下：function r_mm=acf_mue(f,c,tau)r_mm=0。for n=1:length(c), r_mm=r_mm+*c(n)^2*cos(2*pi*f(n)*tau)。endMETHOD=39。jm_j39。N_i=9。f_max=91。sigma_0=1。if METHOD==39。jm_j39。, n=1:N_i1。 f_i_n=f_max*[[cos(pi*n/(2*(N_i1/2))),1]39。,... [cos(pi*n/(2*(N_i1/2))),1]39。]。 c_i_n=2*sigma_0/sqrt(N_i1/2)*[[sin(pi*n/(N_i1)),1/2]39。,... [cos(pi*n/(N_i1)),1/2]39。]。 %K=2。 theta_i_n=zeros(size(f_i_n))。 PHASE=39。none39。else error (39。Method is unknown39。)endif PHASE==39。rand39。, theta_i_n=rand(N_i,1)*2*pi。 elseif PHASE==39。perm39。, n=(1:N_i)39。 Z=rand(size(n))。 [dummy,I]=sort(Z)。 theta_i_n=2*pi*n(I)/(N_i+1)。end。 if METHOD==39。jm_j39。 subplot(2,3,1) stem([f_i_n(N_i:1:1,1)。f_i_n(:,1)],... 1/4*[c_i_n(N_i:1:1,1)。c_i_n(:,1)].^2) title(39。i=139。) xlabel(39。f (Hz)39。) ylabel(39。PSD39。) subplot(2,3,2) stem([f_i_n(N_i:1:1,2)。f_i_n(:,2)],... 1/4*[c_i_n(N_i:1:1,2)。c_i_n(:,2)].^2) title(39。i=239。) xlabel(39。f (Hz)39。) ylabel(39。PSD39。) tau_max=N_i/(2*f_max)。 tau=linspace(0,1000)。 r_mm=sigma_0^2*besselj(0,2*pi*f_max*tau)。 r_mm_tilde1=acf_mue(f_i_n(:,1),c_i_n(:,1),tau)。 subplot(2,3,4) plot(tau,r_mm,39。r39。,tau,r_mm_tilde1,39。g39。) title(39。i=139。) xlabel(39。tau (s)39。) ylabel(39。ACF39。) r_mm_tilde2=acf_mue(f_i_n(:,2),c_i_n(:,2),tau)。 subplot(2,3,5) plot(tau,r_mm,39。r39。,tau,r_mm_tilde2,39。g39。) title(39。i=239。) xlabel(39。tau (s)39。) ylabel(39。ACF39。) subplot(2,3,3) stem([f_i_n(N_i:1:1,1)。f_i_n(:,1)],... 1/4*[c_i_n(N_i:1:1,1)。c_i_n(:,1)].^2+... 1/4*[c_i_n(N_i:1:1,2)。c_i_n(:,2)].^2) title(39。i=1,239。) xlabel(39。f (Hz)39。) ylabel(39。PSD39。) subplot(2,3,6) plot(tau,2*r_mm,39。r39。,tau,r_mm_tilde1+r_mm_tilde2,39。g39。) title(39。i=1,239。) xlabel(39。tau (s)39。) ylabel(39。ACF39。)end