除小數三階Delta Sigma 頻譜問題

ja1688888

* O# D( A9 b' P& v0 u" e剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
7 R6 o$ S* ~6 S7 d手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
  g! s  d5 m: f& y+ T結果卻不像其他論文上的noise shaping
2.JPG
以下為matlab code
7 z( R9 ?' Z5 D/ G0 m: e# N
, r& A! z$ P( H& T' I/ k2 |- F/ {clear, close all
%% initialize parameters
samplerate=250; % in Hz
N=64; % data length
sig=load('c:\xshell\d1.txt');
t=[1:64];
- H4 b  O' V. R" d6 Vplot(sig(:,1))           
sig_freq=fft(sig);f
( M5 n# b0 q* o1 x0 |' ~6 h' W) mPS=abs(sig_freq).^2;
PS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);

figure,
; _& h; V) M) J9 Q8 ?* B, i* Q$ K( nsubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
! |( y8 u3 J! b# `. @( |( H  Dxlabel('Time (s)')

/ [# f: Z# H" x* V2 Ysubplot(1,2,2),plot(faxis,PS(1:N))  
5 G# ?9 ^. E* `! cxlim([faxis(1) faxis(end)])
" \0 Q, j. C" r% O1 }title('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
: [6 q/ y# H+ ]. Jylabel('Power Spectrum')

" i, q2 ~5 U' Z6 t5 f
請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝