有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
因為我在用這個軟體在跑SDM
那我想要請問一下
# N3 L9 B  M3 g, d! ?我要查一個東西
是功率頻譜密度
. S. n3 `( }9 l" w5 D: Z2 c那如果我要用Simulink下去看的話
請問一下Block 要打什麼
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
8 u* ~7 w$ x- _. d" w3 Z請各位教我一下~~
謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
, i* m0 D, `! k' w! R: y5 D3 k都是寫Code的方式的?
有人可以分享一下經驗的嗎?

hcjian

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

" v) C6 e6 ]1 L要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
7 [% M7 r/ Y+ h3 F, {2 `Peroidogram
; c1 p3 q9 \! h" Z4 @Modified periodogram
1 z6 v+ W  o- V! E% ?) jBartlett method
0 x% e+ d9 `( }( }& ]Welch method
' K2 j, W3 G4 D& i; a4 VBlackman-Tukey method

* @4 A" M& y& w; I  H0 g( n+ [; G若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
n1  : starting index, x(n1)
3 I9 l1 A% t+ c' ~5 X1 A- r: Mn2  : ending index, x(n2)

function Px = per(x,n1,n2)
x = x(;
if nargin == 1
& o+ z% T2 E& X& g5 E/ B5 M    n1 = 1;  n2 = length(x);  end;
/ D# P7 B; S+ T! y7 F9 d1 x6 D$ r3 sPx = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
Px(1) = Px(2);

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
; N2 v8 v$ P+ B+ ^) |- m那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
' C5 i. C& j9 I- i! S那是不是要先把這些量化的值先做FFT
3 u' K5 b! V/ x" {+ g! S那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.

% g) e: n7 G# g* d

xyzabc0930

恩~~謝謝kevin所提供的~~
7 O# m  @9 j6 M# i: K不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
那我把量化所產生出來的值(1,-1,1,1,-1.......)
% q1 \# j: `4 [) @( ^那請問一下我要怎麼轉換成左邊的圖
4 O- ?  Z% L/ `% s是經由FFT的嗎?
' c, P# \' m  q; `' [算出來的嗎?
可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
6 H- u8 E; N1 qspec=fft(v.*ds_hann(N))/(N/4);
plot(linspace(0,1,N/2),
+ t; N4 E6 d: h  a& idbv(spec(1:N/2)))
axis([0 1 -120 0]);
. H( {0 V4 @. E; J0 S5 G% agrid on;
+ h% M! d4 G2 \# N3 b3 _ylabel('dB')

從左邊這張圖可以去算出屬於它的SNR,參考附件.
, q& p; e8 }& Z- y/ T3 Q+ v2 d

xyzabc0930

請問kevin一下
那N指的是量化器的位元數?
還是什麼意思?
' q. P# R; M7 G4 |9 S再麻煩回答一下
# B1 q2 Z/ }8 K謝謝你喔

xyzabc0930

再問一個問題一下~~~
+ v5 ]4 O7 z0 B; K+ [5 D6 b那如果我是用
# ^( d7 z0 _8 K2 U+ w" @Hann-windowed fft.
求出PSD之後
那再使用附件所附的檔案
再去求SNR是嗎?

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
% C! z% F: }1 {" x/ N
* g2 W9 A, ], D) c  cHann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

謝謝您的分享，有助於學習PLL modeling，
7 k9 a# }% l; V8 o謝謝您

BTTB

最近也正在學習matlab解析SDM的SNR
4 m6 s* U6 g0 K0 O) T, h" }一直搞不出來，慢慢研究中

s13571357

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

mm770910

最近也正在學習matlab解析PSD
7 E' W+ w* n1 F% j8 H: u8 d+ B希望有進展