matlab编的龙格库塔法解微分方程，求改错

最新推荐文章于 2024-04-03 15:51:35 发布

原创

最新推荐文章于 2024-04-03 15:51:35 发布 · 2.5k 阅读

3 ·

CC 4.0 BY-SA版权

文章标签：

#MATLAB #龙格库塔法 #解微分方程

使用MATLAB的ode15s函数基于龙格库塔法求解微分方程，文章中给出了具体的数学模型、程序代码及变量定义，但在运行时遇到错误，作者寻求帮助修改程序。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

@[TOmatlab编的龙格库塔法解微分方程，求改错

数学模型：
在这里插入图片描述
主程序数据：
R1=0.019; R2=0.022; R3=0.044; R4=0.062;% 一级杆直径，二级杆直径，泵径,油管径
L1=1100; L2=500; %L1碳纤维杆的长度，L2钢杆长度
L=5; %单元节长
h1=1600;h2=1047; %泵挂深度，动液面
hc=h1-h2; %泵沉没度
S1=1900; S2=7800; S3=850; %碳纤维的密度，钢杆密度,井液密度
E1=1.110^11;
E2=2.110^11; %钢杆弹性模量
Ps=4.510^6; %沉没压力
Pd=13.510^6; %排出压力
ke=410^6; c0=10;%抽油杆悬挂系统的等效弹簧常数,和阻尼
A1=pi/4R1^2; A2=pi/4R2^2; A3=pi/4R3^2; A4=pi/4*R4^2; %碳纤维杆截面积，钢杆截面积,抽油泵柱塞横截面积,油管截面积
c_=2;
n=64; %总单元节数
n1=44; n2=20;
L_0=1;%余隙容积高度
Rqy=0.8; %吸入口气液比
T=10;
%
x0=zeros(128,1); %初始条件
t1=0;t2=10;n=201;%,初始点，结束点，求解间隔
[t,x]=ode15s(@longgefa,linspace(t1,t2,n),x0);%求解过程
函数：
function [dx] = longgefa(t,x)