求助优化一个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/07/02 00:16:22
![求助优化一个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](/uploads/image/z/1340355-3-5.jpg?t=%E6%B1%82%E5%8A%A9%E4%BC%98%E5%8C%96%E4%B8%80%E4%B8%AAmatlab%E7%A8%8B%E5%BA%8F%2Cfunction+kuttasym+tt%3D0.001kd%3D6%2A10%5E%28-4%29%2A%5B42.62+20.16%3B20.16+26.88%5Dkt%3D6%2A10%5E%28-4%29%2Ft%2A2%2A%5B132.2+9.6%3B9.6+12.8%5Dktt%3D6%2A10%5E%28-4%29%2F%28t-0.001%29%2A2%2A%5B132.2+9.6%3B9.6+12.8%5Df%3D%5B-1000+0%5D%27u1%3Dinv%28kd%2Bkt%29%2Afu1%3Dinv%28kd%2Bkt%29%2A%28f%2Bkt%2Au1%29t%3D0.002kd%3D6%2A10%5E%28-4%29%2A%5B4)
求助优化一个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