有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
' _3 V" o% m* j2 m& v因為我在用這個軟體在跑SDM
: z) f; S7 J! V$ I3 G那我想要請問一下
; I% i6 X1 M# ^6 h$ u我要查一個東西
是功率頻譜密度
/ S( N% f! g" @) u; x那如果我要用Simulink下去看的話
請問一下Block 要打什麼
1 _1 d" ~9 h+ W5 Q8 e9 j: r9 z才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
請各位教我一下~~
0 {- I5 s# r  V, [謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
6 ^6 H% {1 i* Q4 [8 k3 b6 V: _都是寫Code的方式的?
7 B* {3 x" G; m6 o' k: V有人可以分享一下經驗的嗎?

hcjian

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

$ k5 S; i! @" d+ y要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
Peroidogram
Modified periodogram
$ J, F: C& p6 q9 iBartlett method
  D6 D; G5 C$ D4 zWelch method
Blackman-Tukey method
+ S1 c" ^# w6 r' H/ L6 G  ^3 t+ G
' H8 j: j9 H9 H1 ^若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
- A9 |, g6 Q$ A+ h( P! b# z! k& _, Zx   : input sequence
: H2 R3 o/ I5 y7 ?+ ^0 C* v9 i# mn1  : starting index, x(n1)
n2  : ending index, x(n2)

! ?# G: j) i" @: N) t2 O8 T) Nfunction Px = per(x,n1,n2)
x = x(;
if nargin == 1
    n1 = 1;  n2 = length(x);  end;
/ q. t: z& B) v: Y6 LPx = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
Px(1) = Px(2);

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
" F; i1 l% D5 e, b5 _4 q那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
那是不是要先把這些量化的值先做FFT
; [0 G  p+ z  D- J4 G7 g& _* k那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.

; N7 J; ^' o8 `4 T& G

xyzabc0930

恩~~謝謝kevin所提供的~~
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
那我把量化所產生出來的值(1,-1,1,1,-1.......)
那請問一下我要怎麼轉換成左邊的圖
2 }" B, G$ \- {/ q" c1 d是經由FFT的嗎?
$ q- F& ~) @# U1 y( v, M9 O- D算出來的嗎?
可是不是也可以從左邊這張圖去算出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)))
axis([0 1 -120 0]);
! L' H# D6 T+ w0 M0 Egrid on;
ylabel('dB')

6 t2 q# \9 k0 D' h! C從左邊這張圖可以去算出屬於它的SNR,參考附件.

xyzabc0930

請問kevin一下
: x) J' w) L7 ?) l) O8 Q那N指的是量化器的位元數?
" T. E* R* G( C7 o5 Z5 w0 G& {還是什麼意思?
再麻煩回答一下
謝謝你喔

xyzabc0930

再問一個問題一下~~~
那如果我是用
0 M* u" s- m. C( L% F& g' j, vHann-windowed fft.
求出PSD之後
* q5 x# \- _' A' x  }. c) ?那再使用附件所附的檔案
; ]  U) Z- x" H( |/ M再去求SNR是嗎?

kevin

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

3 y' u5 Q4 N6 D# G5 UHann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

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

BTTB

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

s13571357

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

mm770910

最近也正在學習matlab解析PSD
/ p# b! g2 x$ ~2 |' s0 N希望有進展