Chip123 科技應用創新平台

標題: 有關Simulink [打印本頁]

作者: xyzabc0930    時間: 2008-5-17 11:02 PM
標題: 有關Simulink
我想要請問一下各位有用過Simulink的先進! x& Y% P3 H, l$ H) T& j
因為我在用這個軟體在跑SDM5 ]( t. x4 G& h5 Y7 n/ }# e
那我想要請問一下
! _( @8 [, _% Y1 D& E我要查一個東西
: k% ^4 l8 o3 X; W5 S* u是功率頻譜密度
7 I2 c: J8 S- r. M$ U那如果我要用Simulink下去看的話/ t1 X6 f1 u! a9 i. d
請問一下Block 要打什麼/ Z" o3 f* A8 p: ?5 X
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
. w$ v$ |& V4 Q; r請各位教我一下~~9 n, `1 A* z# F5 U' W
謝謝
作者: xyzabc0930    時間: 2008-5-19 01:46 PM
還是大家要用Simulink來看功率頻譜密度; v$ G6 s% S7 `3 Z$ n
都是寫Code的方式的?
7 p% M! X, o) g7 A( w; e有人可以分享一下經驗的嗎?
作者: hcjian    時間: 2008-5-20 03:24 PM
通常是用periodogram,但是它只有2的冪次方點數可以選,所以你跑出來的FFT會有其他的tone出現,
# A- V8 }8 y/ T: o不是很準的(本來用FFT的話,點數沒有取對,就不會準,而且取不同段,跑出來差很多)。- w( c4 J/ j( v6 Z9 S) X9 j# }

- J- B- C: b* t( ~) }- [8 c' [要算SNR的話,最好還是自己寫。
作者: kevin    時間: 2008-5-21 11:20 AM
power spectrum estimation 有幾種方式,
4 R& Q! s, ]' cPeroidogram
8 \% t4 R' C3 C  aModified periodogram
2 Z, g5 E' @( A$ D  C7 Q/ ABartlett method
7 E6 L3 ~: g. N) e$ [Welch method3 p3 H# M0 d7 l8 R2 y: H, x
Blackman-Tukey method
, a& R9 g. X9 {% O5 h: p5 B
, Z) Q9 o" Q( S: S- \3 _# F% |4 t若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
  G; e2 w: r/ D1 Q6 bx   : input sequence
+ U: @9 G8 l) u. `7 fn1  : starting index, x(n1), Q3 O: _3 j" }' {% `" }  m
n2  : ending index, x(n2)2 I" F4 W  ?7 G3 f# G3 w$ }& X: _

9 ]3 H" C& ^) h8 n7 _! d- \function Px = per(x,n1,n2)
6 j& k, q) x+ P; x6 Fx = x(;4 x5 B+ B# ~- a* L* r
if nargin == 11 \8 b; Z, E) L- ]1 q' m
    n1 = 1;  n2 = length(x);  end;3 x6 G1 k% u9 L3 g  Y, h
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
1 X8 i2 H8 b- `+ P+ C! [Px(1) = Px(2);
作者: xyzabc0930    時間: 2008-5-21 02:48 PM
我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
0 o, j0 n7 K) i* C" j& R那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
* y7 S$ _, D* d6 W3 w那是不是要先把這些量化的值先做FFT/ _' V: V' T  u' X2 _6 P
那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?
作者: kevin    時間: 2008-5-21 09:32 PM
power Spectrum 跟 SNR 是不同的圖.
2 J7 ]4 p5 B4 K
9 p# F& Q& @$ x1 k( Y) U[attach]4135[/attach]
作者: xyzabc0930    時間: 2008-5-22 12:31 AM
恩~~謝謝kevin所提供的~~
3 |- Y) w( F( r3 k2 W: I不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
4 J4 y; M3 {& z4 I那我把量化所產生出來的值(1,-1,1,1,-1.......)4 ^3 ?* e2 `; B8 i1 p; l# \) {4 t$ u
那請問一下我要怎麼轉換成左邊的圖9 [. o) I* K; t
是經由FFT的嗎?9 d& ?8 `, B' d. E9 s' X. F  i" e& B
算出來的嗎?
$ D( Z+ C0 k9 {* B/ X, e可是不是也可以從左邊這張圖去算出SNR是多少的嗎?
作者: kevin    時間: 2008-5-22 12:02 PM
左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
( Z5 f" K+ E8 I0 S( V, Tspec=fft(v.*ds_hann(N))/(N/4);+ g+ x5 T0 I. Y* z8 T0 x
plot(linspace(0,1,N/2),
6 `/ @1 y0 k* Bdbv(spec(1:N/2)))
" i& Z- {) c. u& f% d( maxis([0 1 -120 0]);$ O. G3 s3 |5 {4 c3 b! H- W
grid on;
, @: X0 u5 T# g8 c* |ylabel('dB')( }& V' s: u- Y* N& ]1 B3 b
! {- a! j6 ~7 z3 r4 {) F
從左邊這張圖可以去算出屬於它的SNR,參考附件.( ?" L1 \4 O% M6 e! |! l5 _0 A
[attach]4138[/attach]
作者: xyzabc0930    時間: 2008-5-22 01:56 PM
請問kevin一下
/ w7 k5 w. T' \0 s7 A; T那N指的是量化器的位元數?) H  h, F  |" T1 \1 R
還是什麼意思?2 U2 c( l7 y+ e" _% D4 y7 i3 h
再麻煩回答一下$ x. i( Z  O, h+ Q* l2 x
謝謝你喔
作者: xyzabc0930    時間: 2008-5-22 02:18 PM
再問一個問題一下~~~. e# `) o. m7 N' }8 F2 p' C3 N& o
那如果我是用, r9 Y* |6 b. ~& @/ U: f- N2 h
Hann-windowed fft.
$ _4 e, h- x; v2 q6 W1 J求出PSD之後$ D% m& X3 b$ ?0 ]& h- ]5 V
那再使用附件所附的檔案
: W: H# {) p! ]& K" r6 }3 R. _$ t再去求SNR是嗎?
作者: kevin    時間: 2008-5-22 04:47 PM
N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
6 K! G  L% a. b- i% \7 x[attach]4144[/attach], W* u) L% f- t; U5 h6 F
Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.
作者: ainge    時間: 2009-8-19 06:47 PM
謝謝您的分享,有助於學習PLL modeling,; X, ~' Q* V% r$ H
謝謝您
作者: BTTB    時間: 2011-10-25 03:23 PM
最近也正在學習matlab解析SDM的SNR3 x5 G# V2 x. K4 s4 S
一直搞不出來,慢慢研究中
作者: s13571357    時間: 2012-1-16 10:23 AM
最近也正在學習matlab解析PSD/ p: p$ X( J+ G( J9 r  A
一直搞不出來,慢慢研究中
作者: mm770910    時間: 2012-4-12 03:19 PM
最近也正在學習matlab解析PSD
" M% H5 m% x0 i6 u% B. f希望有進展




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