用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/04/28 08:45:58
用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.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等求解.我用后面的求解过,可以解,但是解比较奇怪,是常数.这与你的方程有关系

用matlab画,dv/dt=-v+v.^1.75 ,plot(v,t) 用matlab画dv/dt=2*g-9*u/p*v+1.152*(p/u).^(-0.25)*R.^(-1.25)*v.^1.75,plot(v,t)其中u=210000;p=760;R=0.001;g=9.8; 求 解微分方程组的matlab程序方程组是:dv/dt=-c*n*p/ta*v.^2*cx+g*sin(th)d(th)/dt=g*cos(th)/vdx/dt=v*cos(th)dy/dt=v*sin(th)dp/dt=p/(r*ta)*v*sin(th)已知:c,n,cx,g,r,ta 初值有:p,v,th,x,y,用符号表示,最后给定,如何编写matl 如何解阻力微分方程dv/dt=g-kv/m 为什么 dv/dt=v*dv/dy MATLAB解微分方程dv/(p/v-kv^2)=dt,常数是p,k,变量是v 用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 证明受阻力a=dv/dt=-kv^2 怎么推 v dv/dt =d(v^2/2)/dt 求一阶微分方程dv/dt=g-uv-kv^2的解其中v为所求函数u,v为常数,g为重力加速度. 大一物理中a =dv/dt 已知C,p,m都是常数,求解微分方程:dv/dt=g- (CpAv^2)/m 求出速度v的表达式,以及高度h(好)的表达式?已知C,p,m都是常数,求解微分方程:dv/dt=g- (CpAv^2)/m 求出速度v的表达式,以及高度h(v=dh/dt)的表 用Matlab求常微分方程的解析解dv/dt=(G-F-bv)/mG=527.436*0.4536*9.8=2344.6 F=470.327*0.4536*9.8=2090.7 m=527.436*0.4536=239.24 b=0.08*0.4536*9.8/0.3048=1.1667参考答案是v(t)=217.6224*(-exp(-0.0048776*t)+1) 要有完整的程序代 dv/dt=(dv*dx)/(dt*dx)的详细过程, 请帮我把这个式子g-(kv^2)/m=dv/dt化成v关于t的式子?g,k,m是常数.g-(v^2)k/m=dv/dt这样看的可能更清楚一楼的好象不对吧 解道积分的题g*sin f(t)dt=dv求f(t)9.8*=乘 这个微分方程怎么解,m(dv/dt)=mg-Cr²v²,C,r,g,m是常数 加速度微分方程疑惑a=dv/dt=d(dR/dt)/dt=d^2R/dt^2.中:最后的dt^2指的是(dt)^2还是d* t^2?