有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
& b- y5 J1 S: V因為我在用這個軟體在跑SDM
" t  j; h2 @! X3 T! V7 n' C那我想要請問一下
我要查一個東西
6 f9 U! A+ }, j- H* P是功率頻譜密度
8 i0 |0 z6 X& K: o  D& M! y那如果我要用Simulink下去看的話
請問一下Block 要打什麼
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
' j- O) ~2 g0 |% `1 q請各位教我一下~~
* `+ U+ V1 y/ e# O5 m謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
# C6 B  n. j9 M0 j" Y) k6 g+ v都是寫Code的方式的?
9 f2 [" a0 P2 a) ^5 A有人可以分享一下經驗的嗎?

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
; c: P( Y! X2 M+ q6 D! X不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。
( R- B4 ^7 m  r, v& C0 V
( E: G& h3 N& `要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
Peroidogram
Modified periodogram
Bartlett method
Welch method
3 H2 `& }* p0 Y7 w/ kBlackman-Tukey method

5 y5 U. e. A2 u5 g+ V0 E若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
" v4 r3 g, X# ~  z$ k# x  qx   : input sequence
n1  : starting index, x(n1)
n2  : ending index, x(n2)

function Px = per(x,n1,n2)
x = x(;
1 G6 c& i# N, J5 X- Gif nargin == 1
$ t) w# `9 ^0 B% b    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.......等等這種值丫
5 @) D; Z' ], }7 w" g, I那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
那是不是要先把這些量化的值先做FFT
  {% R" C: a" m! U" o' d那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.
4 I! _% Y! m  N  z" [/ S

xyzabc0930

恩~~謝謝kevin所提供的~~
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
那我把量化所產生出來的值(1,-1,1,1,-1.......)
. X7 Y! o  K. F那請問一下我要怎麼轉換成左邊的圖
' Z4 o& ~+ T! Y) j4 {是經由FFT的嗎?
算出來的嗎?
6 y" T6 t/ R# W0 [可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
spec=fft(v.*ds_hann(N))/(N/4);
) S/ L- C& Q5 B% ~% J6 \2 lplot(linspace(0,1,N/2),
9 o; `+ Q& F5 a* [dbv(spec(1:N/2)))
/ C+ v7 d9 U" G, _axis([0 1 -120 0]);
2 L' A0 T4 N. rgrid on;
ylabel('dB')
8 X* D0 G8 U. l7 V
從左邊這張圖可以去算出屬於它的SNR,參考附件.
# V8 _" C5 T& R2 S1 E

xyzabc0930

請問kevin一下
' q0 I2 t/ r& j! j那N指的是量化器的位元數?
還是什麼意思?
再麻煩回答一下
謝謝你喔

xyzabc0930

再問一個問題一下~~~
那如果我是用
. m' f" {* w& x( z5 W5 |2 sHann-windowed fft.
求出PSD之後
) S7 i* k3 Y- \" y, d/ T; `, r那再使用附件所附的檔案
1 H! b! M' i% h再去求SNR是嗎?

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
- X0 G; o$ [+ ?- R  }; Z
Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

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

BTTB

最近也正在學習matlab解析SDM的SNR
. @7 |2 Z' \+ v% }8 K% r* x一直搞不出來，慢慢研究中

s13571357

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

mm770910

最近也正在學習matlab解析PSD
; |4 ]' V- i) M( \( L  ]9 `, e希望有進展