除小數三階Delta Sigma 頻譜問題

ja1688888

! U1 ]& V" ^- Y5 I5 h5 d4 U
剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
; ~7 @$ A: Q5 E4 n# `手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
結果卻不像其他論文上的noise shaping
- {/ S' p# j4 b: P3 Z# L' x# ]7 ~# q2.JPG
以下為matlab code

/ ~! U3 P7 _# P6 T7 M( pclear, close all
/ e1 H2 [. ]5 `" V/ _3 S5 f%% initialize parameters
samplerate=250; % in Hz
N=64; % data length
sig=load('c:\xshell\d1.txt');
: h. a# z6 p( u+ ~' h; v8 Kt=[1:64];
plot(sig(:,1))           
sig_freq=fft(sig);f
PS=abs(sig_freq).^2;
PS=PS/max(PS);
9 G0 d! i; A* b5 @! _0 ~# p) Mfaxis=samplerate/2*linspace(0,1,N);
3 \6 s+ _9 C* \" c; h' _: }) c' Q
2 B! z7 f2 S8 p9 S4 \$ D) ]figure,
+ @  x3 S; M6 }- _6 c6 |subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
, \% X. Q& S0 [4 {3 g& l. i! N! n5 @xlabel('Time (s)')

subplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
! B. {( |- G. {" Q, ztitle('Spectral Analysis (FFT)')
/ G& F2 r1 ]4 c/ Bxlabel('Frequency (Hz)')
  \3 c; O# H, q3 F/ U" xylabel('Power Spectrum')


請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
6 s" w3 k) F: c9 n$ G3 J( e