作者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