Chip123 科技應用創新平台
標題:
有關Simulink
[打印本頁]
作者:
xyzabc0930
時間:
2008-5-17 11:02 PM
標題:
有關Simulink
我想要請問一下各位有用過Simulink的先進
' @+ q: Z2 G ~8 B% a
因為我在用這個軟體在跑SDM
$ r+ b- E" \$ Z( Z& v' H N
那我想要請問一下
8 Q5 ~* p; z! G3 e/ V7 h
我要查一個東西
: ^2 f$ v0 i+ L
是功率頻譜密度
; O. I! {3 Q7 V1 y9 c
那如果我要用Simulink下去看的話
! D; q. o2 T: U2 Y! T
請問一下Block 要打什麼
! ]. o3 T2 X7 H& n, h8 ^
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
# ~. |2 p- o4 m8 j- [
請各位教我一下~~
) F8 \6 |( f3 H3 B0 ~: k) r
謝謝
作者:
xyzabc0930
時間:
2008-5-19 01:46 PM
還是大家要用Simulink來看功率頻譜密度
2 Q6 n9 c4 L5 A7 W7 t
都是寫Code的方式的?
7 G: Y, ~- L% x8 b* U6 w$ t1 s
有人可以分享一下經驗的嗎?
作者:
hcjian
時間:
2008-5-20 03:24 PM
通常是用periodogram,但是它只有2的冪次方點數可以選,所以你跑出來的FFT會有其他的tone出現,
9 s( k0 Z8 s/ b) e) G+ e- g
不是很準的(本來用FFT的話,點數沒有取對,就不會準,而且取不同段,跑出來差很多)。
' e: B r# B8 J
! @7 m+ d; G+ k2 E
要算SNR的話,最好還是自己寫。
作者:
kevin
時間:
2008-5-21 11:20 AM
power spectrum estimation 有幾種方式,
7 M- |& X" u. [$ e
Peroidogram
- N$ g; J# C! `8 G- w* J
Modified periodogram
' M" l0 r% c' G) I/ j' s# J
Bartlett method
& P: X& B5 R: v& {0 R
Welch method
4 K% R% }6 s. l) p7 s2 ]
Blackman-Tukey method
3 ~* Y. J% c. e ?9 J5 w
6 i. Y @( q8 d8 O
若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
9 a+ i- y/ ~0 q# D+ ]( z
x : input sequence
# c9 ], k! ~" ^) X4 ?3 L% @+ \& x
n1 : starting index, x(n1)
- F7 g$ i2 t% H
n2 : ending index, x(n2)
8 a: Y4 L- ^4 W* D
# i5 W& O- {# `8 b. q
function Px = per(x,n1,n2)
( ~ x# T4 _8 i: C
x = x(
;
; I* q: ^) i5 C& \' C' b
if nargin == 1
: V: _% k b d
n1 = 1; n2 = length(x); end;
' y+ A( W/ j9 z& @2 w! L' Z
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
5 s1 `$ x4 m o6 d* e9 k/ H
Px(1) = Px(2);
作者:
xyzabc0930
時間:
2008-5-21 02:48 PM
我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
! v4 B0 U$ I6 p. k6 Z
那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
' l7 r% e5 |4 a
那是不是要先把這些量化的值先做FFT
* T8 E/ u- ^5 r, O$ ?0 k
那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?
作者:
kevin
時間:
2008-5-21 09:32 PM
power Spectrum 跟 SNR 是不同的圖.
0 j0 Z& l$ ~" s2 a, n: s
9 G! Q1 [0 Z/ @+ t2 a
[attach]4135[/attach]
作者:
xyzabc0930
時間:
2008-5-22 12:31 AM
恩~~謝謝kevin所提供的~~
* _) j$ u, g m) S% C, [8 B' q
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
+ \9 Y# i3 c5 V+ e2 `3 s$ Z( `3 R
那我把量化所產生出來的值(1,-1,1,1,-1.......)
, u$ ?- G0 `! l/ o5 u0 l- Z7 O, D; _
那請問一下我要怎麼轉換成左邊的圖
2 d/ o2 M e# C! o/ `) |
是經由FFT的嗎?
8 G4 I# p8 T& T: L+ b' t3 W" Z! X
算出來的嗎?
/ I2 ~, g" ?' ?$ @# |" X
可是不是也可以從左邊這張圖去算出SNR是多少的嗎?
作者:
kevin
時間:
2008-5-22 12:02 PM
左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
4 E) N6 X' D; R4 t* u# \
spec=fft(v.*ds_hann(N))/(N/4);
v4 y; F; e) i) m" N. v
plot(linspace(0,1,N/2),
; c# e) i4 P0 p" M8 V
dbv(spec(1:N/2)))
% g4 g" w6 b+ f% L7 _; q: W+ ^
axis([0 1 -120 0]);
, r! Y+ L. L4 x! U1 w1 Z
grid on;
$ z; T. l% A3 v9 L/ Q. N* s
ylabel('dB')
0 n/ D/ n( H: @+ c+ |
& n# l9 ]& i8 E
從左邊這張圖可以去算出屬於它的SNR,參考附件.
6 Y3 \0 y0 B/ t9 s+ e6 O6 f, h4 F" o
[attach]4138[/attach]
作者:
xyzabc0930
時間:
2008-5-22 01:56 PM
請問kevin一下
5 l( Q4 ?0 m3 b8 N
那N指的是量化器的位元數?
1 [8 k6 ]9 u5 R* m/ O1 Z
還是什麼意思?
7 |; a, T6 Q! M
再麻煩回答一下
+ i2 p- t' v$ n! i& |
謝謝你喔
作者:
xyzabc0930
時間:
2008-5-22 02:18 PM
再問一個問題一下~~~
7 Z, |6 u2 ]* H
那如果我是用
- r' T! Y! w% S: c( [
Hann-windowed fft.
1 K# S5 @6 q: I5 }
求出PSD之後
( ~. I0 ?% L# Q7 S
那再使用附件所附的檔案
8 B H2 N" J7 Y& Q* T) p
再去求SNR是嗎?
作者:
kevin
時間:
2008-5-22 04:47 PM
N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
( H6 s7 y6 e( e
[attach]4144[/attach]
( L. r! T5 ^, @) m" e
Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.
作者:
ainge
時間:
2009-8-19 06:47 PM
謝謝您的分享,有助於學習PLL modeling,
. `' z5 G7 m) U3 |: `) j
謝謝您
作者:
BTTB
時間:
2011-10-25 03:23 PM
最近也正在學習matlab解析SDM的SNR
1 j( b! a4 I! w. R3 p
一直搞不出來,慢慢研究中
作者:
s13571357
時間:
2012-1-16 10:23 AM
最近也正在學習matlab解析PSD
% n7 X$ t9 g4 i3 b3 t9 V
一直搞不出來,慢慢研究中
作者:
mm770910
時間:
2012-4-12 03:19 PM
最近也正在學習matlab解析PSD
' m0 ?; M$ R- h& K9 j. e
希望有進展
歡迎光臨 Chip123 科技應用創新平台 (http://www.chip123.com.tw/)
Powered by Discuz! X3.2