除小數三階Delta Sigma 頻譜問題

ja1688888

剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
  f5 y+ q7 O" x& G! d4 m$ w# j  x手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
結果卻不像其他論文上的noise shaping
2.JPG
0 h* _; l$ i3 a! r以下為matlab code
) V& ~1 P* h8 q) v, w8 ]- I
( P% l* t! X6 m0 ^& Gclear, close all
% P# Q. ^3 y" w1 F3 m%% initialize parameters
samplerate=250; % in Hz
N=64; % data length
sig=load('c:\xshell\d1.txt');
& L- R5 z. u% m1 X" `t=[1:64];
+ }$ F6 L3 Y+ }6 y. z) kplot(sig(:,1))           
sig_freq=fft(sig);f
PS=abs(sig_freq).^2;
PS=PS/max(PS);
+ [! k) ~- k: Z5 ^( wfaxis=samplerate/2*linspace(0,1,N);
4 `  H8 L, x! b, ~) I
figure,
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))  
xlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
4 j' X! X& h  Vylabel('Power Spectrum')
( E( _0 g# s1 B+ x8 t6 u3 T9 D6 T
' a8 z1 a0 h/ s; w
請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝