【高等数学】求用matlab得出某物理题中微分方程的数值解

来源：学生作业帮编辑：拍题作业网作业帮分类：数学作业时间：2024/04/29 07:47:44

程序如下：
function dx=DOF2_m(t,x)
dx=zeros(4,1);
global k
%-------------------
L2=9;
L0=10;
E=2.05e11;
I=0.063528;
rou=7.85e3;
A=0.114976;
mc=1.2e4;
md=3e4;
mj=6.2e5;
g=9.8;
m1=rou*A;
k=24*E*I/(L2^3);
m=2*m1*L2*156/420+mj;
M=m+mc;
%----------------------------------------------
if t<=3
    Q=3e4;
else
    Q=0;
end
%---------------------------------------------------------------
dx(1)=x(2);
dx(2)=(-md*L0*dx(4)*cos(x(3))+md*L0*(x(4)^2)*sin(x(3))-k*x(1)+Q)/(M+md);
dx(3)=x(4);
dx(4)=(-dx(2)*cos(x(3))-g*sin(x(3)))/L0;
end

运行以下文件求解：
global k
[t,y]=ode45('DOF2_m',[030],[0 0 0 0]);
Fx=-k*y(:,1)/2;
subplot(3,1,1)
Sta=y(:,3)*180/pi;
plot(t,Sta,'m')
grid on
xlabel('t/s','fontsize',18);
ylabel('¦È/rad','fontsize',18);

subplot(3,1,2)
plot(t,y(:,1))
grid on
xlabel('t/s','fontsize',18);
ylabel('Ux/m','fontsize',18);

subplot(3,1,3)
plot(t,Fx(:,1))
grid on
xlabel('t/s','fontsize',18);
ylabel('Fx/N','fontsize',18);

运行没有问题,我把m文件和结果图都发给你,你看看,

求大神帮忙用matlab求解微分方程的数值解请用Matlab求解下面二阶常微分方程的数值解用matlab求微分方程 y’=-y2+x,y(0)=1的数值解, MATLAB 求解微分方程数值解常微分方程数值解Matlab 求大神支招,我想用ode45来求解一个微分方程的数值解,但是MATLAB老是说有错,怎么办? matlab 求微分方程初值问题的符号解,并与数值解进行比较用matlab求解二阶微分方程数值解,程序出现错误,求大神指点如何处理matlab微分方程数值解,具体就是解微分方程组,得出一系列的（x,y,t）如何作sqrt（x^2+y^2)=z 高等数学求微分方程的通解 MATLAB 常微分方程数值解出错 matlab求解泛函微分方程数值解.