除小數三階Delta Sigma 頻譜問題

ja1688888

, c$ j+ V: B& m1 N
8 J$ I  x7 H) h: m0 e+ B0 v* o剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
8 P( H6 _7 t1 v, v" L) x手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
結果卻不像其他論文上的noise shaping
* A# {! |- [( s' c/ ^2.JPG
以下為matlab code

clear, close all
9 t; q6 I0 _4 ^+ V. y1 {%% initialize parameters
, f/ @8 r* O. L  wsamplerate=250; % in Hz
7 i6 K' |3 r3 }N=64; % data length
2 N1 d6 O5 v0 C- O4 ^1 }sig=load('c:\xshell\d1.txt');
# R5 x* c7 U0 b4 Ht=[1:64];
9 o9 e5 c  L- h' n" y0 S8 bplot(sig(:,1))           
sig_freq=fft(sig);f
PS=abs(sig_freq).^2;
PS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);

6 C/ ~' X9 n/ A# Nfigure,
0 O3 T  i! Q  o3 {2 D# |' \8 L# P# S4 hsubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
, r0 W( A: @: j, l& T3 E9 L" L$ nxlabel('Time (s)')
. c( ~$ K, j  f% }
- w' }  `! t" n) Xsubplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
  P7 C# J4 e7 [: {" l+ {ylabel('Power Spectrum')
2 e4 J$ g% d- g! n

請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
/ x) R! I7 e- p' e