同伦问题与同伦算法

  • 同伦问题
    • 据我所知,这篇博客是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
  • 对于t\in [0,1],假设H(x,t)=0x(t)
    • 如果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 数 B=\Pi_{i=1}^nd_i
    • Bezout 数定义了全部的解曲线
  • Bezout 定理
    • P(x) 孤立解的个数小于等于 Bezout 数

退化

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

同伦算法的计算机实现

  • 私信交流

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/94493.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

在项目中如何解除idea和Git的绑定

在项目中如何解除idea和Git的绑定 1、点击File--->Settings...(CtrlAltS)--->Version Control--->Directory Mappings--->点击取消Git的注册根路径: 2、回到idea界面就没有Git了: 3、给这个项目初始化 这样就可以重新绑定远程仓库了&#x…

部署piwigo网页 通过cpolar分享本地电脑上的图片

通过cpolar分享本地电脑上有趣的照片:发布piwigo网页 文章目录 通过cpolar分享本地电脑上有趣的照片:发布piwigo网页前言1. 设定一条内网穿透数据隧道2. 与piwigo网站绑定3. 在创建隧道界面填写关键信息4. 隧道创建完成 总结 前言 首先在本地电脑上部署…

13、Redis与Mysql数据双写一致性

Redis与Mysql数据双写一致性 1、 复习+面试题 先动Mysql,再动Redis,两害相衡取其轻,避免Redis业务key突然消失,多线程请求集火打满mysql 动,写操作,先更新数据库,在删除缓存。回写的…

实验篇——亚细胞定位

实验篇——亚细胞定位 文章目录 前言一、亚细胞定位的在线网站1. UniProt2. WoLFPSORT3. BUSCA4. TargetP-2.0 二、代码实现1. 基于UniProt(不会)2. 基于WoLFPSORT后续(已完善,有关代码放置于[python爬虫学习(一&#…

【大数据Hive】hive 事务表使用详解

目录 一、前言 二、Hive事务背景知识 hive事务实现原理 hive事务原理之 —— delta文件夹命名格式 _orc_acid_version 说明 bucket_00000 合并器(Compactor) 二、Hive事务使用限制 参数设置 客户端参数设置 客户端参数设置 三、Hive事务使用操作演示 操作步骤 客…

Windows系统Git安装教程(详细Git安装过程)

获取Git安装程序 到Git官网下载,网站地址:https://git-scm.com/downloads,如下图: 因为我们是用Windows系统上的浏览器访问的,Git官网自动之别到了我使用的操作系统,所以右侧直接显示下载使用Windows系统的…

gromacs教程练习1

gromacs能在win上运行,还是个开源的软件,这都很值得入手学习 记录下gromacs教程的练习情况: Lysozyme in water 水中的溶菌酶,嗯,估计就是把蛋白处理后放在显试溶剂里跑MD这个模拟。 1、文件的准备: 1、…

【仿写tomcat】四、解析http请求信息,响应给前端,HttpServletRequest、HttpServletResponse的简单实现

思考 在解析请求之前我们要思考一个问题,我们解析的是其中的哪些内容? 对于最基本的实现,当然是请求类型,请求的url以及请求参数,我们可以根据请求的类型作出对应的处理,通过url在我们的mapstore中找到se…

【数据结构OJ题】合并两个有序链表

原题链接:https://leetcode.cn/problems/merge-two-sorted-lists/description/ 目录 1. 题目描述 2. 思路分析 3. 代码实现 1. 题目描述 2. 思路分析 可以先创建一个空链表,然后依次从两个有序链表中选取最小的进行尾插操作。(有点类似双…

基于docker搭建pytest自动化测试环境(docker+pytest+jenkins+allure)

pytest搭建自动化测试环境(dockerpytestjenkinsallure) 这里我以ubuntu18为例 如果有docker环境,可以直接拉取我打包好的镜像docker pull ziyigun/jenkins:v1.0 1 搭建Docker 1.1 安装docker # 配置docker安装环境 sudo apt-get install ap…

Android Studio实现解析HTML获取图片URL将图片保存到本地

目录 效果activity_main.xmlMainActivityImageItemImageAdapter 效果 项目本来是要做成图片保存到手机然后读取数据后瀑布流展示&#xff0c;但是有问题&#xff0c;目前只能做到保存到手机 activity_main.xml <?xml version"1.0" encoding"utf-8"?…

C++写文件,直接写入结构体

C写文件&#xff0c;直接写入结构体 以前写文件都是写入字符串或者二进制再或者就是一些配置文件&#xff0c;今天介绍一下直接写入结构体&#xff0c;可以在软件参数较多的时候直接进行读写&#xff0c;直接将整个结构体写入和读取&#xff0c;看代码&#xff1a; #include&…

文本图片怎么转Excel?分享一些好用的方法

在处理数据时&#xff0c;Excel 是一个非常强大的工具&#xff0c;但有时候需要将文本和图片转换为 Excel 格式&#xff0c;这可能会让人感到困惑。在本文中&#xff0c;我们将介绍一些好用的方法&#xff0c;以便您能够轻松地将文本和图片转换成 Excel 格式。 将文本图片为Exc…

vue3 videojs实现播放器,动态更改src

一、背景 vue3下载第三方插件videojs&#xff0c;达到播放器的效果&#xff0c;并且点击事件能够动态更改播放器的src。实现思路&#xff1a; 场景一&#xff1a;只有一个播放器&#xff0c;当点击事件&#xff0c;直接赋值&#xff0c;动态更改封装好的组件的src参数&#xff…

线程|线程的使用、四种实现方式

1.线程的实现方式 1.用户级线程 开销小&#xff0c;用户空间就可以创建多个。缺点是&#xff1a;内核无法感知用户级多个线程的存在&#xff0c;把其当作只有一个线程&#xff0c;所以只会提供一个处理器。 2.内核级线程 相对于用户级开销稍微大一点&#xff0c;可以利用多…

无涯教程-Perl - setgrent函数

描述 此功能将枚举设置(或重置)到组条目集的开头。该函数应在第一次调用getgrent之前调用。 语法 以下是此函数的简单语法- setgrent返回值 此函数不返回任何值。 例 以下是显示其基本用法的示例代码- #!/usr/bin/perl -wwhile( ($name,$passwd,$gid,$members)getgrent…

c语言每日一练(9)

前言&#xff1a;每日一练系列&#xff0c;每一期都包含5道选择题&#xff0c;2道编程题&#xff0c;博主会尽可能详细地进行讲解&#xff0c;令初学者也能听的清晰。每日一练系列会持续更新&#xff0c;暑假时三天之内必有一更&#xff0c;到了开学之后&#xff0c;将看学业情…

QT的设计器介绍

设计器介绍 Qt制作 UI 界面&#xff0c;一般可以通过UI制作工具QtDesigner和纯代码编写两种方式来实现。纯代码实现暂时在这里不阐述了在后续布局章节详细说明&#xff0c;QtDesigner已经继承到开发环境中&#xff0c;在工程中直接双击ui文件就可以直接在QtDesigner设计器中打…

unity发布WebGL遇到的坑(持续更新)

1、unity默认字体在网页中不会显示 解决方法&#xff1a;自己新导入一个字体&#xff0c;使用导入的字体 2、之前打过包并运行过&#xff0c;后面又在unity中进行了修改&#xff0c;重新打包&#xff0c;运行发现还是修改之前的效果&#xff0c;虽然是新包&#xff0c; 解决方…

Windows上使用dump文件调试

dump文件 dump文件记录当前程序运行某一时刻的信息&#xff0c;包括内存&#xff0c;线程&#xff0c;线程栈&#xff0c;变量等等&#xff0c;相当于调试程序时运行到某个断点上&#xff0c;把程序运行的信息记录下来。可以通过Windbg打开dump&#xff0c;查看程序运行的变量…