Chip123 科技應用創新平台
標題:
有關Simulink
[打印本頁]
作者:
xyzabc0930
時間:
2008-5-17 11:02 PM
標題:
有關Simulink
我想要請問一下各位有用過Simulink的先進
! x& Y% P3 H, l$ H) T& j
因為我在用這個軟體在跑SDM
5 ]( 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, ]' c
Peroidogram
8 \% t4 R' C3 C a
Modified periodogram
2 Z, g5 E' @( A$ D C7 Q/ A
Bartlett method
7 E6 L3 ~: g. N) e$ [
Welch method
3 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 b
x : input sequence
+ U: @9 G8 l) u. `7 f
n1 : 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 F
x = x(
;
4 x5 B+ B# ~- a* L* r
if nargin == 1
1 \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, T
spec=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* B
dbv(spec(1:N/2)))
" i& Z- {) c. u& f% d( m
axis([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的SNR
3 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
;4 x5 B+ B# ~- a* L* r