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)
|