除小數三階Delta Sigma 頻譜問題

ja1688888

/ H* n, @% c! L剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
  }6 w* w1 {( R9 z$ B手key了前64個輸出到txt檔給matlab做頻譜分析，
; |/ i& J0 c9 z# C; Z1.JPG
2 C+ K3 Q2 J1 v% u/ n# O7 d結果卻不像其他論文上的noise shaping
* z: ^& M1 ~5 X; k; T0 k2.JPG
# Z. d0 B4 v1 D1 ^- {: N8 h" z% ]以下為matlab code

0 H4 M, y% {) Iclear, close all
%% initialize parameters
samplerate=250; % in Hz
9 V- |5 k+ C0 gN=64; % data length
sig=load('c:\xshell\d1.txt');
7 w; t, p$ |  g4 o  h5 It=[1:64];
plot(sig(:,1))           
sig_freq=fft(sig);f
) y6 U. D& k& q  K9 \PS=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 ^+ r" _1 h9 g4 mxlabel('Time (s)')
+ _4 k. ^0 X( }+ k1 Q$ {
subplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
- M: ^0 C# c3 s4 x3 Wtitle('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
ylabel('Power Spectrum')

* u1 s+ p: R! F* o
請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
; l9 g9 _$ G+ A4 k; G) E