likeshu1103
2011-05-09, 13:28
t=0:10^-12:40000*10^-12; dt=10^(-12);N=40000;
l=1;xs=0.2;
a=2.5*10^(-3);b=9.345*10^(-3);
H=5.76*10^7; %对应文献符号:H=σ
erel=2.5;e0=8.85*10^(-12);u0=4*pi*10^(-7);
L=u0/2*pi*log(b/a);C=(2*pi*e0*erel)/(log(b/a));R=(1/pi*H)*(1/a+1/2*b*0.6*10^(-3));
RL1=100;RL2=10;
x=xs;
G=(1/2*pi)*sqrt(u0/H)*(1/a+1/b); %对应文献符号:G=ξ
v=1/sqrt(L*C);
a1=(G/2)*sqrt(C/L);a2=(R/2)*sqrt(C/L);
for k=0:N
dk=dt*k;
if dk<x/v
g1(k)=0;
else
g1(k)=exp(-a2.*x).*(a1.*x/(2.*sqrt(pi.*(k.*dt-x/v)^3)).*exp(-(a1.*x)^2/(4.*(k.*dt-x/v)));
end
end
有哪位大侠知道的指出来,不胜感激。
l=1;xs=0.2;
a=2.5*10^(-3);b=9.345*10^(-3);
H=5.76*10^7; %对应文献符号:H=σ
erel=2.5;e0=8.85*10^(-12);u0=4*pi*10^(-7);
L=u0/2*pi*log(b/a);C=(2*pi*e0*erel)/(log(b/a));R=(1/pi*H)*(1/a+1/2*b*0.6*10^(-3));
RL1=100;RL2=10;
x=xs;
G=(1/2*pi)*sqrt(u0/H)*(1/a+1/b); %对应文献符号:G=ξ
v=1/sqrt(L*C);
a1=(G/2)*sqrt(C/L);a2=(R/2)*sqrt(C/L);
for k=0:N
dk=dt*k;
if dk<x/v
g1(k)=0;
else
g1(k)=exp(-a2.*x).*(a1.*x/(2.*sqrt(pi.*(k.*dt-x/v)^3)).*exp(-(a1.*x)^2/(4.*(k.*dt-x/v)));
end
end
有哪位大侠知道的指出来,不胜感激。