用matlab画dv/dt=2*g-4.5*u*v/p/R/R+0.969*(p/u).^(-0.25)*R.^(-1.25)*v.^1.75clear;clcu=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
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/17 08:27:47
用matlab画dv/dt=2*g-4.5*u*v/p/R/R+0.969*(p/u).^(-0.25)*R.^(-1.25)*v.^1.75clear;clcu=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
用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 ,帮忙看下程序,其中范围和初值也是我随便写的,不知道是不是这个的原因
用matlab画dv/dt=2*g-4.5*u*v/p/R/R+0.969*(p/u).^(-0.25)*R.^(-1.25)*v.^1.75clear;clcu=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等求解.我用后面的求解过,可以解,但是解比较奇怪,是常数.这与你的方程有关系