R=normrnd(50,1,100,1),h=normrnd(25*1e-3,1.3*1e-3,100,1),d=normrnd(4.5*1e-3,1.1*1e-3,100,
1);
Rs=R,RL=R,Rne=R,Rfe=R,h1=h,h2=h,r=0.4*1e-3,r1=r,r2=r,uo=4*pi*1e-7,Eo=1/(36*pi)*1e-9,Er=3,E
e=2/3,l=1.25,t=0.6*1e-3,t1=t,t2=t;
lG=uo/(2*pi)*log(2*h1/r1),lR=uo/(2*pi)*log(2*h2/r2),lm=uo/(4*pi)*log(1+4.*h1*h2/d.^2);
L=[lG lm;lm lR];
c1=1/(2*pi*Eo)*(1/Er*log(1/r1)+Ee*log(1/(r1+t1))-log(1/(2*h1))),c2=1/(4*pi*Eo)*log(1+4*h1*h2
/d^2),c3=c2,c4=1/(2*pi*Eo)*(1/Er*log(1/r2)+Ee*log(1/(r2+t2))-log(1/(2*h2)));
C=[c1 c2;c3 c4]^-1;
Z=i*w*L,Y=i*w*C;
y=(Z*Y)^(1/2);
TI^-1*Y*Z*TI=y^2,Zc=Z*TI*y^-1*TI^-1;
Q11=Zc*cosh(l*(Y*Z)^(1/2))*Zc^-1,Q12=-sinh(l*(Z*Y)^(1/2))*Zc,Q21=-sinh(l*(Y*Z)^(1/2))*Zc^-1,
Q22=Zc^-1*cosh(l*(Z*Y)^(1/2))*Zc;
T11=-Q12^-1*Q11,T12=-Q12^-1,T21=Q21-Q22*Q12^-1*Q11,T22=Q22*Q12^-1;
[k11 k12;k21 k22]=[Rs 0;0 Rne]*[T11+T12([1/RL 0;0 1/Rfe]-T22)^-1*T21]+[1];
Hnext=1/(k12-k11*k21^-1*k22);
am=20*log10(abs(Hnext));
plot(w/(2*pi),am);
set(gca,'XScale','log');
title('均值');
评论0