[問題] 信號處裡的問題..

[cfleu (世事難預料)]

※ [本文轉錄自 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