有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
因為我在用這個軟體在跑SDM
/ |' f2 ]/ G; m8 p( ~4 ?$ k# s那我想要請問一下
我要查一個東西
; l4 T. r9 Q2 d8 k; p% q是功率頻譜密度
' q: b0 ]* ]( U: u8 y6 G% G( N" K7 g9 j那如果我要用Simulink下去看的話
請問一下Block 要打什麼
: F. v( M9 Y9 m" [才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
請各位教我一下~~
謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
都是寫Code的方式的?
  b4 J# G$ {# U) a有人可以分享一下經驗的嗎?

hcjian

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

# R6 O! {* X+ J& ]要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
Peroidogram
Modified periodogram
Bartlett method
2 N* M/ T% m- C+ U+ R' j1 d$ I' eWelch method
Blackman-Tukey method
( X: X6 F1 O9 U6 j" U
0 @. I* ~3 w5 C9 r% v/ l6 E若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
n1  : starting index, x(n1)
n2  : ending index, x(n2)

0 J) k0 Z* [, \( s2 O( I; ~function Px = per(x,n1,n2)
+ n' v9 t0 u* w- J: Sx = x(;
if nargin == 1
5 @7 L# _. {. V8 |    n1 = 1;  n2 = length(x);  end;
/ x* @- a- K1 L& oPx = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
Px(1) = Px(2);

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
那是不是要先把這些量化的值先做FFT
- ?. ^9 T8 V7 p$ n# v9 i9 l% a' d那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.

7 ]& ?# V. L+ u3 x0 v) L

xyzabc0930

恩~~謝謝kevin所提供的~~
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
那我把量化所產生出來的值(1,-1,1,1,-1.......)
那請問一下我要怎麼轉換成左邊的圖
是經由FFT的嗎?
& a2 j2 _  A0 w1 G算出來的嗎?
& ~! p1 O; m( L0 X  a可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
2 f! X5 j. H4 m# w+ E, }spec=fft(v.*ds_hann(N))/(N/4);
! t; E& l5 T  p1 Jplot(linspace(0,1,N/2),
dbv(spec(1:N/2)))
2 W  Q* C9 L( S7 g7 Qaxis([0 1 -120 0]);
* r( y8 I8 t! D* e) h5 lgrid on;
ylabel('dB')
' X" _' b( c3 C. _" b+ ?
從左邊這張圖可以去算出屬於它的SNR,參考附件.
, V) d8 e& p, ?6 ]4 p, R- Q

xyzabc0930

請問kevin一下
* u; A% N: b5 d) d# b" ]& j- \那N指的是量化器的位元數?
還是什麼意思?
/ U) K$ A, p" {! y再麻煩回答一下
+ F$ t# B$ t1 W8 H  @7 M3 r9 |謝謝你喔

xyzabc0930

再問一個問題一下~~~
那如果我是用
8 q" D, r, z' cHann-windowed fft.
4 K( d# O6 c, l; \! L* a# K  m" w# E求出PSD之後
那再使用附件所附的檔案
3 z: V! z3 T7 a. v$ k6 v, c再去求SNR是嗎?

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
  ~5 D9 ^4 k; g2 _7 t
Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

謝謝您的分享，有助於學習PLL modeling，
; O& c" X+ `7 `5 h9 ]( q" U, F8 ^謝謝您

BTTB

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

s13571357

最近也正在學習matlab解析PSD
2 s8 d7 Y+ h' h) p+ D. p- ~一直搞不出來，慢慢研究中

mm770910

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