不同方式解調下FSK、PSK、DPSK的誤碼率比較

1、摘要一般來說，數(shù)字解調與模擬調制的基本原理相同，但是數(shù)字信號有離散取值的特點。因此數(shù)字調制技術有兩種方法：1利用模擬調制的方 法去實現(xiàn)數(shù)字式調制，即數(shù)字調制看成是模擬調制的一個特例，把數(shù)字 基帶信號當做模擬信號的特殊情況處理；2利用數(shù)字信號的離散取值特 點通過開關鍵控載波，從而實現(xiàn)數(shù)字調制。這種方法通常稱為鍵控法， 比如對載波的振幅，頻率和相位進行鍵控，便可獲得振幅鍵控(Ampolitude shift keying, ASK )、頻移鍵控( Frequency shift keying,FSK )、和相移鍵控(Phase shift keying， PSK 三種基本的數(shù)字調制方式。這次為期一

2、周的通信傳輸課程設計的實習，就是通過MATLA編程仿真，來更好的理解FSK PSK DPSK的調制和解調過程。在 這次的實習中，主要是應用MATLAB進行編程仿真并顯示結果。仿真的是FSK的相干、非相干和過零解調，PSK的相干解調及DPSK的相干和差 分解調。并比較相同調制后的信號不同的解調方式和不同調制后的信號 相同的解調方式。關鍵字：頻移鍵控，相移鍵控，誤碼率，信噪比ABSTRACTGenerally speaking, digital demodulation and analog modulation of the basic principles of the same, but t

3、he digital signal has the characteristics of discrete values. Therefore, digital modulation techniques, there are two methods: one using analog modulation methods to achieve digital modulation, digital modulation that is seen as a special case of analog modulation to digital baseband signal as an an

4、alog signal a special case; 2 using digital signal characteristics of discrete values by keying switch carrier in order to achieve digital modulation. This method is usually referred to as keying method, such as the right carrier amplitude, frequency and the phase shift keying, can receive amplitude

5、 shift keying (Ampolitude shift keying, ASK), frequency shift keying (Frequency shift keying, FSK), and the phase shift keying (Phase shift keying, PSK) digital modulation of three basic ways. The week-long internship curriculum design communication transmission is through MATLAB simulation program,

6、 to a better understanding of FSK, PSK, DPSK modulation and demodulation process. In this attachment, the main is the application programming MATLAB simulation and displays the results. Simulation is the FSK coherent, non-coherent and zero-crossing demodulation, PSK and DPSK coherent demodulation of

7、 coherent and differential demodulation. And compare the same modulated signal demodulation in different ways and different modulated signal demodulation the same way.Keywords: Frequency shift keying, phase shift keying, the bit error rate, signal to noise ratio3目錄摘要 .2ABSTRACT 3緒論 .5第一章 FSK .61. FS

8、K產生原理2. FSK信號的解調3. FSK相關仿真第二章 PSK .81. PSK產生原理2. PSK信號的與解調3. PSK相關仿真第三章DPSK 101. DPSK產生的原理2. DPSK信號的解調3. DPSK相關仿真第四章FSK、PSKDPSK誤碼率的比較 12第五章 結論和心得體會 .15參考文獻 .16附錄 .175緒論數(shù)字信號的傳輸方式分為基帶傳輸和帶通傳輸然而,實際中的大多數(shù)信號(如無線信道)因具有帶通特征而不能直接傳送基帶信號，然而用數(shù)字基帶信號往往具有 豐富的低頻分量為了使數(shù)字信號在通信中傳輸，必須用數(shù)字基帶信號對載波進行調 制,以使信號與信道當然特征相匹配這種用數(shù)字基帶

9、信號控制載波，把數(shù)字基帶信 號轉變?yōu)閿?shù)字帶通信號(已調制號)的過程為數(shù)字調制(digital modulation).在接收端通過解調器把帶通信號還原成數(shù)字基帶信號的過程稱為數(shù)字解調(digitaldemodulation).通常把包括調制和人解調過程的數(shù)字傳輸系統(tǒng)叫做數(shù)字帶通傳輸系 統(tǒng)。為了與“基帶” 一詞相對應，帶通傳輸也也稱為頻率傳輸，又因為是借助于正 弦載波的幅度、頻率和相位來傳遞數(shù)字基帶信號的，所以帶通傳輸也叫載波傳輸。 因此，熟知數(shù)字調制、解調的過程是非常重要的。為了更好的理解FSK PSK DPSK的調制和解調過程，此次試驗選擇了通過 MATLAB仿真FSK PSK DPSK勺調

10、制、 解調過程，并通過實驗來比較它們的誤碼率和信噪比的關系。本次課程設計，我們小組采用分工合作的方式完成。關于課題的選擇是由我們三人共同商討后確定的。在程序設計中，唐一文同學負責FSK的調制相關和非相關解調，陳婭靚同學負責FSK調制的過零解調和PSK的調制相關解調，劉榮華同學負 責DPSK的調制相關和差分解調，并且各自負責相應設計報告的編寫。 最后大家一起 共同完成設計報告的整理。第一章FSK1.2FSK產生原理.頻移鍵控是利用載波的頻率變化來傳遞數(shù)字信息。在2FSK中，載波的頻率隨二 進制信號在f1和f2倆個頻點間變化，也就是說，一個2FSK信號可以看成是兩個不 同載頻的2ASK信號的疊加，

11、因此，2FSK信號的表達式為：S2FSK (t) = S(t)cos(ctnTsJCos代)1Cos(E 2舛 g(it) nTs) cosn=A cos國 0t +。tS 2FSK =' a n g (t - nT s ) COS ，1 t ' a n g (t - nTn66s(t1tn) S(n) S發(fā))q”s( 2t=2兀fcos ；0t 'nF0，cosC 1t1,2 _發(fā)PP“2時f2,an 為an1,n)an的反碼62FSK信號的產生方法主要有兩種。 一種可以采用模擬調頻電路來實現(xiàn)； 另一種 可以采用鍵控法來實現(xiàn)，即在二進制基帶矩形脈沖序列的控制下通過開關

12、電路對兩 個不同的獨立頻率源進行選通，使其在每一個碼元Ts期間輸出fl或f2兩個載波之 一。這兩種方法產生2FSK信號的差異在于：有調頻法產生的 2FSK信號在相鄰碼元 之間的相位是連續(xù)變化的，而鍵控法產生的2FSK信號，是由電子開關在兩個獨立的 頻率源之間轉換形成，故相鄰碼元之間的相位不一定連續(xù)。2.FSK信號的解調2FSK信號的常用解調方法是采用如同所示的非相干解調、相干解調和過零檢 測法。相干解調和非相干解調的原理是：將 2FSK信號分解為上下兩路2ASK信號分 別進行解調，然后進行抽樣判決。這里的抽樣判決時直接比較兩路信號抽樣值的大 小，可以不設置專門的門限。判決規(guī)則應與調制規(guī)則相呼應

13、， 調制時若規(guī)定“ 1”符 號對應載波頻率fl，則接收時上支路的樣值較大，應判為“ T;反之，判為“ 0”。非相干解調相干解調過零檢測法過零檢測的原理基于2FSK信號的過零點數(shù)隨不同頻率而異，通過過零點數(shù)目 的多少，從而區(qū)分兩個不同頻率的信號碼元。 如上圖中，2FSK信號經過限幅、微分、整流后形成與頻率變化相對應的尖脈沖序列，這些尖脈沖的密集程度反映了信號的 頻率高低，尖脈沖的個數(shù)就是信號過零點的數(shù)。把這些尖脈沖變換成較寬的矩形脈 沖，以增大其直流分量，該直流分量的大小和信號頻率的高低成正比。然后經過低 通濾波器取出此直流分量，這樣就完成了頻率一幅度轉換，從而根據(jù)直流分量幅度 上的區(qū)別還原出數(shù)

14、字信號“ 1”和“ 0”。3.FSK相關仿真FSK的相干和非相干解調、FSK的過零檢測流程圖如如所示9#第二章PSK1. PSK產生原理相移鍵控是利用載波的相位變換來傳遞數(shù)字信息，而振幅和頻率保持不變。在 2PSK中，通常用初始相位0和n分別表示二進制“ T和“ 0”。因此，2PSK信號的 時域表達式為：S B ( t)=na n g (t-nT s.)cosctan不同于ASK、FSKf+11 , PPan二I-11,1 - 1P-PCos 3 ct (0 相) 發(fā) 1 ”寸 二 eo (t) = y)一般為幅度是g(t) 般偽幅度是為T °1、寬度為矩形|的矩形脈發(fā)中0”時BPS

15、K (t)二 土二cOsPs0(!t = 士 cos(電t ctcos侔 $計“)£ 0 弋 0、只2PSK言號的調制原理如圖所示相乘法產生數(shù)字選相器載波發(fā)生e。( t)2.PSK信號的與解調2PSK信號的解調原理通常采用相干解調法，解調器原理如圖所示相干解調3.PSK相關仿真PSK的相干解調流程圖如圖所示:第三章DPSK1.DPSK產生的原理在PSK信號中，相位變換是以未調載波的相位變換作為參考基準的。 相干解調時， 由于載波回復中相位有0，模糊性，導致解調過程中出現(xiàn)“反向工作”現(xiàn)在，恢復 出的數(shù)字信號“ T, “0”倒置。為克服此缺點，提出了 2DPSK調制方式。2DPSK是利用

16、前后相鄰碼元的載波相對相位變換傳遞數(shù)字信息，又稱相對相移鍵控。假設 W為當前碼元與前一碼元的載波相位，可定義一種數(shù)字信息與屮之間的關系為： W = 0表示數(shù)字信息“ 0”180表示數(shù)值信息“ 1”數(shù)字信息與屮之間的關系也可以相反定義。于是可以將一組二進制數(shù)字信息與其對應的 2DPSK言息的載波相位關系示例如 下：二進制數(shù)字信息：1 1 0 1 0 0 1 1 02DPSK言號相位：(0) n 0 0 n n n 0 n n(n ) 0 n n 0 0 0 n 0 0也就是說，2DPSK信號的相位并不直接代表基帶信號，而前后碼元相對相位的 差才唯一決定信息符號。2DPSK信號的產生方法：先對二進

17、制數(shù)字基帶信號進行差分編碼，即轉換為相 對碼，然后再根據(jù)相對碼進行絕對調相，從而產生二進制差分相移鍵控。差分碼的編碼規(guī)則：bn = bn-1 二 an式中：二為模2加；bn-1為bn的前一碼元，最初的bn-1可任意設定。將產生的相 對碼進行絕對調相即產生二進制差分相移鍵2. DPSK信號的解調其解調原理是：對2DPSK言號進行相干解調，恢復出相對碼，在經碼反變換變 成絕對碼，從而恢復出發(fā)送的二進制信息。差分編碼的逆過程稱為差分譯碼,即：a n= b n-1 二 bnDPSK言號的解調分方式為兩種，一種為相干解調，一種為差分解調。相干解調：差分解調：用這種方法不需要專門的相干載波，只需有收到的2

18、DPSK言號延時一 個碼元間隔T,然后與2DPSK言號本身相乘。相乘器起著相位比較的作用，相乘結 果反映了前后碼元的相位差，經低通濾波器后再進行抽樣判決，即可直接恢復出原 始數(shù)字信息。3. DPSK相關仿真DPSK勺相干和差分解調流程如圖所示:第四章FSK、PSK DPSK誤碼率的比較誤碼率是指接收的碼元數(shù)在傳輸總碼元數(shù)中所占的比例，即:誤碼率錯誤碼元數(shù) 一傳輸總碼元數(shù)FSK相干解調:誤碼PeFSK非相干解調：誤碼率pe =亠-"22PSK相干解調：誤碼率 Pe = erfc ( , r2DPSKf干解調:誤碼率 Pe = erfc ( rDPSK1分非相干解調：誤碼率 pe =丄e

19、 "2r為信噪比。誤碼率是衡量一個數(shù)字通信系統(tǒng)性能的重要指標。 在信道高斯白噪聲的干擾下， 各種二進制數(shù)字調制系統(tǒng)的誤碼率取決于解調器輸入信噪比，而誤碼率表達式的形 式則取決于解調方式。對于所有的數(shù)字調制系統(tǒng)誤碼率與信噪比的關系的圖表來看，所有的曲線呈減 函數(shù)的下降曲線，即隨著信噪比的增大，誤碼率降低。橫向比較來看，對于同一種調制方式，當信噪比相同時，采用相干解調方式的 誤碼率低于非相干解調方式的誤碼率，即2FSK的相干解調的誤碼率低于非相干解調 的誤碼率，由于過零解調步驟復雜，2FSK非相干解調的誤碼率低于過零解調的誤碼 率；縱向比較來看，對2PSK,2DPSK,2FS三種調制方式

20、，若采用同一種解調方式（相 干解調或非相干解調），則2PSK的誤碼率最低，2DSPK勺誤碼率次之，2FSK的誤碼 率最高。當信噪比一定時，誤碼率由低到高依次是：2PSK的相干解調，2DPSK勺相干解 調，2DPSK勺差分解調，2FSK的相干解調，2FSK的非相干解調。FSK相干解調誤碼率：13相干解調后誤碼率統(tǒng)計0EbNo (dB)REB dna RESFSK非相干解調誤碼率:非相干解調后誤碼率統(tǒng)計REB dna BtsFSK過零解調的誤碼率:15BtB dna RES過零解調后誤碼率統(tǒng)計101010丨Theoretical selTheoretical BERSimulated SERSim

21、ulated BER:4t?Jcr3.5100.51.52.5EbNo (dB)4.5PSK相干解調誤碼率統(tǒng)計:PSK相干解調后誤碼率統(tǒng)計REB dna RLEDPSK相干解調誤碼率統(tǒng)計:17相干解調后誤碼率統(tǒng)計DPSK的差分解調誤碼率統(tǒng)計:REB dna RESDPSK差分解調后誤碼率統(tǒng)計REB dna RLE第五章結論和心得體會結論：衡量一個數(shù)字通信系統(tǒng)性能的好壞的指標有多種，但主要是誤碼率性能,頻帶利用率，對信道的適應能力和設配的的復雜程度等。在本次用MATLAB語言設計程序比較不同調制方式和相應的不同的解調方法下的誤碼率大小的比較課程設計 中，其中2DPSK，2PSK，2FSK三種調制

22、方式在高斯白噪聲的干擾下，2PSK性能最好，2DPSK次之，2FSK最差；對同一種調制方法，相干解調的誤碼率要比非相干解調的誤碼率低。心得體會：在為期一周的通信原理課程設計中 ,我們小組選擇了基帶信號調制 和解調后誤碼率的比較作為課題。在本次實踐中，我們小組三人分工明確，相互合 作。在第一天確定課程設計的題目之后，便去圖書館、網(wǎng)上查找相應的資料，但結 果讓人失望。書上和網(wǎng)上關于解調的 MATLAB 源代碼幾乎沒有。于是，在開始做 的時候，我們對照通信原理上的不同解調方式的模塊，自己寫程序，分步運行 仿真，然后再將各個模塊合在一起，形成一個完整的調制解調。讓我們意外的是， 在第四天的上午，所有的

23、程序都已經編寫出來，仿真結果也與理論相符。在這次的 課程設計中， 要感謝老師的幫助指導， 期間在編寫程序的時間， 遇到了很多的問題， 幾次都想放棄，但老師耐心、逐步的指導，最終讓我們小組完成了課程設計的全部 任務。同時，這次的課程設計對我們的影響也是很大的， 看似很難或是不懂的問題， 只要自己去動手做，一定是可以做出來的。這也讓我們知道了，在以后碰到了難題 的時候，與其逃避，不如直接面對，不管怎樣，只要去想，只要去做，結果總是能 出來的。此次為期一周的課程設計， 是感覺收獲最大的一次。 除了能完成所有的設計任 務，同時還熟悉了通信原理上相關的內容。希望在以后的學習中，還能有同樣 的機會，將我們

24、所學的東西運用在實際中。參考文獻【1】鄧華 .MATLAB 通信系統(tǒng)仿真詳解 .中國電力出版社 .2003. 【2】樊昌信.通信原理 （第六版）.國防工業(yè)出版社 .2008【3】現(xiàn)代通信系統(tǒng)（ MATLAB 版） .電子工業(yè)出版社 .2003【4】 J. D. Gibson. Principles of Digital and Analog Communication, Macmillan Publishing Company.New York.1989【5】薛定宇.基于MATLAB/Semolina的系統(tǒng)仿真技術與應用.清華大學出版 社.2003附：%FS碉制解調MATLA源代碼fun ct

25、io n FSKFc=10; % 載頻 Fs=100; % 系統(tǒng)采樣頻率Fd=1; % 碼速率 N=Fs/Fd;df=10;numSymb=25;進行仿真的信息代碼個數(shù) M=2;% 進制數(shù)SNRpBit=60;%言噪比 SNR=SNRpBit/log2(M);seed=12345 54321;numPlot=25;%產生 25 個二進制隨機碼 x=randsrc(numSymb,1,0:M-1);% 產生 25個二進制隨機碼 figure(1)stem(0:numPlot-1,x(1:numPlot),'bx');title(' 二進制隨機序列 ') xlabe

26、l('Time');ylabel('Amplitude'); %調制 y=dmod(x,Fc,Fd,Fs,'fsk',M,df);numModPlot=numPlot*Fs;t=0:numModPlot-1./Fs;figure(2) plot(t,y(1:length(t),'b-');axis(min(t) max(t) -1.5 1.5); title(' 調制后的信號 ') xlabel('Time');ylabel('Amplitude'); %在已調信號中加入高斯白噪聲

27、randn('state',seed(2);在已調信y=awgn(y,SNR-10*log10(0.5)-10*log10(N),'measured','dB');% 號中加入高斯白噪聲figure(3)plot(t,y(1:length(t),'b-');%畫出經過信道的實際信號axis(min(t) max(t) -1.5 1.5);title(' 加入高斯白噪聲后的已調信號 ') xlabel('Time');ylabel('Amplitude'); %相干解調z1=ddemo

28、d(y,Fc,Fd,Fs,'fsk',M,df);獅輸出波形的相干M元頻移鍵控解調figure(4) stem(0:numPlot-1,x(1:numPlot),'bx'); hold on;stem(0:numPlot-1,z1(1:numPlot),'ro');hold off;axis(0 numPlot -0.5 1.5);title(' 相干解調后的信號原序列比較 ')legend(' 原輸入二進制隨機序列 ',' 相干解調后的信號 ') xlabel('Time');yl

29、abel('Amplitude');%非相干解調z2=ddemod(y,Fc,Fd,Fs,'fsk/noncoh',M,df);獅輸出波形的非相干M元頻移鍵控解調figure(5) stem(0:numPlot-1,x(1:numPlot),'bx'); hold on;stem(0:numPlot-1,z2(1:numPlot),'ro');hold off;axis(0 numPlot -0.5 1.5);title(' 非相干解調后的信號 ')legend(' 原輸入二進制隨機序列 ',

30、9; 非相干解調后的信號 ') xlabel('Time');ylabel('Amplitude');%誤碼率統(tǒng)計errorSym ratioSym=symerr(x,z1);figure(6) simbasebandex(0:1:5); title(' 相干解調后誤碼率統(tǒng)計 ')errorSym ratioSym=symerr(x,z2);figure(7) simbasebandex(0:1:5);title(' 非相干解調后誤碼率統(tǒng)計 ')%FS調制過零解調MATLA源代碼 function FSKFc=10; %

31、載頻Fs=100; % 系統(tǒng)采樣頻率Fd=1; % 碼速率N=Fs/Fd;df=10;numSymb=25;進行仿真的信息代碼個數(shù) M=2;% 進制數(shù)SNRpBit=60;%信噪比 SNR=SNRpBit/log2(M);seed=12345 54321;numPlot=25;%產生 25 個二進制隨機碼x=randsrc(numSymb,1,0:M-1);% 產生 25個二進制隨機碼 figure(1)stem(0:numPlot-1,x(1:numPlot),'bx');title(' 二進制隨機序列 ') xlabel('Time'); y

32、label('Amplitude');%調制 y=dmod(x,Fc,Fd,Fs,'fsk',M,df); numModPlot=numPlot*Fs; t=0:numModPlot-1./Fs;figure(2) plot(t,y(1:length(t),'b-'); axis(min(t) max(t) -1.5 1.5); title(' 調制后的信號 ') xlabel('Time'); ylabel('Amplitude');%在已調信號中加入高斯白噪聲 randn('state&

33、#39;,seed(2);在已調信y=awgn(y,SNR-10*log10(0.5)-10*log10(N),'measured','dB');% 號中加入高斯白噪聲figure(3)plot(t,y(1:length(t),'b-');%畫出經過信道的實際信號axis(min(t) max(t) -1.5 1.5);title(' 加入高斯白噪聲后的已調信號 ') xlabel('Time');ylabel('Amplitude');%限幅for t=1:2500 if y(t)>0 y(

34、t)=1;elsey(t)=0;endend figure(4) stem(y) title(' 限幅后的信號 ') xlabel(' 點數(shù) '); ylabel('Amplitude');%差分for t=1:2499if abs(y(t+1)-y(t)=1y(t)=1;else y(t)=0;end end figure(5) stem(y) title(' 差分后的信號 ') xlabel(' 點數(shù) '); ylabel('Amplitude');%低通濾波器 B,A=butter(15,0.

35、15); z1=filter(B,A,y);figure(6) plot(z1) title(' 通過低通濾波器的信號 ') xlabel(' 點數(shù) '); ylabel('Amplitude');%抽樣判決 a=zeros(1,25) for n=1:25a(n)=z1(100*n-50);endfor n=1:25if a(n)>0.2a(n)=1;else a(n)=0;end end figure(7)stem(0:numPlot-1,a(1:numPlot),'bx') title(' 解調后的信號 

36、9;) xlabel('Time');ylabel('Amplitude');%誤碼率統(tǒng)計errorSym ratioSym=symerr(x,a'); figure(8)simbasebandex(0:1:5);title(' 過零解調后誤碼率統(tǒng)計 ')%PS調制和相干解調MATLA源代碼function PSK Fc=10; % 載頻 Fs=100; % 系統(tǒng)采樣頻率 Fd=1;% 碼速率N=Fs/Fd;df=10; numSymb=25;進行仿真的信息代碼個數(shù) M=2;% 進制數(shù)SNRpBit=60;%言噪比 SNR=SNRpBit

37、/log2(M);seed=12345 54321; numPlot=25;%產生 25 個二進制隨機碼x=randsrc(numSymb,1,0:M-1);% 產生 25個二進制隨機碼 figure(1)stem(0:numPlot-1,x(1:numPlot),'bx'); title(' 二進制隨機序列 ') xlabel('Time');ylabel('Amplitude');%調制 y=dmod(x,Fc,Fd,Fs,'psk',M,df); numModPlot=numPlot*Fs; t=0:numM

38、odPlot-1./Fs;figure(2) plot(t,y(1:length(t),'b-'); axis(min(t) max(t) -1.5 1.5); title(' 調制后的信號 ') xlabel('Time');ylabel('Amplitude');%在已調信號中加入高斯白噪聲 randn('state',seed(2);在已調信y=awgn(y,SNR-10*log10(0.5)-10*log10(N),'measured','dB');% 號中加入高斯白噪聲fi

39、gure(3)plot(t,y(1:length(t),'b-');% 畫出經過信道的實際信號 axis(min(t) max(t) -1.5 1.5);title(' 加入高斯白噪聲后的已調信號 ') xlabel('Time');ylabel('Amplitude');%相干解調z1=ddemod(y,Fc,Fd,Fs,'psk',M,df);獅輸出波形的相干M元頻移鍵控解調figure(4) stem(0:numPlot-1,x(1:numPlot),'bx'); hold on;stem(0

40、:numPlot-1,z1(1:numPlot),'ro');hold off;axis(0 numPlot -0.5 1.5);title(' 相干解調后的信號原序列比較 ')legend(' 原輸入二進制隨機序列 ',' 相干解調后的信號 ') xlabel('Time');ylabel('Amplitude');%誤碼率統(tǒng)計errorSym ratioSym=symerr(x,z1);figure(6) simbasebandex(0:1:5); title(' 相干解調后誤碼率統(tǒng)計

41、')%DPS調制和相干解調MATLA源代碼function DPSKFc=10; % 載頻Fs=100; % 系統(tǒng)采樣頻率Fd=1;% 碼速率N=Fs/Fd;df=10;numSymb=25;進行仿真的信息代碼個數(shù) M=2;% 進制數(shù)SNRpBit=60;%言噪比SNR=SNRpBit/log2(M);seed=12345 54321; numPlot=25;%產生 25 個二進制隨機碼 x=randsrc(numSymb,1,0:M-1);% 產生 25 個二進制隨機碼 figure(1)stem(0:numPlot-1,x(1:numPlot),'bx');titl

42、e(' 二進制隨機序列 ') xlabel('Time'); ylabel('Amplitude');%將 25 個絕對碼轉換為相對碼 b=zeros(1,26);b(2)=xor(b(1),x(1)for n=2:25 b(n+1)=xor(x(n),b(n);endfor n=1:25 x(n)=b(n+1);endfigure(2) stem(0:numPlot-1,x(1:numPlot),'bx'); title(' 二進制隨機序列相對碼 ') xlabel('Time');ylabel(

43、'Amplitude');%調制 y=dmod(x,Fc,Fd,Fs,'psk',M,df); numModPlot=numPlot*Fs; t=0:numModPlot-1./Fs;figure(3) plot(t,y(1:length(t),'b-'); axis(min(t) 2 -1.5 1.5); title(' 調制后的信號 ') xlabel('Time');ylabel('Amplitude');%在已調信號中加入高斯白噪聲 randn('state',seed(2)

44、;在已調信y=awgn(y,SNR-10*log10(0.5)-10*log10(N),'measured','dB');% 號中加入高斯白噪聲figure(4)plot(t,y(1:length(t),'b-');% 畫出經過信道的實際信號 axis(min(t) 2 -1.5 1.5);title(' 加入高斯白噪聲后的已調信號 ') xlabel('Time');ylabel('Amplitude');%相干解調 z1=ddemod(y,Fc,Fd,Fs,'psk',M,df)

45、; %將恢復出的相對碼變換為絕對碼 a=zeros(1,26);a(1)=xor(b(1),0); for n=2:25a(n)=xor(b(n),b(n-1)end for n=1:25 z1(n)=a(n+1);end獅輸出波形的相干M元頻移鍵控解調 figure(5) stem(0:numPlot-1,z1(1:numPlot),'ro'); axis(0 numPlot -0.5 1.5); title(' 相干解調后的信號的絕對碼 ') xlabel('Time');ylabel('Amplitude');%誤碼率統(tǒng)計e

46、rrorSym ratioSym=symerr(x,z1); figure(6)simbasebandex(0:1:5); title(' 相干解調后誤碼率統(tǒng)計 ')%DPS調制差分解調MATLA源代碼function DPSKFc=10; % 載頻Fs=100; % 系統(tǒng)采樣頻率Fd=1;% 碼速率N=Fs/Fd;df=10;numSymb=25;進行仿真的信息代碼個數(shù)M=2;% 進制數(shù)SNRpBit=60;%言噪比SNR=SNRpBit/log2(M);seed=12345 54321;numPlot=25;%產生 25 個二進制隨機碼x=randsrc(numSymb,1

47、,0:M-1);% 產生 25 個二進制隨機碼 figure(1)stem(0:numPlot-1,x(1:numPlot),'bx');title(' 二進制隨機序列 ')xlabel('Time'); ylabel('Amplitude');%將 25 個絕對碼轉換為 26 個相對碼 b=zeros(1,26);for n=1:25b(n+1)=xor(x(n),b(n);end figure(2) stem(0:25,b(1:26),'bx'); title(' 二進制隨機序列相對碼 ') x

48、label('Time');ylabel('Amplitude');%調制 y=dmod(b,Fc,Fd,Fs,'psk',M,df); numModPlot=26*Fs;t=0:numModPlot-1./Fs; figure(3)plot(t,y(1:length(t),'b-');axis(min(t) max(t) -1.5 1.5);title(' 調制后的信號 ') xlabel('Time');ylabel('Amplitude');在已調信%在已調信號中加入高斯白噪聲

49、 randn('state',seed(2); y=awgn(y,SNR-10*log10(0.5)-10*log10(N),'measured','dB');% 號中加入高斯白噪聲figure(4)stem(y)plot(t,y(1:length(t),'b-');% 畫出經過信道的實際信號 axis(min(t) max(t) -1.5 1.5);title(' 加入高斯白噪聲后的已調信號 ') xlabel('Time');ylabel('Amplitude');%差分解調帶通濾波器B,A=butter(10,0.2,'high'); a=filter(B,A,y);wc=0.3,0.9; B,A=butter(10,wc,'stop'); b=filter(B,A