matlab課程設(shè)計(jì)--ssb信號(hào)的仿真分析

1、<p><b>　　課程設(shè)計(jì)任務(wù)書</b></p><p>　　學(xué)生姓名： 專業(yè)班級(jí)： </p><p>　　指導(dǎo)教師： 工作單位： </p><p>　　題 目: SSB信號(hào)的仿真分析</

2、p><p><b>　　初始條件：</b></p><p><b>　　MATLAB軟件</b></p><p>　　數(shù)字信號(hào)處理與圖像處理基礎(chǔ)知識(shí)</p><p>　　要求完成的主要任務(wù):</p><p>　　調(diào)制信號(hào)：分別為300Hz正弦信號(hào)和三角波信號(hào)；載波頻率：30kHz

3、；解調(diào)方式：同步解調(diào)；</p><p>　　要求：畫出以下三種情況下調(diào)制信號(hào)、已調(diào)信號(hào)、解調(diào)信號(hào)的波形、頻譜以及解調(diào)器輸入輸出信噪比的關(guān)系曲線；</p><p>　　1）調(diào)制信號(hào)幅度=0.8×載波幅度；2）調(diào)制信號(hào)幅度=載波幅度；</p><p>　　3）調(diào)制信號(hào)幅度=1.5×載波幅度。</p><p><b>

4、　　時(shí)間安排</b></p><p>　　第17周，安排任務(wù)（鑒主5樓實(shí)驗(yàn)室）</p><p>　　第17-18周，仿真設(shè)計(jì)（鑒主5樓實(shí)驗(yàn)室）</p><p>　　第19周，完成（答辯，提交報(bào)告，演示） </p><p>　　指導(dǎo)教師簽名：

5、 年 月 日</p><p>　　系主任（或責(zé)任教師）簽名： 年 月 日</p><p><b>　　目 錄</b></p><p><b>　　摘要1</b></p><p>　

6、　Abstract2</p><p>　　1 SSB調(diào)制與解調(diào)原理3</p><p>　　1.1SSB調(diào)制原理3</p><p>　　1.2 SSB解調(diào)原理與抗噪性能4</p><p>　　2 SSB調(diào)制解調(diào)分析的MATLAB實(shí)現(xiàn)5</p><p>　　3 SSB調(diào)制的實(shí)現(xiàn)6</p><

7、p>　　3.1 調(diào)制信號(hào)為正弦信號(hào)6</p><p>　　3.1.1 調(diào)制信號(hào)幅度=0.8×載波幅度8</p><p>　　3.1.2 調(diào)制信號(hào)幅度=載波幅度10</p><p>　　3.1.3 調(diào)制信號(hào)幅度=1.5×載波幅度11</p><p>　　3.2 調(diào)制信號(hào)為三角波信號(hào)13</p>

8、<p>　　3.2.1 調(diào)制信號(hào)幅度=0.8×載波幅度15</p><p>　　3.2.2 調(diào)制信號(hào)幅度=載波幅度17</p><p>　　3.2.3 調(diào)制信號(hào)幅度=1.5×載波幅度19</p><p><b>　　4 心得體會(huì)22</b></p><p><b>　　5

9、參考文獻(xiàn)23</b></p><p><b>　　摘要</b></p><p>　　MATLAB軟件廣泛用于數(shù)字信號(hào)分析，系統(tǒng)識(shí)別，時(shí)序分析與建模，神經(jīng)網(wǎng)絡(luò)、動(dòng)態(tài)仿真等方面有著廣泛的應(yīng)用。新版的MATLAB增強(qiáng)了圖形處理功能，并在WINDOWS環(huán)境下運(yùn)行?，F(xiàn)今，MATLAB的發(fā)展已大大超出了“矩陣實(shí)驗(yàn)室”的范圍，它的配備了涉及到自動(dòng)控制、信息處理、計(jì)算機(jī)

10、仿真等種類繁多的工具箱（Tool Box），這些工具箱有數(shù)理統(tǒng)計(jì)、信號(hào)處理、系統(tǒng)辨識(shí)、最優(yōu)化、穩(wěn)健等等。本次課程設(shè)計(jì)主要利用MATLAB集成環(huán)境下的M文件，編寫程序來實(shí)現(xiàn)SSB解調(diào)，分別利用300HZ正弦波和矩形波，對(duì)30KHZ正弦波進(jìn)行調(diào)制，觀察調(diào)制信號(hào)、已調(diào)信號(hào)和解調(diào)信號(hào)的波形和頻譜分布，并在解調(diào)時(shí)引入高斯白噪聲，對(duì)解調(diào)前后信號(hào)進(jìn)行信噪比的對(duì)比分析，估計(jì)SSB調(diào)制解調(diào)系統(tǒng)的性能。</p><p><b

11、>　　Abstract</b></p><p>　　MATLAB software is widely used in digital signal analysis, system identification, time series analysis and modeling, neural networks, dynamic simulation have a wide range of

12、 applications. The new version of MATLAB enhanced graphics processing functions, and in WINDOWS environment. Today, MATLAB development has gone far beyond the "matrix laboratory" the scope, it is equipped with

13、a related to automatic control, information processing, computer simulation, such as a wide variety of toolbox (Tool </p><p>　　1 SSB調(diào)制與解調(diào)原理</p><p>　　1.1SSB調(diào)制原理</p><p>　　SSB調(diào)制屬于幅度調(diào)制。

14、幅度調(diào)制是用調(diào)制信號(hào)去控制高頻載波的振幅，使其按調(diào)制信號(hào)的規(guī)律而變化的過程。</p><p>　　與標(biāo)準(zhǔn)幅度調(diào)制相比，單邊帶調(diào)制（SSB）對(duì)于頻譜和輸出功率的利用率更高。盡管很少用于數(shù)據(jù)傳送，SSB仍廣泛地用于HF和VHF低端的語(yǔ)音通訊。雙邊帶調(diào)制信號(hào)包含有兩個(gè)完全相同的基帶信號(hào)，即上、下邊帶。由于兩個(gè)邊帶含的信息相同，因而從信息傳輸角度考慮，傳送一個(gè)邊帶同樣可以達(dá)到信息傳輸?shù)哪康?本設(shè)計(jì)只考慮上邊帶信號(hào)。單邊帶

15、調(diào)制，就是通過某種辦法，只傳送一個(gè)邊帶的調(diào)制方法。</p><p>　　設(shè)調(diào)制信號(hào)為單頻信號(hào)f(t)=Amcosωmt，載波為c(t)=cosωct，則調(diào)制后的雙邊帶時(shí)域波形為：</p><p>　　SDSB(t)＝ Amcosωmt cos t＝ [Amcos(ωc+ωm)t+ Amcos(ωc-ωm)t] /2

16、 保留上邊帶，波形為：</p><p>　　SUSB(t)＝[Amcos(ωc+ωm)t]/2＝Am(cosωctcosωmt-sinωctsinωmt) /2

17、 保留下邊帶，波形為：</p><p>　　SLSB(t)＝[Amcos(ωc-ωm)t]/2＝Am(cosωctcosωmt+sinωctsinωmt) /2 上兩式中的

18、第一項(xiàng)與調(diào)制信號(hào)和載波信號(hào)的乘積成正比，稱為同相分量；而第二項(xiàng)的乘積則是調(diào)制信號(hào)與載波信號(hào)分別移相90°后相乘的結(jié)果，稱為正交分量。</p><p>　　SSB調(diào)制框圖如下:</p><p>　　圖1 SSB調(diào)制框圖</p><p>　　1.2 SSB解調(diào)原理與抗噪性能</p><p>　　解調(diào)是調(diào)制的逆過程，其作用是從接收的已調(diào)信

19、號(hào)中恢復(fù)原基帶信號(hào)（即調(diào)制信號(hào)）。解調(diào)的方法可分為兩類：相干解調(diào)和非相干解調(diào)（包絡(luò)檢波）。</p><p>　　相干解調(diào)，也稱同步檢波，為了無失真地恢復(fù)原基帶信號(hào)，接收端必須提供一個(gè)與接收的已調(diào)載波嚴(yán)格同步（同頻同相）的本地載波（稱為相干載波），它與接受的已調(diào)信號(hào)相乘后，經(jīng)低通濾波器取出低頻分量，即可得到原始的基帶調(diào)制信號(hào)。</p><p>　　包絡(luò)檢波器就是直接從已調(diào)波的幅度中提取原調(diào)制

20、信號(hào)，通常由半波或全波整流器和低通濾波器組成。</p><p>　　由于SSB信號(hào)是抑制載波的已調(diào)信號(hào)，它的包絡(luò)不再與調(diào)制信號(hào)的變化規(guī)律一致，因而不能采用簡(jiǎn)單的包絡(luò)檢波來恢復(fù)調(diào)制信號(hào)。SSB信號(hào)解調(diào)時(shí)需采用相干解調(diào)。</p><p>　　SSB相干解調(diào)性能分析模型如圖3所示：</p><p>　　圖2 SSB相干解調(diào)</p><p>　　解調(diào)

21、器的輸入信噪比為：</p><p>　　解調(diào)器的輸出信噪比為：</p><p><b>　　制度增益為：</b></p><p>　　這是因?yàn)樵赟SB系統(tǒng)中，信號(hào)和噪聲有相同的表示形式，所以相干解調(diào)過程中，信號(hào)和噪聲中的正交分量均被抑制掉，故信噪比沒有改善。 </p><p>　　2 SSB調(diào)制解調(diào)分析的MATLA

22、B實(shí)現(xiàn)</p><p>　　三角波函數(shù)sawtooth：調(diào)用格式為x = sawtooth(t, width).功能：產(chǎn)生一個(gè)周期為2π、幅度在-1到+1之間的周期性三角波信號(hào)。其中width表示最大幅度出現(xiàn)的位置：即在一個(gè)周期內(nèi)，信號(hào)從t=0到width×2π時(shí)函數(shù)值從-1到+1線性增加，而從width×2π到2π又是從+1到-1線性下降。width取值在0 ~ 1之間。 若x = sawt

23、ooth(Ωt, width)，則對(duì)應(yīng)的周期為2π/Ω。</p><p>　　信號(hào)DSB調(diào)制采用MATLAB函數(shù)modulate實(shí)現(xiàn)，其函數(shù)格式為：</p><p>　　Y = MODULATE(X,Fc,Fs,METHOD,OPT)</p><p>　　X為基帶調(diào)制信號(hào)，F(xiàn)c為載波頻率，F(xiàn)s為抽樣頻率，METHOD為調(diào)制方式選擇，SSB調(diào)制時(shí)為’am’，OPT在S

24、SB調(diào)制時(shí)可不選，F(xiàn)s需滿足Fs > 2*Fc + BW，BW為調(diào)制信號(hào)帶寬。</p><p>　　SSB信號(hào)解調(diào)采用MATLAB函數(shù)demod實(shí)現(xiàn)，其函數(shù)使用格式為：</p><p>　　X = DEMOD(Y,Fc,Fs,METHOD,OPT)</p><p>　　Y為SSB已調(diào)信號(hào)，F(xiàn)c為載波頻率，F(xiàn)s為抽樣頻率，METHOD為解調(diào)方式選擇，SSB解調(diào)時(shí)

25、為’am’，OPT在SSB調(diào)制時(shí)可不選。</p><p>　　觀察信號(hào)頻譜需對(duì)信號(hào)進(jìn)行傅里葉變換，采用MATLAB函數(shù)fft實(shí)現(xiàn)，其函數(shù)常使用格式為：Y=FFT(X,N)，X為時(shí)域函數(shù)，N為傅里葉變換點(diǎn)數(shù)選擇，一般取值。頻域變換后，對(duì)頻域函數(shù)取模，格式：Y1=ABS(Y)，再進(jìn)行頻率轉(zhuǎn)換，轉(zhuǎn)換方法：f=(0:length(Y)-1)’*Fs/length(Y)</p><p>　　3 SS

26、B調(diào)制的實(shí)現(xiàn)</p><p>　　3.1 調(diào)制信號(hào)為正弦信號(hào)</p><p><b>　　程序：</b></p><p>　　Fs=100000; %抽樣頻率fs=100000;</p><p>　　t=[0:1/Fs:0.01];</p><p>　　Fc=30000;

27、 %載波頻率</p><p><b>　　a=0.8</b></p><p>　　m=a*cos(300*2*pi*t); %調(diào)制信號(hào)</p><p><b>　　X=fft(m);</b></p><p>　　X=abs(X(1:length(X)/2+1)

28、); %調(diào)制信號(hào)頻譜</p><p>　　frqX=(0:length(X)-1)*Fs/length(X)/2</p><p>　　S = modulate(m,Fc,Fs,'amssb'); %對(duì)信號(hào)進(jìn)行調(diào)制</p><p><b>　　Y=fft(S);</b></p><p>　　Y=abs

29、(Y(1:length(Y)/2+1)); </p><p>　　frqY=(0:length(Y)-1)*Fs/length(Y)/2 ;%已調(diào)信號(hào)頻譜</p><p>　　set(gcf,'color','w')</p><p><b>　　figure(1)</b></p><

30、;p>　　subplot(221) %繪制曲線</p><p><b>　　plot(t,m)</b></p><p>　　xlabel('調(diào)制信號(hào)波形')</p><p>　　subplot(222)</p><p>　　plot(frqX,X)</p&g

31、t;<p>　　axis([0 3000 0 max(X)])</p><p>　　xlabel('調(diào)制信號(hào)頻譜')</p><p>　　subplot(223)</p><p><b>　　plot(t,S)</b></p><p>　　xlabel('已調(diào)信號(hào)波形')&l

32、t;/p><p>　　subplot(224)</p><p>　　plot(frqY,Y)</p><p>　　axis([0 60000 0 max(Y)])</p><p>　　xlabel('已調(diào)信號(hào)頻譜')</p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

33、%%%%%%%%%%</p><p>　　sn=awgn(S,4); %加入高斯白噪聲</p><p>　　sn1=awgn(S,10);</p><p>　　sn2=awgn(S,15);</p><p>　　sn3=awgn(S,20);</p><p>　　sn4=awgn(S,25);<

34、;/p><p>　　Y1=demod(S,Fc,Fs,'amssb'); %無噪聲已調(diào)信號(hào)解調(diào)</p><p>　　YYN=demod(sn,Fc,Fs,'amssb'); %加噪聲已調(diào)信號(hào)解調(diào)</p><p>　　YYN1=demod(sn1,Fc,Fs,'amssb');</p><

35、p>　　YYN2=demod(sn2,Fc,Fs,'amssb');</p><p>　　YYN3=demod(sn3,Fc,Fs,'amssb');</p><p>　　YYN4=demod(sn4,Fc,Fs,'amssb');</p><p>　　J1=fft(sn);</p><p&g

36、t;　　J1=abs(J1(1:length(J1)/2+1))</p><p>　　frqJ1=(0:length(J1)-1)*Fs/length(J1)/2; %加噪聲后已調(diào)信號(hào)頻譜</p><p>　　J2=fft(YYN)</p><p>　　J2=abs(J2(1:length(J2)/2+1))</p><p>　　frqJ2=

37、(0:length(J2)-1)*Fs/length(J2)/2 %加噪聲后解調(diào)信號(hào)頻譜</p><p>　　set(gcf,'color','w')</p><p><b>　　figure(2)</b></p><p>　　subplot(221)</p><p>　　plot(

38、t,YYN);</p><p>　　xlabel('加噪聲解調(diào)信號(hào)波形')</p><p>　　subplot(222)</p><p>　　plot(frqJ2,J2);</p><p>　　axis([0 3000 0 max(J2)])</p><p>　　xlabel('加噪聲解調(diào)信號(hào)頻

39、譜')</p><p>　　subplot(223)</p><p>　　plot(t,Y1)</p><p>　　xlabel('無噪聲解調(diào)信號(hào)波形')</p><p>　　subplot(224)</p><p>　　plot(frqJ1,J1)</p><p>　　

40、xlabel('無噪聲解調(diào)信號(hào)頻譜')</p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p>　　dyi=sn-S; %高斯白噪聲</p><p>　　s_ni=var(S)/var(dyi); %輸入信噪比</p

41、><p>　　dyo=YYN-Y1; %解調(diào)后噪聲</p><p>　　s_no=var(Y1)/var(dyo); %輸出信噪比 </p><p>　　dyi1=sn1-S; </p><p>　　s_ni1=var(S)/var(dy

42、i1); </p><p>　　dyo1=YYN1-Y1; </p><p>　　s_no1=var(Y1)/var(dyo1); </p><p>　　dyi2=sn2-S; </p><p>　　s_ni2=var(S)/var(dyi2);

43、 </p><p>　　dyo2=YYN2-Y1; </p><p>　　s_no2=var(Y1)/var(dyo2); </p><p>　　dyi3=sn3-S; </p><p>　　s_ni3=var(S)/var(dyi3

44、); </p><p>　　dyo3=YYN3-Y1; </p><p>　　s_no3=var(Y1)/var(dyo3); </p><p>　　dyi4=sn4-S; </p><p>　　s_ni4=var(S)/var(dy

45、i4); </p><p>　　dyo4=YYN4-Y1; </p><p>　　s_no4=var(Y1)/var(dyo4); </p><p>　　in=[s_ni,s_ni1,s_ni2,s_ni3,s_ni4]; </p><p>　

46、　out=[s_no,s_no1,s_no2,s_no3,s_no4];</p><p>　　set(gcf,'color','w');</p><p>　　figure(3);</p><p>　　plot(in,out,'*') </p><p><b>　　hold on<

47、;/b></p><p>　　plot(in,out)</p><p>　　xlabel('輸入信噪比')</p><p>　　ylabel('輸出信噪比')</p><p>　　3.1.1 調(diào)制信號(hào)幅度=0.8×載波幅度</p><p>　　調(diào)用程序，程序中a=0.8。

48、</p><p>　　調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜如圖3所示：</p><p>　　圖3 調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜圖</p><p>　　解調(diào)信號(hào)的波形、頻譜如圖4所示：</p><p>　　圖4 有噪聲、無噪聲的解調(diào)信號(hào)波形和頻譜</p><p>　　輸入輸出信噪比關(guān)系曲線如圖5所示：</p>

49、<p>　　圖5 輸入輸出信噪比關(guān)系曲線</p><p>　　3.1.2 調(diào)制信號(hào)幅度=載波幅度</p><p>　　調(diào)用程序，程序中a=1。</p><p>　　調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜如圖6所示：</p><p>　　圖6 調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜圖</p><p>　　解調(diào)信號(hào)的波形、頻譜

50、如圖7所示：</p><p>　　圖7 有噪聲、無噪聲的解調(diào)信號(hào)波形和頻譜</p><p>　　輸入輸出信噪比關(guān)系曲線如圖8所示： </p><p>　　圖8 輸入輸出信噪比關(guān)系曲線</p><p>　　3.1.3 調(diào)制信號(hào)幅度=1.5×載波幅度</p><p>　　調(diào)用程序，程序中a=1.5。</p&

51、gt;<p>　　調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜如圖9所示：</p><p>　　圖9 調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜圖</p><p>　　解調(diào)信號(hào)的波形、頻譜如圖10所示：</p><p>　　圖10 有噪聲、無噪聲的解調(diào)信號(hào)波形和頻譜</p><p>　　輸入輸出信噪比關(guān)系曲線如圖11所示： </p>&l

52、t;p>　　圖11 輸入輸出信噪比關(guān)系曲線</p><p>　　3.2 調(diào)制信號(hào)為三角波信號(hào)</p><p><b>　　程序：</b></p><p>　　Fs=100000; %抽樣頻率fs=100000;</p><p>　　t=[0:1/Fs:0.01];</p>&

53、lt;p>　　Fc=30000; %載波頻率</p><p><b>　　a=0.8</b></p><p>　　m=a*sawtooth(300*2*pi*t); %調(diào)制信號(hào) </p><p><b>　　X=fft(m);</b></p><p>　

54、　X=abs(X(1:length(X)/2+1)); %調(diào)制信號(hào)頻譜</p><p>　　frqX=(0:length(X)-1)*Fs/length(X)/2</p><p>　　S = modulate(m,Fc,Fs,'amssb'); %對(duì)信號(hào)進(jìn)行調(diào)制</p><p><b>　　Y=fft(S);</b><

55、;/p><p>　　Y=abs(Y(1:length(Y)/2+1)); </p><p>　　frqY=(0:length(Y)-1)*Fs/length(Y)/2 ;%已調(diào)信號(hào)頻譜</p><p>　　set(gcf,'color','w')</p><p><b>　　figure(1)

56、</b></p><p>　　subplot(221) %繪制曲線</p><p><b>　　plot(t,m)</b></p><p>　　xlabel('調(diào)制信號(hào)波形')</p><p>　　subplot(222)</p><p&

57、gt;　　plot(frqX,X)</p><p>　　axis([0 3000 0 max(X)])</p><p>　　xlabel('調(diào)制信號(hào)頻譜')</p><p>　　subplot(223)</p><p><b>　　plot(t,S)</b></p><p>　　x

58、label('已調(diào)信號(hào)波形')</p><p>　　subplot(224)</p><p>　　plot(frqY,Y)</p><p>　　axis([0 60000 0 max(Y)])</p><p>　　xlabel('已調(diào)信號(hào)頻譜')</p><p>　　%%%%%%%%%%

59、%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p>　　sn=awgn(S,4); %加入高斯白噪聲</p><p>　　sn1=awgn(S,10);</p><p>　　sn2=awgn(S,15);</p><p>　　sn3=awgn(S,20);&

60、lt;/p><p>　　sn4=awgn(S,25);</p><p>　　Y1=demod(S,Fc,Fs,'amssb'); %無噪聲已調(diào)信號(hào)解調(diào)</p><p>　　YYN=demod(sn,Fc,Fs,'amssb'); %加噪聲已調(diào)信號(hào)解調(diào)</p><p>　　YYN1=demod(sn1

61、,Fc,Fs,'amssb');</p><p>　　YYN2=demod(sn2,Fc,Fs,'amssb');</p><p>　　YYN3=demod(sn3,Fc,Fs,'amssb');</p><p>　　YYN4=demod(sn4,Fc,Fs,'amssb');</p>&

62、lt;p>　　J1=fft(sn);</p><p>　　J1=abs(J1(1:length(J1)/2+1))</p><p>　　frqJ1=(0:length(J1)-1)*Fs/length(J1)/2; %加噪聲后已調(diào)信號(hào)頻譜</p><p>　　J2=fft(YYN)</p><p>　　J2=abs(J2(1:lengt

63、h(J2)/2+1))</p><p>　　frqJ2=(0:length(J2)-1)*Fs/length(J2)/2 %加噪聲后解調(diào)信號(hào)頻譜</p><p>　　set(gcf,'color','w')</p><p><b>　　figure(2)</b></p><p>　　

64、subplot(221)</p><p>　　plot(t,YYN);</p><p>　　xlabel('加噪聲解調(diào)信號(hào)波形')</p><p>　　subplot(222)</p><p>　　plot(frqJ2,J2);</p><p>　　axis([0 3000 0 max(J2)])<

65、;/p><p>　　xlabel('加噪聲解調(diào)信號(hào)頻譜')</p><p>　　subplot(223)</p><p>　　plot(t,Y1)</p><p>　　xlabel('無噪聲解調(diào)信號(hào)波形')</p><p>　　subplot(224)</p><p>

66、;　　plot(frqJ1,J1)</p><p>　　xlabel('無噪聲解調(diào)信號(hào)頻譜')</p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p>　　dyi=sn-S; %高斯白噪聲</p&g

67、t;<p>　　s_ni=var(S)/var(dyi); %輸入信噪比</p><p>　　dyo=YYN-Y1; %解調(diào)后噪聲</p><p>　　s_no=var(Y1)/var(dyo); %輸出信噪比 </p><p>　　dyi1=sn1-S;

68、 </p><p>　　s_ni1=var(S)/var(dyi1); </p><p>　　dyo1=YYN1-Y1; </p><p>　　s_no1=var(Y1)/var(dyo1); </p><p>　　dyi2=sn2-S;

69、 </p><p>　　s_ni2=var(S)/var(dyi2); </p><p>　　dyo2=YYN2-Y1; </p><p>　　s_no2=var(Y1)/var(dyo2); </p><p>　　dyi3=sn3-S;

70、 </p><p>　　s_ni3=var(S)/var(dyi3); </p><p>　　dyo3=YYN3-Y1; </p><p>　　s_no3=var(Y1)/var(dyo3); </p><p>　　dyi4=sn4-S;

71、 </p><p>　　s_ni4=var(S)/var(dyi4); </p><p>　　dyo4=YYN4-Y1; </p><p>　　s_no4=var(Y1)/var(dyo4); </p><p>　　in=[s_ni,s_n

72、i1,s_ni2,s_ni3,s_ni4]; </p><p>　　out=[s_no,s_no1,s_no2,s_no3,s_no4];</p><p>　　set(gcf,'color','w');</p><p>　　figure(3);</p><p>　　plot(in,o

73、ut,'*') </p><p><b>　　hold on</b></p><p>　　plot(in,out)</p><p>　　xlabel('輸入信噪比')</p><p>　　ylabel('輸出信噪比')</p><p>　　3.2.

74、1 調(diào)制信號(hào)幅度=0.8×載波幅度</p><p>　　調(diào)用程序，程序中a=0.8。</p><p>　　調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜如圖12所示：</p><p>　　圖12 調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜圖</p><p>　　解調(diào)信號(hào)的波形、頻譜如圖13所示：</p><p>　　圖13 有噪聲、無

75、噪聲的解調(diào)信號(hào)波形和頻譜</p><p>　　輸入輸出信噪比關(guān)系曲線如圖14所示： </p><p>　　圖14 輸入輸出信噪比關(guān)系曲線</p><p>　　3.2.2 調(diào)制信號(hào)幅度=載波幅度</p><p>　　調(diào)用程序，程序中a=1。</p><p>　　調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜如圖15所示：</p&g

76、t;<p>　　圖15 調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜圖</p><p>　　解調(diào)信號(hào)的波形、頻譜如圖16所示：</p><p>　　圖16 有噪聲、無噪聲的解調(diào)信號(hào)波形和頻譜</p><p>　　輸入輸出信噪比關(guān)系曲線如圖17所示： </p><p>　　圖17 輸入輸出信噪比關(guān)系曲線</p><p>

77、;　　3.2.3 調(diào)制信號(hào)幅度=1.5×載波幅度</p><p>　　調(diào)用程序，程序中a=1.5。</p><p>　　調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜如圖18所示：</p><p>　　圖18 調(diào)制信號(hào)、已調(diào)信號(hào)的波形、頻譜圖</p><p>　　解調(diào)信號(hào)的波形、頻譜如圖19所示：</p><p>　　圖19

78、 有噪聲、無噪聲的解調(diào)信號(hào)波形和頻譜</p><p>　　輸入輸出信噪比關(guān)系曲線如圖20所示： </p><p>　　圖20 輸入輸出信噪比關(guān)系曲線</p><p><b>　　4 心得體會(huì) </b></p><p>　　在此次MATLAB應(yīng)用課程設(shè)計(jì)中我感覺自己收獲很大，無論是查閱資料的能力還是對(duì)MATLAB這個(gè)軟件

79、的掌握都是一個(gè)很大的進(jìn)步。通過對(duì)程序的設(shè)計(jì)，我進(jìn)一步熟悉了MATLAB開發(fā)環(huán)境，對(duì)MATLAB的一些基本操作和應(yīng)用有了更深入的了解。如：有要求的正弦信號(hào)的產(chǎn)生，基本圖形的繪制和各種的函數(shù)的使用等。</p><p>　　剛接到這個(gè)題目真的感到有點(diǎn)束手無策，因?yàn)橐郧爸皇菃渭兊膹臅旧蠈W(xué)習(xí)數(shù)字信號(hào)處理、通信原理的知識(shí)，而這次卻要用MATLAB這個(gè)不熟悉的軟件實(shí)現(xiàn)通信原理中SSB信號(hào)的調(diào)制與解調(diào)。但是，通過不斷查閱相關(guān)資

80、料，許多問題都迎刃而解了。這次設(shè)計(jì)使我對(duì)數(shù)字信號(hào)處理和通信原理課本上學(xué)到的知識(shí)點(diǎn)有了更深入的理解和掌握。比如對(duì)信號(hào)的調(diào)制和解調(diào)過程有了更深層次的理解，學(xué)會(huì)了如何使用MATLAB對(duì)信號(hào)進(jìn)行SSB調(diào)制和解調(diào)，了解了低通濾波器的MATLAB設(shè)計(jì)方法。還有很重要的一點(diǎn)是，我學(xué)會(huì)了如何安排設(shè)計(jì)所需的時(shí)間及合理利用網(wǎng)絡(luò)資</p><p>　　源等普遍實(shí)用的學(xué)習(xí)方法，通過和和同學(xué)探討，拓寬了我的眼界，學(xué)習(xí)了別人好的設(shè)計(jì)思路和設(shè)

81、計(jì)方法等。</p><p><b>　　5 參考文獻(xiàn)</b></p><p>　　[1]劉敏 魏玲 編著.《MATLAB通信仿真與應(yīng)用》.國(guó)防工業(yè)出版社.2007.12 </p><p>　　[2]王立寧 樂光新 詹菲編著.《MATLAB與通信仿真》.人民郵電出版社.2000.01 </p>