Chip123 科技應用創新平台
標題:
有關Simulink
[打印本頁]
作者:
xyzabc0930
時間:
2008-5-17 11:02 PM
標題:
有關Simulink
我想要請問一下各位有用過Simulink的先進
, k: m* i# b, {, t/ T
因為我在用這個軟體在跑SDM
4 T2 s" p. k( r3 } X3 O, f
那我想要請問一下
6 P3 H" N# F1 k. c
我要查一個東西
6 p) A& ]. ^7 U6 n5 t
是功率頻譜密度
! Y( y% w* A F3 t" D. i
那如果我要用Simulink下去看的話
! Y- J" p$ o2 |# H
請問一下Block 要打什麼
! H4 o v4 `, d% M1 e
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
6 \0 G# }3 B* Z, P8 w4 S9 ~/ C; Z
請各位教我一下~~
% Q2 N/ K+ O5 x
謝謝
作者:
xyzabc0930
時間:
2008-5-19 01:46 PM
還是大家要用Simulink來看功率頻譜密度
5 e) R. R. M, T q
都是寫Code的方式的?
[8 c @- z! ` k$ K; X4 Y
有人可以分享一下經驗的嗎?
作者:
hcjian
時間:
2008-5-20 03:24 PM
通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
! s" b* l8 S2 C
不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。
n/ |0 W9 Z; B
X7 S- K; X8 z
要算SNR的話，最好還是自己寫。
作者:
kevin
時間:
2008-5-21 11:20 AM
power spectrum estimation 有幾種方式,
* f [* C9 l( P1 H" t$ c
Peroidogram
( j; M3 `$ ~0 {% u3 [; |5 n, F
Modified periodogram
( ~% h- I- c6 ?$ I. `1 y
Bartlett method
0 {/ ~: L3 M9 v M4 N" f
Welch method
; S; [9 n# C0 \
Blackman-Tukey method
: X& Z4 Y( P& M& w( C3 h: N
1 `; m6 }. B6 s- H
若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
1 K0 @2 L1 {, \2 V* `$ K
x : input sequence
( T7 S# U: d: a% F- L+ g$ [
n1 : starting index, x(n1)
9 {; q/ v( \0 T: K) Q3 \: i
n2 : ending index, x(n2)
6 ^# }1 N# E. U6 c, i& \% b( r* V2 b
4 L) Z& ^6 k! T1 F
function Px = per(x,n1,n2)
! c- {4 B! O% Y/ X* A$ |0 x* u+ t
x = x(
;
3 B# [9 X+ g, x
if nargin == 1
* v$ w0 E& P p; v: y
n1 = 1; n2 = length(x); end;
; K* E8 s- [! A. t8 ?4 R
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
3 h1 A+ t5 j" @
Px(1) = Px(2);
作者:
xyzabc0930
時間:
2008-5-21 02:48 PM
我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
+ C. G$ d7 P4 r* p
那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
) N m& H8 d3 c
那是不是要先把這些量化的值先做FFT
' h( v' W- o" D) h
那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?
作者:
kevin
時間:
2008-5-21 09:32 PM
power Spectrum 跟 SNR 是不同的圖.
) M( q A: ]) P& |; | b1 j& |* o
) K4 W8 |& \/ B0 r
[attach]4135[/attach]
作者:
xyzabc0930
時間:
2008-5-22 12:31 AM
恩~~謝謝kevin所提供的~~
$ L( Q+ l+ x9 ^$ f. f
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
& k- g- t! G7 V" C/ A2 ^( ?9 l; {
那我把量化所產生出來的值(1,-1,1,1,-1.......)
1 k+ \( k8 X8 B) t( i( B! E+ B
那請問一下我要怎麼轉換成左邊的圖
2 w! Z1 I0 n; I0 j4 z5 X
是經由FFT的嗎?
& d+ M% v; D; a! |: R2 X& E4 }: O
算出來的嗎?
5 c" k: O" ^- ~2 M; R
可是不是也可以從左邊這張圖去算出SNR是多少的嗎?
作者:
kevin
時間:
2008-5-22 12:02 PM
左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
# @6 p/ {, P2 ]6 S6 N( \; M7 u. n
spec=fft(v.*ds_hann(N))/(N/4);
7 P# C: e/ s" ?
plot(linspace(0,1,N/2),
3 A0 D! `: t4 M/ D& E
dbv(spec(1:N/2)))
: b+ O& C$ h4 J9 c1 E
axis([0 1 -120 0]);
8 P* n/ o* N/ }! s
grid on;
6 N( x: S5 b& A1 W% `5 ^, e
ylabel('dB')
8 e* e P5 @. c2 Q( W: l
$ O, g7 Z1 U, A1 B& a
從左邊這張圖可以去算出屬於它的SNR,參考附件.
P( a- ^! A2 x; h# {2 ?: y
[attach]4138[/attach]
作者:
xyzabc0930
時間:
2008-5-22 01:56 PM
請問kevin一下
1 r* W# t- ^, U+ p' ^! g: H
那N指的是量化器的位元數?
% t0 F) J& ]& \0 a I7 b e
還是什麼意思?
$ B1 P6 L3 Z+ V. v
再麻煩回答一下
|1 V- s) t* F |- e$ k
謝謝你喔
作者:
xyzabc0930
時間:
2008-5-22 02:18 PM
再問一個問題一下~~~
3 D7 q$ z& e; l
那如果我是用
. U0 s' H; n7 w% t
Hann-windowed fft.
: c3 l* s, h$ `4 j6 D- k
求出PSD之後
4 E! B, g3 L% X1 [9 J3 h- c$ T; h! e
那再使用附件所附的檔案
& U/ h6 y) o0 C3 A+ d1 b
再去求SNR是嗎?
作者:
kevin
時間:
2008-5-22 04:47 PM
N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
: Z: f/ k6 c% ]" K% _% y
[attach]4144[/attach]
}, s0 v" t! s3 `6 l( o0 r
Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.
作者:
ainge
時間:
2009-8-19 06:47 PM
謝謝您的分享，有助於學習PLL modeling，
+ a/ ^/ [4 M- \1 j, t! z
謝謝您
作者:
BTTB
時間:
2011-10-25 03:23 PM
最近也正在學習matlab解析SDM的SNR
+ }* I N) j5 }- j* @; w
一直搞不出來，慢慢研究中
作者:
s13571357
時間:
2012-1-16 10:23 AM
最近也正在學習matlab解析PSD
) p n0 q9 d! F! ~$ q' U
一直搞不出來，慢慢研究中
作者:
mm770910
時間:
2012-4-12 03:19 PM
最近也正在學習matlab解析PSD
/ I5 j0 ~; V+ X9 V+ {7 ^
希望有進展
歡迎光臨 Chip123 科技應用創新平台 (http://www.chip123.com.tw/)
Powered by Discuz! X3.2
;