有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
1 E" }" a! x! m- I8 Y8 W2 J因為我在用這個軟體在跑SDM
那我想要請問一下
我要查一個東西
是功率頻譜密度
$ T8 a$ L* h- q那如果我要用Simulink下去看的話
請問一下Block 要打什麼
/ L" c1 U" }( z6 }9 I7 s才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
請各位教我一下~~
. o& g6 D2 ?% s* p0 H# p1 n謝謝

xyzabc0930

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

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。

要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
& ?% b2 ?2 W- P. ^6 q- H  IPeroidogram
* }8 a: p# h3 k& e) {( R$ dModified periodogram
Bartlett method
Welch method
" C- B4 @6 q0 OBlackman-Tukey method
; d7 p7 r5 m$ V4 p! m# }
若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
n1  : starting index, x(n1)
( I. t# [, D, D0 j# R- j2 ]) v8 Cn2  : ending index, x(n2)
- L+ r) q3 X5 \; X: d
# }) T) i, ~) @0 Hfunction Px = per(x,n1,n2)
x = x(;
if nargin == 1
    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.......等等這種值丫
那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
那是不是要先把這些量化的值先做FFT
那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.

; o: b+ m1 _$ h' h1 o

xyzabc0930

恩~~謝謝kevin所提供的~~
2 A) z& d5 E% p6 d不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
那我把量化所產生出來的值(1,-1,1,1,-1.......)
那請問一下我要怎麼轉換成左邊的圖
; ~0 z4 F9 V# Q: B是經由FFT的嗎?
! L" y$ Q4 S% k3 t5 k9 H算出來的嗎?
; _. S' H( L9 `7 l可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
! i+ P8 L  H+ |% k: ~spec=fft(v.*ds_hann(N))/(N/4);
, C. l. s: n. F# R( o1 Cplot(linspace(0,1,N/2),
dbv(spec(1:N/2)))
axis([0 1 -120 0]);
grid on;
ylabel('dB')
' O8 ]% j2 S7 i# l9 R; v
從左邊這張圖可以去算出屬於它的SNR,參考附件.

xyzabc0930

請問kevin一下
那N指的是量化器的位元數?
/ M1 m: v0 ?% Q( P$ C# [7 ]還是什麼意思?
再麻煩回答一下
$ q6 f# s: X: X/ [' U謝謝你喔

xyzabc0930

再問一個問題一下~~~
那如果我是用
+ l( T- i# t2 u7 b2 b0 KHann-windowed fft.
1 J& i7 h* r! @2 r6 H' e求出PSD之後
那再使用附件所附的檔案
2 q7 M1 ?- g+ F/ D2 Q! f* G再去求SNR是嗎?

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
! s9 i" ^( K1 R; n7 o+ R
# m/ K$ s" T# B* E0 {" uHann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

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

BTTB

最近也正在學習matlab解析SDM的SNR
2 K' g: V' K" w" y  a  d4 d一直搞不出來，慢慢研究中

s13571357

最近也正在學習matlab解析PSD
$ {( ^* l5 \, c: d一直搞不出來，慢慢研究中

mm770910

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