關於SDM SNR的求法

jsp0520

我在simulink中建了一個簡單的二階的SDM

http://www.wretch.cc/album/show.php?i=bleezebub&b=14&f=1513912854&p=0
' P- n! a8 f3 H! q
5 j$ h! \5 B0 c* _輸入頻率是500Hz OSR為32

; ~$ ]9 z3 i4 @+ e( V想請問板上的前輩

%calculate SNR from simulink output
  K7 w. e- `7 C$ r, Y- S. sclear;
& h' z5 J; l8 t: Z' Yload simout.mat -ascii;
3 k/ P/ A5 A! r# N" N: y) a; p
! l2 r- S; n9 D  a  I7 v: _3 cosr=32;
: s7 c! g3 I- j) o% Z' W
( ]3 b, ?7 X3 @- \. da=0;
for a=1:64001
    x(1,a)=simout(a,1);
end
6 Q0 y% f5 s  p: o7 g+ b+ z7 q

$ [$ F4 T/ J/ `$ jN=length(x)-1;
y=ds_hann(N);
s=x(2:N+1).*y;
3 Q  U8 c/ q8 q9 p" k& {# psigfft=fft(s)./(N/4);
snr_simulink=calculateSNR(sigfft(1:ceil(N/(2*osr))),500,1)

以上這段計算SNR的code有錯嗎

因為在我改變輸入頻率或OSR的時候
, J$ w( _. ~2 D& E2 T! n
* r# }6 U0 ~+ c; r有時候會出現
7 c6 P) M2 O1 e* D" Z8 \/ i; I
( `: B/ ^9 Q& h2 s- t1 X. Tsnr_simulink = -Inf

0 I9 G. N( J- `3 h請問這是甚麼情況呢
9 I& _/ N! Z( f/ _+ d8 k, s
想請高手指點