求助优化一个matlab程序,function kuttasym tt=0.001kd=6*10^(-4)*[42.62 20.16;20.16 26.88]kt=6*10^(-4)/t*2*[132.2 9.6;9.6 12.8]ktt=6*10^(-4)/(t-0.001)*2*[132.2 9.6;9.6 12.8]f=[-1000 0]'u1=inv(kd+kt)*fu1=inv(kd+kt)*(f+kt*u1)t=0.002kd=6*10^(-4)*[4
来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/19 03:40:52
求助优化一个matlab程序,function kuttasym tt=0.001kd=6*10^(-4)*[42.62 20.16;20.16 26.88]kt=6*10^(-4)/t*2*[132.2 9.6;9.6 12.8]ktt=6*10^(-4)/(t-0.001)*2*[132.2 9.6;9.6 12.8]f=[-1000 0]'u1=inv(kd+kt)*fu1=inv(kd+kt)*(f+kt*u1)t=0.002kd=6*10^(-4)*[4
求助优化一个matlab程序,
function kutta
sym t
t=0.001
kd=6*10^(-4)*[42.62 20.16;20.16 26.88]
kt=6*10^(-4)/t*2*[132.2 9.6;9.6 12.8]
ktt=6*10^(-4)/(t-0.001)*2*[132.2 9.6;9.6 12.8]
f=[-1000 0]'
u1=inv(kd+kt)*f
u1=inv(kd+kt)*(f+kt*u1)
t=0.002
kd=6*10^(-4)*[42.62 20.16;20.16 26.88]
kt=6*10^(-4)/t*2*[132.2 9.6;9.6 12.8]
ktt=6*10^(-4)/(t-0.001)*2*[132.2 9.6;9.6 12.8]
u2=inv(kd+kt)*f
u2=inv(kd+kt)*(f+kt*u2+ktt*u1)
按此循环.
求助优化一个matlab程序,function kuttasym tt=0.001kd=6*10^(-4)*[42.62 20.16;20.16 26.88]kt=6*10^(-4)/t*2*[132.2 9.6;9.6 12.8]ktt=6*10^(-4)/(t-0.001)*2*[132.2 9.6;9.6 12.8]f=[-1000 0]'u1=inv(kd+kt)*fu1=inv(kd+kt)*(f+kt*u1)t=0.002kd=6*10^(-4)*[4
function kutta
k = 1;
f=[-1000 0]'
kd=6*10^(-4)*[42.62 20.16;20.16 26.88]
for t= 0.001 : 0.001: 0.005
kt=6*10^(-4)/t*2*[132.2 9.6;9.6 12.8]
ktt=6*10^(-4)/(t-0.001)*2*[132.2 9.6;9.6 12.8]
u(k)=inv(kd+kt)*f
u(k)=inv(kd+kt)*(f+kt*u(k))
k = k+1;
end