除小數三階Delta Sigma 頻譜問題

ja1688888

: ?/ Z/ V2 W+ A! f7 a# I
  |. V1 R. B5 P9 g) c! U' n剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
手key了前64個輸出到txt檔給matlab做頻譜分析，
1.JPG
結果卻不像其他論文上的noise shaping
4 R- |  Z" r" I1 B0 Q# D( l2.JPG
# E  x3 n+ L# x8 x+ L以下為matlab code

3 |* Y9 D/ W. y+ Q% R- T8 W! @clear, close all
) s0 K9 y# j7 N# I%% initialize parameters
samplerate=250; % in Hz
4 c. t3 i- J# Q( n9 jN=64; % data length
sig=load('c:\xshell\d1.txt');
2 s5 u% q4 |& S6 at=[1:64];
plot(sig(:,1))           
' ~9 c! t, \( H" X, d, A+ hsig_freq=fft(sig);f
PS=abs(sig_freq).^2;
PS=PS/max(PS);
faxis=samplerate/2*linspace(0,1,N);

- ?& M/ [: X1 N8 R3 hfigure,
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
title(123)
xlabel('Time (s)')

subplot(1,2,2),plot(faxis,PS(1:N))  
xlim([faxis(1) faxis(end)])
* w; L! `+ `7 ^7 etitle('Spectral Analysis (FFT)')
" a7 Y% h" D6 M; i" |xlabel('Frequency (Hz)')
ylabel('Power Spectrum')
( `+ X7 E7 b! J9 o4 N& \% e

請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝