有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
因為我在用這個軟體在跑SDM
8 _2 [0 {2 ]# \( f% m& {那我想要請問一下
我要查一個東西
$ I. y6 f: X5 v4 t/ j2 _& f是功率頻譜密度
那如果我要用Simulink下去看的話
1 M# m  a' W+ S) C請問一下Block 要打什麼
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
請各位教我一下~~
: }( i  j6 n3 [# \# Q/ I謝謝

xyzabc0930

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

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。
/ k& w( B$ @8 I- V* b
要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
/ w2 |' n3 a4 z  n( x# K# TPeroidogram
Modified periodogram
Bartlett method
Welch method
) i5 K) p/ `3 I" k$ q: [) X) }Blackman-Tukey method

若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
n1  : starting index, x(n1)
3 F* G2 v6 L5 ^, {$ D) sn2  : ending index, x(n2)
' k/ K1 Q" ?0 C+ F
% [! i5 ]+ u$ r& efunction Px = per(x,n1,n2)
! _9 K6 U1 {4 N/ nx = x(;
if nargin == 1
    n1 = 1;  n2 = length(x);  end;
$ ]# S5 x7 b' o( qPx = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
Px(1) = Px(2);

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
$ n7 U" n  V) Z. ?/ \那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
那是不是要先把這些量化的值先做FFT
, d" ~$ f3 q( {1 J! W! r那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.
1 m9 O+ c( X! L# N- \
) {- V) f( n+ f9 U4 J

xyzabc0930

恩~~謝謝kevin所提供的~~
$ [$ @! x; c5 v6 M! T不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
9 _. u1 t! z* X7 E  _" l5 E1 C7 p那我把量化所產生出來的值(1,-1,1,1,-1.......)
那請問一下我要怎麼轉換成左邊的圖
是經由FFT的嗎?
算出來的嗎?
) C, c# J  P/ c可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
spec=fft(v.*ds_hann(N))/(N/4);
plot(linspace(0,1,N/2),
dbv(spec(1:N/2)))
! u7 `4 K+ o% i9 I* v' ~3 I' h6 naxis([0 1 -120 0]);
grid on;
5 {' V& ^' Y4 ]- Eylabel('dB')

- y' h& C" g% b* j( \從左邊這張圖可以去算出屬於它的SNR,參考附件.

xyzabc0930

請問kevin一下
那N指的是量化器的位元數?
還是什麼意思?
再麻煩回答一下
0 O) C( S  e1 \2 n3 X謝謝你喔

xyzabc0930

再問一個問題一下~~~
6 _$ J0 u5 k. @- V7 Y& l那如果我是用
Hann-windowed fft.
( [4 X! O9 @4 c求出PSD之後
3 T% X* ^5 y' v4 Q; g6 T那再使用附件所附的檔案
: j8 M, B9 S- A& ]7 ~1 T再去求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
希望有進展