【高等数学】求用matlab得出某物理题中微分方程的数值解程序如下: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=3e
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/23 22:12:39
![【高等数学】求用matlab得出某物理题中微分方程的数值解程序如下: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=3e](/uploads/image/z/10926854-62-4.jpg?t=%E3%80%90%E9%AB%98%E7%AD%89%E6%95%B0%E5%AD%A6%E3%80%91%E6%B1%82%E7%94%A8matlab%E5%BE%97%E5%87%BA%E6%9F%90%E7%89%A9%E7%90%86%E9%A2%98%E4%B8%AD%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%9A%84%E6%95%B0%E5%80%BC%E8%A7%A3%E7%A8%8B%E5%BA%8F%E5%A6%82%E4%B8%8B%EF%BC%9Afunction+dx%3DDOF2_m%28t%2Cx%29dx%3Dzeros%284%2C1%29%3Bglobal+k%25-------------------L2%3D9%3B+++++++++L0%3D10%3B+++++++++E%3D2.05e11%3BI%3D0.063528%3B+++rou%3D7.85e3%3B+++A%3D0.114976%3B+++mc%3D1.2e4%3Bmd%3D3e)
【高等数学】求用matlab得出某物理题中微分方程的数值解程序如下: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=3e
【高等数学】求用matlab得出某物理题中微分方程的数值解
程序如下:
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);
【高等数学】求用matlab得出某物理题中微分方程的数值解程序如下: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=3e
运行没有问题,我把m文件和结果图都发给你,你看看,
你好 ,请问你是问这程序哪里有错吗