Chip123 科技應用創新平台

標題: 除小數三階Delta Sigma 頻譜問題 [打印本頁]
作者: ja1688888    時間: 2021-9-17 11:22 PM
標題: 除小數三階Delta Sigma 頻譜問題
本帖最後由 ja1688888 於 2021-9-17 11:24 PM 編輯
- Z3 r& F- h; ~* |6 e: e
! u+ W" ]& M* N剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
, C+ J! @/ t: n6 V6 {手key了前64個輸出到txt檔給matlab做頻譜分析，4 O6 v. |% O/ f  e# N3 c
1.JPG" O6 h0 W* [8 D) F& E6 {0 f, L
結果卻不像其他論文上的noise shaping6 m0 d+ K; x3 B5 i! F
2.JPG: i2 P# R$ @& @4 h; ^
以下為matlab code% f# X# M: J, N# B" G; j

* _/ V. {7 p* e6 e+ |( q6 i! X" Pclear, close all
7 K4 o8 z( s1 R5 Y; Z%% initialize parameters- D& y1 F4 r1 D: k" i
samplerate=250; % in Hz1 d( U+ y3 O% {  r
N=64; % data length, A+ B' s# Q: t/ z) G+ ~9 B) s: ]
sig=load('c:\xshell\d1.txt');
* [" ]* y, Q  c+ G: W$ _t=[1:64];
! t2 k& U9 C* W; i! B, u6 Eplot(sig(:,1))           
; Z: Q0 z+ V6 z2 W1 Q" wsig_freq=fft(sig);f2 n: Q0 C. `3 N8 e3 f, S+ _. s
PS=abs(sig_freq).^2;
& E# z' @3 ]9 ?8 |! H. oPS=PS/max(PS);
& g% s  P5 V& I/ R+ s. `faxis=samplerate/2*linspace(0,1,N);) E- F& l0 @8 g5 j% T9 X' y) W

$ ^9 |  u6 |" p2 Q# Y5 Bfigure,
; i" a, M) a- csubplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])6 b# D. W+ D' G% c2 W
title(123)$ _% X# \* s* F9 f. [: C# E4 A
xlabel('Time (s)')0 p8 E7 y0 v& M2 u  x
1 g$ m2 q. d: g% L: [& ^) ?
subplot(1,2,2),plot(faxis,PS(1:N))  
+ I0 x% {0 V5 E& Oxlim([faxis(1) faxis(end)])/ B0 @, S3 B% \, O4 Q$ y+ L
title('Spectral Analysis (FFT)'). _6 [) u# l& E
xlabel('Frequency (Hz)')
" w: d* E9 g, ]; qylabel('Power Spectrum')' R: a. R) Y2 i8 B, a  o! p" X0 h
; o5 s& ^* ^  u/ |" _
! O- b1 P8 C; G- _. y2 {2 k
請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝( z) z! O6 C. b- A( s' E4 A- j




歡迎光臨 Chip123 科技應用創新平台 (http://www.chip123.com/) Powered by Discuz! X3.2