Matlab求微分方程的数值解

您所在的位置：网站首页 › matlab求二阶微分方程的解析解 › Matlab求微分方程的数值解

Matlab求微分方程的数值解

2024-07-09 09:32:02| 来源: 网络整理| 查看: 265

注：首先计算微分方程的解析解，如果发现没有解析解，再用数值解

一、Matlab中求微分方程的数值解函数

[x,y]=solver('f',ts,x0,options)

1）x代表自变量

2）y代表函数值

3）solver代表求解函数，常用的为ode45->函数图像变化较为平稳，ode15s->函数图像存在突变（一个运行时间长时，换另一个）

4）ts=[t0,t1]，代表自变量函数的初值和终值，也可以指定以一个向量

5）x0表示函数的初始值

6）options（可选）：options=odeset('reltol','rt','abstol',at)，rt，at分别为设定的相对误差，默认为10e-3，10e-6

7）f表示由待解微分方程编写的m函数文件名

二、一阶微分方程的数值解

1、求解步骤：

1）化为标准型（等号左边为微分）如下图

2）带入函数求解

2、例子

例1（求方程）：

题目： $y-y{^'}=2x$ ,初始值为y(0)=3,x范围为[0,2]

解：

1）化为标准型： $y{^'}=y-2x$

2）matlab求解如下：

2.1）新建一个m文件，用来编辑f函数

function dy = f(x,y) % dy是标准化后左边的式子，x是式子右边的自变量，y是因变量 dy = y - 2*x; end

2.2）新建一个m文件用来计算结果

[x y] = ode45('f',[0,2],3); %x是自变量，y是因变量 figure(1) plot(x,y,'*-') axis([0 2 -inf inf]) grid on

3）运行结果

例2：（求方程组）

题目：

$\left\{\begin{matrix} y{_1^'}=y{_2}y{_3}\\ y{_2^'}=-y{_1}y{_3}\\ y{_3^'}=-0.51y{_1}y{_2} \end{matrix}\right.$ ，自变量 $x\in[0.4\pi]$ ，初始值 $y{_1}(0)=0,y{_2}(0)=y{_3}(0)=1$

解：

1）化为标准型（题目已经为标准型）

2）matlab代码

2.1）新建一个m文件，用来编辑f函数

function dy = f(x,y) % dy是标准化后左边的式子，x是式子右边的自变量，y是因变量 dy = zeros(3,1); dy(1) = y(2)*y(3); dy(2) = -y(1)*y(3); dy(3) = -0.51*y(1)*y(2); end

2.2）新建一个m文件用来计算结果

[x, y] = ode45('f', [0 4*pi], [0 1 1]); % 这里的y是一个有3列的矩阵哦！x是自变量，y是因变量 plot(x, y(:,1), 'o', x, y(:,2), '*', x, y(:,3), '+') legend('y1','y2','y3') % 加上标注 axis([0, 4*pi, -inf, +inf]) % 设置横坐标范围为0-4pi，纵坐标范围不需要设置，写成-inf到+inf

3）结果

三、高阶微分方程的数值解

1、步骤

1）标准化

高阶方程转化为一阶方程组合

2）按一阶求解

2、例子

1）标准化

2）matlab求解

【本文地址】

公司简介

联系我们

今日新闻

点击排行

实验室常用的仪器、试剂和: 说到实验室常用到的东西，主要就分为仪器、试剂和耗

不用再找了，全球10大实验: 01、赛默飞世尔科技（热电）Thermo Fisher Scientif

三代水柜的量产巅峰T-72坦: 作者：寞寒最近，西边闹腾挺大，本来小寞以为忙完这

通风柜跟实验室通风系统有: 说到通风柜跟实验室通风，不少人都纠结二者到底是不

集消毒杀菌、烘干收纳为一: 厨房是家里细菌较多的地方，潮湿的环境、没有完全密

实验室设备之全钢实验台如: 全钢实验台是实验室家具中较为重要的家具之一，很多

图片新闻

实验室药品柜的特性有哪些: 实验室药品柜是实验室家具的重要组成部分之一，主要

小学科学实验中有哪些教学: 计算机计算器一般打孔器打气筒仪器车显微镜

实验室各种仪器原理动图讲: 1.紫外分光光谱UV分析原理：吸收紫外光能量，引起分

高中化学常见仪器及实验装: 1、可加热仪器：2、计量仪器：（1）仪器A的名称：量

微生物操作主要设备和器具: 今天盘点一下微生物操作主要设备和器具，别嫌我啰嗦

浅谈通风柜使用基本常识: 　众所周知，通风柜功能中最主要的就是排气功能。在

Matlab求微分方程的数值解

Matlab求微分方程的数值解

今日新闻

点击排行

推荐新闻

图片新闻

专题文章