早教吧作业答案频道 -->数学-->
请高手用MATLAB帮忙解下微分方程组教下：Dy(1)=y(2);Dy(2)=y(3)^2*u*A*N^2/(2*m*(L-y(1))^2)-C/m*y(2)-k*(l+y(1))/m-(f+F)/m;Dy(3)=(L-y(1))/(u*A*N^2)*(V-y(3)*y(2)*u*A*N^2/(L-y(1))^2-R*y(3));初步带入数据如下：m=0.6;N=500;V=24;C=0.5
题目详情
请高手用MATLAB帮忙解下微分方程组教下：
Dy(1)=y(2);
Dy(2)=y(3)^2*u*A*N^2/(2*m*(L-y(1))^2)-C/m*y(2)-k*(l+y(1))/m-(f+F)/m;
Dy(3)=(L-y(1))/(u*A*N^2)*(V-y(3)*y(2)*u*A*N^2/(L-y(1))^2-R*y(3));
初步带入数据如下：
m=0.6;
N=500;
V=24;
C=0.5;
k=500;
l=0.0095;
f=16;
F=7;
L=0.0149;
r=0.034;
u=1.256e-6;
R=1;
A=pi*r^2;
时间t=10s，初始值y=【0 0 0】
Dy(1)=y(2);
Dy(2)=y(3)^2*u*A*N^2/(2*m*(L-y(1))^2)-C/m*y(2)-k*(l+y(1))/m-(f+F)/m;
Dy(3)=(L-y(1))/(u*A*N^2)*(V-y(3)*y(2)*u*A*N^2/(L-y(1))^2-R*y(3));
初步带入数据如下：
m=0.6;
N=500;
V=24;
C=0.5;
k=500;
l=0.0095;
f=16;
F=7;
L=0.0149;
r=0.034;
u=1.256e-6;
R=1;
A=pi*r^2;
时间t=10s，初始值y=【0 0 0】
▼优质解答
答案和解析
clear;clc
[T,Y]=ode45('nobk',[0 5],[0;0;0]) %微分方程有问题,求出来的都是0
plot(T,Y(:,1))
---------------------
function dy=nobk(t,y)
m=0.6;
N=500;
V=24;
C=0.5;
k=500;
l=0.0095;
f=16;
F=7;
L=0.0149;
r=0.034;
u=1.256e-6;
R=1;
A=pi*r^2;
dy=[0;0;0];
Dy(1)=y(2);
Dy(2)=y(3)^2*u*A*N^2/(2*m*(L-y(1))^2)-C/m*y(2)-k*(l+y(1))/m-(f+F)/m;
Dy(3)=(L-y(1))/(u*A*N^2)*(V-y(3)*y(2)*u*A*N^2/(L-y(1))^2-R*y(3));
[T,Y]=ode45('nobk',[0 5],[0;0;0]) %微分方程有问题,求出来的都是0
plot(T,Y(:,1))
---------------------
function dy=nobk(t,y)
m=0.6;
N=500;
V=24;
C=0.5;
k=500;
l=0.0095;
f=16;
F=7;
L=0.0149;
r=0.034;
u=1.256e-6;
R=1;
A=pi*r^2;
dy=[0;0;0];
Dy(1)=y(2);
Dy(2)=y(3)^2*u*A*N^2/(2*m*(L-y(1))^2)-C/m*y(2)-k*(l+y(1))/m-(f+F)/m;
Dy(3)=(L-y(1))/(u*A*N^2)*(V-y(3)*y(2)*u*A*N^2/(L-y(1))^2-R*y(3));
看了 请高手用MATLAB帮忙解下...的网友还看了以下：
偏微分方程∂^2*u/∂x^2+∂^2*u/∂y^2=2(x^2+y^2)的数值解法..定解条件u 2020-05-13 …
函数(x,y)=x^4+y^4-2(x-y)^2,先求偏导,得到（0,0）（根2,-根2）（-根2 2020-05-13 …
不等式组x-y+5≥0,x+y≥0,x≤3表示的平面区域是什么形状?若u=(x+2)^2+(y+1 2020-05-16 …
椭圆方程式的题椭圆方程x^2/a^2+y^2/b^2=t怎么用y=f(u),x=f(u)表示如果椭 2020-05-16 …
试用三元一次方程组的解法解四元一次方程组X+Y+Z+U=4,X+Y-Z-U=-2,X-Y+Z-U= 2020-05-21 …
请高手用MATLAB帮忙解下微分方程组教下：Dy(1)=y(2);Dy(2)=y(3)^2*u*A 2020-07-21 …
验证方程组F(x，y，u，v)＝u2+v2−x2−y＝0G(x，y，u，v)＝−u+v−xy+1＝ 2020-07-31 …
设X~U(0,1),U(0,1),且X与Y相互独立,求关于z的二次方程Xz^2+z+Y=0有实根的概 2020-11-01 …
Uxx+2Uxy-3Uyy=0,x属于R,y>0,U(x,0)=3X^2,Uy(x,0)=0求解U( 2020-11-01 …
matlab解微分方程组出错?>>[x,y,z,v,w]=dsolve('Dx=-2*(A+A1)* 2020-12-12 …