二阶边值问题的数值解matlab,二阶线性微分方程边值问题的MATLAB求解
2012 年 3 月第 26 卷 第 1 期 阴 山 学 刊 YINSHAN ACADEMIC JOURNAL Mar. 2012 Vo1. 26 No. 1 二阶线性微分方程边值问题的 MATLAB 求解* 云 文 在 ( 包头师范学院 数学科学学院,内蒙古 包头 014030) 摘 要: 本文给出二阶线性微分方程边值问题数值算法的 MATLAB 实现,并举例进行了求解仿真及与解 析解的精度比较。 关键词: 边值问题; 初值问题; MATLAB 中图分类号: O175 文献标识码: A 文章编号:1004 -1869( 2012) 01 -0023 -02 微分方程数值解中,侧重研究初值问题,即已知 x0 对其他时刻状态变量值进行求解的方法。在实际问题中,经常会遇到这样的问题: 已知部分状态在t = 0时刻的值,还知道部分状态在时刻t0 = tf 的值,这类问题即所谓边值问题。而在 MATLAB 语言中边值问题也是 ode45( ) 类函数无法直接求解的一类问题。本文采用将边值问题转化为初值问题的方法,给出二阶线性微分方程的边值问题的计算机求解。 1 边值问题的数学描述 二阶线性微分方程的边值问题的数学描述: ¨ y( x) + p( x) y( x) + q( x) y( x) = f( x) ( 1) 其中 p( x) 、q( x) 和 f( x) 均为给定函数。假设在区间[a,b]上研究该方程的解,且已知在这两个边界点 上满足边界条件 y( a) = γa,y( b) = γb ( 2) 2 数值求解方法 由于不能直接获得在初始时刻的各个变量的值,因此求解初值问题的通常算法在解边值问题时是不能直接使用的。边值问题数值解法的基本思想是找出能够满足式( 2) 边值的相应初值 y( 0) 和 y( 0) ,然后再利用初值问题算法来求解这一初值问题。该算法也称为打靶算法( shooting method) 。 算法步骤: 求出下面方程初值问题的数值解 y1( b) ¨ y1( x) + p( x) y1( x) + q( x) y1( x) = 0, y1( a) = 1, y1( a) = 0 求出下面方程初值问题的数值解 y2( b) ¨ y2( x) + p( x) y2( x) + q( x) y2( x) = 0, y2( a) = 0, y2( a) = 1 求出下面方程初值问题的数值解 yp( b) ¨ yp( x) + p( x) yp( x) + q( x) yp( x) = 0, yp( a) = 0, yp( a) = 1 若 y2( b) ≠ 0,则计算 m = γb - γay1( b) - yp( b) y2( b) 求出下面方程初值问题的数值解,则 y( x) 即为原边值问题的数值解 ¨ y( x) + p( x) y( x) + q( x) y( x) = f( x) , y( a) = γa, y( a) = m 3 上面算法的 MATLAB 实现 求解时应该首先得出对应的一阶微分方程组模型,即设x1 = y,x2 = y,则得出式( 1) 对应的方程组为: x1 = x2 x2 = - q( x) x1 - p( x) x2 + f( x{ ) 则上面算法的 MATLAB 实现为 Function [t,y] = shooting( f1,f2,tspan,x0f,varargin) 32 * 收稿日期:2011 -11 -04 基金项目: 内蒙古自治区自然科学基金项
二阶边值问题的数值解matlab,二阶线性微分方程边值问题的MATLAB求解相关推荐
- 迭代函数系统matlab实现,线性随机IFS迭代Matlab示例程序.doc
线性随机IFS迭代Matlab示例程序 线性随机IFS迭代 Matlab示例程序 Tree Leaf Mountain %% 线性随机IFS迭代[iterated function systems迭代 ...
- matlab 一元线性微分,怎样用matlab求一元线性函数极值
clear; %%%%% 清变量 clc; %%%%% 清屏幕 syms x; %%%% 定变量 f = sqrt(5^2+(8-x)^2)*28.8+sqrt(15^2+(5+x)^2)*7.2; ...
- MATLAB求线性代数的参数范围,MATLAB科学计算04(线性代数问题求解一)
文章目录 特殊矩阵 矩阵的基本概念 求解线性方程组 直接求解 判定求解 特殊矩阵 零矩阵.1矩阵及单位矩阵 生成nxn方阵: A=zeros(n), B=ones(n), C=eye(n) 生成mxn ...
- matlab二阶非线性微分方程求解,二阶非线性常微分方程的打靶法.ppt
二阶非线性常微分方程的打靶法 二阶非线性常微分方程的打靶法 计算思路 主要分为以下五步: 给定容许误差ε,迭代初始值γ1,对k=1,2,...做: (1)用四阶Runge-Kutta 方法求解初值问, ...
- matlab二阶阻尼震荡衰减,二阶欠阻尼电路的零响应输入-MATLAB课程设计.doc
二阶欠阻尼电路的零响应输入-MATLAB课程设计.doc 摘要 MATLAB具有强大的运算和作图功能,为电路分析引入计算机数值方法提供了技术保证.MATLAB中的变量与常量都是矩阵,其元素可以使复数和 ...
- 用matlab求微分方程系数,用Matlab软件求解微分方程的解析解和数值解.pdf
用Matlab软件求解微分方程的解析解和数值解.pdf Matlab软件求解微分方程 的解析解和数值解 数学与信息科学学院 孔祥庆 数学建模实验项目2 (1) 一.实验名称: Matlab软件求解微分 ...
- 2021-01-07 matlab数值分析 常微分方程初边值问题数值解 标准龙格库塔四阶四段公式 欧拉法
matlab数值分析 常微分方程初边值问题数值解 标准龙格库塔四阶四段公式 欧拉法 1.标准龙格库塔四阶四段公式 function y=rk4(fun,a,b,y0,n) h=(b-a)/n; y(1 ...
- 智能软开关 配电网重构matlab 二阶锥 以33节点为研究对象,编制配电网故障重构模型
智能软开关 配电网重构matlab 二阶锥 编程方法:matlab+yalmip(cplex为求解器) 基本内容:以33节点为研究对象,编制配电网故障重构模型,采用图论知识保证配电网的连通性和辐射性, ...
- 二阶锁频环辅助三阶锁相环 原理及matlab实现
接收接收信号时,需要我们对信号进行捕获,跟踪,通过载波对信号进行跟踪. 跟踪:通过卫星信号的载波频率与码相位 跟踪环路:由载波跟载波环路,与码跟踪环路两部分构成,分别用来跟踪信号的中的杂波与伪码. 载 ...
最新文章
- 关于某日访问次数最多的IP的topK问题的三种解法
- codevs1032
- 雅虎书面问题——最大的问题汇总
- Python语言学习之字母D开头函数使用集锦:del、dict使用方法之详细攻略
- 如何使用oracle查询,oracle 表查询
- Linux命令之 -- export 设置/显示系统环境变量
- C++学习——string
- 初识大数据(四. 大数据与人工智能的关系)
- Virtualbox 2.1突发性错误解决办法(也许是BUG)
- aws s3 獲取所有文件_Url从Amazon S3获取文件
- 构造方法(设计一个Fan类来表示一个风扇)
- fastreport 打印有 demo version字样_用户增长留存的数字化管理实操, 11大行业demo重构你的运营方法论...
- 前端小白-HTML简介
- windows关闭445、3389端口
- 孤岛惊魂5 for Android,孤岛惊魂5手机版
- 并行调用框架封装实战-completionService
- [leetcode] 4 寻找两个有序数组的中位数(二分+递归查找第K小数)(重要)
- 概率论考点总结类型26 统计中的卡方分布 t分布F分布
- android中Service使用startService
- Android集成极光推送踩坑(二)升级篇