如何使用MATLAB解非线性方程组

2021-11-03 科技 236阅读

用solve函数。

例如:

x^2+y^3=10

x^3-y^2=1

其中x,y为方程组的未知量

在Matlab的命名窗口中输入:

syms x y

[x y]=solve('x^2+y=10','x^2-y^2=1','x','y')
输出计算结果为:
x =

(37^(1/2)/2 + 21/2)^(1/2)

(21/2 - 37^(1/2)/2)^(1/2)

-(21/2 - 1/2*37^(1/2))^(1/2)

-(1/2*37^(1/2) + 21/2)^(1/2)

y =

- 37^(1/2)/2 - 1/2

37^(1/2)/2 - 1/2

37^(1/2)/2 - 1/2

- 37^(1/2)/2 - 1/2

扩展资料:

注意事项

solve 是基本的用于符号解方程的内置函数,返回类型为符号变量矩阵(m×nm×n sym)。当无法符号求解时,抛出警告并输出一个数值解。

基本形式为:solve(eqn, var, Name, Val);

% eqn为符号表达式/符号变量/符号表达式的函数句柄, var为未知量; Name为附加要求,Val为其值可以用solve解一维方程。对于多项式,solve可以返回其所有值。

func1 = @(x)x^3 - 20*x^2 - 25*x + 500;    % 创建函数句柄.句柄中的变量不属符号变量,不需要定义。

syms x exp1;    % 定义符号变量 x, exp1;

exp1 = x^3 - 20*x^2 - 25*x + 500;    % 符号表达式,包含符号变量. 符号变量必须先有上一行定义。

solve(exp1 == 0, x)    % 命令行输入a,传入一个包含符号表达式的等式,x为所要求的变量

solve(exp1, x)    % 命令行输入b,传入一个符号表达式,函数默认求其零点

solve(func1(x), x)    % 命令行输入c,传入参数func1(x)等价于传入了符号表达式,和输入b完全一样

solve(func1(x) == 0, x)    % 命令行输入d,这句话和a完全一样

solve(func1, x)    % 命令行输入e,传入参数func1,这是一个函数句柄,函数默认求其零

ans =    % 命令行输出,三个解,为3*1的符号向量。对以上五种输入输出都完全一样

-5

5

20

对于不可符号求解的函数零点/方程解,solve抛出警告并返回一个数值解:

exp1 = atan(x) - x - 1;    % 不可符号求零点的表达式

solve(exp1 == 0, x)    % 命令行输入

% 命令行输出:

警告: Cannot solve symbolically. Returning a numeric approximation instead.

ans =

-2.132267725272885131625420696936


声明:你问我答网所有作品(图文、音视频)均由用户自行上传分享,仅供网友学习交流。若您的权利被侵害,请联系fangmu6661024@163.com