有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
3 ]$ C9 c5 H/ @( ], u: f因為我在用這個軟體在跑SDM
$ m5 B: ~+ Y* k# _7 h! {那我想要請問一下
& B$ m" `7 d* x6 ?) a我要查一個東西
2 {) Q9 p& i5 d* k' U) R- U, c4 Y1 x是功率頻譜密度
) r4 E+ M& \& |" q: i+ F9 u& _, \那如果我要用Simulink下去看的話
' O, y5 N! I( \' W請問一下Block 要打什麼
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
9 f/ g  h$ V- Y. Y  O請各位教我一下~~
! ?& z  ?7 t3 e! e: c謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
8 A  g" a3 w( V; |& F. K% P都是寫Code的方式的?
. l, i, |: v$ V' R, p% M  u有人可以分享一下經驗的嗎?

hcjian

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

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

kevin

power spectrum estimation 有幾種方式,
Peroidogram
8 V9 N+ L- K$ k- zModified periodogram
Bartlett method
Welch method
4 R, z" ^+ f1 s/ Q  kBlackman-Tukey method

若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
, J3 V/ A& ]# Z$ nx   : input sequence
n1  : starting index, x(n1)
n2  : ending index, x(n2)
% b; `% m' `5 |, h9 x# ]
function Px = per(x,n1,n2)
x = x(;
  C! v% k6 j9 K) v% g3 s  Vif nargin == 1
2 G/ R# l" [  Z+ J    n1 = 1;  n2 = length(x);  end;
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
Px(1) = Px(2);

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
7 y$ C% [8 y% g* U) f4 {那是不是要先把這些量化的值先做FFT
2 c% S, f, R' ?0 Z  _那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.
1 o( O6 H  a+ Z+ v

xyzabc0930

恩~~謝謝kevin所提供的~~
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
那我把量化所產生出來的值(1,-1,1,1,-1.......)
4 X" ]# E# H1 ^  S1 w3 ]. z& g5 E那請問一下我要怎麼轉換成左邊的圖
是經由FFT的嗎?
算出來的嗎?
7 N7 ^, k+ T+ f8 Z% E可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
9 K3 L4 Z5 @0 m5 |1 L5 V. mspec=fft(v.*ds_hann(N))/(N/4);
: t# Q  Q" H# }8 V2 C8 V3 hplot(linspace(0,1,N/2),
7 ^9 T4 U' u5 s- udbv(spec(1:N/2)))
* S. D/ F8 ]4 _axis([0 1 -120 0]);
( @0 g8 `' [  O0 T, T( g0 |grid on;
ylabel('dB')
( `. ?8 u& O5 [  t: N9 r6 l
從左邊這張圖可以去算出屬於它的SNR,參考附件.

xyzabc0930

請問kevin一下
那N指的是量化器的位元數?
還是什麼意思?
再麻煩回答一下
' L! N1 B$ U5 {6 W5 u% t謝謝你喔

xyzabc0930

再問一個問題一下~~~
/ X6 I" R" r( z$ n6 R$ T) _那如果我是用
( t2 n7 `, x* ]: B- G6 EHann-windowed fft.
" T. G# O, Z9 b7 L/ |7 n$ e1 o求出PSD之後
  i9 M5 |! I* d2 a% [9 N* v那再使用附件所附的檔案
再去求SNR是嗎?

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.

* K3 S: s8 R6 b) I" W3 ~Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

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

BTTB

最近也正在學習matlab解析SDM的SNR
2 R# Q4 Y# e) Q6 {* [3 p/ P一直搞不出來，慢慢研究中

s13571357

最近也正在學習matlab解析PSD
; j4 N# P5 S) C9 D一直搞不出來，慢慢研究中

mm770910

最近也正在學習matlab解析PSD
希望有進展