除小數三階Delta Sigma 頻譜問題

ja1688888

0 ]% D/ s, k. Q0 }剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
- E4 ~7 W. X7 R' ~1 g0 A結果卻不像其他論文上的noise shaping
2 b  S" Y) g: I! R7 d2.JPG
以下為matlab code

5 b6 n' c2 p8 a) Iclear, close all
3 b) N* e$ Z$ a9 @%% initialize parameters
3 Q, \4 o0 k; E9 g) b+ G4 zsamplerate=250; % in Hz
N=64; % data length
9 x0 l; x5 I; k5 V" E: L( Dsig=load('c:\xshell\d1.txt');
& }* e0 B( t8 k  H% k4 ^t=[1:64];
plot(sig(:,1))           
sig_freq=fft(sig);f
+ p( g: r' A8 l1 f; k' ?, i5 @5 RPS=abs(sig_freq).^2;
PS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);

figure,
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
6 W( C  M; N% x/ M. v8 x$ Wxlabel('Time (s)')

2 k6 L/ {& b/ Hsubplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
3 M- i: Q; O: G$ l3 ~title('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
ylabel('Power Spectrum')
  W; n" e; ~9 X* Q( [8 l3 A8 b

7 E6 X! l9 R+ f1 S9 x( O+ p請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
4 G: j- n, h& y* H( R. f