中聾）夫普

信息科學與工程學院

課程設計報告書

課題：通信原理

班級：

學號：

姓名：

指導教師：郭麗梅

-0-0年十月

目錄

令目錄—一....——-———......---

----------------------------3

■■■■■■■■■■，*

令一、課程設計目的?一一......一一一一

■■一■■■■■■?Q1

令二、課程設計題目.一—一一一------

-......-------------------2

■■■■■■J

令三、課程設計實驗原理--------

……-———3

令四、建立模型描述一一一..........一一

———6

令五、模塊功能分析及源程序代碼…一……

-6

令六、倜試過程及結論—------——---

--——22

令七、調試分析一......................

————-34

令八、基于MATLAB的MASK、MFSK、

MPSK調制拓展(M=4;M=8)——-

-——34

令九、心得體會一------——....一

■■---■-■-■-------■---■■■■An■■

《十、參考文獻——-——--一

........................—41

摘要：

通信原理是通信工程專業相稱重要的學科，對日后就業和科研有重大的意義,

通過MATLAB,我們可以清楚地理解通信原理中難以理解的一面，對理論的知

識加以深化。

關鍵字：

MATLAB通信原理GUI序列頻譜相位相干非相干

一、課程設計題目

1）、應用MATLAB編制信號生成程序，并對信號進行HDB3編碼和譯碼。

2）、課程設計需要運用MATLAB編程實現2ASK,2FSK,2PSK,2DPSK

調制解調過程，并且輸出其源碼，調制后碼元以及解調后碼元的波形。

二、課程設計目的

通信原理是一門以算法為核心,理論和實踐性較強的學科。是電子信息工

程、通信工程專業、電子信息科學與技術專業的一門重要的專業技術基礎課。通

信原理課程設計是在學習完通信原理的相關理論后，進行的綜合性訓練課程，其

目的是：

1.使學生進一步鞏固通信原理的基本概念、理論、分析方法和實現方法；

2.增強學生應用Ma11ab語言編寫數字信號解決的應用程序及分析、解決實

際問題的能力；

3.全面復習課程所學理論知識，鞏固所學知識重點和難點，將理論與實踐很好

地結合起來；

4.提高綜合運用所學知識獨立分析和解決問題的能力；

5.純熟使用一種高級語言進行編程實現。

三、課程設計實驗原理

1.HDB3編碼解碼原理

HDB3碼：三階高密度雙極性碼。

HDB3碼與二進制序列的關系：

(1)二進制信號序列中的“0”碼在HDB3碼中仍編為“0”碼,二進制信

號中“1”碼，在HDB3碼中應交替地成+1和-1碼，但序列中出現四個連“0”

碼時應按特殊規律編碼(引入傳號交替反轉碼的“破壞點”V碼)；

(2)二進制序列中四個連“0”按以下規則編碼:信碼中出現四個連“0”

碼時，要將這四個連“0”碼用000V或B00V取代節來代替(B和V也是“1”

碼，可正、可負)。這兩個取代節選取原則是，使任意兩個相鄰v脈沖間的傳號

數為奇數時選用000V取代節，偶數時則選用B00V取代節。

2.二進制數字調制技術原理

數字信號的傳輸方式分為基帶傳輸和帶通傳輸，在實際應用中，大多數信

道具有帶通特性而不能直接傳輸基帶信號。為了使數字信號在帶通信道中傳輸，

必須使用數字基帶信號對載波進行調制，以使信號與信道的特性相匹配。這種用

數字基帶信號控制載波,把數字基帶信號變換為數字帶通信號的過程稱為數字調

制。通常使用鍵控法來實現數字調制，比如對載波的振幅、頻率和相位進行鍵控。

(1)2ASK信號的產生方法通常有兩種：模擬調制和鍵控法。解調有相干

解調和非相干解調。P=1時f(t)=Acoswt;p=0時f(t)=0;其功率譜密度是

基帶信號功率譜的線性搬移。

解調原理圖

定時脈沖

圖12ASK非相干解調接受系統

圖22ASK相干解調接受系統

(2)一個2FSK信號可以當作是兩個不同載波的2ASK信號的疊加。其解

調和解調方法和ASK差不多。2FSK信號的頻譜可以當作是f1和f2的兩個2

ASK頻譜的組合。

解調原理圖

圖32FSK相干解調接受系統

圖42FSK非相干解調接受系統

(3)2PSK以載波的相位變化作為參考基準的，當基帶信號為0時相位相

對于初始相位為0,當基帶信號為1時相對于初始相位為180°。

(4)2DPSK調制原理方框圖如下圖：

t

圖5間接法信號調制器原理方框圖

2DPSK信號的解調，重要有兩種方法，即相位比較法和相干解調法。相干解調法

原理方框圖如下圖:

圖6相干解調法原理方框圖

四、建立模型描述

用MATLAB實現二進制振幅鍵控（2ASK）的調制和解調

在這里用MATLAB函數來對這個二進制振幅鍵控（2ASK）來實現調制與解調

的仿真。二進制振幅鍵信號可以表達完畢一個單極性矩形脈沖序列與一個正弦型

載波的乘積。通常它的調制方法有兩種，即模擬幅度調制方法和鍵控方法，在

MATLAB里我們采用模擬幅度調制的方法,解調采用相干解調（包絡檢波法）的方

式。我們用SOURCE函數來產生一個原始二進制基帶信號，即一個單矩形脈沖序列。

以askM。du函數來進行模擬幅度調制,得到一個已調2ASK信號,并用此函數進

行此2ASK信號的頻譜分析。然后用gussian函數加入加性高斯白噪聲，再用demoA

SK函數進行想干解調并分別輸出各點的輸出波形，最后通過抽樣判決后得出輸

出波形。同時我們用CheekRatePe函數來得出誤碼率，最后運營主函數ASK_m

ain可以看到各種波形。用MATLAB函數來對2FSK、2PSK、2DPSK實現調

制與解調的仿真與2ASK的調制與解調過程的仿真類似。

五、模塊功能分析及源程序代碼

5.1HDB3碼的編碼和譯碼

globa1x

xn=x;

yn=xn；

num=0;%計數器初始化

fork=1:1ength(xn)

ifxn(k)==l%1的計數器

num=num+1;

ifnum/2==fix(num/2)/奇數個一時輸出一1,進行記性交替

yn(k)=l;

eIse

yn(k)=-l；

end

end

end

%HDB3編碼

num=0;

yh=yn;

sign=0；

v=zeros(1,1ength(yn));

B=zeros(1,length(yn))；%B脈沖位置記錄

fork=l:length(yn)

ifyn(k)==0

num=num+1；%連0個數計數

ifnum==4%假如4連0

num=O；%計數器清零

yh(k)=l*yh(k-4)*最后一個0改變與前一個非零符號同極性

v(k)=yh(k);

ifyh(k)==sign%假如當前V與前一個極性相同

yh(k)=-l*yh(k);%則當前V符號極性反轉，

yh(k-3)=yh(k);%添加B脈沖

B(k-3)=yh(k)；%8脈沖位置

v(k)=yh(k);%V脈沖位置

yh(k+1:length(yn))=—1*yh(k+1:length(yn));

end

sign二yh(k)；%計算前一個V符號的極性

end

eIse

nuni=O;%當前輸入為0，則計數器清零

end

end

re=[xn，,yn，,yh),v，,B，];

%HDB3解碼

input=yh;

decode=input;

sign=0;為極性標志初始化

fork=1:length(yh)

ifinput(k)"=0

ifsign==yh(k)%假如當前碼與前一?個非零碼的極性相同

decode(k-3:k)=[0000];

end

sign=input(k);%極性標志

end

end

decode=abs(decode)；

subplot(3,1,1);stairs(xn)；axis([l1ength(xn)-22]);

grid;ylabel('xn');

title。HDB3碼的編碼前圖形');

subplot(3,1,2);stairs(yh)；axis([1length(xn)-22]);

grid;ylabel('HDB3codes');

title(,HDB3碼的編碼后(解碼前)圖形’)；

subplot(3,1,3)；stairs(decode);axis([1length(xn)-22]);

grid;ylabe1('decodedxn');

titleCHDB3碼的解碼后圖形');

codel=yh;

decode1=decode；

5.2ASK的調制與解調

5.2.1source函數

functionsendSigna1=source(n,N)

sendSignal=randint(1,n)

bit=[];

fori=l:length(sendSignal)

ifsendSigna1(i)=0

bit1=zeros(1,N);

e1se

bitl=ones(1,N);

end

bit=[bit,bit1];

end

figure(1)

plot(l:length(bit),bit),title('發送端二進制波形'),gridon；

axis([0,N*length(sendSignal),-2,2]);

end

5.2.2askModu函數源程序

functiontransmittedSignal=askModu(signal,bitRate,fc,N)

%signal=[l0101001]；

%bitRate=l000000;

%fc=1000000;

%N=32；

t=linspace(0,1/bitRate,N);

c=sin(2*pi*t*fc);

transmittec!Signal=[];

fori=l:length(signal)

transmittedSignal=[transmittedSigna1,signa1(i)*c];

end

figure(2)

plot(1:length(transmittedSigna1),transmittedSigna1);tit1e('ASK調

制波形');gridon;

figure(3)

m=0:length(transmittedSigna1)-1;

F=fft(transmittedSignal)；

plot(m,abs(rea1(F))),title('ASK仿真頻譜分析');

gridon;

%figure(4)

%plot(m,imag(F))；title。ASKfrequency-domainanalysisimag,)；

%gridon;

End

5.2.3gussian函數源程序

functionsignal=gussian(transmittedSignal,noise)

signal=sqrt(2)*transmittedSignal;

signal=awgn(signa1,noise);

figure(5)

plot(1：length(signal),signal);

tit1e(z包含噪聲的波形')，gridon;

end

5.2.4CheckRatepe函數源程序

functionPeWrong=CheckRatePe(signalLsigna12,s)

rights=0;

wrongs=0;

forki=l:s—2

if(signall(ki)==signal2(ki))

rights=rights+1；

else

wrongs=wrongs+1;

end

end

PeWrong=wrongs/(wrongs+rights);

End

5.2.5demoASK函數源程序

functionbitstream=demoASK(receivedSignal,bitRate,fc,n,N)

loadnum

signa11=receivedSignal;

signa12=abs(signal1);%整流

signal3=filter(numl,1,signa12);%LPF,包絡檢波”

IN=fix(1ength(numl)/2);%延遲時間

bitstream=口;

LL=fc/bitRate*N;

i=IN+LL/2;

whi1e(i<=1ength(signa13))外判決

bitstream=[bitstrearn,signa13(i)>=0.51;

i=i+LL;

end

figure(6)

subplot(3,1,1);

plot(l:length(signa11),signal1);title('接受端波形(包含噪聲)')；gridon;

subplot(3,1,2);

p1ot(1：length(signa12),signal2)；title(，整流之后的波形')；gridon;

subplot(3,1,3);

plot(l:length(signa13),signa13);tit1e('LPF濾波后的包絡波形')；grid

on;

bit=[];

fori=l:1ength(bitstream)

ifbitstream(i)=0

bit1=zeros(l,N);

else

bit1=ones(1,N);

end

bit=[bit,bit1];

end

figure(7)

plot(bit),titleC接受端二進制波形'),gridon;

axis([0,N*1ength(bitstream),-2.5,2.5]);

end

5.2.6ASK_main函數源程序

closeal1

c1earal1

%(

ti=0；

fpeask=[]；

startn=-6;

endn=18；

forti=startn:endn

n=1000;

%}

n=16;

fc=1000000;%fc>=bitRatefc/bitRate為每個包含sin周期個數

bitRate=1000000；

N=50；

%noise=ti;

noise=l0;

signa1=source(n,N);

transmittedSignal=askModu(signa1,bitRate,fc,N);

signall=gussian(transmittedSignal,noise);

configueSignal=demoASK(signall,bitRate,fc,n,N);

%(

configueSigna1；

P=CheckRatePe(signal,configueSignal,n)

fpeask=[fpeask,P]；

end

figure(8);

semilogy(startn：length(fpeask)+startn-l,fpeask)；

gridon;

title('ASK—誤碼率');

xlabelCr/dB');

ylabe1CPeASK‘);

savePeRate.matfpeask

%}

5.3FSK的調制與解調

5.3.1source函數

functionsendSignal=source(n,N)

sendSignal=randint(1,n)

bit=[];

fori=l:length(sendSigna1)

ifsendSignal(i)==0

bit1=zeros(1,N);

else

bit1=ones(1,N)；

end

bit=[bit,bitl]；

end

figure(1)

plot(bit),titleC發送端二進制波形'),gridon；

axis([0,N*length(sendSigna1),-2.5,2.5]);

end

5.3.2FskModu函數源程序

functiontransmittedSignal=fskModu(signal,bitRate,fl,f2,N)

t=1inspace(0,1/bitRate,N);

cl=sin(2*pi*t*fl);

c2=sin(2*pi*t*f2);

transmittedSignal=[];

fori=l:length(signal)

ifsignal(i)==l

transmittedSigna1=[transmittedSignal,c1];

else

transmittedSignal=[transmittedSigna1,c2];

end

end

figure(2)

plot(1:1ength(transmittedSignal),transmittedSignal);titleCFSK

調制波形')；gridon;

figure(3)

m=0:1ength(transmittedSigna1)-1;

F=fft(transmittedSignal);

P1ot(m,abs(rea1(F))),tit1e('FSK仿真頻譜分析')；

gridon;

end

5.3.3gussian函數源程序

functionsignal=gussian(transmittedSignal,noise)

signal=sqrt(2)*transmittedSigna1;

signal=awgn(signa1,noise);

figure(4)

P1ot(l:length(signal),signal),title(z包含噪聲的波形')；

gridon;

end

5.3.4CheckRatepe函數源程序

functionPeWrong=CheckRatePe(signal1,signa12,s)

rights=0；

wrongs=0；

forki=1:s-2

if(signa11(ki)==signal2(ki))

rights=rights+1;

else

wrongs=wrongs+1；

end

end

PeWrong=wrongs/(wrongs+rights);

end

5.3.5demoFSK函數源程序

functionbitstream=demoFSK(receivedSigna1,bitRate,f1,f2,N)

1oadFSKnum

signa11=receivedSignal;

signa12=fiIter(gaotong,1,signall);為通過HPF,得到高頻分量

signa13=abs(signa12);/整流

signal3=filter(lowpass,1,signal3);%通過LPF,形成包絡

bitstream=[];

IN1=fix(length(lowpass)/2)+fix(length(gaotong)

/2);%延遲時間

bitstreaml=[];

LL=N;%每個bil的抽樣點數

i=IN1+LL/2；

while(i<=length(signal3))%判決

bitstreaml=[bitstream1,signa13(i)>=0.5];

i=i+LL;

end

bitstreaml

figure(5)

subplot(3,1,1);

plot(1：1ength(signall),signall)接受端波形(包含噪聲)’)；grido

n;

subplot(3,1,2)；

plot(l:length(signa12),signa12);title。通過HPF得到的高頻分量波形');

gridon；

subplot(3,1,3);

Plot(1:length(sigoal3),signal3)；titleC通過LPF后的包絡波形，);gr

idon;

signa14=fi1ter(daitong,1,signall)；先通過BPF,得到低頻分量

signa15=abs(signa14);%整流

signal5=fi1ter(1owpass,1,signal5);%通過LPF,形成包絡

IN2=fix(1ength(1owpass)/2)+fix(length(daitong)/

2);%延遲時間

bitstream2=口；

LL=N;%每個bit的抽樣點數

i=IN2+LL/2;

while(i<=length(signal5))%判決

bitstream2=[bitstream2,signal5(i)>=0.5];

i=i+LL;

end

bitstream2

figure(6)

subplot(3,1,1);

p1ot(l:length(signall),signall);tit1e('接受端波形(包含噪聲)')；grid

on;

subplot(3,1,2);

plot(1：length(signal4),signa14)；title。通過BPF得的低頻分量波形')；

gridon;

subp1ot(3,1,3);

p1ot(l:length(signal5),signal5);title('通過LPF后的包絡波形');gridon；

fori=1：min(length(bitstream1),length(bitstream2))%判決

if(bitstreaml(i)>bitstream2(i))

bitstream(i)=l;

else

bitstream(i)=0;

end

end

bitstream

bit二口；外接受端波形

fori=1：length(bitstream)

ifbitstream(i)==0

bit1=zeros(1,N);

else

bitl=ones(1,N)；

end

bit=[bit,bit1];

end

figure(7)

p1ot(bit),title('接受端波形(解調后波形)’)，gridon;

axis([0,N*length(bitstream),-2.5,2.5]);

end

5.3.6FSK_main函數源程序

closea11

c1earall

%(

ti=0；

fpefsk=[];

startn=-6;

endn=18；

forti=startn:endn

n=1000;

%}

n=16；

f1=18000000;

f2=6000000;

bitRate=l000000；

N=50;

%noise=ti;

noise=10;

signal=FSKsource(n,N);

transmittedSignal=fskModu(signa1,bitRate,f1,f2,N);

signa11=FSKgussian(transmittedSigna1,noise);

configueSigna1二demoFSK(signal1,bitRate,f1,f2,N);

%(

configueSigna1;

P=CheckRatePe(signal,configueSignal,n)

fpefsk=[fpefsk,P];

end

figure(8);

semilogy(startn:1ength(fpefsk)+startn-1,fpefsk);

gridon;

title—BitErrorRateOfFSK');

x1abe1('r/dB1);

y1abeICPeFSK*);

1oadPeRate

savePeRate.matfpefskfpeask

%}

5.4PSK的調制與解調

5.4.1source函數

functionsendSigna1=source(n,N)

sendSignal=randint(1,n)

bit=[];

fori=l:1ength(sendSigna1)

ifsendSigna1(i)==0

bit1=zeros(1,N);

e1se

bit1=ones(1,N);

end

bit=[bit,bitl];

end

figure(1)

P1ot(bit),titie('發送端二進制波形')，gridon;

axis([0,N*1ength(sendSignal),-2.5,2.5]);

end

5.4.2PSKModu函數源程序

functiontransmittedSigna1=bpskModu(signal,bitRate,fc,N)

t=linspace(0,1/bitRate,N);

cl=sin(2*pi*t*fc)；

c2=sin(2*pi*t*fc+pi);

transmittedSignal=[1;

fori=l:1ength(signal)

ifsigna1(i)==1

transmittedSignal=[transmittedSignal,cl];

else

transmittedSignal=[transmittedSignal,c2]；

end

end

figure(2)

plot(1:1ength(transmittedSignal),transmittedSignal)；title('BPSK')；gr

idon;

figure(3)

m=0:length(transmittedSignal)-1;

F=fft(transmittedSigna1);

plot(m,abs(real(F))),title。2PSK仿真頻譜分析')；

gridon;

end

5.4.3gussian函數源程序

functionsigna1=BPSKgussian(transmittedSignal,noise)

signa1=sqrt(2)*transmittedSignal;

signal=awgn(signal,noise);

figure(4)

plot(1:length(signal),signal),gridon;

title('包含噪聲的波形')

end

5.4.4CheckRatepe函數源程序

functionPeWrong=BPSKCheckRatePe(signa11,signal2,s)

rights=0;

wrongs=0;

forki=1：s-2

if(signa11(ki)二二signa12(ki))

rights=rights+1;

else

wrongs=wrongs+1;

end

end

PeWrong=wrongs/(wrongs+rights);

End

5.4.5demoPSK函數源程序

functionbitstream=demoBPSK(receivedSigna1,bitRate,fc,n,N)

loadnum%讀取num中存儲的低通濾波用的數據

signa11=receivedSigna1;

t=linspace(0,1/bitRate,N);

c=sin(2*pi*t*fc)；

signal二口；

fori=1:n

signal=［signal,c］；

end

signa12=signall.*signal;為乘同頻同相sin

signa13=fiIter(numl,1,signal2)；%LPF,包絡檢波

IN=fix(1ength(numl)/2)；%延遲時間

bitstream=1］；

LL=fc/bitRate*N；

i=IN+LL/2;

while(i<=1ength(signal3))%判決

bitstream=[bitstream,signal3(i)>=0];

i=i+LL;

end

figure(5)

subplot(3,1,1);

p1ot(l:length(signal1),signa11);title('接受端波形(包含噪聲)');gridon；

subplot(3,1,2);

plot(1：length(signal2),signa12)；

丫126。1('相干解調’)；

tit1e(z乘同頻同相sin之后的波形');gridon;

subplot(3,1,3)；

p1ot(1：1ength(signal3),signal3);title。LPF濾波后的包絡波形')；grido

n;

bit=[]；

fori=1:length(bitstream)

ifbitstream(i)==O

bit1=zeros(1,N);

else

bitl=ones(1,N);

end

bit=[bit,bit1];

end

figure(6)

plot(bit)；titlec接受端二進波形(解調后波形)')；gridon;

axis([0,N*1ength(bitstream),-2.5,2.5]);

end

5.4.6PSK_main函數源程序

c1osea11

c1earall

%(

ti=0;

fpepsk=口；

startn=-30;

endn=—6;

forti=startn:endn

n=l000；

%}

n=16;

fc=1000000；

bitRate=1000000;

N=50；

%noise=ti;

noise=10;

signal=BPSKsource(n,N);

transmittedSignal=bpskModu(signa1,bitRate,fc,N);

signa11=BPSKgussian(transmittedSignal,noise);

configueSignal=demoBPSK(signail,bitRate,fc,n,N)；

%(

configueSignal;

P=CheckRatePe(signa.l,configueSigna1,n)

fpepsk=[fpepsk,P]；

end

figure(8);

semi1ogy(startn:1ength(fpepsk)+startn-1,fpepsk);

gridon;

title(，BitErrorRateOfPSK');

x1abe1('r/dB');

ylabel('PePSK');

%loadPeRate

%savePeRate.matfpepskfpeaskfpefsk

savePeRatep.matfpepsk

%}

5.5DPSK的調制與解調

fs=30000;

Time_Hold_0n=0.1；

NumUnit=fs*TimeHoldOn;

High_Leve1=ones(1,Num_Unit)；

LowLevel=zeros(1,Num_Unit)；

w=300;

A=1;

Sign_Set=[0,1,1,0,1,0,0,1]

Lenth_0f_Sign=length(Sign_Set);

st=zeros(1,Nuui_Unit*Lenth_Of_Sign)；

sign_orign=zeros(1,Num_Unit*Lenth_Of_Sign)；

sign_resu1t=zeros(1,Num_Unit*Lenth_0f_Sign);

t=0:1/fs:TimeHo1d_0n*Lenth_0f_Sign—1/fs;

forI=1:Lenth_Of_Sign

ifSign_Set⑴==1

sign_orign((1-1)*Num_Unit+1：I*Num_Unit)=High_Leve1;

else

sign_orign((I-1)*Num_Unit+1：I*Num_Unit)=Low_Leve1；

end

end

for1=1:Lenth_0f_Sign

ifSign_Set(I)==1

st((1-1)*NumUnit+1:T*Num_Unit)=A*cos(2*pi*w*

t((1-1)*Num_Unit+1:I*Num_Unit)+(pi/2))；

else

st((I-l)*Num_Unit+1：I*Num_Unit)=A*cos(2*pi*

w*t((I-1)*Num_Unit+1：I*Num_Unit))；

end

end

figure

subplot(2,1,1)

plot(t,sign_orign);

axis([0,Time_Hold_On*(LenthOf_Sign+1),-(A/2),A+(A/

2)]);

title（'原始信號’）；

grid

subplot(2,1,2);

p1ot(t,st)；

axis([0,Time_Hold_On*(Lenth_Of_Sign+1),-3*(A/2),3*(A/

2)]);

title('調制后的信號’);

grid

dt=st.*cos(2*pi*w*t);

figure

subplot(2,1,1)

plot(t,dt);

axis([0,Time_Hold_0n*(Lenth_OfSign+1),—3*(A/2),3*(A/

2)])；

title('相乘后的波形’);

grid

[N,Wn]=buttord(2*pi*50,2*pi*150,3,25/s');%臨界頻率采用角頻率表達

[b,a]=butter(N,Wn,s');

[bz,az]=impinvar(b,a,fs)；%映射為數字的

dt=fiIter(bz,az,dt)；

subplot(2,1,2)

plot(t,dt);

axis([0,Time_Ho1d_On*(Lenth_0f_Sign+1),—3*(A/2),3*(A/

2)]);

title('低通濾波后的波形');

grid

for1=1：Lenth_Of_Sign

ifdt((2*1-1)*Num_Unit/2)<0.25

sign_resuit((I-1)*Num_Unit+1：I*Num_Unit)=High_Level;

else

sign_resu1t((1-1)*Num_Unit+1:I*Num_Unit)=Low_Level;

end

end

figure

plot(t,signresult)；

axis([0,Time_Hold_0n*(Lenth_Of_Sign+1),-3*(A/2),3*(A/2)]);

title('逆碼變換后的波形’);

grid

六、調試過程及結論

6.1基于MATLAB的2ASK調制解調仿真過程及結論

6.1.1仿真過程中的各點波形

圖7發送端二進制波形

Figure2匚］叵］區|

FileEditViewInsertToolsDesktopWindowHelp

D■口昌6⑥要口回■0

圖8已調2ASK波形

圖92ASK頻譜分析

圖10加入高斯白噪聲后的2ASK波形

圖11解調過程中各點的輸出波形

圖13解調后的波形

6.2基于MATLAB的2FSK調制解調仿真過程及結論

圖14發送端二進制波形

圖15已調2FSK波形

圖16FSK仿真頻譜分析圖

圖17

圖18

圖19解調后的波形

6.3基于MATLAB的2PSK調制解調仿真過程及結論

圖20接受端二進制波形

圖212PSK調制后的波形

圖222PSK仿真頻譜分析圖

圖23加入高斯白噪聲后的2PSK波形

Figure5-□X

FileEditViewInsertToolsDesktopWindowHelp

口方0昌4戲0㈠⑤要□日■口

接收端波形（包含噪聲:

!IIIIII

5illIIII

0

-5iiiiiii

0100200300400500600700800

乘同頻同相sin之后的波形

5

%

犍

0

H-

要

-5

0100200300400500600700800

LPF濾波后的包絡波形

1

0

-1

0100200300400500600700800

圖24

T.口X

FileEditViewInsertToolsDesktopWindowHelpi

口自。與Q煢口回回國

圖25接受端二進制波形

6.4基于MATLAB的2DPSK調制解調仿真過程及結論

圖26

圖27

圖28

6.5基于MATLAB的HDB3碼的編碼與譯碼仿真過程

及結論