Chip123 科技應用創新平台
標題:
除小數三階Delta Sigma 頻譜問題
[打印本頁]
作者:
ja1688888
時間:
2021-9-17 11:22 PM
標題:
除小數三階Delta Sigma 頻譜問題
本帖最後由 ja1688888 於 2021-9-17 11:24 PM 編輯
/ k4 ~/ t o1 C4 p; A0 D1 d- y
, ^/ _9 w, p8 C
剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
+ a$ V P8 u- J- b( v% R* v
手key了前64個輸出到txt檔給matlab做頻譜分析，
5 w5 F% U# b" A6 {- R8 J
1.JPG
2 Y9 h- ^0 k4 S/ b' h( X% Q
結果卻不像其他論文上的noise shaping
/ h2 {1 A! P9 P6 u( ]% {" w
2.JPG
, q; ~$ w2 S9 O+ S
以下為matlab code
: H2 D/ O& \" D
; J/ A7 R* w# v( ?: b/ _( `
clear, close all
* f" j, B4 N7 Y& L' j( Z5 l& ^
%% initialize parameters
* ]" Q/ c- D! G
samplerate=250; % in Hz
' E ]. V8 M4 K& f
N=64; % data length
+ Z* e, q: Y$ u: g* |
sig=load('c:\xshell\d1.txt');
1 N4 [3 |& l8 _4 b
t=[1:64];
' Y, M6 K' B* _: N4 @' ]0 `
plot(sig(:,1))
M4 a$ N, D& ]+ a
sig_freq=fft(sig);f
& P# g3 F/ f5 @! C+ O
PS=abs(sig_freq).^2;
! [/ { G- a0 O# P1 {. B
PS=PS/max(PS);
' _! ~8 \) e7 X- L9 F" y
faxis=samplerate/2*linspace(0,1,N);
$ F9 {6 e) Y- ~$ y9 q' w
8 A& V! z# y$ T k3 ^
figure,
U" Q. @& }% r9 T. t
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])
5 c- {: C# L/ u3 A( I
title(123)
1 ?$ |! T& K1 _& s! R1 [4 W
xlabel('Time (s)')
; l3 m* \$ M L3 F# ^
" O, f; ^: R$ Y+ X+ ^
subplot(1,2,2),plot(faxis,PS(1:N))
3 r' n( \& T" C1 o Y/ F
xlim([faxis(1) faxis(end)])
4 l5 Y4 R) l$ v* t
title('Spectral Analysis (FFT)')
: D8 e% _$ T4 ^9 {% a- W
xlabel('Frequency (Hz)')
3 z- k6 d, w, O8 V- ]
ylabel('Power Spectrum')
8 |& v5 {% C L0 s
, k8 Y. x0 ]3 \4 g3 T6 B) x
% M; u3 T5 I+ H6 q4 k- p8 B2 r' }6 v
請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝
/ s" a1 v- L( Y$ s$ |# @2 h# z4 R3 J
歡迎光臨 Chip123 科技應用創新平台 (http://free.vireal.world/chip123_website/innoingbbs/)
Powered by Discuz! X3.2