matlab代码clear all;
clc;
syms r1 r2 c1 c2 s
pz=@(z1,z2) z1*z2/(z1+z2)
sz=@(z1,z2) z1+z2
hs11=pz(1/s/c1,sz(1/s/c2,r2))/sz(r1,pz(1/s/c1,sz(1/s/c2,r2)))
hs12=r2/sz(1/s/c2,r2)
hs1=hs11*hs12
hs1=collect(hs1,s)
Z1=[r1+1/s/c1,1/s/c1;1/s/c1,1/s/c1]
Y1=inv(Z1)
T1=[-Y1(2,2)/Y1(2,1),-1/Y1(2,1);-Y1(2,2)*Y1(1,1)/Y1(2,1)+Y1(1,2),-Y1(1,1)/Y1(2,1)]
Z2=[r2+1/s/c2,r2;r2,r2]
Y2=inv(Z2)
T2=[-Y2(2,2)/Y2(2,1),-1/Y2(2,1);-Y2(2,2)*Y2(1,1)/Y2(2,1)+Y2(1,2),-Y2(1,1)/Y2(2,1)]
T=T1*T2
hs2=1/T(1,1)
hs2=collect(hs2,s)
|