除小數三階Delta Sigma 頻譜問題

ja1688888

7 J, ^0 _  Q- ~7 ~3 X# w9 f剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
( r6 y# q, F3 [. O5 A& T手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
5 ~- P! c6 g+ r# z結果卻不像其他論文上的noise shaping
2.JPG
以下為matlab code
+ u! S' c, ?+ S6 |
clear, close all
%% initialize parameters
3 \! d! m/ s; Q( T- _( u8 W: }samplerate=250; % in Hz
, f2 H: Y0 |. R; iN=64; % data length
( S* m- [$ ^5 C% k' ssig=load('c:\xshell\d1.txt');
$ z0 k3 [# C# |! I9 Qt=[1:64];
plot(sig(:,1))           
0 o5 b% F1 O8 J; U6 N7 v  Y0 psig_freq=fft(sig);f
PS=abs(sig_freq).^2;
1 U8 {8 @$ N6 M* n, zPS=PS/max(PS);
! j& M. L" N" s7 w/ e7 sfaxis=samplerate/2*linspace(0,1,N);

figure,
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
1 F8 _$ ^  R+ y9 k1 K9 k7 c+ Q) ztitle(123)
7 T3 A4 a, v. w% ~- m" Vxlabel('Time (s)')
- e( X, p' r7 S$ ~$ `0 j
subplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
8 y+ T, P/ x1 j( R$ P" z! Fylabel('Power Spectrum')
+ H) B" Y9 p( R3 w
3 m7 s' o$ H' o4 d% g5 J0 \
請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
  S6 g1 m8 z0 \! B