除小數三階Delta Sigma 頻譜問題

ja1688888

% L- h' t6 _* \% T
剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
5 u$ ~0 Z' X6 N0 q手key了前64個輸出到txt檔給matlab做頻譜分析，
5 A$ J0 F6 x5 e- Y/ x8 h& v; G1.JPG
& J9 v' @% y2 g5 C$ I9 r" @結果卻不像其他論文上的noise shaping
2.JPG
以下為matlab code
: g( A8 q5 N$ ?$ `7 z; A4 G
clear, close all
%% initialize parameters
samplerate=250; % in Hz
N=64; % data length
sig=load('c:\xshell\d1.txt');
4 W; D2 ?3 d# Ft=[1:64];
plot(sig(:,1))           
sig_freq=fft(sig);f
; [, P6 w" N& l$ W4 g3 {PS=abs(sig_freq).^2;
. i7 ]% ~1 x4 b! v6 \PS=PS/max(PS);
% ^" N; K3 W* B* X; Xfaxis=samplerate/2*linspace(0,1,N);

figure,
( V# M* o2 O$ G0 x, C6 o, Jsubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
# Z. I2 d9 t% c- d0 q/ B+ [xlabel('Time (s)')

9 A$ I3 u) A4 Osubplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
' i* E. x6 o. Z' ^( E1 ptitle('Spectral Analysis (FFT)')
) r" [& S! i& _5 [3 b0 {# pxlabel('Frequency (Hz)')
' g% Q1 q7 z3 f- h- y" M% w) s5 s7 {3 Kylabel('Power Spectrum')

" s! f! ?  k7 o
3 D. ^7 [/ x2 V# t3 M: ?請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
2 u+ L, j! k% n) g) Z