除小數三階Delta Sigma 頻譜問題

ja1688888

8 \; a  k- S" ~. ]" Q剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
8 C. V0 X+ U) W" F" R. L結果卻不像其他論文上的noise shaping
4 R% N# F! K0 T. I, `3 R" O& N2.JPG
以下為matlab code
, A, A. z* y+ H& v
clear, close all
$ D5 x- y3 e. @3 P+ p%% initialize parameters
samplerate=250; % in Hz
N=64; % data length
+ ?& x! s# y- e" s% u2 b; b; Bsig=load('c:\xshell\d1.txt');
0 N9 u5 }" K# a: o* \; qt=[1:64];
plot(sig(:,1))           
sig_freq=fft(sig);f
& a. L2 g% t0 E6 r. V! y* ePS=abs(sig_freq).^2;
. i6 T+ _, h" R7 X+ t* C2 ?PS=PS/max(PS);
, U8 W9 ]5 K6 Hfaxis=samplerate/2*linspace(0,1,N);

0 V: {: z2 v, V" C+ Xfigure,
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
xlabel('Time (s)')

subplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
; C* _" D1 h3 ?xlabel('Frequency (Hz)')
& m3 G. h" [# fylabel('Power Spectrum')

7 j7 `9 T. g0 K# f2 `
+ t4 Q0 S) K+ f請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
( l. k& O$ t0 @" s5 c) L