有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
因為我在用這個軟體在跑SDM
那我想要請問一下
# p9 v+ u0 F* u6 _8 A0 Y5 f8 Q' {8 `我要查一個東西
是功率頻譜密度
那如果我要用Simulink下去看的話
5 L4 V) c0 _4 C" I4 g2 {請問一下Block 要打什麼
" Z/ y: ^) f, c. F9 C才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
請各位教我一下~~
謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
: C) o3 i& @/ @都是寫Code的方式的?
& h  @% g. @( o* Z" j; q. G' l有人可以分享一下經驗的嗎?

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
! h4 f2 U2 b) f2 y: F5 ^& n3 O不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。

1 k0 f! E0 w4 l. A要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
Peroidogram
1 \" W, V& D. E% o9 N! ^. Q1 YModified periodogram
Bartlett method
Welch method
. K% ?) B6 ~" ]% h2 f; `. bBlackman-Tukey method

, _. X, `# p6 M# @) L若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
n1  : starting index, x(n1)
, t9 y' Z7 t9 B- @n2  : ending index, x(n2)

3 z2 g! E: b9 ~$ y+ jfunction Px = per(x,n1,n2)
x = x(;
" g, J) ~2 G2 e/ j3 e& Hif nargin == 1
% A" Z- z2 i, d* _% `    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.......等等這種值丫
9 z0 w5 a, y: B6 m: ]# I那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
那是不是要先把這些量化的值先做FFT
" Q6 ?( ]# q! I5 w/ u. d/ q那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.
, W2 s/ n) e; @
, h9 d6 e) u' [, _2 G

xyzabc0930

恩~~謝謝kevin所提供的~~
; e+ L$ z- Y, d1 B2 J不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
那我把量化所產生出來的值(1,-1,1,1,-1.......)
那請問一下我要怎麼轉換成左邊的圖
是經由FFT的嗎?
算出來的嗎?
可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
spec=fft(v.*ds_hann(N))/(N/4);
' j! [+ x4 E( W$ Bplot(linspace(0,1,N/2),
dbv(spec(1:N/2)))
4 h- Q' ~- L! n: Z  Iaxis([0 1 -120 0]);
grid on;
ylabel('dB')

, |; N. i9 a/ v) R7 y從左邊這張圖可以去算出屬於它的SNR,參考附件.
8 \5 h: M1 @: Q

xyzabc0930

請問kevin一下
) q6 M& ~& f# g# g1 G6 p9 ?那N指的是量化器的位元數?
還是什麼意思?
7 [2 E" t* }0 k! a再麻煩回答一下
1 p9 R+ q' x7 J9 M. o2 j( Y謝謝你喔

xyzabc0930

再問一個問題一下~~~
; f& K) @  M% d! P- s那如果我是用
Hann-windowed fft.
6 O1 Q( C0 `7 }$ H# M求出PSD之後
那再使用附件所附的檔案
再去求SNR是嗎?

kevin

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

Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

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

BTTB

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

s13571357

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

mm770910

最近也正在學習matlab解析PSD
% u7 V$ a# I0 g; J希望有進展