除小數三階Delta Sigma 頻譜問題

ja1688888

剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
" X: M2 ~9 J% d& X; u手key了前64個輸出到txt檔給matlab做頻譜分析，
9 a% o' S" q  }" y1.JPG
' T* y* @& _. G* `9 y結果卻不像其他論文上的noise shaping
0 a; I# n" f+ \; s) T7 F1 Z2.JPG
$ @: G# N* r. u& w' n7 G6 L以下為matlab code

clear, close all
+ Z( [3 Q9 y1 n5 v  p5 M%% initialize parameters
8 g/ l7 A) q$ ^; u9 c: `* r8 e) L7 fsamplerate=250; % in Hz
N=64; % data length
4 p4 b7 b2 S7 L( l# Z* Msig=load('c:\xshell\d1.txt');
' W2 i, Y5 z) M0 yt=[1:64];
, b9 [0 u* h! J* Uplot(sig(:,1))           
; G3 s) c' Q# K3 y; Hsig_freq=fft(sig);f
PS=abs(sig_freq).^2;
- _9 h2 p+ ?0 n5 nPS=PS/max(PS);
9 n: z# [% J  h4 Nfaxis=samplerate/2*linspace(0,1,N);
* }# K& ]) U& v, }: i3 b# l8 v0 ?
figure,
4 i: \  f5 I$ Jsubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
: K% F1 t4 V5 Vtitle(123)
( v+ B% i. y  H; W7 Zxlabel('Time (s)')
( g4 G( O  `' M3 q% S  b
subplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
* l* c% h' o' Y+ d' X9 ixlabel('Frequency (Hz)')
' G# y" i) k2 ~6 K; ]" _ylabel('Power Spectrum')
. q/ E- d! f$ W8 t

- Y* ?2 H4 \. f請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝