1、<p><b>　　課程設計任務書</b></p><p>　　學生姓名： 專業(yè)班級： </p><p>　　指導教師： 工作單位： </p><p>　　題 目: SSB信號的仿真分析</

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

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

4、　　時間安排</b></p><p>　　第17周，安排任務（鑒主5樓實驗室）</p><p>　　第17-18周，仿真設計（鑒主5樓實驗室）</p><p>　　第19周，完成（答辯，提交報告，演示） </p><p>　　指導教師簽名：

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

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

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

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

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

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

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

15、調制，就是通過某種辦法，只傳送一個邊帶的調制方法。</p><p>　　設調制信號為單頻信號f(t)=Amcosωmt，載波為c(t)=cosωct，則調制后的雙邊帶時域波形為：</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、第一項與調制信號和載波信號的乘積成正比，稱為同相分量；而第二項的乘積則是調制信號與載波信號分別移相90°后相乘的結果，稱為正交分量。</p><p>　　SSB調制框圖如下:</p><p>　　圖1 SSB調制框圖</p><p>　　1.2 SSB解調原理與抗噪性能</p><p>　　解調是調制的逆過程，其作用是從接收的已調信

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

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

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

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

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

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

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

26、B調制的實現(xiàn)</p><p>　　3.1 調制信號為正弦信號</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); %調制信號</p><p><b>　　X=fft(m);</b></p><p>　　X=abs(X(1:length(X)/2+1)

28、); %調制信號頻譜</p><p>　　frqX=(0:length(X)-1)*Fs/length(X)/2</p><p>　　S = modulate(m,Fc,Fs,'amssb'); %對信號進行調制</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 ;%已調信號頻譜</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('調制信號波形')</p><p>　　subplot(222)</p><p>　　plot(frqX,X)</p&g

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

32、t;/p><p>　　subplot(224)</p><p>　　plot(frqY,Y)</p><p>　　axis([0 60000 0 max(Y)])</p><p>　　xlabel('已調信號頻譜')</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'); %無噪聲已調信號解調</p><p>　　YYN=demod(sn,Fc,Fs,'amssb'); %加噪聲已調信號解調</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; %加噪聲后已調信號頻譜</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 %加噪聲后解調信號頻譜</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('加噪聲解調信號波形')</p><p>　　subplot(222)</p><p>　　plot(frqJ2,J2);</p><p>　　axis([0 3000 0 max(J2)])</p><p>　　xlabel('加噪聲解調信號頻

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

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

41、><p>　　dyo=YYN-Y1; %解調后噪聲</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 調制信號幅度=0.8×載波幅度</p><p>　　調用程序，程序中a=0.8。

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

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

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

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

52、t;p>　　圖11 輸入輸出信噪比關系曲線</p><p>　　3.2 調制信號為三角波信號</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); %調制信號 </p><p><b>　　X=fft(m);</b></p><p>　

54、　X=abs(X(1:length(X)/2+1)); %調制信號頻譜</p><p>　　frqX=(0:length(X)-1)*Fs/length(X)/2</p><p>　　S = modulate(m,Fc,Fs,'amssb'); %對信號進行調制</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 ;%已調信號頻譜</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('調制信號波形')</p><p>　　subplot(222)</p><p&

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

58、label('已調信號波形')</p><p>　　subplot(224)</p><p>　　plot(frqY,Y)</p><p>　　axis([0 60000 0 max(Y)])</p><p>　　xlabel('已調信號頻譜')</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'); %無噪聲已調信號解調</p><p>　　YYN=demod(sn,Fc,Fs,'amssb'); %加噪聲已調信號解調</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; %加噪聲后已調信號頻譜</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 %加噪聲后解調信號頻譜</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('加噪聲解調信號波形')</p><p>　　subplot(222)</p><p>　　plot(frqJ2,J2);</p><p>　　axis([0 3000 0 max(J2)])<

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

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

67、t;<p>　　s_ni=var(S)/var(dyi); %輸入信噪比</p><p>　　dyo=YYN-Y1; %解調后噪聲</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 調制信號幅度=0.8×載波幅度</p><p>　　調用程序，程序中a=0.8。</p><p>　　調制信號、已調信號的波形、頻譜如圖12所示：</p><p>　　圖12 調制信號、已調信號的波形、頻譜圖</p><p>　　解調信號的波形、頻譜如圖13所示：</p><p>　　圖13 有噪聲、無

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

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

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

78、 有噪聲、無噪聲的解調信號波形和頻譜</p><p>　　輸入輸出信噪比關系曲線如圖20所示： </p><p>　　圖20 輸入輸出信噪比關系曲線</p><p><b>　　4 心得體會 </b></p><p>　　在此次MATLAB應用課程設計中我感覺自己收獲很大，無論是查閱資料的能力還是對MATLAB這個軟件

79、的掌握都是一個很大的進步。通過對程序的設計，我進一步熟悉了MATLAB開發(fā)環(huán)境，對MATLAB的一些基本操作和應用有了更深入的了解。如：有要求的正弦信號的產生，基本圖形的繪制和各種的函數的使用等。</p><p>　　剛接到這個題目真的感到有點束手無策，因為以前只是單純的從書本上學習數字信號處理、通信原理的知識，而這次卻要用MATLAB這個不熟悉的軟件實現(xiàn)通信原理中SSB信號的調制與解調。但是，通過不斷查閱相關資

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

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