除小數三階Delta Sigma 頻譜問題

ja1688888

剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
- c+ J4 W. y3 d2 n/ O手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
" t) k4 f& H" o9 e" v, l2 L結果卻不像其他論文上的noise shaping
- w) M( P4 J/ H) j& Y2.JPG
: g% J0 e) G; ?以下為matlab code
% |, b" M9 [( F! ^/ f- P
! a$ n& h6 m7 qclear, close all
( u- h% h. ?/ a) S%% initialize parameters
7 f/ @, c2 M4 z- c1 c; _; \5 E  S! qsamplerate=250; % in Hz
) u7 L1 s$ v- w% {6 P! H2 DN=64; % data length
* Y5 _: v0 s0 c, w8 f/ S  |sig=load('c:\xshell\d1.txt');
% V/ I7 v8 E! h2 L$ ^$ s& t  e( Ht=[1:64];
( X* T5 D% F* d# xplot(sig(:,1))           
sig_freq=fft(sig);f
PS=abs(sig_freq).^2;
* h$ Y* a$ Q0 q, g8 SPS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);
. O5 w4 h6 X! X1 m9 U
3 I$ _$ x0 r! Ffigure,
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
xlabel('Time (s)')

subplot(1,2,2),plot(faxis,PS(1:N))  
8 Q+ y# J; H! cxlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
, o1 D( g& o! ^xlabel('Frequency (Hz)')
# a; A- J  ?$ Hylabel('Power Spectrum')
% K9 q5 E: _0 ?
5 {, v# q+ ^. P1 X6 X5 X: x( X
, K6 v/ r2 e  ~" H請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
+ u" e9 u0 W+ u1 r$ i6 f3 S