有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
( C$ I) T  Q% k- z) `' ]因為我在用這個軟體在跑SDM
那我想要請問一下
我要查一個東西
是功率頻譜密度
" F. Q4 P: q2 r% _# ^" }# V那如果我要用Simulink下去看的話
& {5 e' @4 U0 v4 Z5 T& T請問一下Block 要打什麼
; A( ]: }1 H6 b# k* K7 u" ?/ ?, m才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
請各位教我一下~~
6 K3 j% Q% K$ `* \0 i2 [3 Q1 S謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
都是寫Code的方式的?
1 J" [' K3 P) u) f# D( Y有人可以分享一下經驗的嗎?

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。
$ r8 I6 _8 E2 P) k' ?7 d% j0 a) Y
要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
8 ~' C, d) ?9 d8 Y% U* `; K* h4 qPeroidogram
9 i4 R( y+ O- F+ E$ v2 `Modified periodogram
: e; a( m" M: {* N) ?. v) wBartlett method
Welch method
Blackman-Tukey method
6 |+ c; `' k! E: F! k- X
1 v* {. s7 U7 Y( v& O若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
0 K( e0 U8 z$ y2 M; s* _. ux   : input sequence
n1  : starting index, x(n1)
n2  : ending index, x(n2)

function Px = per(x,n1,n2)
3 C! U, t$ J& Ox = x(;
if nargin == 1
    n1 = 1;  n2 = length(x);  end;
9 A8 E7 \! |+ \& N8 K( r5 cPx = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
Px(1) = Px(2);

xyzabc0930

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

kevin

power Spectrum 跟 SNR 是不同的圖.

' Y9 K* u5 ~# }8 l/ \  P

xyzabc0930

恩~~謝謝kevin所提供的~~
$ H4 l- ~8 I/ c9 }1 V# l3 t不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
. L2 S+ F" T/ T; W( Y, ]2 i/ W5 n那我把量化所產生出來的值(1,-1,1,1,-1.......)
那請問一下我要怎麼轉換成左邊的圖
+ R4 m6 o# P! I1 X$ N) B是經由FFT的嗎?
; `: `9 w5 R! A2 [0 e算出來的嗎?
# F5 x2 d* W$ k; F- ]! {5 Y( k可是不是也可以從左邊這張圖去算出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),
( K  p# G. i4 ^dbv(spec(1:N/2)))
axis([0 1 -120 0]);
grid on;
ylabel('dB')

: h2 G' r3 K% f' O" u  @從左邊這張圖可以去算出屬於它的SNR,參考附件.

xyzabc0930

請問kevin一下
- F/ d/ B" f& c+ z那N指的是量化器的位元數?
還是什麼意思?
再麻煩回答一下
謝謝你喔

xyzabc0930

再問一個問題一下~~~
那如果我是用
Hann-windowed fft.
求出PSD之後
1 F# Z" @; [$ N2 u$ V: F: _$ A7 t那再使用附件所附的檔案
再去求SNR是嗎?

kevin

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

! W0 K) ]! U+ I+ i/ sHann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

謝謝您的分享，有助於學習PLL modeling，
8 t3 d: `- D3 z謝謝您

BTTB

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

s13571357

最近也正在學習matlab解析PSD
, `" e4 S- G* L0 h8 h一直搞不出來，慢慢研究中

mm770910

最近也正在學習matlab解析PSD
7 p- }( r& A; q4 b& p" Z; V希望有進展