[MATLAB毕业设计] 请问这段代码错误在哪 - MATLAB爱好者论坛-LabFans.com

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
有哪位大侠知道的指出来，不胜感激。

2011-05-09, 13:28	#1
likeshu1103 初级会员注册日期: 2011-04-11 年龄: 36 帖子: 8 声望力: 0	请问这段代码错误在哪 t=0:10^-12:4000010^-12; dt=10^(-12);N=40000; l=1;xs=0.2; a=2.510^(-3);b=9.34510^(-3); H=5.7610^7; %对应文献符号：H=σ erel=2.5;e0=8.8510^(-12);u0=4pi10^(-7); L=u0/2pilog(b/a);C=(2pie0erel)/(log(b/a));R=(1/piH)(1/a+1/2b0.610^(-3)); RL1=100;RL2=10; x=xs; G=(1/2pi)sqrt(u0/H)(1/a+1/b); %对应文献符号：G=ξ v=1/sqrt(LC); a1=(G/2)sqrt(C/L);a2=(R/2)sqrt(C/L); for k=0:N dk=dtk; 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 有哪位大侠知道的指出来，不胜感激。