除小數三階Delta Sigma 頻譜問題

ja1688888

; r" B/ {" ~' t, s
: G6 |" R2 f( f8 d8 |剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
% m8 [# S& u$ ~8 y+ |# L% f: q手key了前64個輸出到txt檔給matlab做頻譜分析，
4 ?2 p4 N3 y/ F0 \. V: Y1.JPG
結果卻不像其他論文上的noise shaping
- l$ a9 A8 U- w; Q" A% h" _2.JPG
以下為matlab code

! U  n4 M' y5 J8 B* X# V! Zclear, close all
- Z% a+ V! x7 Z# e%% initialize parameters
samplerate=250; % in Hz
3 D1 H! X) g9 }% s. X" w# K+ N+ ~, z! aN=64; % data length
sig=load('c:\xshell\d1.txt');
. z9 N4 k% b/ n1 nt=[1:64];
8 x; m) |8 F$ q* {6 `+ ?plot(sig(:,1))           
sig_freq=fft(sig);f
$ U1 S9 B/ W6 GPS=abs(sig_freq).^2;
5 b' h& I8 s* L& t& f2 C9 QPS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);
/ s8 c# C, Z4 B# r* _% `# W2 _
, g( c# `, f7 J1 M0 e, k, m5 \9 Ufigure,
+ k% H+ V$ q( e4 |5 Nsubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
* Q& ]# {5 B4 C$ d& Ztitle(123)
xlabel('Time (s)')
& e4 Y5 w( x0 ~6 p. n! s. {
7 a6 w. {! a0 I: bsubplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
title('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
ylabel('Power Spectrum')

5 `# d- u6 L2 l. ^# z4 }+ n6 m- K
) l3 F: R  E4 {: j請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝