除小數三階Delta Sigma 頻譜問題

ja1688888

剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
手key了前64個輸出到txt檔給matlab做頻譜分析，
; ~3 N1 s9 S" o  ^1.JPG
% j2 Q$ G" h' u1 q/ S$ u3 w2 T結果卻不像其他論文上的noise shaping
2.JPG
- e" c/ u0 c. }: j& N/ O. w4 M* J6 V4 x以下為matlab code

3 D) h. b4 Q# r. Fclear, close all
%% initialize parameters
samplerate=250; % in Hz
N=64; % data length
sig=load('c:\xshell\d1.txt');
t=[1:64];
& g, \& q2 e: V, Fplot(sig(:,1))           
" G; h! R0 M1 {. ~2 {; s) wsig_freq=fft(sig);f
PS=abs(sig_freq).^2;
6 W. D# A% ]4 V7 Z5 n2 \' rPS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);

figure,
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
1 M" K' }4 k2 K' x: ?% p; E+ o5 `title(123)
# x6 M; A+ y9 o2 D# |# i# Y$ Rxlabel('Time (s)')
3 L( i# h. {; n& s
0 t; K) F+ w; V6 H+ i: msubplot(1,2,2),plot(faxis,PS(1:N))  
& w& R  j, M  _3 ?xlim([faxis(1) faxis(end)])
7 q1 _& o) |. ]- C' Ltitle('Spectral Analysis (FFT)')
. y- l/ C8 o5 p; h/ axlabel('Frequency (Hz)')
; i8 B& i* `  a; u; x# C1 P- bylabel('Power Spectrum')


請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝