【正文】 %dff = df/8。 sf = fft(st)。 %f=N/2*df+df/2:df:N/2*dfdf/2。 df = 1/T。s Fourier %Translation %Input is the time and the signal vectors,the length of time must greater %than 2 %Output is the frequency and the signal spectrum dt = t(2)t(1)。 out = reshape(out,1,M*N)。 out = zeros(M,N)。 axis([2 2 1 1])。f39。)。title(39。)。xlabel(39。信道幅頻特性 39。 axis([2 2 0 2])。LineWidth39。)。xlabel(39。輸出信號幅度譜 39。 axis([2 2 0 60])。LineWidth39。 subplot(224)。輸入信號幅度譜 39。 axis([2 2 0 60])。LineWidth39。 subplot(222)。t39。)。title(39。,2)。 subplot(223) plot(t,yt1,39。輸入信號 39。 axis([20 40 0 ])。LineWidth39。 [f yf1] = T2F(t,yt1)。 hf = m(1)*exp(j*2*pi*f*tao(1)) +m(2)* exp(j*2*pi*f*tao(2))+m(3)* 齊魯工業(yè)大學(xué) 2020 屆本科畢業(yè)生畢業(yè)設(shè)計(jì)（論文） exp(j*2*pi*f*tao(3))。 %3 徑信道 m=[ ]。 %對序列間隔插入 N_sample1 個 0 st1 = conv(data,gt1)。 %NRZ 非歸零波形 d = ( sign( randn(1,N) ) +1 )/2。 dLen = length(t)。 %抽樣時間間隔 N = 1000。 N_sample = 8。 close all。 %f = NN/2*dff:dff:NN/2*dffdff。 %mf = sinc(ff*T)。 %dff = df/8。 sf = fft(st)。 %f=N/2*df+df/2:df:N/2*dfdf/2。 df = 1/T。s Fourier %Translation %Input is the time and the signal vectors,the length of time must greater %than 2 %Output is the frequency and the signal spectrum dt = t(2)t(1)。)。title(39。)。 ylabel(39。f39。 plot(ff,abs( rf ) )。 axis([20 20 0 5])。輸入單頻信號頻譜 39。)。 ylabel(39。f39。 plot(ff,abs( sf ) )。 [ff rf]=T2F(t,rt)。)。 title(39。s(t)39。)。 xlabel(39。 axis([0 2 ])。輸入單頻信號 39。)。 ylabel(39。t39。 %將信號經(jīng)過 20 徑的結(jié)果相加 figure(1) subplot(211) plot(t, st)。 s(i, : ) = cos(2*pi*f*t + fai(i,: ) + fai0(i) )。 %路徑的初始相位 st = cos(2*pi*f*t)。 %徑數(shù) 齊魯工業(yè)大學(xué) 2020 屆本科畢業(yè)生畢業(yè)設(shè)計(jì)（論文） taof = *rand(1,L)。 t = 0:dt:10。 f = 10。 %多徑時變 clear all。 sff = ifftshift(sf)。 %t=T/2:dt:T/2dt。 N = length(sf)。 Fmx = ( f(end)f(1) +df)。 function [t st]=F2T(f,sf) %This function calculate the time signal using ifft function for the input %signal39。 %sf = conv(mf,sf) %NN = length(f)+length(ff)1。 %ff = 50*df:dff:50*df。 sf = T/N*fftshift(sf)。 齊魯工業(yè)大學(xué) 2020 屆本科畢業(yè)生畢業(yè)設(shè)計(jì)（論文） f=N/2*df:df:N/2*dfdf。 N = length(st)。 T=t(end)。 function [f,sf]= T2F(t,st) %This is a function using the FFT function to calculate a signal39。 out(1,:) = d。 函數(shù) function [out]=sigexpand(d,M)%將輸入的序列擴(kuò)展成間隔為 M1個 0的序列 N = length(d)。t39。grid on。 plot(t4,real(yt4))。)。xlabel(39。)。 title(39。幅頻、相頻特性 39。 subplot(223) plotyy(f,abs(hf4),f,angle(hf4)/pi)。 axis([0,20, ])。經(jīng)過信道后的輸出信號 39。 plot(t3,real(yt3))。grid on。相頻失真、群時延無失真信道 39。)。ylabel(39。 齊魯工業(yè)大學(xué) 2020 屆本科畢業(yè)生畢業(yè)設(shè)計(jì)（論文） figure(2) subplot(221)。t39。grid on。) subplot(224) plot(t2,real(yt2))。 xlabel(39。)。 title(39。幅頻、相頻特性 39。 subplot(223) plotyy(f,abs(hf2),f,angle(hf2)/pi)。 axis([0,20, ])。經(jīng)過信道后的輸出信號 39。 subplot(222) plot(t1,real(yt1) )。)。 title(39。幅頻、相頻特性 39。 figure(1) subplot(221) plotyy(f,abs(hf1),f,angle(hf1)/pi)。 yf4 = xf.*hf4。 [t3,yt3]=F2T(f,yf3)。 hf3(f1) = exp( j*pi*f(f1)j*pi )。 f1 = find(f0)。 [t2,yt2] = F2T(f,yf2)。 %幅頻失真信道 齊魯工業(yè)大學(xué) 2020 屆本科畢業(yè)生畢業(yè)設(shè)計(jì)（論文） hf2 = sinc(f).*exp(j*pi*f)。 yf1 = xf.*hf1。 %無失真信道 [f,xf] = T2F(t,xt)。 st3 = conv(data,gt3)。 d = 2*d1。 %對序列間隔插入 N_sample1 個 0 st1 = conv(data,gt1)。 d = ( sign( randn(1,N) ) +1 )/2。 mt3 = sinc((t5)/Ts)。 %NRZ 非歸零波形 gt2 = ones(1,N_sample/2)。 %碼元數(shù) t = 0:dt:(N*N_sample1)*dt。 %每個碼元的抽樣點(diǎn)數(shù) dt = Ts/N_sample。 Ts=1。 最后衷心感謝電氣學(xué)院的所有給予我指導(dǎo)和幫助的領(lǐng)導(dǎo)、老師和同學(xué)，衷心的感謝你們 ！ 齊魯工業(yè)大學(xué) 2020 屆本科畢業(yè)生畢業(yè)設(shè)計(jì)（論文） 附 錄 主程序： %信道失真示意 clear all。在論文完成之際，首先衷心感謝林老師對我的深切關(guān)懷和諄諄教導(dǎo)。 [5] 鄧華 .MATLAB 通信仿真及應(yīng)用實(shí)例詳解 .北京：人民郵電出版社， :294~318 [6] Andrea Goldsmith、 Stanford University. WIRELESS COMMUNICATIONS. Cambridge University [7] , . MultipleSymbol Differential Detection of MPSK[J].IEEETransactions on :300～ 308. [8] 曹志剛、錢亞生 .現(xiàn)代通信原理 .北京：清華大學(xué)出版社， 1992 [9] 求是科技 .MATLAB 從入門到精通 .人民郵電出版社 .2020 [10] 郭文斌、桑林《通信原理：基于 MATLAB 的計(jì)算機(jī)仿真》，北京郵電大學(xué)出版社， 。希望經(jīng)過這次設(shè)計(jì)，自己能夠取長補(bǔ)短，為以后的學(xué)習(xí)和工作積累寶貴的經(jīng)驗(yàn)。在實(shí)際的仿真過程當(dāng)中，第一次得到的結(jié)果往往與理論值有很大的誤差，經(jīng)過多次更改程序可以縮小誤差使得仿真結(jié)果更貼近理論值。經(jīng)過一系列前期準(zhǔn)備，在搜集了更多相關(guān)知識之后，結(jié)合所學(xué)利用 MATLAB 軟件編程仿真得到仿真圖。然后再次逐個進(jìn)行講解討論，通過 MATLAB 仿真，對連續(xù)信道的模型有了更加清楚地了解和認(rèn)識。通常可用信道的時延擴(kuò)展 來表示信道的多徑擴(kuò)展情況，多徑時延擴(kuò)展的倒數(shù)稱為信道的相干帶寬m21???B，設(shè)輸入信號的碼元間隔為sT ，當(dāng) Bs 》 1 時，信號的衰減是平坦的；反之，信號的衰減是頻率選擇性的（需要強(qiáng)調(diào)：多徑非時變信道時恒參信道，在輸出端不存在衰落現(xiàn)象）。當(dāng)輸入信號的帶寬遠(yuǎn)小 于信道帶寬時（第一個零點(diǎn)帶寬），則信道對輸入信號的所有頻率分量的衰減幾乎相同，這種情況下，信號經(jīng)歷平坦性衰減，如圖齊魯工業(yè)大學(xué) 2020 屆本科畢業(yè)生畢業(yè)設(shè)計(jì)（論文） 46 所示。 （ 1）用 MATLAB 畫出信道的幅頻響應(yīng)特性和相頻響應(yīng)特性； （ 2）設(shè)信道輸入信號為)()( sn n nTtgatb ?? ? ，其中 ? T sttg ??? 010)( 其他 ，sT=1，畫出輸出信號波形； （ 3）同（ 2）相同形式的輸入信號，但sT =8，畫出輸出信號波形。當(dāng)改齊魯工業(yè)大學(xué) 2020 屆本科畢業(yè)生畢業(yè)設(shè)計(jì)（論文） 變多徑數(shù)為 2，且時延變化頻率參數(shù) taof=100 時，此時相對于快衰落情況，接收包絡(luò)快速起伏，但由于徑數(shù)不多，接收信號包絡(luò)起伏具有明顯的周期性，見下圖 44。當(dāng) taof= 時，輸出信號額波形圖如圖 43 所示，圖中示意了長時觀察和 短時觀察時輸出信號的包絡(luò)，（ a）圖時間范圍為 0~10s，（ b）圖時間范圍為 0~1000s。 齊魯工業(yè)大學(xué) 2020 屆本科畢業(yè)生畢業(yè)設(shè)計(jì)（論文） 0 0 . 2 0 . 4 0 . 6 0 . 8 1 1 . 2 1 . 4 1 . 6 1 . 8 2101ts(t)輸入單頻信號0 1 2 3 4 5 6 7 8 9 1021012ts(t)經(jīng)過 20 徑后接收信號 2 0 1 5 1 0 5 0 5 10 15 20012345fs(f)輸入單頻信號頻譜 2 0 1 5 1 0 5 0 5 10 15 20012345fr(f)多徑信道輸出信號頻譜 42 時變多徑信道輸出信號頻譜 從圖中可以看到，單頻信號經(jīng)過 20 徑時變信道后，輸出信號的包絡(luò)隨時間隨機(jī)起伏，輸出信號的頻譜從沖激譜變成一個窄帶頻譜。 多徑信道模型仿真 單頻信號經(jīng)過多徑時變信道 以下示意了一個幅度為 頻率為 10HZ 的單頻信號經(jīng)過 20 條路徑傳輸?shù)玫降牟ㄐ渭捌漕l譜，這 20 條路徑的衰減相同，但時延的大小是 隨時間變化的，每徑時延的變化規(guī)律為正弦型，變化的頻率從 0~2Hz 隨機(jī)均勻抽取?？梢钥吹?，數(shù)字信號無論經(jīng)過幅頻失真還是相頻失真的信道時，對信號接收都有影響。 齊魯工業(yè)大學(xué) 2020 屆本科畢業(yè)生畢業(yè)設(shè)計(jì)（論文） （ 1）fjefH ???)( （無失真信道） （ 2）fjef ffH ??? ?? sin)( （幅度失真信道） （ 3）????? ??? ??? 01)fjππ(e 01)fjππ(eH(f) （相位失真信道） （ 4）)1( 2)( ???? ffjefH ? 編寫主程序，在 MATLAB 窗口中生成 m 文件，運(yùn)行程序得到數(shù)字信號經(jīng)過不同失真信道時的輸出信號畸變情況，如圖 41 所示。除內(nèi)部函數(shù)以外，所有 MATLAB 的核心文件和工具箱文件都是可讀可改的源文件。 MATLAB 包含兩個部分：核心部分和各種可選工具箱。由于MATLAB 的程序不用編譯等預(yù)處理，也不生成可執(zhí)行文件，程序?yàn)榻忉寛?zhí)行，所以速度較慢。 MATLAB 還具有較強(qiáng)的編輯圖形界面的能力。 程序的可移植性很好，基本上不做修改就可以在各種型號的計(jì)算機(jī)和操作系統(tǒng)上運(yùn)行。 MATLAB 既具有結(jié)構(gòu)化的控制語