有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
. T& n3 u3 ]* Q7 x7 k' C, C0 o& F" M因為我在用這個軟體在跑SDM
2 Y& v" q* s( P" q7 E* ~3 E! i那我想要請問一下
我要查一個東西
2 B  e; J4 f' P" d是功率頻譜密度
那如果我要用Simulink下去看的話
請問一下Block 要打什麼
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
請各位教我一下~~
- P: c( V: h; t+ z4 H/ }. ?謝謝

mm770910

最近也正在學習matlab解析PSD
6 C! N7 D/ ~9 g# v希望有進展

s13571357

最近也正在學習matlab解析PSD
一直搞不出來，慢慢研究中

BTTB

最近也正在學習matlab解析SDM的SNR
6 e. v) b' ]2 G) T; P& I' w. m一直搞不出來，慢慢研究中

ainge

謝謝您的分享，有助於學習PLL modeling，
謝謝您

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
$ w' _" C7 i( e7 c  ^1 z, w1 v! e
Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

xyzabc0930

再問一個問題一下~~~
) G# p' Z6 g/ ?6 `& y那如果我是用
Hann-windowed fft.
) ~7 K% c, Q4 Y" c" T" v求出PSD之後
那再使用附件所附的檔案
再去求SNR是嗎?

xyzabc0930

請問kevin一下
4 j$ v: ~- y/ I6 w4 k% j: ^4 e那N指的是量化器的位元數?
還是什麼意思?
$ G: c$ j7 F( x% V' H9 |% u2 @; Z0 @再麻煩回答一下
& p6 L9 G: n3 Q8 F謝謝你喔

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
9 @/ s! l7 f5 x% p) ~7 z1 \/ Pspec=fft(v.*ds_hann(N))/(N/4);
  G( g/ z6 M# ^0 q9 R, ~, cplot(linspace(0,1,N/2),
8 {# `+ K/ k- f$ K5 D) O, N+ zdbv(spec(1:N/2)))
" e6 N6 P7 s$ y: S) g8 iaxis([0 1 -120 0]);
1 p2 r( y, s4 ?: lgrid on;
ylabel('dB')
" T" M" e* u; F* Q1 S' }+ C
  U! \: a9 \0 n2 x從左邊這張圖可以去算出屬於它的SNR,參考附件.
( ?5 T# y. ~0 [3 `: X! M: z$ c

xyzabc0930

恩~~謝謝kevin所提供的~~
; W! q' s; f5 J# Y不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
那我把量化所產生出來的值(1,-1,1,1,-1.......)
那請問一下我要怎麼轉換成左邊的圖
是經由FFT的嗎?
算出來的嗎?
, d, I7 k" n  A" w可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.
$ F+ m- \; h' A+ _! B! x, L9 r- w) u
3 v" J7 z0 R. K0 V

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
那是不是要先把這些量化的值先做FFT
那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power spectrum estimation 有幾種方式,
Peroidogram
* P1 z5 W7 I; s1 Q3 e* YModified periodogram
3 I% }6 }9 p1 m) h& E( Y" ]" x; M& OBartlett method
Welch method
0 u6 |3 a. u, C- j% uBlackman-Tukey method

若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
8 k+ Z& U1 [- S; M4 Yx   : input sequence
) I! y% y! J3 V& a" T7 D  z8 [n1  : starting index, x(n1)
* ?# J1 d( L( W8 Gn2  : ending index, x(n2)
4 \/ U9 w8 [  J. J! s* m# b8 s" l
function Px = per(x,n1,n2)
% \; m: d' B/ [# L; e: u' ox = x(;
8 f; f/ z) B1 Z$ s7 J& T" aif nargin == 1
    n1 = 1;  n2 = length(x);  end;
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
3 X) R- C) v' CPx(1) = Px(2);

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。
0 f; o* d/ l" R) {6 a; C3 H
5 k* _/ @, w* x要算SNR的話，最好還是自己寫。

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
$ |" j( i0 l' @, V) f" v& Y  O, Q都是寫Code的方式的?
# }, c8 z3 \, J& N( J9 o有人可以分享一下經驗的嗎?