Chip123 科技應用創新平台
標題:
有關Simulink
[打印本頁]
作者:
xyzabc0930
時間:
2008-5-17 11:02 PM
標題:
有關Simulink
我想要請問一下各位有用過Simulink的先進
1 Y( [" i. A- u y
因為我在用這個軟體在跑SDM
; @% S! V$ i0 m R
那我想要請問一下
9 |# \; g2 ~% Y2 G5 T/ E/ o
我要查一個東西
) D# k+ L, h' |9 {
是功率頻譜密度
. E3 X6 f" N# z5 o! z9 r
那如果我要用Simulink下去看的話
+ `$ N8 ~& M" b; J7 X9 N
請問一下Block 要打什麼
1 `5 e" U2 X, u9 h+ a- z) i
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
] y% s! S5 L4 t2 g+ t
請各位教我一下~~
/ s4 ~( c w% C
謝謝
作者:
xyzabc0930
時間:
2008-5-19 01:46 PM
還是大家要用Simulink來看功率頻譜密度
$ } ? p7 l- N
都是寫Code的方式的?
' ` C( r% J& x! T
有人可以分享一下經驗的嗎?
作者:
hcjian
時間:
2008-5-20 03:24 PM
通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
; p0 j/ i/ `0 P5 s3 ]9 B p
不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。
$ J% y: h* o: {) p
2 ]! x4 g, v0 b0 r, t9 I& g" m
要算SNR的話，最好還是自己寫。
作者:
kevin
時間:
2008-5-21 11:20 AM
power spectrum estimation 有幾種方式,
# F7 {! f5 e$ |5 P1 ^" _3 Z" k
Peroidogram
# J8 m7 R$ Z+ t: m' E8 A
Modified periodogram
* [ f D: s& ?! v- O
Bartlett method
% O$ R3 x$ p1 P, ]1 i
Welch method
/ O2 k, j2 m9 ?5 x4 P) _5 e% f4 M
Blackman-Tukey method
% v# m6 w, K4 l6 o, v
# S: z! K9 B1 J& h/ T
若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
- G% L6 R' J' z5 c. f* F
x : input sequence
7 w3 m* H0 g) z+ K0 V; D6 x
n1 : starting index, x(n1)
# x8 S$ W, U+ \5 F5 _- P7 d: Y. P
n2 : ending index, x(n2)
- a* m# Z- R! Z, Q5 L: T% ]1 N: ?
* U6 c5 b5 f) ~; h. m
function Px = per(x,n1,n2)
0 j2 T% k# [$ F' }2 G# J5 d
x = x(
;
& J9 T' G6 ^9 q+ j! {) J( Q
if nargin == 1
) }8 D' V( C" a4 B, R
n1 = 1; n2 = length(x); end;
3 a' t7 n$ ~4 u$ C0 t/ h$ U
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
; y1 ^0 l* H: q% W" V/ ]
Px(1) = Px(2);
作者:
xyzabc0930
時間:
2008-5-21 02:48 PM
我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
9 S* S! K# ~, b4 |9 }6 A
那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
7 j8 H N6 H, O6 j5 L2 W) b
那是不是要先把這些量化的值先做FFT
4 o# y! o! }6 q; {6 r( h
那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?
作者:
kevin
時間:
2008-5-21 09:32 PM
power Spectrum 跟 SNR 是不同的圖.
, C4 [5 _- P7 F1 _* n
: D* {) O5 L: U" g
[attach]4135[/attach]
作者:
xyzabc0930
時間:
2008-5-22 12:31 AM
恩~~謝謝kevin所提供的~~
+ t7 ?1 a- v6 z
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
/ v2 m6 G( {) @* w
那我把量化所產生出來的值(1,-1,1,1,-1.......)
q) Z1 I) z8 ]
那請問一下我要怎麼轉換成左邊的圖
4 I# p [1 I6 S, j+ B5 O9 z# ]& i
是經由FFT的嗎?
! y. t( Z. {( x% K% r: {
算出來的嗎?
2 n' ]& z# Q" S; w3 k
可是不是也可以從左邊這張圖去算出SNR是多少的嗎?
作者:
kevin
時間:
2008-5-22 12:02 PM
左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
; E. \0 v+ F' }& G! ?7 y( o
spec=fft(v.*ds_hann(N))/(N/4);
. G) T4 t$ n! ~7 s+ N+ x/ j5 d( p
plot(linspace(0,1,N/2),
" \* h# Z" q, _: D1 z
dbv(spec(1:N/2)))
# f8 z! z7 b. F& t& I, F
axis([0 1 -120 0]);
$ O' a7 w3 X. i
grid on;
% u* M+ D L# h+ K4 {: C3 W8 l( ?! J
ylabel('dB')
1 ^# n* R' u+ ~* a) {
/ `8 m4 N$ ?3 a
從左邊這張圖可以去算出屬於它的SNR,參考附件.
' d: ? |% N0 `9 F0 {' E5 I
[attach]4138[/attach]
作者:
xyzabc0930
時間:
2008-5-22 01:56 PM
請問kevin一下
+ n( t# i% q8 X! A
那N指的是量化器的位元數?
/ l7 P7 C! I. J+ j, b; k' T4 u
還是什麼意思?
+ R" T( v+ ^$ A W$ K4 {8 ~: N
再麻煩回答一下
* A8 P' Q0 }8 O4 [; h, k5 _( f
謝謝你喔
作者:
xyzabc0930
時間:
2008-5-22 02:18 PM
再問一個問題一下~~~
$ r" P- e- Z" Z, t' V
那如果我是用
' }) E% r( ~. j k7 T+ }5 D, E
Hann-windowed fft.
7 `4 w# z4 h2 w2 {- W" s5 r0 }, t
求出PSD之後
3 Q" U3 R% z& ^" ~
那再使用附件所附的檔案
- ]9 U/ y4 p0 P$ x8 t+ L
再去求SNR是嗎?
作者:
kevin
時間:
2008-5-22 04:47 PM
N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
5 k+ _% W3 I) o) U4 F
[attach]4144[/attach]
9 ] c% S1 G0 B2 ?+ m7 s
Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.
作者:
ainge
時間:
2009-8-19 06:47 PM
謝謝您的分享，有助於學習PLL modeling，
* I6 y( {% L$ E( \
謝謝您
作者:
BTTB
時間:
2011-10-25 03:23 PM
最近也正在學習matlab解析SDM的SNR
+ L* E' H# z, [& ?& e4 J
一直搞不出來，慢慢研究中
作者:
s13571357
時間:
2012-1-16 10:23 AM
最近也正在學習matlab解析PSD
8 B4 {+ |4 i+ H9 b
一直搞不出來，慢慢研究中
作者:
mm770910
時間:
2012-4-12 03:19 PM
最近也正在學習matlab解析PSD
6 U2 n/ Z) Y1 J1 S6 I9 Z$ M
希望有進展
歡迎光臨 Chip123 科技應用創新平台 (http://www.chip123.com.tw/)
Powered by Discuz! X3.2
;