有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
因為我在用這個軟體在跑SDM
那我想要請問一下
" F7 v; i% h; z, x我要查一個東西
8 `* T4 Y. K* ?. \4 ?: v! i# r2 H0 S' M是功率頻譜密度
那如果我要用Simulink下去看的話
- r% ~+ L9 z6 l+ K% A. g請問一下Block 要打什麼
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
0 w  u+ X1 R9 ?7 N# e請各位教我一下~~
6 g3 k7 E" t, |1 d5 {謝謝

xyzabc0930

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

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
9 g5 @5 c, Z; n1 A! a4 z不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。

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

kevin

power spectrum estimation 有幾種方式,
7 b" H# x( ]# ~2 d, B. HPeroidogram
. M7 {2 Y$ x( V4 @Modified periodogram
Bartlett method
Welch method
; G' ?. t. {$ ~$ zBlackman-Tukey method
# I6 s& f3 u2 j; P
* ^1 ]; U+ b1 a$ R# ^1 n若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
/ B: m$ W0 P+ {) L1 w4 Vx   : input sequence
& M! l8 d0 ?$ g' Yn1  : starting index, x(n1)
! }1 m. A6 V  Wn2  : ending index, x(n2)

function Px = per(x,n1,n2)
x = x(;
if nargin == 1
    n1 = 1;  n2 = length(x);  end;
' m* c$ `  {5 \/ d. \Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
( r3 b4 e* P- T" o1 O, Q. [4 `Px(1) = Px(2);

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
, r; K6 t( r7 j& e8 |那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
那是不是要先把這些量化的值先做FFT
/ C8 g  A0 x* b1 D6 Z那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.
6 v; V' K: x" }6 N% y
. ^7 K9 k  z3 H; y8 s; I

xyzabc0930

恩~~謝謝kevin所提供的~~
! C' t! \/ t! V' w% n不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
' Y2 g  g1 ?* F4 f那我把量化所產生出來的值(1,-1,1,1,-1.......)
, m4 @- V# v2 b* Z那請問一下我要怎麼轉換成左邊的圖
是經由FFT的嗎?
算出來的嗎?
可是不是也可以從左邊這張圖去算出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),
+ H. E8 K2 Y+ _' Edbv(spec(1:N/2)))
% h) j- B' Q' h% X( S1 Z( H3 F' _axis([0 1 -120 0]);
1 R; o9 d* g5 T4 q+ \  wgrid on;
" u( _" p, E) U' Qylabel('dB')
7 |' h! c* }: B# V
* l0 Z" R& F% D7 W- k從左邊這張圖可以去算出屬於它的SNR,參考附件.

xyzabc0930

請問kevin一下
那N指的是量化器的位元數?
還是什麼意思?
再麻煩回答一下
; z4 o. w$ q# o! c0 Z, R1 n謝謝你喔

xyzabc0930

再問一個問題一下~~~
/ W2 ~6 K( B* X2 [$ V1 n8 a那如果我是用
Hann-windowed fft.
4 D: @0 m& v1 D. l* Z求出PSD之後
$ B: D' Z- t9 ~" Z* {% E5 X- Z那再使用附件所附的檔案
* a8 e* l. V6 g再去求SNR是嗎?

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
; [" B# K7 ?2 }
( F+ I* D  N2 h5 K* LHann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

謝謝您的分享，有助於學習PLL modeling，
. U0 N- C7 Z0 g* U" z" ?謝謝您

BTTB

最近也正在學習matlab解析SDM的SNR
' ]5 k3 |. b! X3 z- s; B一直搞不出來，慢慢研究中

s13571357

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

mm770910

最近也正在學習matlab解析PSD
5 ]# y, E7 [% b6 ]  G  C2 ]希望有進展