除小數三階Delta Sigma 頻譜問題

ja1688888

, _! l% I8 T  x7 a
2 K6 T4 ^; ?) U. F$ O剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
結果卻不像其他論文上的noise shaping
* }; {. y1 V' @+ Q; V2.JPG
% `5 r/ G1 `* [; P1 l# y以下為matlab code

9 O4 T8 e# z) Jclear, close all
) E7 F9 R+ i7 X! d: R3 H: d%% initialize parameters
. R0 E6 A' v- I5 G( g  |+ Usamplerate=250; % in Hz
" x. O0 U% j# ~- lN=64; % data length
9 ], e( `0 P' G* D; R3 O! csig=load('c:\xshell\d1.txt');
* N& \- a- P# p3 l) X5 at=[1:64];
plot(sig(:,1))           
sig_freq=fft(sig);f
PS=abs(sig_freq).^2;
PS=PS/max(PS);
! G, a& K" E6 q1 T- W- `faxis=samplerate/2*linspace(0,1,N);
$ J; i0 p  t: e3 P( F. d
; t5 c. k( R, m0 D* rfigure,
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
xlabel('Time (s)')
- G: U. p3 e, @/ t
6 Z' u& l1 o2 R# {6 r; r- }subplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
ylabel('Power Spectrum')

3 E! h9 E- z% ~) Z5 `
, i& K' G6 |% _  {請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝