除小數三階Delta Sigma 頻譜問題

ja1688888

8 x2 |6 z& ]- g5 r4 V3 k0 o剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
3 e) @% \% @6 s' T; Q" c手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
5 p& G5 i$ t8 A. Z結果卻不像其他論文上的noise shaping
! {' V& D# @  Y9 l2.JPG
0 @& R: n$ ^: r# D7 G以下為matlab code
% T- H9 E, A$ w6 M
clear, close all
' ~  Q" s9 Q1 Z2 q, F%% initialize parameters
/ r" U) N4 e) d: G' [" Tsamplerate=250; % in Hz
N=64; % data length
sig=load('c:\xshell\d1.txt');
t=[1:64];
plot(sig(:,1))           
sig_freq=fft(sig);f
) Q/ X8 _' T# O0 `+ u5 F! fPS=abs(sig_freq).^2;
* `9 j4 P, G( k  c" APS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);
- u# j' j7 v5 k$ [7 j. X
figure,
0 Z, @8 o0 B' N# s, ksubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
0 J/ I! U' W( Z2 X2 A# @7 yxlabel('Time (s)')

subplot(1,2,2),plot(faxis,PS(1:N))  
6 F! X- T& \; G* rxlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
ylabel('Power Spectrum')
( H* k2 W% |. J% R  C& B
/ N" X1 y9 u8 s0 @
請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
8 {, \4 I) o* T+ J! Z( a' i& `  [