有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
( d4 p3 ~( V" `) @( t因為我在用這個軟體在跑SDM
那我想要請問一下
我要查一個東西
# f) ~) L9 }3 s* o, g是功率頻譜密度
那如果我要用Simulink下去看的話
& a) j7 y) E5 k請問一下Block 要打什麼
8 K6 d- P! G7 t" w9 T7 a# B; t才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
請各位教我一下~~
謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
) U7 {9 w$ m; Y, F# t都是寫Code的方式的?
/ E+ ]5 s: [' c有人可以分享一下經驗的嗎?

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
2 b* ~+ V  D" G5 t9 h" G! R不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。
4 \$ P  E- J) e) j$ x
3 ^8 K4 U( ~) Q+ }要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
Peroidogram
/ X2 A- ^) B3 N/ r: A' [Modified periodogram
Bartlett method
Welch method
Blackman-Tukey method

  H6 [  w) Q( R- r0 A* R若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
# d* ^4 w$ k# H- ?. g9 ]x   : input sequence
n1  : starting index, x(n1)
n2  : ending index, x(n2)

8 O# z, E: I( a* X6 Bfunction Px = per(x,n1,n2)
, \$ i, C2 J: N' l. L) hx = x(;
) K+ i& }5 v1 D5 h. Q& ?if nargin == 1
3 @! Z0 j8 j4 \( `1 X# j1 |* n    n1 = 1;  n2 = length(x);  end;
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
0 d# }2 K0 l* ~Px(1) = Px(2);

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
5 e7 Z0 f$ y6 u9 N那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
那是不是要先把這些量化的值先做FFT
那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.

xyzabc0930

恩~~謝謝kevin所提供的~~
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
% @1 Z' S6 ~) C& S$ ^那我把量化所產生出來的值(1,-1,1,1,-1.......)
那請問一下我要怎麼轉換成左邊的圖
是經由FFT的嗎?
算出來的嗎?
可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
& [- G; n+ G! [4 K- r, kspec=fft(v.*ds_hann(N))/(N/4);
plot(linspace(0,1,N/2),
6 A. o( [! U2 k1 B' E+ Adbv(spec(1:N/2)))
axis([0 1 -120 0]);
, P$ m0 M" z8 S" B  _6 u% L1 F4 O9 ogrid on;
" B) G! c) ~$ E9 Q3 E) xylabel('dB')
' n) @! `+ N" _$ v% v) L! [* f
從左邊這張圖可以去算出屬於它的SNR,參考附件.

xyzabc0930

請問kevin一下
6 }, `; `. \+ i1 ]' [) `7 Q那N指的是量化器的位元數?
還是什麼意思?
再麻煩回答一下
& v0 O/ [% f$ b9 x謝謝你喔

xyzabc0930

再問一個問題一下~~~
8 A: l5 p7 v$ v. a# l4 u: Y那如果我是用
Hann-windowed fft.
' i2 _' q/ ^, ^求出PSD之後
* L& K$ x- \3 c  P' E! I那再使用附件所附的檔案
再去求SNR是嗎?

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
( }' ~5 M' V) Y# f4 F
& a; i$ c) z- E1 q2 h: ~' `Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

謝謝您的分享，有助於學習PLL modeling，
2 f: y: f4 X( Z6 A) o$ |謝謝您

BTTB

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

s13571357

最近也正在學習matlab解析PSD
. P; \+ r6 U" m4 R1 K, f4 d一直搞不出來，慢慢研究中

mm770910

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