除小數三階Delta Sigma 頻譜問題

ja1688888

剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
結果卻不像其他論文上的noise shaping
2.JPG
( f  `( ^, E% s% T, y# K/ J  X0 m以下為matlab code

7 J# {# ?: c6 L) U3 a/ t  Fclear, close all
: a6 u5 `9 p1 ^2 v- N%% initialize parameters
samplerate=250; % in Hz
N=64; % data length
# f: U9 ^1 H- Tsig=load('c:\xshell\d1.txt');
0 O" f7 [1 f5 C0 d- b+ F3 Lt=[1:64];
) P" o+ G; C0 Hplot(sig(:,1))           
sig_freq=fft(sig);f
PS=abs(sig_freq).^2;
PS=PS/max(PS);
4 h8 L  c" s, j, R$ `1 mfaxis=samplerate/2*linspace(0,1,N);

" W1 @* @4 u$ N" mfigure,
. o+ J. Q- d4 E; ?/ [+ m+ T, K  asubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
xlabel('Time (s)')
! }% ~; J  j8 a  o" [! o, Z
subplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
' k  m  X8 O. Z  W9 Rylabel('Power Spectrum')
% i2 x7 X5 v) d2 D0 _
8 }1 {: H& b4 Y9 {/ D' O
5 L3 ?  a: J& v# Q- Y1 X請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝