非线性微分方程 matlab,用ODE45求解Matlab中微分方程的非线性系统

我正在尝试使用ODE45 MATLAB解决非线性微分方程系统，我做了很多次成功，但这次我得到以下错误，我真的不知道有什么不对，我很困惑 . 这是代码 .

%%这是错误：

Subscript indices must either be real positive integers or logicals.

Error in non_L_ss (line 6)

(-Fk*(ds0+x(3)-x(1))+Fk*ds0-Fc(x(4)-x(2)))/ms +Fa/ms ] ;

Error in odearguments (line 87)

f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ode45 (line 115)

odearguments(FcnHandlesUsed,solver_name,ode, tspan,y0,optio varargin);

Error in solve (line 50)

[t X]=ode45(@non_L_ss,t_span,IC);

%%方程式在函数中定义：

function dX=non_L_ss(t,x)

global Fk Fc kt Fa q ds0 ms mu

dX=[ x(2);

(Fk*(ds0+x(3)-x(1))-Fk*ds0+Fc*(x(4)-x(2))-kt*x(1))/mu-Fa/m-kt*q/mu ;

x(4);

(-Fk*(ds0+x(3)-x(1))+Fk*ds0-Fc(x(4)-x(2)))/ms +Fa/ms ] ;

end

%%，这里函数被ODE45调用来解决：

clear

clc

global Fk Fc kt Fa q ds0 ms mu qdot v2

mu = 100 ;

ms = 1242 ;

k1s = 80000 ;

k2s = 32000 ;

kt = 405000 ;

c1s = 4000 ;

c2s = 1600 ;

v = 20 ;

Gq = 256e-6 ;

ds0 = 0.1538 ;

a = 1 ;

b = 0.001 ;

n0 = 0.1 ;

f0 = 0.011*v ;

w = 0.5 ;

Fa = 2000 ;

q = 0.05 ;

xs = 0.1 ;

xu = 0.1 ;

dxs = 0.1 ;

dxu = 0.2 ;

Fk = k1s+k2s*(ds0+xs-xu).^2 ;

if dxs >= dxu

Fc = c1s ;

elseif dxs < dxu

Fc = c2s ;

end

t_span=[0 1];

IC=[2 3 2 2];

[t X]=ode45(@non_L_ss,t_span,IC);

非线性微分方程 matlab,用ODE45求解Matlab中微分方程的非线性系统相关推荐

Matlab中用ode45求解速率方程，一直显示输入参量太多，什么原因呢？
脉冲泵浦时的m文件为: function dy = rate_eq(t,y,T0,N_T,d) P_in=80; R=0.8; sigma = 8e-25; %铒离子受激发射截面 sigma_y=1. ...
【matlab】ode45求解二阶微分方程，绘制曲线图 | 使用函数句柄的方法
朋友问题: 有微分方程如下: md2ydt2+dydtexp(t)−y2=5m \frac{d^2y}{dt^2} + \frac{dy}{dt} exp(t) - y^2 = 5mdt2d2y+d ...
matlab ode45设置步长,MATLAB中用ode45求解微分方程，如何设置最大步长？
如果你用过simulink里的ode45配置,我觉得你就会发现高赞就是在扯淡,你给的时间序列只是采样点,根本不是设置步长用的.我最近偷懒不想用simulink就研究了下ode45的函数配置项,在mat ...
齐次弦振动方程的matlab解法,ode45求解振动微分方程
function dy=order22(t,y);%用积分形式定义函数: dy=zeros(16,1); mc=854.97;mp=745.26;ms=332.69;m1=2135.1;m2=94.2 ...
机器学习（MACHINE LEARNING）MATLAB中微分方程的求解
文章目录 1 MATLAB之极限.积分.微分 2 matlab中微分方程的求解 2.1 一阶微分方程 2.2 求解二阶线性微分方程是指含有未知函数及其导数的关系式.解微分方程就是找出未知函数.微分方 ...
【Matlab 控制】微分方程 ode45() 求解并绘制曲线
Matlab 微分方程 ode45 求解并绘制曲线 2. 用 ode45() 求解 2.1 ode45() 函数用法 2.2 示例:求解一阶微分方程 2.2.1 Matlab 代码如下 2.2.2 代 ...
用Matlab求解高等数学中的问题（不定积分，定积分，多元微积分，微分方程，级数问题）
用Matlab求解高等数学中的问题一.求解不定积分利用matlab符号工具箱中的求积函数int, 可求函数的不定积分,int函数的调用格式如下: int(S) int(S,v) 说明: (1)in ...
matlab欧拉方程求解微分方程并和ode45对比结果
1.内容简介 matlab欧拉方程求解微分方程并和ode45对比结果 2.内容说明略 3.仿真分析 clc close all clear %% ode45方法 y0 = [8.5;2;1];%初始 ...
MATLAB实战应用案例：欧拉法、改进欧拉法、ode45求解微分方程实例
前言 ode45求解 clc clear 以下是我为大家准备的几个精品专栏,喜欢的小伙伴可自行订阅,你的支持就是我不断更新的动力哟! MATLAB-30天带你从入门到精通 MATLAB深入理解高级教程 ...
matlab程序中ode45,关于matlab中ode45的问题
关于matlab中ode45的问题关注:261 答案:2 mip版解决时间 2021-02-02 07:51 提问者森很绿却致人迷途 2021-02-02 04:45 function wlc ...

非线性微分方程 matlab,用ODE45求解Matlab中微分方程的非线性系统

非线性微分方程 matlab,用ODE45求解Matlab中微分方程的非线性系统相关推荐

最新文章

热门文章