Chip123 科技應用創新平台

標題: 除小數三階Delta Sigma 頻譜問題 [打印本頁]
作者: ja1688888    時間: 2021-9-17 11:22 PM
標題: 除小數三階Delta Sigma 頻譜問題
本帖最後由 ja1688888 於 2021-9-17 11:24 PM 編輯
8 b+ ?1 D  G# n& B1 i
3 `, u. W" l) w+ o. H( h% {剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
( ?$ `: ~8 k8 D5 m手key了前64個輸出到txt檔給matlab做頻譜分析，  |% x$ p) m8 u
1.JPG
; Z1 y4 ]& o6 I! e; W* p# v結果卻不像其他論文上的noise shaping
$ K: j7 `0 M% l( E% m7 ?2.JPG
6 @0 r* |1 C) ^5 X; c8 H7 ~- B以下為matlab code9 W: e: F2 @$ M5 [) ]0 ]/ y3 K  l

5 [9 h/ w. ~3 ^5 G. @9 L1 @4 E( dclear, close all# u6 I2 R5 Z$ ?# ]5 V
%% initialize parameters1 E+ t, H5 L" h' R
samplerate=250; % in Hz" Y" k9 W$ T+ S) Q7 ^, b
N=64; % data length8 z- C5 ^; O7 M, p
sig=load('c:\xshell\d1.txt');
+ Q5 L, j" L9 M, j5 t" {* ?. ft=[1:64]; ; o4 R1 R2 h7 r2 G( ^1 r$ g
plot(sig(:,1))           
6 m; M6 A6 `! [! v3 ?sig_freq=fft(sig);f: l5 g0 A: H1 x
PS=abs(sig_freq).^2;6 d; |! y4 x' K
PS=PS/max(PS);+ a' I1 |8 W; ?1 v% R
faxis=samplerate/2*linspace(0,1,N);
& P9 e" e+ \. v! }( R& e, v
" }# D9 n4 ]: b# f+ Afigure,3 f8 ], @! e: ?+ D: \1 g3 x
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])! c5 Q  W. w* n* Q  i0 p9 V" W' Q
title(123)
- ?! b# N' L# {8 T4 ?xlabel('Time (s)')7 }. w' s& y/ U  C) c! L
! j2 j% I  \1 J/ l8 f3 T& c* V
subplot(1,2,2),plot(faxis,PS(1:N))  1 M0 Q9 z+ J3 C3 x
xlim([faxis(1) faxis(end)])
. u0 |! ^/ D5 n7 V$ ~  k9 A2 M7 S% Gtitle('Spectral Analysis (FFT)')* m, E. i* a. V3 {( Z' B
xlabel('Frequency (Hz)')
, M# x; p$ `/ ]6 c0 H5 \ylabel('Power Spectrum')) w' K$ h& W, ^( n

2 M* y# a* j& S0 n) U0 b" a
3 R4 S3 h) m5 t% t3 l請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝+ B! ^& i, `1 g/ e. h  L. a




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