有關Simulink

xyzabc0930

我想要請問一下各位有用過Simulink的先進
因為我在用這個軟體在跑SDM
那我想要請問一下
$ b+ j3 D; [) |5 ]. \; x$ ~我要查一個東西
+ u! v3 R& S' m, J* }是功率頻譜密度
那如果我要用Simulink下去看的話
請問一下Block 要打什麼
& h8 _6 I; C% c3 J3 |才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
; |' Q6 O9 r3 x6 d$ Q' y* [/ v% u' U請各位教我一下~~
謝謝

xyzabc0930

還是大家要用Simulink來看功率頻譜密度
都是寫Code的方式的?
8 {8 x$ F# r7 e; f) F# W7 ]有人可以分享一下經驗的嗎?

hcjian

通常是用periodogram，但是它只有2的冪次方點數可以選，所以你跑出來的FFT會有其他的tone出現，
3 p0 t# g- L/ X1 w# K不是很準的（本來用FFT的話，點數沒有取對，就不會準，而且取不同段，跑出來差很多）。
# J1 ~. u& a# O3 {  Z& |9 G% W6 x
要算SNR的話，最好還是自己寫。

kevin

power spectrum estimation 有幾種方式,
Peroidogram
Modified periodogram
& V" h+ m" @9 ]Bartlett method
Welch method
3 H) x5 p. e5 a" g$ `3 }Blackman-Tukey method

若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias
x   : input sequence
n1  : starting index, x(n1)
& x( s0 B, `2 M2 ]% |* Wn2  : ending index, x(n2)

function Px = per(x,n1,n2)
x = x(;
if nargin == 1
    n1 = 1;  n2 = length(x);  end;
Px = abs(fft(x(n1:n2),1024)).^2/(n2-n1+1);
3 ]5 n/ f1 b3 _' r: P7 EPx(1) = Px(2);

xyzabc0930

我想要問的是因為在Simulink量化輸出的值是類似-1,1,1,-1,-1.......等等這種值丫
+ i( ?3 {  W, o- ^+ m那如果我要把這些值變成可以知道SNR(就是可能要先畫成PSD)
% `9 U0 A6 J+ d/ r, Z, Y那是不是要先把這些量化的值先做FFT
+ v1 q% ^5 ?. N! w5 b8 J* ^" x3 A那我可以請問各位~~如果我FFT做完了之後~~那就可以看到了嗎?

kevin

power Spectrum 跟 SNR 是不同的圖.
& B9 I* w$ ~6 n) U. w1 @. }
# @7 E  ]( |) ?

xyzabc0930

恩~~謝謝kevin所提供的~~
不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話
那我把量化所產生出來的值(1,-1,1,1,-1.......)
7 I4 z0 N3 V- e6 ]那請問一下我要怎麼轉換成左邊的圖
3 `( }! g7 Q5 ~1 B: ^' s是經由FFT的嗎?
- C) q) X# g8 |算出來的嗎?
可是不是也可以從左邊這張圖去算出SNR是多少的嗎?

kevin

左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
spec=fft(v.*ds_hann(N))/(N/4);
. S0 ~) S* s. {1 N8 Qplot(linspace(0,1,N/2),
dbv(spec(1:N/2)))
axis([0 1 -120 0]);
grid on;
ylabel('dB')
* k! ~5 x. ~+ {0 r  h: K
從左邊這張圖可以去算出屬於它的SNR,參考附件.
% ?% Q! o3 T0 |* f- \9 e9 Q

xyzabc0930

請問kevin一下
那N指的是量化器的位元數?
還是什麼意思?
$ P0 e) c2 k% b5 }* K4 k再麻煩回答一下
謝謝你喔

xyzabc0930

再問一個問題一下~~~
! F) o6 Y5 f$ Q$ v6 \那如果我是用
Hann-windowed fft.
求出PSD之後
那再使用附件所附的檔案
, X6 v+ b, u6 \7 Z1 h再去求SNR是嗎?

kevin

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

) w0 I) M9 G& y9 S3 C  l* BHann-windowed fft 之後再使用SNR附件所附的檔案去求SNR.

ainge

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

BTTB

最近也正在學習matlab解析SDM的SNR
% z/ X, |# e; l2 _* m6 S- b" M一直搞不出來，慢慢研究中

s13571357

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

mm770910

最近也正在學習matlab解析PSD
希望有進展