除小數三階Delta Sigma 頻譜問題

ja1688888

. k8 ^  ^/ n, W+ ]1 k# [
剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
; P, m0 K/ S' h1 V9 ~5 L手key了前64個輸出到txt檔給matlab做頻譜分析，
+ _* B- _6 R( v3 F: y1.JPG
結果卻不像其他論文上的noise shaping
- p1 \5 @( l( {9 C0 m2.JPG
4 O% l& a  l8 H7 v以下為matlab code
  h9 A% N3 ~2 Y% R6 ^( ~* V  Q
$ C5 z3 {* p" J8 w; oclear, close all
6 Z$ @% [: i0 K/ @4 m%% initialize parameters
  F8 H$ u0 N6 ~! D5 H$ ?samplerate=250; % in Hz
N=64; % data length
sig=load('c:\xshell\d1.txt');
t=[1:64];
plot(sig(:,1))           
sig_freq=fft(sig);f
2 w8 P9 M8 M" _7 ?+ S3 W9 e& uPS=abs(sig_freq).^2;
9 O( j5 b- W4 x! q9 IPS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);
" W! z& A" B% U
figure,
$ ?7 ?7 o, N3 u$ d" P2 `3 Q( fsubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
$ @6 \& f# k0 p  f0 C* Vtitle(123)
* a$ \$ `6 {+ M9 _. |. z# rxlabel('Time (s)')
3 t, r; H& }. q  Z4 q. @  J
, d7 F0 `" o8 e/ Q( H; U# Jsubplot(1,2,2),plot(faxis,PS(1:N))  
, ^' k, m9 u& A6 t: S1 ixlim([faxis(1) faxis(end)])
- A: n! s) H6 Y& I1 k" Rtitle('Spectral Analysis (FFT)')
xlabel('Frequency (Hz)')
ylabel('Power Spectrum')


  L4 ]* T- G5 b; v' n請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝