%For?a?given?omega?the?function?computes?the?transform?function
%V=exact（omega）
%V-->transfer?function
%omega->2*pi*f
%the?rest?parameters?must?be?changed?inside?the?function

clear?all;
clc;

R=0;
G=.1;
L=1;
C=1;
Rs=2;
Rl=1;
Cl=1/（4*pi）;
h?=?0.125;

omega?=?-16*pi:0.05*pi:16*pi;
number?=?max（size（omega））;

for?k=1:number
????gama=sqrt（（R+i*L*omega（k））*（G+i*C*omega（k）））;
????Zo=sqrt（（R+i*omega（k）*L）/（G+i*omega（k）*C））;
????Zl=Rl/（1+i*omega（k）*Cl*Rl）;
????A=（Zl-Zo）/（Zl+Zo）;
????exact_T（k）=（1+A）/（（1+Rs/Zo）*exp（gama*h）+A*（1-Rs/Zo）*exp（-gama*h））;
end

figure;
plot（omegaabs（exact_T）‘b‘）;?%axis（[-16*pi?16*pi?0?.5]）;
title（‘Exact?Transfer?function‘）;
xlabel（‘\omega‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????670??2007-11-15?05:28??analy.m
?????文件?????????230??2007-11-15?15:25??cal_rms.m
?????文件?????????143??2007-11-15?04:25??clfftf.m
?????文件?????????391??2007-11-15?16:26??err.m
?????文件?????????752??2007-11-23?02:04??impedance.m
?????文件?????????191??2007-11-15?03:52??source1.m
?????文件?????????231??2007-11-15?03:52??source2.m
?????文件????????2516??2007-11-25?04:26??timeUpdt.m
?????文件?????????460??2007-10-28?00:44??trans.m