作业帮用matlab画dv/dt=2*g-4.5*u*v/p/R/R+0.969*(p/u).^(-0.25)*R.^(-1.
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/27 04:01:05
用matlab画dv/dt=2*g-4.5*u*v/p/R/R+0.969*(p/u).^(-0.25)*R.^(-1.25)*v.^1.75
clear;clc
u=210000;
p=760;
R=0.001;
g=9.8;
fs=@(t,v)(2*g-4.5*u*v/p/R/R+0.969*(p/u).^(-0.25)*R.^(-1.25)*v.^1.75);
[t,v]=ode45(fs,[0,1],0.1);
plot(t,v);
请教问什么总是显示busy ,帮忙看下程序,其中范围和初值也是我随便写的,不知道是不是这个的原因
clear;clc
u=210000;
p=760;
R=0.001;
g=9.8;
fs=@(t,v)(2*g-4.5*u*v/p/R/R+0.969*(p/u).^(-0.25)*R.^(-1.25)*v.^1.75);
[t,v]=ode45(fs,[0,1],0.1);
plot(t,v);
请教问什么总是显示busy ,帮忙看下程序,其中范围和初值也是我随便写的,不知道是不是这个的原因
首先希望你检查你的方程,因为你的方程中的数量级差别很大：2*g和4.5*u*v/p/R/R相差10个数量级左右（尽管v初值是0.1）.这也使得方程可能要采用其他求解方式,也就是说4-5变阶的龙哥库塔法可能无法求解（比如刚性的情况）.可以采用ode23s,ode15s等求解.我用后面的求解过,可以解,但是解比较奇怪,是常数.这与你的方程有关系