除小數三階Delta Sigma 頻譜問題

ja1688888

1 F7 X6 ~2 V* x* E
/ @% Q& [2 S& T) A剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
: E3 H" F* V- [( F+ r手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
( P  A. A, i( P, O結果卻不像其他論文上的noise shaping
2.JPG
以下為matlab code

clear, close all
8 F0 A2 `0 q2 ]2 f%% initialize parameters
# G! v; o" _: P$ B/ ssamplerate=250; % in Hz
N=64; % data length
! F* J3 B' n% S: ksig=load('c:\xshell\d1.txt');
t=[1:64];
+ o7 |1 D) r3 m. p% Z3 rplot(sig(:,1))           
sig_freq=fft(sig);f
, d7 s( {1 [% J2 ?; EPS=abs(sig_freq).^2;
& X8 i8 m; C7 ~PS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);
: C- D0 m4 l- l' Q5 O
figure,
9 f( P8 A6 i+ m* \subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
xlabel('Time (s)')
6 T* K2 Z! R3 R; ^( G
subplot(1,2,2),plot(faxis,PS(1:N))  
5 A% E- [1 p3 J+ ]" b, txlim([faxis(1) faxis(end)])
6 Q6 p% Z3 |1 C+ ititle('Spectral Analysis (FFT)')
4 |  V. |) H) S7 U8 ]4 Bxlabel('Frequency (Hz)')
( R2 v* D8 ]( \2 ?& A4 d$ z, Nylabel('Power Spectrum')
" v* D1 o  Z4 ^- v- ~

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