有關Simulink

kevin

power spectrum estimation 有幾種方式,
1 b# e2 ]) v+ v7 xPeroidogram
; Q8 k. ^+ D  z+ v$ v1 sModified periodogram
Bartlett method
* Z9 p% l/ c* qWelch method
Blackman-Tukey method
/ i5 |" h; m8 Z% w, |
若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
n1  : starting index, x(n1)
n2  : ending index, x(n2)

function Px = per(x,n1,n2)
x = x(;
/ C2 q# t7 [/ }1 xif nargin == 1
- k  k& K- s, _0 `    n1 = 1;  n2 = length(x);  end;
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
Px(1) = Px(2);

kevin

power Spectrum 跟 SNR 是不同的圖.

: g& Y2 S$ y6 n5 A

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
! N* s! k9 W4 g" E& k8 ~6 Qspec=fft(v.*ds_hann(N))/(N/4);
plot(linspace(0,1,N/2),
# G( U: a# o. L  p5 M) ~" jdbv(spec(1:N/2)))
( @# Y& p/ I9 n7 _' Caxis([0 1 -120 0]);
grid on;
. g7 o1 Z1 l& s& Xylabel('dB')
1 y- H# O2 j2 }* T, G& N: R
! J' ?' F' Z. q0 E! ?; F從左邊這張圖可以去算出屬於它的SNR,參考附件.

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.

Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.