除小數三階Delta Sigma 頻譜問題

ja1688888

+ F4 @4 F# ^1 n! t. V
. D0 h( Z& v" v7 r' b剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
: I5 I7 l  u  @. @  H' B, e2 {$ y手key了前64個輸出到txt檔給matlab做頻譜分析，
# X6 X# v. x4 S( h: G1.JPG
- W1 h5 O3 n1 k. q% a結果卻不像其他論文上的noise shaping
2.JPG
. w  e; u" S: H+ Q! p以下為matlab code

0 _1 I+ n, `3 p# ?clear, close all
. l. t8 D4 r; {+ @  a%% initialize parameters
samplerate=250; % in Hz
. W& b; Z  r: x0 u! j) IN=64; % data length
! f9 R2 l( [) y+ D, `' B6 D9 H# ysig=load('c:\xshell\d1.txt');
t=[1:64];
plot(sig(:,1))           
sig_freq=fft(sig);f
PS=abs(sig_freq).^2;
PS=PS/max(PS);
. W% B3 f( H+ k7 X0 u& X$ rfaxis=samplerate/2*linspace(0,1,N);
$ o2 a0 d& e* v+ x9 z8 O' l
% E4 [/ ?9 s8 k; K! T& rfigure,
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
/ n: H, ^- u5 ttitle(123)
xlabel('Time (s)')
/ g% X& @) s6 a
! A! Z' `( u0 U1 hsubplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
6 w: S; h+ L- }& ^: p% g- \title('Spectral Analysis (FFT)')
" {; c1 l* E0 p8 k! dxlabel('Frequency (Hz)')
ylabel('Power Spectrum')

+ O. H- f( U& Z% g, B' g8 x- j
: S' q) s; J8 ]' U  D0 n, U" n請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝