Chip123 科技應用創新平台
標題:
除小數三階Delta Sigma 頻譜問題
[打印本頁]
作者:
ja1688888
時間:
2021-9-17 11:22 PM
標題:
除小數三階Delta Sigma 頻譜問題
本帖最後由 ja1688888 於 2021-9-17 11:24 PM 編輯
) O% z+ `6 y) S
$ X* l% q# s/ L, @
剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
$ }6 ]! N7 _# Y+ {1 l$ p0 s6 ?* n
手key了前64個輸出到txt檔給matlab做頻譜分析，
5 }: t7 C2 c$ W7 G4 q z; u) _
1.JPG
8 @9 ]. y( {) [, P+ k% F
結果卻不像其他論文上的noise shaping
9 M$ O3 y/ m7 |2 h- D8 C# R7 ?: J
2.JPG
3 O6 V$ u( j0 E' x0 U
以下為matlab code
' ?4 W% x' s: B. Z+ K! ]& U
. r3 ~: r/ `: X9 |. g# @1 ?$ Z
clear, close all
' e! Z( w9 {. O- n
%% initialize parameters
" ~* B0 P( A. U+ Y, t* e
samplerate=250; % in Hz
3 S' q0 Z' }1 m4 x, H: G
N=64; % data length
* y7 o6 D0 }" b& ^! e5 l# F
sig=load('c:\xshell\d1.txt');
6 o1 D$ g/ \2 q l/ c: @
t=[1:64];
) X$ G0 o V/ h9 w
plot(sig(:,1))
/ M f0 K0 |: n, K
sig_freq=fft(sig);f
/ U' v, r6 a& ?0 N. v/ {
PS=abs(sig_freq).^2;
. |3 G* v+ } T' @0 H, w
PS=PS/max(PS);
8 j) U) N \) A1 `
faxis=samplerate/2*linspace(0,1,N);
! G. Y( _% {" k" ]" o. i+ q% N
0 {, H& F8 O2 m
figure,
8 f2 v2 _4 i) K% b; I6 r
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
# ^: _4 k: y6 `8 P9 j
title(123)
7 M2 h/ i4 ?( S6 i. B/ N3 X
xlabel('Time (s)')
: B! [/ h! I8 Z- D7 m, j
: T1 F& U1 A7 G
subplot(1,2,2),plot(faxis,PS(1:N))
& v" U/ u( h" @* J
xlim([faxis(1) faxis(end)])
7 s0 C4 Q3 C- [% e3 S6 D1 p
title('Spectral Analysis (FFT)')
8 p+ b/ Z+ I: [; p; P: r
xlabel('Frequency (Hz)')
$ @6 O7 x/ D* G. t
ylabel('Power Spectrum')
1 y# t) L7 W5 Q$ x u5 i5 F9 B
# Y! [+ Y& R* W% G9 Z k
" v0 i/ e: u" | ~; Q0 u! z' C* S
請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
0 K8 K1 h% k; X2 z: S C4 q
歡迎光臨 Chip123 科技應用創新平台 (http://chip123.com/)
Powered by Discuz! X3.2