通信系統(tǒng)課群綜合訓(xùn)練與設(shè)計說明書

word格式完美整理/課程設(shè)計任務(wù)書題目:通信系統(tǒng)課群綜合訓(xùn)練與設(shè)計初始條件：MATLAB軟件,電腦,通信原理知識要求完成的主要任務(wù):1、利用仿真軟件〔如Matlab或SystemView,或硬件實驗系統(tǒng)平臺上設(shè)計完成一個典型的通信系統(tǒng)2、學(xué)生要完成整個系統(tǒng)各環(huán)節(jié)以及整個系統(tǒng)的仿真,最終在接收端或者精確或者近似地再現(xiàn)輸入〔信源,計算失真度,并且分析原因。時間安排：序號設(shè)計內(nèi)容所用時間1根據(jù)設(shè)計任務(wù),分析電路原理,確定實驗方案2天2根據(jù)實驗條件進行電路的測試,并對結(jié)果進行分析7天3撰寫課程設(shè)計報告1天合計2周指導(dǎo)教師簽名：2013年1月11日系主任〔或責(zé)任教師簽名：2013年1月11日摘要在當(dāng)今信息化社會中,隨著通信技術(shù)的飛速發(fā)展,相關(guān)的通信理論、技術(shù)也得到了飛速發(fā)展。數(shù)字通信已成為信息傳輸?shù)闹匾侄?數(shù)字通信的新設(shè)備不斷涌現(xiàn),全球數(shù)字化已成為當(dāng)今社會發(fā)展的主要潮流,人們越來越離不開數(shù)字通信,越來越期望了解和掌握數(shù)字通信技術(shù)。數(shù)字通信是迅速發(fā)展的信息技術(shù),目前將科研的發(fā)展前沿與基礎(chǔ)的通信理論有機結(jié)合闡述的著作還很少。Matlab是MathWorks公司推出的一套高性能的數(shù)值計算和可視化軟件,其強大的計算和圖形功能使其在科學(xué)計算和工程領(lǐng)域贏得了眾多的用戶。經(jīng)過許多專家、工程師在自己相關(guān)領(lǐng)域的擴充,Matlab成為了一個多領(lǐng)域、多學(xué)科、多功能的優(yōu)秀科技應(yīng)用軟件,從而被廣泛地應(yīng)用于各領(lǐng)域和學(xué)科的研究與仿真。Matlab是一種編程語言和可視化工具,它能對數(shù)據(jù)以圖形的方式顯示出來,使數(shù)據(jù)間的關(guān)系明了Matlab所包含的應(yīng)用工具箱的功能非常豐富,提供的Simulink是一個用來對動態(tài)系統(tǒng)進行建模、仿真和分析的軟件包,支持線性和非線性系統(tǒng),能夠在連續(xù)時間域、離散時間域或者兩者的混合時間域里進行建模,其中通信系統(tǒng)工具箱中包含了對通信系統(tǒng)進行分析和仿真所需的信源編碼、糾錯編碼、信道、調(diào)制解調(diào)以及其他所用的庫函數(shù)和模塊。關(guān)鍵詞：增量調(diào)制,Miller碼,循環(huán)碼,PSK,衰落信道AbstractIntoday'sinformationsociety,withtherapiddevelopmentofcommunicationtechnology,communicationtheory,technologyhasalsobeenarapiddevelopment.Digitalcommunicationhasbecomeanimportantmeansofinformationtransmission,newdigitalcommunicationdevicesemergingglobaldigitaltechnologyhasbecomethedevelopmenttrendintoday'ssociety,peopleincreasinglyinseparablefromdigitalcommunicationsincreasinglyexpectedtounderstandandmasterdigitalcommunicationtechnology.Digitalcommunicationsistherapiddevelopmentofinformationtechnology,thedevelopmentofcutting-edgeresearchbasedcommunicationtheorycombineelaboratewritingsrarely.MatlabMathWorksInc.launchedasetofhigh-performancenumericalcomputationandvisualizationsoftware,powerfulcomputingandgraphicscapabilitiesmakeitwonalargenumberofusersinthefieldofscientificcomputingandengineering.Manyexperts,engineersintheirownfieldexpansion,Matlabhasbecomeamulti-field,multi-disciplinary,multi-functionalexcellenceintechnologyandapplicationsoftware,whichiswidelyusedinvariousfieldsanddisciplinesandsimulation.Matlabisaprogramminglanguageandvisualizationtools,datacanbedisplayedingraphicalformtherelationshipbetweenthedatatounderstandtheMatlabtheapplicationtoolboxfunctionscontainedveryrich,providedSimulinkisusedtodynamicallysystemmodeling,simulationandanalysissoftwarepackagethatsupportslinearandnonlinearsystems,continuous-timedomain,thediscrete-timedomain,oramixofbothtimedomainmodeling,communicationsystemtoolboxcontainscommunicationsystemanalysisandsimulationrequiredsourcecoding,errorcorrectioncoding,channel,modulationanddemodulation,andotherlibraryfunctionsandmodules.Keywords:deltamodulation,Millercodes,cycliccodes,PSKfadingchannel1.設(shè)計任務(wù)利用Matlab設(shè)計完成一個典型的通信系統(tǒng)〔如圖1所示。圖1典型的通信系統(tǒng)其中,信源為正弦信號,數(shù)字化方式為增量調(diào)制,基帶碼為Miller碼,信道碼為循環(huán)碼,調(diào)制方式為PSK,信道類型為衰落信道。2.實驗電路原理分析2.1增量調(diào)制原理增量調(diào)制是由PCM發(fā)展而來的模擬信號數(shù)字化的一種編碼方式,它是PCM的一種特例。增量調(diào)制編碼基本原理是指用一位編碼,這一位碼不是表示信號抽樣值的大小,而是表示抽樣幅度的增量特性,即采用一位二進制數(shù)碼"1"或"0"來表示信號在抽樣時刻的值相對于前一個抽樣時刻的值是增大還是減小,增大則輸出"1"碼,減小則輸出"0"碼。輸出的"1","0"只是表示信號相對于前一個時刻的增減,不表示信號的絕對值。增量調(diào)制最主要的特點就是它所產(chǎn)生的二進制代碼表示模擬信號前后兩個抽樣值的差別<增加、還是減少>而不是代表抽樣值本身的大小,因此把它稱為增量調(diào)制。在增量調(diào)制系統(tǒng)的發(fā)端調(diào)制后的二進制代碼1和0只表示信號這一個抽樣時刻相對于前一個抽樣時刻是增加<用1碼>還是減少<用0碼>。收端譯碼器每收到一個1碼,譯碼器的輸出相對于前一個時刻的值上升一個量化階,而收到一個0碼,譯碼器的輸出相對于前一個時刻的值下降一個量化階。增量調(diào)制<DM>是DPCM的一種簡化形式。在增量調(diào)制方式下,采用1比特量化器,即用1位二進制碼傳輸樣值的增量信息,預(yù)測器是一個單位延遲器,延遲一個采樣時間間隔。預(yù)測濾波器的分子系數(shù)向量是[0,1],分母系數(shù)為1。當(dāng)前樣值與預(yù)測器輸出的前一樣值相比較,如果其差值大于零,則發(fā)1碼,如果小于零則發(fā)0碼。增量調(diào)制系統(tǒng)框圖如圖2所示,其中量化器是一個零值比較器,根據(jù)輸入的電平極性,輸出為δ,預(yù)測器是一個單位延遲器,其輸出為前一個采樣時刻的解碼樣值,編碼器也是一個零值比較器,若其輸入為負(fù)值,則編碼輸出為0,否則輸出為1。解碼器將輸入1,0符號轉(zhuǎn)換為δ,然后與預(yù)測值相加后得出解碼樣值輸出,同時也作為預(yù)測器的輸入。解碼編碼二電平量化輸入樣值enen=δsgn<en>傳輸δsgn<en>xn=xn-1+δsgn<en>解碼編碼二電平量化xn+-++單位延遲預(yù)測輸出+n-1單位延遲單位延遲+預(yù)測輸出解碼樣值輸出單位延遲xn-1預(yù)測輸入xn=xn-1+δsgn<en>圖2增量調(diào)制原理框圖2.2Miller碼原理Miller碼也稱延遲調(diào)制碼,是一種變形雙向碼。其編碼規(guī)則：對原始符號"1"碼元起始不躍變,中心點出現(xiàn)躍變來表示,即用10或01表示。對原始符號"0"則分成單個"0"還是連續(xù)"0"予以不同處理；單個"0"時,保持0前的電平不變,即在碼元邊界處電平不躍變,在碼元中間點電平也不躍變；對于連續(xù)"0",則使連續(xù)兩個"0"的邊界處發(fā)生電平躍變。2.3循環(huán)碼原理圖3總原理方框圖循環(huán)碼的編、譯碼系統(tǒng)由下列五部分組成：定時單元、信碼發(fā)生器及顯示部分、編碼器和譯碼器?！?定時單元本單元提供編碼器及譯碼器所需的時序信號。其時鐘重復(fù)頻率〔CP為2MHZ?！?信碼發(fā)生器本單元給編碼器提供一個信號源,手控開關(guān)〔板上CDIN置于+5V時,發(fā)光二極管亮,代表輸出"1"碼元。若開關(guān)置于"0",代表輸出"0"碼元?！?循環(huán)碼編碼器編碼器是本實驗的主要部分。根據(jù)生成多項式,采用5個異或門和D觸發(fā)器組成編碼器。在K1信號的控制下,輸入6位信息碼元CDIN,一方面串行輸入信道,另一方面通過與門送入除法電路進行計算。第6位輸入碼元結(jié)束時,K1信號也為零,在CP脈沖作用下,移位寄存器將計算的結(jié)果〔CDOUT送往信道,即在6位信息碼元后附加了9位監(jiān)督碼,使碼長〔n=K+r為15?！?譯碼方框圖及原理介紹圖4譯碼器方框圖經(jīng)過信道加錯后的信碼,在K1信號的作用下,進入6位移存器,同時另一路進入除法電路進行伴隨式計算,當(dāng)6個信碼全存人移存器時,電子開關(guān)置于"0",此時信碼保存在移存器中,同時另一路已進入除法器的信碼,在CP脈沖的作用下,進入除法電路及正交方程形成網(wǎng)絡(luò)、大數(shù)邏輯判決電路。由于本實驗最小碼間距離d0=6,故最多能糾正兩位錯碼,若錯碼個數(shù)在2個以內(nèi),該系統(tǒng)能自行糾正,糾正后的信碼通過電子開關(guān)進入移存器,并在顯示信號K3的作用下,若發(fā)光二極管亮表示"1"碼,不亮表示"0"碼。此時譯碼信號是并行輸出至顯示部分的,它顯示的信號應(yīng)與信源顯示的一一對應(yīng)〔注意此時信道干擾產(chǎn)生的錯碼只能是1個或2個。假如信道中錯誤個數(shù)已超過該碼糾錯能力〔即超過2,那么譯碼顯示與編碼顯示不能對應(yīng)。正交方程的定義是：假定最高錯誤碼元為e14,其次e13,此類推至e0,即e14在每個方程中均出現(xiàn)一次,而其它錯誤碼元在4個方程組中出現(xiàn)一次,正交方程組如下：由正交形成網(wǎng)絡(luò)輸出至大數(shù)邏輯判決電路,由四個三與非門及四輸入與非門組成。該電路輸出信號通過與門在K2信號作用下,一方面進入除法電路進行伴隨式復(fù)位,另一方面進入異或門,與6位移位寄存器中相應(yīng)的信號相加〔最高位。已糾正的碼經(jīng)過與門又存在相應(yīng)的移位寄存器中。6位移位寄存器將糾正后的正確碼字在顯示脈沖K3作用下并行輸出。在2個錯誤碼元以內(nèi)本電路能自行校正,譯碼器顯示的碼與編碼器顯示的碼應(yīng)一致。如果錯誤碼元超過2個,譯碼電路會產(chǎn)生錯糾現(xiàn)象,譯碼顯示電路顯示的碼與編碼器不一樣,即使一樣也是巧合。2.4PSK原理相位偏移調(diào)制,又稱移相鍵控〔PSK,PhaseShiftKeying是一種利用相位差異的信號來傳送資料的調(diào)制方式。該傳送信號必須為正交信號,其基底更須為單位化信號。一個信號所代表的數(shù)學(xué)公式一般調(diào)制信號的改變部份可分為幅度A〔ASK用、相位〔PSK用及頻率〔FSK用三種。其中PSK即利用相位差異來產(chǎn)生的調(diào)制方式。MPSK通用的傳輸符號之公式。PSK又可稱M-PSK或MPSK,目前有BPSK、QPSK、16PSK、64PSK等等,常用的只有QPSK。而M是代表傳送信號的符號〔symbol種類。符號越多,傳送的比特數(shù)越多,自然在固定時間可傳送越多的資料量〔bps。傳輸量公式。圖5BPSK、QPSK、8PSK及16PSK的BER對SNR圖假設(shè)各MPSK皆在同一能量下傳送,PSK會因為符號種類〔M的提升使比特錯誤率〔BitsErrorRate,BER快速上升。所以在符號數(shù)M大于16后都由QAM來執(zhí)行調(diào)制工作。QPSK如果用格雷碼對映的方式,其BER會和BPSK一樣。所以目前常用的只有QPSK。2.5衰落信道原理瑞利衰落信道〔Rayleighfadingchannel是一種無線電信號傳播環(huán)境的"統(tǒng)計模型<statisticalmodel>"。這種模型假設(shè)信號通過無線信道之后,其信號幅度<amplitude>是隨機的,即"衰落<fading>",并且其包絡(luò)<envelope>服從瑞利分布〔Rayleighdistribution。這一信道模型能夠描述由電離層和對流層反射的短波信道,以及建筑物密集的城市環(huán)境。瑞利衰落只適用于從發(fā)射機到接收機不存在直射信號〔LoS,LineofSight的情況,否則應(yīng)使用萊斯衰落信道〔Riceanfadingchannel作為信道模型。瑞利衰落能有效描述存在能夠大量散射無線電信號的障礙物的無線傳播環(huán)境。若傳播環(huán)境中存在足夠多的散射,則沖激信號到達接收機后表現(xiàn)為大量統(tǒng)計獨立的隨機變量的疊加,根據(jù)中心極限定理,則這一無線信道的沖激響應(yīng)將是一個高斯過程。如果這一散射信道中不存在主要的信號分量,通常這一條件是指不存在直射信號〔LoS,則這一過程的均值為0,且相位服從0到2π的均勻分布。即,信道響應(yīng)的包絡(luò)服從瑞利分布。設(shè)隨機變量R,于是其概率密度函數(shù)為：其中。若信道中存在一主要分量,例如直射信號〔LoS,則信道響應(yīng)的包絡(luò)服從萊斯分布,對應(yīng)的信道模型為萊斯衰落信道。通常將信道增益以等效基帶信號表示,即用一復(fù)數(shù)表示信道的幅度和相位特性。由此瑞利衰落即可由這一復(fù)數(shù)表示,它的實部和虛部服從于零均值的獨立同分布高斯過程。圖6最大多普勒頻移為10Hz的瑞利衰落信道。圖7最大多普勒頻移為100Hz的瑞利衰落信道。圖8瑞利衰落信道的自相關(guān)函數(shù),其多普勒頻移為10Hz。瑞利衰落模型適用于描述建筑物密集的城鎮(zhèn)中心地帶的無線信道。密集的建筑和其他物體使得無線設(shè)備的發(fā)射機和接收機之間沒有直射路徑,而且使得無線信號被衰減、反射、折射、衍射。在曼哈頓的實驗證明,當(dāng)?shù)氐臒o線信道環(huán)境確實接近于瑞利衰落。通過電離層和對流層反射的無線電信道也可以用瑞利衰落來描述,因為大氣中存在的各種粒子能夠?qū)o線信號大量散射。瑞利衰落屬于小尺度的衰落效應(yīng),它總是疊加于如陰影、衰減等大尺度衰落效應(yīng)上。信道衰落的快慢與發(fā)射端和接收端的相對運動速度的大小有關(guān)。相對運動導(dǎo)致接收信號的多普勒頻移。圖中所示即為一固定信號通過單徑的瑞利衰落信道后,在1秒內(nèi)的能量波動,這一瑞利衰落信道的多普勒頻移最大分別為10Hz和100Hz,在GSM1800MHz的載波頻率上,其相應(yīng)的移動速度分別為約6千米每小時和60千米每小時。特別需要注意的是信號的"深衰落"現(xiàn)象,此時信號能量的衰減達到數(shù)千倍,即30～40分貝。3.實驗過程及結(jié)果分析圖9增量調(diào)制圖10密勒碼圖11循環(huán)碼圖12PSK圖13瑞利衰落信道4.總結(jié)在本次課程設(shè)計的學(xué)習(xí)中,讓我受益頗多。一、讓我養(yǎng)成了課前預(yù)習(xí)的好習(xí)慣。一直以來就沒能養(yǎng)成課前預(yù)習(xí)的好習(xí)慣〔雖然一直認(rèn)為課前預(yù)習(xí)是很重要的,但經(jīng)過這一周,讓我深深的懂得課前預(yù)習(xí)的重要。只有在課前進行了認(rèn)真的預(yù)習(xí),才能在課上更好的學(xué)習(xí),收獲的更多、掌握的更多。二、培養(yǎng)了我的動手能力。"實驗就是為了讓你動手做,去探索一些你未知的或是你尚不是深刻理解的東西。"每個步驟我都親自去做,不放棄每次鍛煉的機會。經(jīng)過這兩周,讓我的動手能力有了明顯的提高。三、讓我在探索中求得真知。那些偉大的科學(xué)家之所以偉大就是他們利用實驗證明了他們的偉大。實驗是檢驗理論正確與否的試金石。為了要使你的理論被人接受,你必須用事實〔實驗來證明,讓那些懷疑的人啞口無言。雖說我們的通信原理實驗只是對前人的經(jīng)典實驗的重復(fù),但是對于一個知識尚淺、探索能力還不夠的人來說,這些探索也非一件易事。通信原理實驗都是一些經(jīng)典的給人類帶來了難以想象的便利與財富。對于這些實驗,我在探索中學(xué)習(xí)、在模仿中理解、在實踐中掌握。通信原理實驗讓我慢慢開始"摸著石頭過河"。學(xué)習(xí)就是為了能自我學(xué)習(xí),這正是實驗課的核心,它讓我在探索、自我學(xué)習(xí)中獲得知識。四、教會了我處理數(shù)據(jù)的能力。實驗就有數(shù)據(jù),有數(shù)據(jù)就得處理,這些數(shù)據(jù)處理的是否得當(dāng)將直接影響你的實驗成功與否。經(jīng)過這一周,我學(xué)會了圖像法等處理數(shù)據(jù)的方法,讓我對其它課程的學(xué)習(xí)也是得心應(yīng)手。

經(jīng)過這次課程設(shè)計的學(xué)習(xí),讓我收獲多多。但在這中間,我也發(fā)現(xiàn)了我存在的很多不足。我的動手能力還不夠強,當(dāng)有些實驗需要很強的動手能力時我還不能從容應(yīng)對；我的探索方式還有待改善,當(dāng)面對一些復(fù)雜的實驗時我還不能很快很好的完成；我的數(shù)據(jù)處理能力還得提高,當(dāng)眼前擺著一大堆復(fù)雜數(shù)據(jù)時我處理的方式及能力還不足,不能用最佳的處理手段使實驗誤差減小到最小程度

總之,課程設(shè)計讓我收獲頗豐,同時也讓我發(fā)現(xiàn)了自身的不足。在實驗課上學(xué)得的,我將發(fā)揮到其它中去,也將在今后的學(xué)習(xí)和工作中不斷提高、完善；在此間發(fā)現(xiàn)的不足,我將努力改善,通過學(xué)習(xí)、實踐等方式不斷提高,克服那些不應(yīng)成為學(xué)習(xí)、獲得知識的障礙。在今后的學(xué)習(xí)、工作中有更大的收獲,在不斷地探索中、在無私的學(xué)習(xí)、奉獻中實現(xiàn)自己的人身價值！5.參考文獻[1]劉泉.通信電子線路.XX理工大學(xué)出版社,2004.[2]樊昌信.通信原理.國防工業(yè)出版社,2006.[3]吳大正.信號與線性系統(tǒng)分析.XX電子科技大學(xué)出版社,1998[4]陳懷琛.MATLAB及在電子信息課程中的應(yīng)用.電子工業(yè)出版社,2006.[5]鄭君里.信號與系統(tǒng).高等教育出版社,2000.6.附錄functiony=DMcoding<>Ts=1e-3;t=0:Ts:20*Ts;x=sin<2*pi*50*t>;delta=0.4;%量化階距D<1+length<t>>=0;%預(yù)測器初始狀態(tài)fork=1:length<t>e<k>=x<k>-D<k>;%誤差信號e_q<k>=delta*<2*<e<k>>=0>-1>;%量化器輸出D<k+1>=e_q<k>+D<k>;%延遲器狀態(tài)更新y<k>=<e_q<k>>0>;%編碼輸出endsubplot<2,1,1>;plot<t,x,'-o'>;axis<[020*Ts-22]>;subplot<2,1,2>;stairs<t,y>;axis<[020*Ts-22]>;functionMiller=Millercoding<x>m=length<x>;Miller=zeros<1,2*m>;f=0;d=1;fori=1:mifx<i>==1Miller<2*i-1>=f;Miller<2*i>=not<f>;f=not<f>;d=1;elseifd==1Miller<2*i-1>=f;Miller<2*i>=f;d=0;elsef=not<f>;Miller<2*i-1>=f;Miller<2*i>=f;d=0;endendendTs=1e-3;t=0:Ts:41*Ts;stairs<t,Miller>;axis<[041*Ts-22]>;functionT=Cycliccoding<Input>symsx%求出所有的生成多項式g=[10111];%選擇任意一個作為〔7,4循環(huán)碼的生成多項式r=4;%監(jiān)督為數(shù)m1=poly2sym<Input>;ifm1==0T=[0000000];elsem11=expand<x^r*m1>;%用x^r乘以m1,相當(dāng)于對m1進行左移r位的操作m2=sym2poly<m11>;%將多項式轉(zhuǎn)化為矩陣表示形式len=length<m2>;iflen<7fori=1:lenm2<8-i>=m2<len-i+1>;endfori=1:7-lenm2<i>=0;endend[Q,R]=deconv<m2,g>;%求m2除以g所得的余數(shù)R=abs<R>;fori=1:length<R>R<i>=mod<R<i>,2>;end%轉(zhuǎn)化為模2運算T=R+m2;%T為生成的一個循環(huán)碼字endTs=1e-3;t=0:Ts:44*Ts;stairs<t,T>;axis<[044*Ts-22]>;function[re]=PSKcoding<Input>cp=[];mod1=[];f=4*pi;t=0:2*pi/199:2*pi;forn=1:length<Input>ifInput<n>==0A=zeros<1,200>;%每個值200個點elseifInput<n>==1A=ones<1,200>;endcp=[cpA];%s<t>,碼元寬度200c=cos<f*t>;%載波信號mod1=[mod1c];%與s<t>等長的載波信號變?yōu)榫仃囆问絜ndfigure<1>;subplot<2,1,1>;plot<cp>;gridon;axis<[04000-22]>;title<'二進制信號序列'>;cm=[];mod=[];forn=1:length<Input>ifInput<n>==0B=ones<1,200>;c=cos<f*t>;%載波信號elseifInput<n>==1B=ones<1,200>;c=cos<f*t+pi>;%載波信號endcm=[cmB];%碼元寬度200mod=[modc];%與s<t>等長的載波信號endre=mod;%調(diào)制subplot<2,1,2>;plot<re>;gridon;axis<[04000-22]>;title<'2PSK調(diào)制信號'>;figure<2>;subplot<2,1,1>;plot<abs<fft<cp>>>;axis<[040000800]>;title<'原始信號頻譜'>;subplot<2,1,2>;plot<abs<fft<re>>>;axis<[040000800]>;title<'2PSK信號頻譜'>;function[y_out]=Rayleighchannel<SignalInput>LengthOfSignal=length<SignalInput>;%信號長度<最好大于兩倍fc>fm=51;%最大多普勒頻移fc=510;%載波頻率delay=[03171109173251];power=[0-1-9-10-15-20];%dBy_in=[zeros<1,delay<6>>SignalInput];%為時移補零y_out=zeros<1,LengthOfSignal>;%用于信號輸出fori=1:6f=1:2*fm-1;%通頻帶長度y=0.5./<<1-<<f-fm>/fm>.^2>.^<1/2>>/pi;%多普勒功率譜〔基帶Sf=zeros<1,LengthOfSignal>;Sf1=y;%多普勒濾波器的頻響Sf<fc-fm+1:fc+fm-1>=y;%基帶映射到載波頻率上x1=randn<1,LengthOfSignal>;x2=randn<1,LengthOfSignal>;nc=ifft<fft<x1+i*x2>.*sqrt<Sf>>;%同相分量x3=randn<1,LengthOfSignal>;x4=randn<1,LengthOfSignal>;ns=ifft<fft<x3+i*x4>.*sqrt<Sf>>;%正交分量r0=<real<nc>+j*real<ns>>;%瑞利信號r=abs<r0>;%瑞利信號幅值y_out=y_out+r.*y_in<delay<6>+1-delay<i>:delay<6>+LengthOfSignal-delay<i>>*10^<power<i>/20>;end;figure<3>;subplot<2,1,1>;plot<SignalInput<delay<6>+1:LengthOfSignal>>;%去除時延造成的空白信號title<'SignalInput'>;subplot<2,1,2>;plot<y_out<delay<6>+1:LengthOfSignal>>;title<'SignalOutput'>;%figure<4>;%subplot<2,1,1>;%hist<r,256>;%title<'AmplitudeDistributionOfRayleighSignal'>；瑞利信號的振幅%subplot<2,1,2>;%hist<angle<r0>>;%title<'AngleDistributionOfRayleighSignal'>;瑞利信號的相位%figure<5>;%plot<Sf1>;%title<'TheFrequencyResponseofDopplerFilter'>;多普勒濾波器的頻率響應(yīng)function[output]=PSKdecoding<Input>t=0:2*pi/199:2*pi;f=4*pi;mod1=[];forn=1:length<Input>/200c=cos<f*t>;%載波信號mod1=[mod1c];%與s<t等長的載波信號變?yōu)榫仃囆问絜ndre=mod1.*Input;%同步解調(diào)figure<4>;subplot<2,1,1>;plot<re>;gridon;axis<[04000-22]>;title<'相乘后信號波形'>subplot<2,1,2>;plot<abs<fft<re>>>;axis<[040000400]>;title<'相乘后信號頻譜'>;%低通信號濾波器fp=500;fs=700;rp=3;rs=20;fn=11025;ws=fs/<fn/2>;wp=fp/<fn/2>;%計算歸一化角頻率[n,wn]=buttord<wp,ws,rp,rs>;%計算階數(shù)和截止頻率[b,a]=butter<n,wn>;%計算H〔z>figure<5>;freqz<b,a,1000,11025>;subplot<2,1,1>;axis<[04000-1003]>;title<'LPF幅頻相頻圖'>;jt=filter<b,a,re>;figure<6>;subplot<2,1,1>;plot<jt>;gridonaxis<[04000-22]>;title<'經(jīng)低通濾波器后的信號波形'>subplot<2,1,2>;plot<abs<fft<jt>>>;axis<[040000400]>;title<'經(jīng)低通濾波器后的信號頻譜'>;%抽樣判決form=1:length<Input>ifjt<m><0jt<m>=1;ifmod<m,200>==0output<fix<m/200>+1>=