bvpinit
得出边界值问题求解器的初始估计值
说明
为边界值问题中具有未知值的参数指定初始估计值向量。您可以将此语法用于上述任一输入参量组合。solinit
= bvpinit(___,parameters
)
示例
创建 BVP 的解的初始估计值,用 bvp4c
求解 BVP,然后将解扩展到新域。
得出 BVP 问题的解的良好初始估计值可能是求解问题最困难的部分。BVP 解不一定唯一,因此初始估计值可以是求解器确定要返回的解的决定因素。初始估计值应该满足边界条件,其间的行为应该反映您对此问题的一般期望(解是否振荡,是否为简单的线性函数等)。
假设有以下微分方程:
.
该方程受限于边界条件
.
将方程编码为一阶方程组的函数是
function dydx = bvpfun(x,y) dydx = [y(2) -y(1)]; end
同样,对边界条件进行编码的函数是
function res = bcfun(ya,yb) res = [ya(1) yb(1)]; end
您可以将所需的函数作为局部函数包含在文件末尾(如本处所示),也可以将它们作为单独的命名文件保存在 MATLAB® 路径上的目录中。
通过函数句柄提供初始估计值
您可以合理地预期方程的解是振荡的,因此使用正弦和余弦函数可以很好地得出解及其固定边界点之间导数的初始估计值。
function y = guess(x) y = [sin(x) cos(x)]; end
使用域 中的 10 个等距网格点和初始估计值函数创建一个解结构体。
xmesh = linspace(0,pi,10); solinit = bvpinit(xmesh,@guess);
求解 BVP
使用 ode 函数、边界条件和解估计值调用 bvp4c
。绘制结果。
sol = bvp4c(@bvpfun, @bcfun, solinit);
plot(sol.x,sol.y,'-o')
局部函数
此处列出了 BVP 求解器 bvp4c
为计算解而调用的局部辅助函数。您也可以将这些函数作为它们自己的文件保存在 MATLAB 路径上的目录中。
function dydx = bvpfun(x,y) % equation being solved dydx = [y(2) -y(1)]; end %------------------------------------------- function res = bcfun(ya,yb) % boundary conditions res = [ya(1) yb(1)]; end %------------------------------------------- function y = guess(x) % guess at solution behavior y = [sin(x) cos(x)]; end %-------------------------------------------
在初始区间上求解 BVP,然后使用每个解作为下一区间的初始估计值来迭代扩展区间。
假设使用以下方程
.
如果要表示为一阶方程组,则包含以下两个方程:
,
.
该方程最初在区间 上定义,并受限于边界条件
,
.
将方程编码为一阶方程组的函数是
function dydx = bvpfun(x,y) dydx = [y(2) y(1)]; end
同样,对边界条件进行编码的函数是
function res = bcfun(ya,yb) res = [ya(1) yb(1)-1]; end
您可以将所需的函数作为局部函数包含在文件末尾(如本处所示),也可以将它们作为单独的命名文件保存在 MATLAB 路径上的目录中。
初始估计值
使用指数函数作为解的初始估计值。由于方程有两个解分量,因此需要编写一个 y = guess(x)
形式的初始估计值函数,该函数返回向量。
function y = guess(x) y = [exp(x) exp(x)]; end
包含五个点的网格足以捕获估计值函数的行为。
xmesh = linspace(0,3,5); solinit = bvpinit(xmesh,@guess);
求解方程
在初始区间 中求解方程,并绘制 的结果。
sol = bvp4c(@bvpfun, @bcfun, solinit);
plot(sol.x(1,:),sol.y(1,:),'-o')
扩展区间
现在,使用 bvpinit
来循环扩展积分区间,求解每个新问题并绘图。在每次迭代中,使用外插至新区间 [0 k]
的上一个解 sol
得出初始估计值。在每个新问题中,bvp4c
会在新边界 [0 k]
上强制应用边界条件。
hold on for k = 4:8 solinit = bvpinit(sol,[0 k]); sol = bvp4c(@bvpfun, @bcfun, solinit); plot(sol.x(1,:),sol.y(1,:),'-o') end
此示例显示了连续法的一个简化版本,这是一种很有用的方法,它通过将问题分解成更小的区间或更简单的问题来求解 BVP。有关此方法的更多示例,请参阅:
局部函数
此处列出了 BVP 求解器 bvp4c
为计算解而调用的局部辅助函数。您也可以将这些函数作为它们自己的文件保存在 MATLAB 路径上的目录中。
function dydx = bvpfun(x,y) % equation being solved dydx = [y(2) y(1)]; end %------------------------------------------- function res = bcfun(ya,yb) % boundary conditions res = [ya(1) yb(1)-1]; end %------------------------------------------- function y = guess(x) % guess at solution behavior y = [exp(x) exp(x)]; end %-------------------------------------------
输入参数
初始网格,指定为向量。要在区间 [a,b] 上求解问题,请将 x(1)
指定为 a 并将 x(end)
指定为 b。x
的项必须为升序(如果 a < b)或降序(如果 a > b)。求解器会根据解调整此网格(通过添加、删除和移动网格点),因此诸如 x = linspace(a,b,10)
之类的估计值通常足以满足需求。为了处理比较难的问题,请在解快速变化的位置放置一些网格点。
对于两点边界值问题,
x
的项必须唯一。即,如果 a < b,项必须满足x(1)
<x(2)
< ... <x(end)
。如果 a > b,项必须满足x(1)
>x(2)
> ... >x(end)
。对于多点边界值问题,通过重复
x
中的项,可以在端点 a 和 b 之外,在 [a,b] 中指定应用边界条件的点。例如,假设有以下向量对于此网格,在以下三个点应用边界条件:端点x = [0 0.5 1 1 1.5 2];
0
和2
,以及重复项1
。一般说来,重复项表示 [a,b] 中区域之间的边界点。重复项1
将区间[0 2]
分为两个区域:[0 1]
和[1 2]
。
示例: solinit = bvpinit(linspace(a,b,10),yinit)
数据类型: single
| double
| int8
| int16
| int32
| int64
| uint8
| uint16
| uint32
| uint64
| logical
复数支持: 是
解的初始估计值,指定为向量或函数。
向量 - 对于解的每个分量,
bvpinit
在所有网格点上复制向量的对应元素作为常量估计值。也即,yinit(i)
是x
中所有网格点处解的第i
个分量yinit(i,:)
的常量估计值。函数 - 对于给定网格点而言,估计值函数必须返回一个向量,该向量的元素是解的对应分量的估计值。此函数必须采用以下格式
y = guess(x)
x
是网格点,y
是其长度与解中的分量数相同的向量。例如,如果yinit
是函数,则在每个网格点上,bvpinit
都会调用y(:,j) = guess(x(j))
对于多点边界值问题,估计值函数必须采用以下格式
y = guess(x,k)
y
是区域k
中x
处的解的初始估计值。该函数必须接受输入参量k
,提供此输入参量的目的是为了灵活编写估计值函数。但是,该函数不必使用k
。
示例: solinit = bvpinit(x,[sqrt(3)/2; 0])
示例: solinit = bvpinit(x,@guess)
数据类型: single
| double
| int8
| int16
| int32
| int64
| uint8
| uint16
| uint32
| uint64
| logical
| char
| function_handle
复数支持: 是
未知参数值的初始估计值,指定为标量或向量。
示例: solinit = bvpinit(x, yinit, [0 1 sqrt(2)])
指定三个未知参数的估计值向量。
示例: 键入 edit mat4bvp
可查看具有一个未知参数的 BVP 的示例。
数据类型: single
| double
| int8
| int16
| int32
| int64
| uint8
| uint16
| uint32
| uint64
| logical
复数支持: 是
先前的解,指定为由 bvp4c
或 bvp5c
返回的解结构体。如果 sol
包含参数,这些参数将被复制到 solinit
。
数据类型: struct
扩展功能
此函数完全支持基于线程的环境。有关详细信息,请参阅在基于线程的环境中运行 MATLAB 函数。
版本历史记录
在 R2006a 之前推出
MATLAB Command
You clicked a link that corresponds to this MATLAB command:
Run the command by entering it in the MATLAB Command Window. Web browsers do not support MATLAB commands.
选择网站
选择网站以获取翻译的可用内容,以及查看当地活动和优惠。根据您的位置,我们建议您选择:。
您也可以从以下列表中选择网站:
如何获得最佳网站性能
选择中国网站(中文或英文)以获得最佳网站性能。其他 MathWorks 国家/地区网站并未针对您所在位置的访问进行优化。
美洲
- América Latina (Español)
- Canada (English)
- United States (English)
欧洲
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)