Chip123 科技應用創新平台
標題:
除小數三階Delta Sigma 頻譜問題
[打印本頁]
作者:
ja1688888
時間:
2021-9-17 11:22 PM
標題:
除小數三階Delta Sigma 頻譜問題
本帖最後由 ja1688888 於 2021-9-17 11:24 PM 編輯
3 c9 l2 y3 N, w3 a) T/ v
" [& ~, X9 a( p
剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
; `8 x1 u) g& w" O1 v5 X9 t
手key了前64個輸出到txt檔給matlab做頻譜分析，
; U' O: P' t2 x8 q1 A
1.JPG
- U! y" o2 [% i0 z$ l0 ]
結果卻不像其他論文上的noise shaping
$ K; v4 O- H" V y" v/ N
2.JPG
) s$ c+ Q1 B; m3 U2 y0 R
以下為matlab code
& d) ]; ^( G4 {9 R: @: |/ e
; |1 o2 p, k! M/ b8 O
clear, close all
7 r1 R% D( c- H3 S8 B
%% initialize parameters
0 E& [# O1 z* \9 |6 W4 ~
samplerate=250; % in Hz
4 S5 j+ c( p5 y5 \
N=64; % data length
+ ^: O- V2 v# C9 r0 O7 Y$ ~
sig=load('c:\xshell\d1.txt');
9 |" h' r/ W2 {" {6 O2 _
t=[1:64];
3 z5 i, f! X- c7 G. q$ t( C
plot(sig(:,1))
) T0 D( K9 \: g/ A% z0 `# S
sig_freq=fft(sig);f
, u) b3 e" b% N! o! o
PS=abs(sig_freq).^2;
" a* C1 X/ d' |& M; o' B
PS=PS/max(PS);
+ h. \( O2 W8 @0 B
faxis=samplerate/2*linspace(0,1,N);
; R* p$ w) D+ e ?. i
% s, H& C | K1 l+ n% `0 o0 R! G- O
figure,
6 z6 n, N" H& L7 D% q) P
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
' B4 m J- ?- x& [- {; p0 B
title(123)
: u2 `6 T- q3 I* q5 o
xlabel('Time (s)')
+ E5 \+ }! ~/ [& k) B
" k" O1 L5 e$ e$ w7 M5 W- F
subplot(1,2,2),plot(faxis,PS(1:N))
1 q/ Q* Q/ u( c# ]9 M
xlim([faxis(1) faxis(end)])
% J! ^) k3 p* R
title('Spectral Analysis (FFT)')
! j, _% H& v' o* D$ p
xlabel('Frequency (Hz)')
* U5 r* v' Q% ~( s- c
ylabel('Power Spectrum')
0 T; g- M, {# U0 q: f
. t; W" F: q% A+ @/ G; F& r' a
' ]$ n% l# A( C7 T* c# I: q! t& l
請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
/ m/ Z, h! W- I0 ?0 s
歡迎光臨 Chip123 科技應用創新平台 (http://free.vireal.world/chip123_website/innoingbbs/)
Powered by Discuz! X3.2