有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
因為我在用這個軟體在跑SDM
那我想要請問一下
0 T' g$ V& c+ V# ~我要查一個東西
是功率頻譜密度
那如果我要用Simulink下去看的話
8 d/ C2 r4 r0 x& E) ^請問一下Block 要打什麼
5 a4 Q+ C& r, V% E$ s5 p* }才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
請各位教我一下~~
謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
! J/ d: d8 }' I% U9 ^& l都是寫Code的方式的?
有人可以分享一下經驗的嗎?

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。

要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
Peroidogram
$ G0 s) C, a. i6 fModified periodogram
Bartlett method
Welch method
Blackman-Tukey method

若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
* A: t- I8 q& Ox   : input sequence
& W$ d% W! ~3 w/ \& y! ?n1  : starting index, x(n1)
$ [" I5 ]  ?) u% M% D! N/ h1 K+ P0 en2  : ending index, x(n2)

& V, T  }) G8 I3 ]" ~7 Z; xfunction Px = per(x,n1,n2)
x = x(;
. x( {6 V$ h! Y4 P0 c# z& _4 `( Uif nargin == 1
. n2 d- ~( h- y' m; W) M/ _& M" l% E( S    n1 = 1;  n2 = length(x);  end;
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
7 r4 N' {8 n0 zPx(1) = Px(2);

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
5 Z8 K& ]  h8 p' E那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
3 ~2 \+ R; Z  D+ @' q  q) s那是不是要先把這些量化的值先做FFT
3 ^) x6 Q  K. Z  q! H) |那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.
: M, V0 P4 \4 @5 m% b( p
+ O8 d# q, L! U' _6 J2 u: `6 I

xyzabc0930

恩~~謝謝kevin所提供的~~
- C9 g% N7 q/ y2 T/ z不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
那我把量化所產生出來的值(1,-1,1,1,-1.......)
+ ^* T1 a5 v3 v+ L7 U/ u* H那請問一下我要怎麼轉換成左邊的圖
# w  C# m0 n: H& v是經由FFT的嗎?
% m! @1 @# f/ g1 s) r9 a算出來的嗎?
  F7 A$ V, R6 S+ T9 `( k" F: }可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
spec=fft(v.*ds_hann(N))/(N/4);
; f7 c, N: u; Uplot(linspace(0,1,N/2),
dbv(spec(1:N/2)))
axis([0 1 -120 0]);
! l% B* r1 x) ogrid on;
$ X: _, n( n, W) U# [- c/ G* w2 Qylabel('dB')

5 g& Z! r+ _) q8 f! W, E- R& i; J從左邊這張圖可以去算出屬於它的SNR,參考附件.

xyzabc0930

請問kevin一下
3 |! h5 ~+ T6 f那N指的是量化器的位元數?
1 j" ]& w5 d+ D, s+ t& i! y, Q5 l還是什麼意思?
再麻煩回答一下
' ]9 W5 D" h. ?9 k謝謝你喔

xyzabc0930

再問一個問題一下~~~
1 ]) K& \1 |! h1 i0 K8 G那如果我是用
# ~  g& b, b! O$ Y0 E: rHann-windowed fft.
求出PSD之後
那再使用附件所附的檔案
$ Z) h! L* l" r' o0 y再去求SNR是嗎?

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
$ F; I& X% v( c8 `6 F( g
: c. Y/ n& U4 `Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

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

BTTB

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

s13571357

最近也正在學習matlab解析PSD
- M6 {+ x3 g4 `9 t; i一直搞不出來，慢慢研究中

mm770910

最近也正在學習matlab解析PSD
/ _* @: F1 X1 v+ w/ X; A3 N希望有進展