有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
因為我在用這個軟體在跑SDM
- F+ _8 h' r; k. Z! n$ E2 ^6 n+ z4 b那我想要請問一下
& v( Y; u" e% i2 T& |我要查一個東西
4 H9 d% e1 n: I是功率頻譜密度
2 G& f4 G( L, _' Y" `那如果我要用Simulink下去看的話
請問一下Block 要打什麼
5 s9 A" y8 K7 m0 T7 Z才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
$ w# P+ r. [: G請各位教我一下~~
8 t; }6 b  c4 i0 j0 O' W謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
- [$ s- s  f4 K5 y: C都是寫Code的方式的?
* D9 K6 V6 Z9 f- k! i0 y7 ]% O有人可以分享一下經驗的嗎?

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
) V1 q/ x" Y5 x8 m! y* V不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。

; n, L1 |1 M1 s0 s要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
Peroidogram
Modified periodogram
Bartlett method
Welch method
Blackman-Tukey method
" ?( R+ `# F( G" Y* p# n& n
若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
n1  : starting index, x(n1)
* H( X9 ~6 f! g8 In2  : ending index, x(n2)
. }' R( g6 R" t" V; i9 U1 r
function Px = per(x,n1,n2)
# ]9 t* Y- L/ `/ }  Q0 lx = x(;
. a, g" ?; F. d, nif nargin == 1
    n1 = 1;  n2 = length(x);  end;
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
7 \1 Y2 n" q, n0 T' j& u% x5 \Px(1) = Px(2);

xyzabc0930

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

kevin

power Spectrum 跟 SNR 是不同的圖.
, g6 N9 L! l- \# R" e( d

xyzabc0930

恩~~謝謝kevin所提供的~~
# |" Q$ B* }% m' ?8 o! E不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
' k& C) X4 ]* e) y2 [/ H0 B那我把量化所產生出來的值(1,-1,1,1,-1.......)
# G6 i/ v1 [+ ]0 e- g1 {; Y那請問一下我要怎麼轉換成左邊的圖
5 x* g5 ]+ c' j$ F9 c, G8 e$ E是經由FFT的嗎?
% L# a; f9 W! q8 \) R: v$ ?算出來的嗎?
) }2 E+ t$ [& u可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
spec=fft(v.*ds_hann(N))/(N/4);
+ A; |; S) P8 q8 n- R# l0 `plot(linspace(0,1,N/2),
) w. r0 Z( _5 M1 \! ~% a8 Sdbv(spec(1:N/2)))
axis([0 1 -120 0]);
grid on;
( n% c1 r6 G. `+ d4 B1 Vylabel('dB')
/ D5 [# O1 b; e/ u
從左邊這張圖可以去算出屬於它的SNR,參考附件.

xyzabc0930

請問kevin一下
那N指的是量化器的位元數?
& a6 J' e8 h; G( |; T" X6 y8 D9 r3 H還是什麼意思?
再麻煩回答一下
% Y8 d3 p0 W# N5 L- H* c謝謝你喔

xyzabc0930

再問一個問題一下~~~
% c# q" l, r  o( v那如果我是用
Hann-windowed fft.
$ q! }! ^0 `0 ?% y求出PSD之後
# }- w' k' b* ]- N4 ~那再使用附件所附的檔案
再去求SNR是嗎?

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
  d/ K9 w7 \3 d6 [# y$ V0 \+ F
Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

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

BTTB

最近也正在學習matlab解析SDM的SNR
' l7 q, Z* z* O8 N9 B5 a一直搞不出來，慢慢研究中

s13571357

最近也正在學習matlab解析PSD
1 B( w" g* F& j  x3 R8 O一直搞不出來，慢慢研究中

mm770910

最近也正在學習matlab解析PSD
% m' g8 L7 s" o* _1 T希望有進展