有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
因為我在用這個軟體在跑SDM
那我想要請問一下
9 u1 \9 r; h9 W2 E0 O我要查一個東西
是功率頻譜密度
, k& k4 L6 ^2 J: v3 |! G7 [那如果我要用Simulink下去看的話
% r* S- `6 l( m2 T7 \- D  k, j請問一下Block 要打什麼
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
, a+ q, r9 H: ~6 a% y請各位教我一下~~
謝謝

mm770910

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

s13571357

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

BTTB

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

ainge

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

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
! Q/ z% w% g, Y0 t% `  G( H
" _  I! I9 ?) `/ A* ZHann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

xyzabc0930

再問一個問題一下~~~
那如果我是用
: }* M9 z) d& o( t, B  Q5 K9 X. l3 C! NHann-windowed fft.
9 {# \$ b8 z9 T/ g# \求出PSD之後
那再使用附件所附的檔案
9 ^. ~2 N6 x. X再去求SNR是嗎?

xyzabc0930

請問kevin一下
! z$ |0 T* u% A8 o6 p1 ?7 x那N指的是量化器的位元數?
還是什麼意思?
再麻煩回答一下
9 q. a2 h! |$ K/ z謝謝你喔

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
2 z/ k; P1 S0 R) }spec=fft(v.*ds_hann(N))/(N/4);
plot(linspace(0,1,N/2),
7 V- `; |3 ~+ Odbv(spec(1:N/2)))
  ]9 \) f7 V6 t& K2 l4 ?0 i1 b4 caxis([0 1 -120 0]);
grid on;
ylabel('dB')

從左邊這張圖可以去算出屬於它的SNR,參考附件.
5 ~) P) z* P& C( r+ R& Q  L' F  ?- e

xyzabc0930

恩~~謝謝kevin所提供的~~
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
' g0 e# G; n) r! S0 {  F/ U那我把量化所產生出來的值(1,-1,1,1,-1.......)
那請問一下我要怎麼轉換成左邊的圖
& C2 a4 G+ ?. p" G$ r是經由FFT的嗎?
算出來的嗎?
可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.
( ^/ P! P1 ?$ g; x6 R$ a

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
+ _! ]. A9 n! k, G7 }. S' G; f% F那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
4 r0 }% @* Y; W+ E那是不是要先把這些量化的值先做FFT
! |# a# t/ n& Z4 v2 ^* j那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power spectrum estimation 有幾種方式,
Peroidogram
' u+ \, n( k- o9 lModified periodogram
Bartlett method
0 Q5 P4 ]9 x. ^5 q2 tWelch method
Blackman-Tukey method
" ^8 ], Q% y( y! Q' L7 p- f, y
7 e1 ^' v) X/ ^9 x# D若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
. c1 ^7 Z3 Z. m7 p4 mn1  : starting index, x(n1)
& o9 T! t. r$ k7 g5 [n2  : ending index, x(n2)

function Px = per(x,n1,n2)
x = x(;
: p4 g; }; L1 V) x% U. w- Oif nargin == 1
- [/ x; q2 j) L& [: m1 T* H! q    n1 = 1;  n2 = length(x);  end;
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
Px(1) = Px(2);

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
* X! a4 y2 Z% P不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。
/ n* L0 X& u+ I  y8 b6 n/ p
: _. h; }" V! g5 S2 ~! J& i要算SNR的話，最好還是自己寫。

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
都是寫Code的方式的?
) R1 D( E# ~0 u- z+ R有人可以分享一下經驗的嗎?