除小數三階Delta Sigma 頻譜問題

ja1688888

) C  Z2 y0 }$ m! ?( p" B
$ j; V+ y6 [- {* `剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
結果卻不像其他論文上的noise shaping
2.JPG
以下為matlab code
0 r& e" N. V1 {# M
: c/ y2 _% g2 i" e% kclear, close all
7 }: g# E+ P3 ]; Y; @/ o. ]4 `%% initialize parameters
samplerate=250; % in Hz
  W5 H6 i8 D% s2 p/ S7 g1 L0 ~2 hN=64; % data length
& r# Q2 U1 q0 r# ]5 P& |sig=load('c:\xshell\d1.txt');
! S# n+ C% o, T9 ht=[1:64];
2 ~' j1 }  Y! L: g* ~6 |7 ?plot(sig(:,1))           
8 C- i5 y9 A- a* R* S. Xsig_freq=fft(sig);f
PS=abs(sig_freq).^2;
5 v* n' [! O/ b7 F  g# C% C( G  `# mPS=PS/max(PS);
9 l3 ~3 R& K& d$ \* yfaxis=samplerate/2*linspace(0,1,N);

9 {! I& L/ S; `  r! |+ bfigure,
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
% b# g, ^% C9 B2 Y. Qtitle(123)
5 Q( t( {. w) P" c" lxlabel('Time (s)')

subplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
' G+ X3 u$ r" rylabel('Power Spectrum')
; q1 K- ~# a: Z: Z5 E3 G

' p: H) V7 L' g) O3 ?- \請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
2 N& g' L5 [5 C* \- j  ]) J