程建峰

(甘肅工業(yè)職業(yè)技術學院，甘肅 天水 741025)

?

基于MATLAB的連續(xù)信號的采樣與重構仿真分析

程建峰

(甘肅工業(yè)職業(yè)技術學院，甘肅 天水 741025)

摘要針對連續(xù)信號的采樣與重構問題，利用MATLAB仿真軟件平臺，仿真不同條件下連續(xù)信號的采樣信號時域波形和采樣后信號頻譜、重構信號時域波形和重構后誤差波形圖。通過對采樣周期對采樣頻譜疊加和信號重構精度的影響、以及信號被采樣前后在頻域的變化對比分析，得出在不同采樣頻率的條件下，對應采樣信號的時域、頻域特性以及重構信號與誤差信號也隨之產生變化，連續(xù)信號可以完全恢復過來。

關鍵詞MATLAB；信號采樣；頻譜疊加；信號重構

0引言

連續(xù)信號的數(shù)據(jù)采樣遵循采樣定律[1-2]，滿足ws≥2wm，其中ws為采樣頻率，wm為被測信號的上限截止頻率。但在實際操作中，由于被測信號的復雜性和其他因素的限制，合理選取采樣頻率會存在一定的困難，由此產生的誤差不易發(fā)現(xiàn)。本文利用MATLAB軟件[3-4]仿真平臺再現(xiàn)了采樣信號與頻譜，分析了兩者關系，以此探討影響采樣周期的因素以及信號重構的可行性。

1信號采樣與信號重構

1.1　信號采樣

理想信號采樣器可以看成是一個載波為理想單位脈沖序列δT(t)的幅值調制器，即理想采樣器的輸出信號e*(t)是連續(xù)輸入信號e(t)調制在載波δT(t)上的結果[5]。

用數(shù)學表達式描述理想采樣器調制過程，則有e*(t)=e(t)δT(t)，理想單位脈沖序列δT(t)可以表示為：

式中，δ(t-nT)是出現(xiàn)在時刻t=nT、強度為1的單位脈沖。由于e(t)的數(shù)值僅在采樣瞬時才有意義，同時，假設e(t)=0?t<0，所以e*(t)又可表示為：

1.2　信號重構

設信號f(t)被采樣后形成的采樣信號為fs(t)，信號的重構是指由fs(t)經過內插處理后，恢復出原來信號f(t)的過程。若設f(t)是帶限信號，帶寬為ωm，經采樣后的頻譜為Fs(jω)。設采樣頻率ωs≥2ωm，現(xiàn)選取一個頻率特性為：

顯然，F(xiàn)(jω)=Fs(jω)H(jω)，與之對應的時域表達式為：

f(t)=h(t)*fs(t)，

將h(t)及fs(t)代入得：

上式為f(nTs)求解f(t)的表達式。

通過分析，得到時域采樣定理：一個帶寬為ωm的帶限信號f(t)，可唯一地由它的均勻取樣信號fs(nts)確定，其中，取樣間隔ts<π/ωm,該取樣間隔又稱為奈奎斯特間隔。根據(jù)時域卷積定理，求出信號重構的數(shù)學表達式為：

2仿真與分析

連續(xù)信號是指自變量的取值范圍是連續(xù)的，且對于一切自變量的取值，除了有若干個不連續(xù)點以外，信號都有確定的值與之對應。MATLAB并不能處理連續(xù)信號，而是用等時間間隔點的樣值來近似表示連續(xù)信號。當取樣時間間隔足夠小時，這些離散的樣值就能較好地近似連續(xù)信號。時域對連續(xù)時間信號進行采樣，是給它乘以一個采樣脈沖序列，就可以得到采樣點上的樣本值，信號被采樣前后在頻域的變化，可以通過時域頻域的對應關系分別求得采樣信號的頻譜[8]。

2.1　Sa(t)的臨界采樣及重構

當采樣頻率小于一個連續(xù)的同信號最大頻率的2倍,即ws=2ωm時，稱為臨界采樣。修改門信號寬度和采樣周期等參數(shù)，得到的采樣信號時域和頻域特性，以及重構信號與誤差信號的變化，如圖1所示。

圖1　Sa(t)的臨界采樣及重構圖

為了比較由采樣信號恢復后的信號與原信號的誤差，可以計算出2個信號的絕對誤差。當t選取的數(shù)據(jù)越大，起止的寬度越大。

2.2　Sa(t)的過采樣及重構分析

當采樣頻率大于一個連續(xù)的同信號最大頻率的2倍，即ws>2ωm時，稱為過采樣。在不同采樣頻率的條件下，觀察對應采樣信號的時域和頻域特性，以及重構信號與誤差信號的變化，如圖2所示。

圖2　Sa(t)的過采樣信號、重構信號及兩信號的絕對誤差

從圖2可看出將原始信號分別修改為抽樣函數(shù)Sa(t)、正弦信號sin(20*π*t)+cos(20*π*t)、指數(shù)信號e-2tu(t)時，在不同采樣頻率的條件下，可以觀察到對應采樣信號的時域和頻域特性，以及重構信號與誤差信號的變化[9]。

2.3　Sa(t)的欠采樣及重構分析

當采樣頻率小于一個連續(xù)的同信號最大頻率的2倍，即ws<2ωm時，稱為過采樣。利用頻域濾波的方法修改實驗中的部分程序，完成對采樣信號的重構，結果如圖3所示。

圖3　Sa(t)的欠采樣信號、重構信號及兩信號的絕對誤差

通過圖3可以看出，絕對誤差error已大為增加，其原因是因采樣信號的頻譜混疊，使得在|ω|<ωc區(qū)域內的頻譜相互“干擾”所致。

3結束語

利用MATLAB所生成的臨界采樣、過采樣和欠采樣3種狀態(tài)下的采樣信號、重構信號及兩信號的絕對誤差圖，可以直觀地看到在不同采樣頻率的條件下，對應采樣信號的時域、頻域特性以及重構信號與誤差信號也隨之產生變化，一個連續(xù)時間信號完全可以用該信號在等時間間隔上的瞬時值來表示[10]，并且可以用這些樣本值把信號完全恢復過來。

參考文獻

[1]梁虹，梁潔，陳躍斌．信號與系統(tǒng)分析及其MATLAB實現(xiàn)[M].北京：電子工業(yè)出版社，2002.

[2]周利清，蘇菲． 數(shù)字信號處理基礎[M].北京：北京郵電大學出版社，2005.

[3]朱衡君．MATLAB語言及實踐教程[M].北京：清華大學出版社,北京交通大學出版社，2009.

[4]駱文，趙強．MATLAB在信號與系統(tǒng)中的應用[J].數(shù)字技術與應用，2010(8)：23-25.

[5]陳永勝．基于信號頻域特征實現(xiàn)位定時恢復[J].無線電通信技術，2014，40(2)：94-96.

[6]劉可．通用信號處理模塊功能線程動態(tài)重構技術[J].無線電工程，2014，44(11)：77-80.

[7]余愷，李元實，王智．基于壓縮感知的聲音信號采集方法[J].儀器儀表學報，2012(1)：105-112.

[8]林從客，謝求成．一個時限信號采樣定理及其在數(shù)據(jù)壓縮中的應用[J].南京航空航天大學學報，1997(4)：29-45.

[9]夏輝，靳曉艷．相關運算在信號處理中的應用[J].無線電工程，2014,44(9)：67-69.

[10]艾名舜．多線性調頻信號瞬時頻率估計快速算法[J].無線電通信技術，2014，40(6)：36-39.

程建峰女，(1979—)，講師。主要研究方向：通信和EDA技術。

Simulation and Analysis of Continuous Signal Sampling and

Reconstruction Based on MATLAB

CHENG Jian-feng

(GansuIndustryPolytechnicCollege,TianshuiGansu741025,China)

AbstractTo address the issue of continuous signal sampling and reconstruction,the continuous signal’s sampling signal time domain waveform and signal spectrum after sampling,the reconstruction signal time domain waveform and the error waveform after reconstruction are simulated under different conditions using the MATLAB simulation software platform.Through an analysis on the influence of sampling period on the sampling spectral overlay and signal reconstruction accuracy,and through a comparison between the signals before and after sampling in frequency domain,it is concluded that under different sampling frequencies,the corresponding sampling signal’s time domain,frequency domain characteristics,and the reconstruction signal and the error signal also vary,and the continuous signal can be completely recovered.

Key wordsMATLAB;signal sampling;spectrum overlay;signal reconstruction

作者簡介

基金項目：甘肅工業(yè)職業(yè)技術學院項目(Gzy2011-27)。

收稿日期：2015-03-13

中圖分類號TN911.7

文獻標識碼A

文章編號1003-3106(2015)06-0035-03

doi:10.3969/j.issn.1003-3106.2015.06.10

引用格式：程建峰.基于MATLAB的連續(xù)信號的采樣與重構仿真分析[J].無線電工程，2015，45(6)：35-37.