除小數三階Delta Sigma 頻譜問題

ja1688888

  |" V1 w( Z" ]9 v
剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
結果卻不像其他論文上的noise shaping
2.JPG
  n- g( @6 c" [  Q以下為matlab code
& z& W: R% |. _) o6 U6 _
. j  g6 D1 _3 w' p" e& pclear, close all
%% initialize parameters
- J  O6 A: m+ {6 V: {) W/ c  }samplerate=250; % in Hz
4 Q! r, C, v3 CN=64; % data length
% j8 p$ [$ @. h9 u6 y9 \sig=load('c:\xshell\d1.txt');
t=[1:64];
9 U* k$ s5 n% R/ k9 Zplot(sig(:,1))           
% |; W+ l1 s9 t6 Y8 }* S  p7 O% Esig_freq=fft(sig);f
PS=abs(sig_freq).^2;
PS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);

figure,
% }6 D$ H8 m$ \3 K( `/ r- P/ ?/ Psubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
% _, _! s/ n8 A5 T0 z! ttitle(123)
) p! D; p4 ?3 O8 X( sxlabel('Time (s)')

& S& F2 z' I5 U$ G' ysubplot(1,2,2),plot(faxis,PS(1:N))  
& V) R5 }& t" ]+ Rxlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
# x( o: T0 O5 v7 [1 p9 `0 Axlabel('Frequency (Hz)')
/ @& j8 b# [* Zylabel('Power Spectrum')
/ ^" Y, r  F' e, z+ Q
2 l1 l$ }% }. I" A5 w
- ?! }9 c, i" l* @' \) A; R/ k請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝