除小數三階Delta Sigma 頻譜問題

ja1688888

/ w2 a8 F) n- y1 x% D
剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
手key了前64個輸出到txt檔給matlab做頻譜分析，
# _& J9 F/ r* |) g0 `% T1.JPG
) L$ \. e& v+ G結果卻不像其他論文上的noise shaping
2.JPG
以下為matlab code
- A& @/ ?" R' v: y5 b; P8 l# p
clear, close all
%% initialize parameters
8 \: v, G" C6 L& nsamplerate=250; % in Hz
N=64; % data length
sig=load('c:\xshell\d1.txt');
" j1 F. j7 p( kt=[1:64];
plot(sig(:,1))           
sig_freq=fft(sig);f
PS=abs(sig_freq).^2;
PS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);

8 q0 M5 q( t, M, }: Mfigure,
; P& P/ p2 }/ Z: ~: D+ e& o& f1 bsubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
" z+ x* H2 `8 P8 Hxlabel('Time (s)')
; ?9 y. a( v! t
subplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
( b9 P# z+ ]; V  [: htitle('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
+ V- o' f# ]8 u& q) ~- ^  D! u- `ylabel('Power Spectrum')

- J7 U1 f4 p" [, P* g! t2 q
/ o- v6 m- n1 @# u; t' G* P+ B3 R請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
9 @5 a, Z2 l5 \