Chip123 科技應用創新平台

標題: 除小數三階Delta Sigma 頻譜問題 [打印本頁]
作者: ja1688888    時間: 2021-9-17 11:22 PM
標題: 除小數三階Delta Sigma 頻譜問題
本帖最後由 ja1688888 於 2021-9-17 11:24 PM 編輯 9 W2 S  P! w8 y0 C& G& s! n! d

' Y& h# i9 y! r# g剛學hspice、matlab 不久各位多多指教目前是用hspice 跑出mash1-1-1 delta sigma modulator -3～4的輸出，
. X9 b9 `" R- O* W. y; f% ], r: Y手key了前64個輸出到txt檔給matlab做頻譜分析，
8 @& C. ?5 g7 h1 x9 s) Z1.JPG
) A, J# A% k2 d8 \9 g; m結果卻不像其他論文上的noise shaping% C# f+ @1 E6 G
2.JPG
* K. H* Y2 B. a3 E0 e3 z0 l以下為matlab code5 b! f3 X6 H+ c& S1 I" I

5 U9 i. c& |. _+ V; mclear, close all+ h* g* d8 J5 h" R
%% initialize parameters
  ?. V( W1 I. A3 csamplerate=250; % in Hz* C. @# `( M! A$ s
N=64; % data length
4 [+ J. s* D1 k" f7 y! T# E! @sig=load('c:\xshell\d1.txt');
) l% z4 v. a# v6 ]/ z$ [/ c3 lt=[1:64]; 4 D# X5 j9 \# k; |
plot(sig(:,1))           
' c8 {1 ^8 q& v+ ?. S, ?" K/ K% isig_freq=fft(sig);f2 h2 E& Q% |8 ]- U# d. P
PS=abs(sig_freq).^2;
( ^1 l" n! Z# Z/ h0 QPS=PS/max(PS);
; f0 E' D# G- E" Lfaxis=samplerate/2*linspace(0,1,N);5 X: x% ~7 o* D6 x
/ N3 j  }+ D5 p, T3 h  b
figure,; c' N! b$ _1 u, k+ v9 \6 H  R
subplot(1,2,1),plot(t,sig),xlim([t(1) t(end)])! ~+ n+ p0 Z& s5 G
title(123)
3 j5 O+ S$ k2 X* ^% q. Yxlabel('Time (s)'); a, g/ m5 f1 T7 C( D7 O
* c7 @$ _  b/ }6 U5 J: r
subplot(1,2,2),plot(faxis,PS(1:N))  . [$ G! y: D/ C: @  E# }1 b- E* R
xlim([faxis(1) faxis(end)]); A2 C) W" w2 @/ q) ]5 A+ F$ n
title('Spectral Analysis (FFT)')6 q0 ?' R  ]0 U; {& Y- G
xlabel('Frequency (Hz)')
' ?* K, i" x( F6 c9 ^ylabel('Power Spectrum')
3 R0 L2 m# I. K% l
9 P) v1 E" z6 |/ }/ e# r
9 Z- t4 W0 s3 ]' |4 O- K, P請問有甚麼較正確的方法或語法能跑出noise shaping嗎?謝謝' T$ x0 b9 n/ b7 Y1 b! j




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