除小數三階Delta Sigma 頻譜問題

ja1688888

% R3 i( w: \1 Q% ^' c$ s
剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
- q/ ]: b' D2 d' v" e手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
結果卻不像其他論文上的noise shaping
& q6 q3 J/ R- p' V) H' m2.JPG
以下為matlab code

/ G2 z3 x( S, }clear, close all
0 F  d7 A, C! f# s! ]%% initialize parameters
samplerate=250; % in Hz
# j8 {1 s5 G0 ^0 {; x% PN=64; % data length
sig=load('c:\xshell\d1.txt');
5 ?! U  |+ l0 s& b! s; {t=[1:64];
7 V2 Y1 J3 M6 L% S; ]5 }( x4 P! Bplot(sig(:,1))           
sig_freq=fft(sig);f
PS=abs(sig_freq).^2;
8 C& l. G' B0 I4 @# zPS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);
% ]: X9 u- ^: E( D( J6 N
figure,
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
' a; f$ \7 Z+ s0 M1 Hxlabel('Time (s)')

1 p* e9 V0 e6 b' P6 esubplot(1,2,2),plot(faxis,PS(1:N))  
- k) x1 D+ O6 L; h3 Cxlim([faxis(1) faxis(end)])
9 m4 W8 Z7 J) y% \) l8 S2 ~title('Spectral Analysis (FFT)')
# R0 B+ d9 P; n: i& t% Y$ c- Rxlabel('Frequency (Hz)')
6 y  k0 }, F6 ~/ p% }0 {ylabel('Power Spectrum')

2 T3 x5 O6 c7 `
# l# B+ e9 m, C6 `& s) ^1 L請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝