卷積和及取樣定理

2021-07-11 14:22:45 字數 3687 閱讀 1764

乙個頻帶限制在(0

,fc)

赫茲內的模擬訊號m(

t),若以取樣頻率fs

≥2fc

對模擬訊號m(

t)進行取樣,得到最終的取樣值,則可無混疊失真地恢復原始模擬訊號m(t)。

其中,無混疊失真地恢復原始模擬訊號m(

t)是指被恢復訊號與原始模擬訊號在頻譜上無混疊失真,並不是說被恢復訊號與原始訊號在時域上完全一樣。由於取樣和恢復器件的精度限制以及量化誤差等存在,兩者實際是存在一定誤差或失真的。

奈奎斯特頻率:通常把最低允許的取樣頻率fs

=2fc

稱為奈奎斯特頻率。

取樣定理的幾點小結:

1. 乙個帶限模擬訊號 xa

(t) ,其頻譜的最高頻率為 fc

,以間隔ts

得到取樣訊號x̂ 

a(t)

,當且僅當fs

≥2fc

時,x̂ 

a(t)

才可不失真的恢復xa

(t) 。

2. x̂ a

(t) 的頻譜 x̂ 

a(jω

) 是模擬訊號xa

(t) 的頻譜xa

(jω)

以ω=2

πfs 為週期進行週期延拓得到的。

3. 一般稱fs

/2為摺疊頻率,只要訊號的最高頻率不超過這個頻率,就不會出現頻譜混疊現象,否則,超過fs

/2的頻譜會「摺疊」回來形成混疊現象。

分析連續時間訊號的時域波形及其幅頻特性曲線,原訊號: f(

x)=0.5si

n(2π

∗65t)

+0.8co

s(2π

∗40t)

+0.7co

s(2π

∗30t)

對原訊號進行取樣,得到取樣序列,對不同取樣頻率下的取樣序列進行頻譜分析,由取樣序列恢復出連續時間訊號,畫出其時域波形,對比與原連續時間訊號的時域波形。

matlab**如下:

function

fz = cy

(fy,fs)

%實現取樣頻譜分析繪圖函式

%fy:原訊號函式,fy以字串的格式輸入

%fs:取樣頻率

fs0 = 1e4;

tp = 0.1;

t = [-tp:1/fs0:tp];

k1 = 0:999; k2 = -999:-1;

m1 = length(k1); m2 = length(k2);

f = [fs0*k2/m2,fs0*k1/m1];%設定原訊號的頻率陣列

w = [-2*pi*k2/m2,2*pi*k1/m1];

fx1 = eval(fy);%獲取取樣序列

fx1 = fx1*exp(-j*[1:length(fx1)]'*w);%求原訊號的離散時間傅利葉變換

%畫原訊號

figure

subplot(211),plot(t,fx1,'r')

title('原訊號');

xlabel('t(s)'),ylabel('y(t)');

axis([min(t),max(t),min(fx1),max(fx1)])

grid on

%畫原訊號頻譜

subplot(212),plot(f,abs(fx1),'r');

title('原訊號頻譜');

xlabel('f(hz)'),ylabel('fx1');

axis([-100,100,0,max(abs(fx1))+5]);

grid on

%對訊號取樣

ts = 1/fs;%取樣週期

t1 = -tp:ts:tp;%取樣時間序列

f1 = [fs*k2/m2,fs*k1/m1];%設定取樣訊號的頻率陣列

t = t1;

fz = eval(fy);%獲取取樣序列

fz = fz*exp(-j*[1:length(fz)]'*w);%取樣訊號的離散時間傅利葉變換

%畫取樣序列波形

figure

subplot(211),stem(t,fz,'.');

title('取樣訊號')

xlabel('t(s)'),ylabel('y(t)')

line([min(t),max(t)],[0,0]);

grid on

%畫取樣訊號頻譜

subplot(212),plot(f1,abs(fz),'m');

title('取樣訊號頻譜')

xlabel('f(hz)'),ylabel('fz')

grid on

end

function

fh = hf

(fz,fs)

%訊號的恢復及頻譜函式

%fz:取樣序列

%fs:取樣頻率

t = 1/fs;

dt = t/10; tp = 0.1;

t = -tp:dt:tp; n = -tp/t:tp/t;

tmn = ones(length(n),1)*t-n'*t*ones(1,length(t));%t-nt

fh = fz*sinc(fs*tmn);%由取樣訊號恢復原訊號

k1 = 0:999; k2 = -999:-1;

m1 = length(k1);m2 = length(k2);

w = [-2*pi*k2/m2,2*pi*k1/m1];

fh = fh*exp(-j*[1:length(fh)]'*w);%恢復後的訊號的離散時間傅利葉變換

figure

%恢復訊號波形

subplot(211),plot(t,fh,'g');

st1 = sprintf('由取樣頻率fs=%d',fs);

st2 = '恢復後的訊號';

st =[st1,st2];

title(st);

xlabel('t(s)'),ylabel('y(t)');

axis([min(t),max(t),min(fh),max(fh)])

line([min(t),max(t)],[0,0])

grid on

%畫重構訊號的頻譜

f = [10*fs*k2/m2,10*fs*k1/m1];

subplot(212),plot(f,abs(fh),'g')

title('恢復後訊號的頻譜');

xlabel('f(hz)'),ylabel('fh');

axis([-100

1000 max(abs(fh))+2]);

grid on

end

程式結果如下:

原訊號及原訊號頻譜圖

max 時,為原訊號的欠取樣訊號和恢復,即不滿足取樣定理,出現頻譜混疊現象。

低通取樣和帶通取樣定理

耐奎斯特取樣定理 2 我們假設有乙個 時間連續訊號 x t 的頻帶在 0,f h 之間,以採 樣速率為連續訊號頻率 2 倍 f s 2f h 的採 樣速率對 x t 進行等間隔取樣,得到時間離 散的取樣訊號 x n x nt s 其中 ts 1 f s為取樣間隔 則原始訊號 x t 將被所得到的採 ...

關於取樣定理

取樣定理,又稱夏農取樣定理,奈奎斯特取樣定理,是 資訊理論,特別是通訊與 訊號處理 學科中的乙個重要基本結論 e.t.whittaker 1915 年發表的統計理論 克勞德 夏農 與harry nyquist 都對它作出了重要貢獻。另外 v.a.kotelnikov 也對這個定理做了重要貢獻。取樣是...

夏農取樣定理

編輯夏農取樣定理,又稱奈奎斯特取樣定理,是資訊理論,特別是通訊與訊號處理學科中的乙個重要基本結論。1924年奈奎斯特 nyquist 就推導出在理想低通訊道的最高大碼元傳輸速率的公式 理想低通訊道的最高大碼元傳輸速率 2w log 2n 其中w是理想低通訊道的頻寬,n是電平強度 中文名 夏農取樣定理...