Matlab:解非线性方程组

news/2024/10/17 22:27:13/

1、基于问题求解非线性方程组

例:

exp(-exp(-(x_{1}+x_{2})))=x_{2}(1+x_{1}^{2})

x_{1}cosx_{2}+x_{2}sinx_{1}=\frac{1}{2}

x=optimvar('x',2);            %将x定义为一个二元素优化变量
eq1=exp(-exp(-(x(1)+x(2))))==x(2)*(1+x(1)^2); %创建第一个方程作为优化等式表达式
eq2=x(1)*cos(x(2))+x(2)*sin(x(1))==1/2;       %创建第二个方程作为优化等式表达式
prob=eqnproblem;              %创建一个方程问题
prob.Equations.eq1=eq1;
prob.Equations.eq2=eq2;       %将方程放入问题中
show(prob)                    %检查问题
x0.x=[0 0];                   %从[0,0]点开始求解问题
[sol,fval,exitflag]=solve(prob,x0)
disp(sol.x)                   %查看解点

2、基于问题求解多项式非线性方程组

例:

已知x是2x2的矩阵,x^{3}=\begin{bmatrix} 1 & 2\\ 3& 4 \end{bmatrix}求解

x=optimvar('x',2,2);                %将变量x定义为一个2×2矩阵变量
eqn=x^3 == [1,2;3,4];               %用x定义要求解的方程
prob=eqnproblem('Equations',eqn);   %用方程创建一个问题
x0.x=ones(2);                       %从[1,1;1,1]点开始求解问题
sol=solve(prob,x0)
disp(sol.x)                         %检查解
sol.x^3                             %显示解的立方

3、fsolve解非线性方程组

语法:

x = fsolve(fun,x0)    %从x0 开始,尝试求解方程fun(x)=0(全零数组)
x = fsolve(fun,x0,options)
x = fsolve(problem)
[x,fval] = fsolve(___)
[x,fval,exitflag,output] = fsolve(___)
[x,fval,exitflag,output,jacobian] = fsolve(___)

例:求解下面两个非线性方程组

e^{-e^{-(x_{1}+x_{2})}}=x_{2}(1+x_{1}^{2})

x_{1}cosx_{2}+x_{2}sinx_{1}=\frac{1}{2}

(1)先将其转化成F(x)=0的形式

e^{-e^{-(x_{1}+x_{2})}}-x_{2}(1+x_{1}^{2})=0

x_{1}cosx_{2}+x_{2}sinx_{1}-\frac{1}{2}=0

(2)创建root2d.m

function F = root2d(x)F(1) = exp(-exp(-(x(1)+x(2)))) - x(2)*(1+x(1)^2);
F(2) = x(1)*cos(x(2)) + x(2)*sin(x(1)) - 0.5;

(3)

type root2d.m
fun=@root2d;
x0=[0,0];
x=fsolve(fun,x0)

4、迭代输出,展示方程组的求解过程

例:2x_{1}-x_{2}=e^{-x_{1}}

-x_{1}+2x_{2}=e^{-x_{2}}

(1)先将其转化成F(x)=0的形式

2x_{1}-x_{2}-e^{-x_{1}}=0

-x_{1}+2x_{2}-e^{-x_{2}}=0

(2)

F=@(x)[2*x(1)-x(2)-exp(-x(1));-x(1)+2*x(2)-exp(-x(2))]; %编写一个函数用来算F
x0=[-5;-5];                                             %创建初始点x0[-5;-5]
options=optimoptions('fsolve','Display','iter');        %返回迭代输出
[x,fval]=fsolve(F,x0,options)                           %求解

上接非线性https://blog.csdn.net/weixin_73011353/article/details/135141433icon-default.png?t=N7T8https://blog.csdn.net/weixin_73011353/article/details/135141433


http://www.ppmy.cn/news/1283222.html

相关文章

基于采样的自动驾驶规划算法 - PRM,RRT,RRT*,CL-RRT

本文将讲解PRM,RRT,RRT*自动驾驶规划算法原理,不正之处望读者指正 0 前言 机器人运动规划的基本任务:从开始位置到目标位置的运动 (1)如何躲避构型空间出现的障碍物 (2)如何满足机器…

【QT】qt中多线程的使用

前景: C#中的多线程举例 private Thread threadDownLoadParameter;threadDownLoadParameter new Thread(new ThreadStart(parameterInterface.downParameterLoadData)); //parameterInterface 实例对象,downParameterLoadData 方法threadDownLoadParame…

Linux阶段性项目——FTP服务器

一、FTP服务概述 ​FTP为File Transfer Protocol的缩写,即文件传输协议,是TCP/IP 协议族中的协议之一。FTP是一个用于在计算机网络上在客户端和服务器之间进行文件传输的应用层协议。RFC 959定义了FTP协议规范。是互联网上使用的最广泛的文件传输协议&a…

不浪费时间,昂首资本1分钟如何快速学习MT4价差

不要浪费时间在手工计算上,昂首资本解释一下如何快速学习MT4价差,。 想要在MT4中输入交易时,需要在交易窗口中设置未来交易的参数。在同一个窗口中,可以看到卖价和买价。如果在上面的例子中比较这两个价格,就会发现两…

JVM基础篇---02

为什么需要用户自定义类加载器: 扩展类加载器的功能: Java的默认类加载器主要有三个,分别是引导类加载器、扩展类加载器和应用程序类加载器。其中,引导类加载器和扩展类加载器是由JVM实现的,用户无法修改其行为。而应用…

Ensp DHCP 接口地址池(配置命令)

难易度:* 拓补图: 实验环境:一台路由器,两台交换机,四台pc机 需求: 左边两台PC 自动分配 IP为1.0网段 右边两台PC 自动分配 2.0 网段, 路由器配置命令 左边接口 分配1.0网段地址 [Huawei]…

【三维生成】稀疏重建、Image-to-3D方法(汇总)

系列文章目录 总结一下近5年的三维生成算法,持续更新 文章目录 系列文章目录一、LRM:单图像的大模型重建(2023)摘要1.前言2.Method3.实验 二、SSDNeRF:单阶段Diffusion NeRF的三维生成和重建(ICCV 2023&am…

排样问题 强化学习组合优化

目录 排样问题 上海友图: 2d不规则图像排样 java代码 有一定效果: 2023年像是python 矩形框 算法讲解 无代码 强化学习加组合优化 Max Rects Packer 最大矩形装箱算法 排样问题 2.排样问题(nesting problem)又称为下料问题(cutting and stock prob…