除小數三階Delta Sigma 頻譜問題

ja1688888

$ s- u( ]$ B9 p" {! T
剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
結果卻不像其他論文上的noise shaping
5 U( F; M2 K1 V- L7 u$ g2.JPG
9 R" C' F/ U: C% h9 O. Q9 _; p以下為matlab code
2 `0 C+ k; |; U/ h' U1 ~& q8 }" f* f$ S
clear, close all
4 m7 o  i4 ~/ J4 ~%% initialize parameters
0 R& }- ]3 K8 j/ ^8 D' Wsamplerate=250; % in Hz
N=64; % data length
sig=load('c:\xshell\d1.txt');
t=[1:64];
plot(sig(:,1))           
% x9 ~. R2 y: H2 n3 V, Ssig_freq=fft(sig);f
' x% o7 D! w2 |3 aPS=abs(sig_freq).^2;
PS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);

figure,
% a) i  ~6 S: T0 x% B# n2 @subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
' D% H( {  b! p) ?2 l  W: g) Atitle(123)
) r- N0 H) q0 m' _2 f0 R9 zxlabel('Time (s)')

* x8 g" c  z& P8 j# ~, fsubplot(1,2,2),plot(faxis,PS(1:N))  
7 `6 r% V; w: R3 H: ?# x# Z8 |xlim([faxis(1) faxis(end)])
' C" p2 F3 Q* ~9 E- d7 Y2 ^8 u# ktitle('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
ylabel('Power Spectrum')
" v- Y* {0 z! H6 I& V
# @% w0 u4 c8 L; r& j
2 X& m) u& o" F# l請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
, f$ s4 g+ s& e  q/ j1 X6 {' z