有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
8 g; }) j8 B9 e2 {& @因為我在用這個軟體在跑SDM
% @: o9 n0 c& ^+ B# f! G( }那我想要請問一下
我要查一個東西
8 j+ c5 w6 ]! C4 G( ]是功率頻譜密度
& G4 m! K, d3 u1 H3 H! E* \. E那如果我要用Simulink下去看的話
5 f, H: b+ D; G- Y. p2 \7 H請問一下Block 要打什麼
  {7 H: T* ]6 B+ x) O才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
請各位教我一下~~
1 V9 }( x9 P' v* h, t- _! j謝謝

mm770910

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

s13571357

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

BTTB

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

ainge

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

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
) u0 @/ v) A% x6 @, B
Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

xyzabc0930

再問一個問題一下~~~
% S2 ^6 u5 J# F那如果我是用
- @: s2 V5 i% w+ m( K/ ]2 F1 q& SHann-windowed fft.
2 m5 H$ O" `: B) h  J, O1 P求出PSD之後
" M" O. ~$ C1 W那再使用附件所附的檔案
再去求SNR是嗎?

xyzabc0930

請問kevin一下
, y. \  s# Z5 ^# \0 z% R& I5 m那N指的是量化器的位元數?
1 P& @, T( l# L8 k. N. n還是什麼意思?
再麻煩回答一下
- A  B% @; H4 o6 h$ z) @( J& n4 g% Q謝謝你喔

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
spec=fft(v.*ds_hann(N))/(N/4);
; N5 D& a: ]! }* m6 eplot(linspace(0,1,N/2),
' u& w* e. s! H. Ndbv(spec(1:N/2)))
axis([0 1 -120 0]);
grid on;
  n  e% ?" Q2 _5 ^ylabel('dB')
9 t' m. a( G; e/ k! @; V: H
從左邊這張圖可以去算出屬於它的SNR,參考附件.
* Q  m4 o: k. d( P; Z7 ?

xyzabc0930

恩~~謝謝kevin所提供的~~
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
: V1 r$ @; `+ W: Y+ U- B+ q那我把量化所產生出來的值(1,-1,1,1,-1.......)
0 R; r5 [4 p+ [: |7 B2 J9 q% x那請問一下我要怎麼轉換成左邊的圖
! A) p* ?8 D- v- y/ [' X是經由FFT的嗎?
( g* |" X& Y1 t+ Q算出來的嗎?
- [& T6 z: k4 F可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.

5 X# r8 g. B- d" {! G

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
. g& `  z# B: |* u那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
那是不是要先把這些量化的值先做FFT
7 t3 v- Y% m8 X8 a6 ?% D( T0 E那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power spectrum estimation 有幾種方式,
Peroidogram
Modified periodogram
Bartlett method
/ D/ F  w7 E, F! M1 WWelch method
& @) R% \0 y$ _  RBlackman-Tukey method

若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
) Q4 R2 b4 C4 F% ]# gn1  : starting index, x(n1)
n2  : ending index, x(n2)
1 d! p5 M1 P0 `5 u0 u' z4 [2 X  @. O
  ]- Q1 K$ v: T$ h+ [( y: k2 c6 ofunction Px = per(x,n1,n2)
2 P7 H5 s$ z4 {; s9 Jx = x(;
if nargin == 1
$ Q9 l. d1 V. T; o    n1 = 1;  n2 = length(x);  end;
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
0 b& v+ Q# ^4 s. a8 oPx(1) = Px(2);

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
/ j1 ]9 k0 d  {不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。
* _' u% ]% @2 a! w$ Q
要算SNR的話，最好還是自己寫。

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
都是寫Code的方式的?
有人可以分享一下經驗的嗎?