Matlab求解微分方程组

我们采用ode方法：

（1）求解普通微分方程组：使用ode45方法

1. 创建一个函数文件eq2.m，在函数文件中描述这个解的微分方程组：

%eq2.m文件
%描述微分方程组function dy=eq2(t,y) %说明微分变量是二维的，令y(1)=x,y(2)=y
dy=zeros(2,1);
%微分方程组
dy(1)=5*(1-y(1))/sqrt((1-y(1))^2+(t-y(2))^2);
dy(2)=5*(1-y(2))/sqrt((1-y(1))^2+(t-y(2))^2);
end

2.外部调用ode45函数求解微分方程组：

[t,y]=ode45(@eq2,[0,2],[0,0]);

ode45函数说明：第一个参数是方程的名称，第二个参数是指求解时t的范围，第三组参数是指y中每个元素的初值。

[t,y]=ode45(@eq2,[t1,t2],[y1(0),y2(0)]);

3.绘制响应曲线

 plot(t,y(:,1),'-',t,y(:,2),'*')

y(1)用-绘制，y(2)用*绘制（下面的两条曲线重合了）

Matlab求解微分方程组相关推荐

龙格库塔法matlab求解微分方程组,微分方程组的龙格库塔公式求解matlab版.pdf
微分方程组的龙格库塔公式求解matlab版微分方程组的龙格-库塔公式求解matlab版南京大学王寻 1. 一阶常微分方程组考虑方程组     y'f x,y,z , y x y ...
数学建模学习（29）：matlab求解微分方程组详细讲解，代码+案例讲解，学不会找我！
文章目录前言求解微分方程组求解矩阵微分方程组总结前言上一篇我已经详细讲过求微分方程,这一篇是对上一篇的补充,也就是变得稍微复杂一点,就是要求方程组了,如果你学会了上一篇,那么求解方程组其实 ...
MATLAB怎么解方程解,怎么用MATLAB求解微分方程组并画出解函数图？
!using["XSLSF"]; //使用命名空间XSLSF //数组xArray存放x的值:ti为当前有效值的个数:tmax为ti对应的时间:tmi ...
matlab求解微分方程组固有频率,电机固有频率计算
什么是固有频率? 鉴于本论坛是电机论坛,机械和力学方面相对弱些,因此为了便于大家理解起见,下面的讨论用机电类比的办法来加以说明,因为数学模型是完全一样的! 例如一个电感.电容.电阻可以构成一个振荡电路 ...
龙格库塔法解微分方程组的matlab程序,MATLAB实例源码教程：龙格库塔法求解微分方程组源代码实例.doc...
MATLAB实例源码教程:龙格库塔法求解微分方程组源代码实例.doc MATLAB实例源码教程龙格库塔法求解微分方程组源代码实例题目用经典 Runge-Kutta方法求下列一阶微分方程组的近似解y1 ...
matlab解二阶微分方程组,[微分方程组]急急急！用MATLAB按二阶龙格库塔法求解微分方程组，急用于毕业设计！...
急急急!用MATLAB按二阶龙格库塔法求解微分方程组,急用于毕业设计! 问题补充:今天才发现自己之前做的一点都不对,17号就交论文了,我傻了,急死了!求各位大侠帮帮忙.谢谢!要求解的微分方程如图所示. ...
matlab常系数线性矩阵微分方程组,基于Matlab常系数线性微分方程组的求解
·基础数学· 基于 Matlab 常系数线性微分方程组的求解* 严水仙 (赣南师范大学数学与计算机科学学院,江西赣州 341000) 摘要: 在常微分方程课程教学中,常系数线性微分方程组可以通过 ...
2021-01-13 Matlab求解微分代数方程 (DAE)
Matlab求解微分代数方程 (DAE) 什么是微分代数方程? 微分代数方程是一类微分方程,其中一个或多个因变量导数未出现在方程中.方程中出现的未包含其导数的变量称为代数变量,代数变量的存在意味着不能 ...
matlab解符号方程组,matlab 求解符号方程组
1特殊符号可爱组成的小狗图案求解符号方程组: 特殊符号可爱组成的小狗图案,缺失:matlab求解符号方程组4057/9 ▄██████▄ █████████▄ ███ ▄████▄▄▄▄███ ██ ...

Matlab求解微分方程组

Matlab求解微分方程组相关推荐

最新文章

热门文章