有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
因為我在用這個軟體在跑SDM
8 v  m  A. Q$ @( Y9 c6 J9 J4 H+ A$ ^那我想要請問一下
: `' o* V) @9 X+ e- Q* F; C. U我要查一個東西
1 {1 e" W. u4 l. G1 U/ ^- j是功率頻譜密度
那如果我要用Simulink下去看的話
請問一下Block 要打什麼
8 [* i$ F& q5 F' d0 I: u才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
" {0 J7 K& e8 ?6 p1 |請各位教我一下~~
" L* S. v( r1 ~謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
! C6 Z& q% B, N5 r% z都是寫Code的方式的?
& Q2 J8 d4 m& x1 M; f( \有人可以分享一下經驗的嗎?

hcjian

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

+ G. A/ i5 Z3 `8 M$ u要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
% n- x4 D+ M- R6 qPeroidogram
; m6 m+ ^' P* ]  j' e6 ^% j- XModified periodogram
Bartlett method
Welch method
Blackman-Tukey method

若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
0 s7 u3 ^% n1 k8 t3 n% nn1  : starting index, x(n1)
n2  : ending index, x(n2)
/ x( O- B/ e; H: M8 v0 i
function Px = per(x,n1,n2)
: w; ]9 {) }0 y  Bx = x(;
# o: B# `, p' t) m) m( Q# x1 Uif nargin == 1
6 k" W9 ]0 e8 F7 V    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
% U5 p  S2 m# O/ \那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.

xyzabc0930

恩~~謝謝kevin所提供的~~
$ |5 h7 D4 b1 `- C+ g8 v/ a不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
* ~9 Y1 ~( l3 q& Z& U1 b& e那我把量化所產生出來的值(1,-1,1,1,-1.......)
那請問一下我要怎麼轉換成左邊的圖
4 ]/ d2 B! l8 G# z1 q# T- ~是經由FFT的嗎?
算出來的嗎?
可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
9 S6 ^7 a- H8 [. yspec=fft(v.*ds_hann(N))/(N/4);
plot(linspace(0,1,N/2),
* N$ c. Z# l: t' y3 U( @0 udbv(spec(1:N/2)))
  u3 U1 Q; S' A$ E2 D& @  j$ baxis([0 1 -120 0]);
grid on;
ylabel('dB')
, o6 W! _' T0 I, m' |
3 k$ M+ s3 f( N6 {1 m2 D從左邊這張圖可以去算出屬於它的SNR,參考附件.
' N) D2 t7 J+ @+ `$ q- N/ k- u1 z# O

xyzabc0930

請問kevin一下
. v4 X: S! f- s. ~5 ?那N指的是量化器的位元數?
還是什麼意思?
6 _5 ]  g% |' s+ C再麻煩回答一下
謝謝你喔

xyzabc0930

再問一個問題一下~~~
那如果我是用
Hann-windowed fft.
/ \. I2 |5 _. r6 @" ^. m; a求出PSD之後
那再使用附件所附的檔案
- |* h9 C# T) A再去求SNR是嗎?

kevin

N 是Hann window 的長度,ds_hann(N)參考附件.跟MATLAB 原來的HANN window 不一樣.
6 b9 Q/ r( S$ i! j! G! Z3 [; l
Hann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

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

BTTB

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

s13571357

最近也正在學習matlab解析PSD
) x, T3 ?  D1 E5 b一直搞不出來，慢慢研究中

mm770910

最近也正在學習matlab解析PSD
3 L5 D+ v* ?7 S9 t1 \+ Z希望有進展