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# JBartlett 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+ ]( zx   : input sequence
# c9 ], k! ~" ^) X4 ?3 L% @+ \& xn1  : starting index, x(n1)
- F7 g$ i2 t% Hn2  : 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' bif 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/ HPx(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: s9 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. vplot(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 Zgrid 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