作者cfleu (世事難預料)
看板comm_and_RF
標題[問題] 信號處裡的問題..
時間Mon Oct 23 11:23:36 2006
※ [本文轉錄自 MATLAB 看板]
作者: cfleu (世事難預料) 看板: MATLAB
標題: [問題] 信號處裡的問題..
時間: Sun Oct 22 18:43:57 2006
這是類比信號
n
Σ ρ *exp(-j θ )*exp(-s t )
i=1 i i i
H(s) = --------------------------------
n
(Σ ρ ^2)^0.5
i=1 i
thta=thta+tow
y=(sum(ro.*exp(-j.*thta)))*exp(-j.*w);
x=(sum(ro.^(2))).^(0.5);
h=y/x;
ρ=[0.057662 0.176809 0.407163 0.303585 0.258782 0.061831 0.15034 0.051534
0.185074 0.400967 0.295723 0.350825 0.262909 0.225894 0.170996 0.149723
0.240140 0.116587 0.221155 0.259730]
t=[1.003019 5.422093 0.518650 2.751772 0.602895 1.016585 0.143556 0.153832
3.324866 1.935570 0.429948 3.228872 0.848831 0.073883 0.203952 0.194207
0.924450 1.381320 0.640512 1.368671]
單位：uS
θ=[4.855121 3.419109 5.864470 2.215894 3.758058 5.430202 3.952093 1.093586
5.775198 0.154459 5.928383 3.053023 0.628578 22.128544 1.099463 3.462951
3.664773 2.833799 3.334290 0.393889]
單位：rad
必須要畫出H(s)的大小響應與相位響應圖
希望有高手可以幫忙
我有大概寫出上面那個算式
不過問題出在畫圖
我曉得大小響應與相位響應的指令是abs與agu(?)
m = abs(h);
p = angle(h);
subplot(121);
semilogy(w,m,'linewidth',2);
title('Magnitude');
subplot(122);
plot(w,p*180/pi,'linewidth',2);
title('Phase');
但是不曉得如何完整的顯示全部的週期
不知道有沒有高手願意寫寫看跟我指導一下..
真的是很感激..
ps.
我晚一點會po我寫的程式碼
因為再另外一台電腦..等一下在補上
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 210.58.148.212
※ 編輯: cfleu 來自: 60.248.164.205 (10/23 11:23)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 60.248.164.205