Chip123 科技應用創新平台

 找回密碼
 申請會員

QQ登錄

只需一步,快速開始

Login

用FB帳號登入

搜索
1 2 3 4
查看: 7767|回復: 14
打印 上一主題 下一主題

[問題求助] 有關Simulink

  [複製鏈接]
跳轉到指定樓層
1#
發表於 2008-5-17 23:02:15 | 只看該作者 回帖獎勵 |倒序瀏覽 |閱讀模式
我想要請問一下各位有用過Simulink的先進8 K+ i" Q* r, u; k7 G
因為我在用這個軟體在跑SDM; x2 v9 ~& O( `2 S# {7 Q
那我想要請問一下
# N3 L9 B  M3 g, d! ?我要查一個東西/ ~0 T; a) G, L7 P2 D
是功率頻譜密度
. S. n3 `( }9 l" w5 D: Z2 c那如果我要用Simulink下去看的話2 o2 {& x7 o3 I% F6 j
請問一下Block 要打什麼5 {1 }+ R/ c  m4 k# W9 L9 |4 u% \
才可以看我用Simulink所做出來的SDM跑出來的功率頻譜密度是多少勒
8 u* ~7 w$ x- _. d" w3 Z請各位教我一下~~- a/ U: r8 A; y$ e. g0 U
謝謝
分享到:  QQ好友和群QQ好友和群 QQ空間QQ空間 騰訊微博騰訊微博 騰訊朋友騰訊朋友
收藏收藏 分享分享 頂4 踩 分享分享
2#
 樓主| 發表於 2008-5-19 13:46:45 | 只看該作者
還是大家要用Simulink來看功率頻譜密度
, i* m0 D, `! k' w! R: y5 D3 k都是寫Code的方式的?% @+ d( S, e& }$ E$ C4 w8 j# n8 p
有人可以分享一下經驗的嗎?
3#
發表於 2008-5-20 15:24:18 | 只看該作者
通常是用periodogram,但是它只有2的冪次方點數可以選,所以你跑出來的FFT會有其他的tone出現,3 E1 x* L7 `0 u) V
不是很準的(本來用FFT的話,點數沒有取對,就不會準,而且取不同段,跑出來差很多)。5 S2 w* |' _. ]) K" B( F

" v) C6 e6 ]1 L要算SNR的話,最好還是自己寫。
4#
發表於 2008-5-21 11:20:37 | 只看該作者
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* K; j0 q/ B& S3 D4 g" N* z" }' y

* @4 A" M& y& w; I  H0 g( n+ [; G若用simple Peroidogram m file 方式如下,estimated 跟 actual 會有Bias. v2 r7 S" K8 |: z; P# i) h/ n+ G' a9 H
x   : input sequence" Z9 y% l, Z# {+ W8 n
n1  : starting index, x(n1)
3 I9 l1 A% t+ c' ~5 X1 A- r: Mn2  : ending index, x(n2)2 k- [) Z' v  r0 K$ L0 u0 G
0 b3 o  r# m, Y# `# Q: l: e- p
function Px = per(x,n1,n2)! p2 `9 B0 y1 [% K
x = x(;* n# j% _  A' b9 t6 _0 p2 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);/ ^1 J" Z  Z9 d/ S
Px(1) = Px(2);
5#
 樓主| 發表於 2008-5-21 14:48:21 | 只看該作者
我想要問的是因為在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做完了之後~~那就可以看到了嗎?
6#
發表於 2008-5-21 21:32:36 | 只看該作者
power Spectrum 跟 SNR 是不同的圖., \; b' `; J- t  Y4 E3 Y4 I

% g) e: n7 G# g* d

本帖子中包含更多資源

您需要 登錄 才可以下載或查看,沒有帳號?申請會員

x
7#
 樓主| 發表於 2008-5-22 00:31:39 | 只看該作者
恩~~謝謝kevin所提供的~~
7 O# m  @9 j6 M# i: K不過像左邊這個圖~如果我是用Simulink兜SDM model出來的話" `4 j" z2 h# T, |  i" j# {
那我把量化所產生出來的值(1,-1,1,1,-1.......)
% q1 \# j: `4 [) @( ^那請問一下我要怎麼轉換成左邊的圖
4 O- ?  Z% L/ `% s是經由FFT的嗎?
' c, P# \' m  q; `' [算出來的嗎?& e+ F# s" v2 N- e# w: F
可是不是也可以從左邊這張圖去算出SNR是多少的嗎?
8#
發表於 2008-5-22 12:02:43 | 只看該作者
左邊的圖是把輸出做FFT,假設v是輸出(+1,-1.....),但是要用windowed-FFT,底下例子是用Hann-windowed fft.
6 H- u8 E; N1 qspec=fft(v.*ds_hann(N))/(N/4);' H0 N( f: H7 T+ v1 K; `2 I$ s
plot(linspace(0,1,N/2),
+ t; N4 E6 d: h  a& idbv(spec(1:N/2)))3 A& u0 I$ M" k1 V- Y" ]1 I* A
axis([0 1 -120 0]);
. H( {0 V4 @. E; J0 S5 G% agrid on;
+ h% M! d4 G2 \# N3 b3 _ylabel('dB'), y7 w) Y5 n" j) Y' ]9 T) F
- M( W3 S7 W$ {( g( \
從左邊這張圖可以去算出屬於它的SNR,參考附件.
, q& p; e8 }& Z- y/ T3 Q+ v2 d

本帖子中包含更多資源

您需要 登錄 才可以下載或查看,沒有帳號?申請會員

x
9#
 樓主| 發表於 2008-5-22 13:56:21 | 只看該作者
請問kevin一下, F: n/ e8 @) ?! x
那N指的是量化器的位元數?. q% X% O# R+ ]3 `
還是什麼意思?
' q. P# R; M7 G4 |9 S再麻煩回答一下
# B1 q2 Z/ }8 K謝謝你喔
10#
 樓主| 發表於 2008-5-22 14:18:58 | 只看該作者
再問一個問題一下~~~
+ v5 ]4 O7 z0 B; K+ [5 D6 b那如果我是用
# ^( d7 z0 _8 K2 U+ w" @Hann-windowed fft.4 \6 {, r3 W3 `& s- Y. w# o
求出PSD之後- b) m& h# |8 }+ d$ M
那再使用附件所附的檔案  _$ H1 [* q$ k9 d( l
再去求SNR是嗎?
11#
發表於 2008-5-22 16:47:56 | 只看該作者
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.

本帖子中包含更多資源

您需要 登錄 才可以下載或查看,沒有帳號?申請會員

x
12#
發表於 2009-8-19 18:47:52 | 只看該作者
謝謝您的分享,有助於學習PLL modeling,
7 k9 a# }% l; V8 o謝謝您
13#
發表於 2011-10-25 15:23:28 | 只看該作者
最近也正在學習matlab解析SDM的SNR
4 m6 s* U6 g0 K0 O) T, h" }一直搞不出來,慢慢研究中
14#
發表於 2012-1-16 10:23:05 | 只看該作者
最近也正在學習matlab解析PSD3 y; [% w; h# ^% S9 r3 r
一直搞不出來,慢慢研究中
15#
發表於 2012-4-12 15:19:42 | 只看該作者
最近也正在學習matlab解析PSD
7 E' W+ w* n1 F% j8 H: u8 d+ B希望有進展
您需要登錄後才可以回帖 登錄 | 申請會員

本版積分規則

首頁|手機版|Chip123 科技應用創新平台 |新契機國際商機整合股份有限公司

GMT+8, 2024-11-1 06:58 AM , Processed in 0.205011 second(s), 20 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回復 返回頂部 返回列表