同伦问题与同伦算法

River Chandler

已于 2023-08-16 21:22:42 修改

阅读量418

点赞数

分类专栏： # 动力系统专题文章标签：数学建模 python numpy

于 2023-08-16 00:56:24 首次发布

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/Chandler_river/article/details/132310027

版权

动力系统专题专栏收录该内容

5 篇文章 0 订阅

订阅专栏

同伦问题
- 据我所知，这篇博客是CSDN上少数几篇讲同伦算法的博客之一
考虑
同伦算法的目的
- 扩大初值选取范围
- 解决非线性代数方程组的全部解计算问题

同伦算法中的基本概念

考虑 $F:D\subset R^n\rightarrow R^n,$ 求 $F(x)=0$ 的解 $x^*\in D$ 人为地引入参数t,构造一个函数族 $H(x,t)$ 使得

$H(x,0)=F_0(x),H(x,1)=F(x)$

同时假设 $F_0(x)=0$ 的解已知，从 $t=0$ 出发可以求解 $H(x,t)=0$
对于，假设为
- 如果 $x(t)$ 可以形成一条 $R^n$ 中的光滑曲线，其奇点 $x(0)$ 为 $F_0(x)=0$ 的解，据假设它是已知的，曲线的终点 $x(1)$ 正是我们要求的 $x^*\in D$
- $H(x,t)$ 称为一个同伦其解为同伦曲线
同伦的构造，例如
- $H(x,t)=tF(X)+(1-t)F_0(x)$
同伦曲线的存在与光滑
- 假设存在
- 假设其性质好

同伦算法得到的等价问题

$H(x(t),t) \equiv 0,\forall t\in [0,1]$

$\left\{\begin{matrix} H_x(x(t),t)x'(t)+H_t(x(t),t)=0,\forall t\in[0,1]\\ x'(t)=-[H_x(x(t),t)]^{-1}H_t(x,t),t\in(0,1)\\ x(0)=x_0 \end{matrix}\right.$

代数方程组与同伦方程

假定多项式方程组 $P(x)=0$ 其分量形式

$\begin{matrix} p_1(x_1,x_2,...,x_n)=0\\ ...\\ p_n(x_1,x_2,...,x_n)=0 \end{matrix}$

假定 $Q(x)=0$ 易得到

$Q(x)=(q_1(x),q_2(x),...,q_n(x))^T$

对应的同伦方程

$0=H(x,t)=(1-t)Q(x)+tP(x)$

同伦方程的性质

平凡性
- $Q(x)=0$ 的解已知
光滑性
- 每一条同伦曲线都是t的单值函数
可达性
- $P(x)=0$ 的任意一个孤立解都有从 $Q(x)=0$ 出发的同伦方程的解曲线达到
小结
- 跟踪所有的解曲线即可得到 $P(x)=0$ 的全部解

Q(x) 的一种取法

$\begin{matrix} q_1(x)=a_1x_1^{d_1}-b_1\\ ...\\ q_n(x)=a_nx_n^{d_n}-b_n \end{matrix},\,\,\,d_i=max\sum_{i=1}^n\alpha_i$

$\alpha_i$ 是多项式 $p_i$ 某一单项式包含x_i 的幂次（一个p_i 一个d_i）

Bezout 定理

Bezout 数
- Bezout 数定义了全部的解曲线
Bezout 定理
- P(x) 孤立解的个数小于等于 Bezout 数

退化

非线性代数方程组的解的个数小于Bezout 数则为退化的
- 绝大多数实际问题的解都是远远小于Bezout 数的

同伦算法的计算机实现

略
私信交流

关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
打赏
0
评论
同伦问题与同伦算法

同伦问题
复制链接

扫一扫

专栏目录

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

打赏作者

River Chandler 谢谢,我会更努力学习工作的!!

¥1 ¥2 ¥4 ¥6 ¥10 ¥20

扫码支付：¥1

获取中

扫码支付

您的余额不足，请更换扫码支付或充值

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。