關於SDM SNR的求法

jsp0520

我在simulink中建了一個簡單的二階的SDM
, @2 Q2 ^/ B& n
http://www.wretch.cc/album/show.php?i=bleezebub&b=14&f=1513912854&p=0

輸入頻率是500Hz OSR為32

' t1 @& i, j4 Q, _# x. K想請問板上的前輩

6 L3 T% R4 F  y8 L" D. J( V- z6 }%calculate SNR from simulink output
clear;
load simout.mat -ascii;
" V+ D  V4 x  X) j5 H: F  Y' j' X2 B" {
3 z# n: e! U" g1 Posr=32;
9 s& X- Y1 B& o& Q, N  r- H
' \' L/ g. |  {a=0;
for a=1:64001
9 c  T/ a* o; f4 w9 a    x(1,a)=simout(a,1);
& W: f. z! U+ s( b7 Xend
) R3 {2 T" R$ m* Y2 `' G

9 }. q9 ]" m! u2 ]N=length(x)-1;
) k! M) p/ A8 M7 C1 X# Ty=ds_hann(N);
1 v$ r; w) S7 s) A7 As=x(2:N+1).*y;
! Y+ ^, z$ ^  V( P+ Zsigfft=fft(s)./(N/4);
snr_simulink=calculateSNR(sigfft(1:ceil(N/(2*osr))),500,1)

' ?$ J& c$ E4 N- h  E以上這段計算SNR的code有錯嗎
4 j: n0 T! [5 B. s
因為在我改變輸入頻率或OSR的時候
# p( W2 |* D% \: C
+ p5 m' c2 w. f有時候會出現
9 s# b: E- Z- Z
3 B1 |  o% Q- S  L; V; Y" ^! A6 Rsnr_simulink = -Inf
8 v! Q& h. T' }6 D
5 }- p) Y% I7 U請問這是甚麼情況呢
# y( l) E5 T& Z# K& h, d% q/ d
想請高手指點