matlab解非线性方程组?matlab语言常用算法程序集
matlab解非线性方程组?
用solve函数。例如:x^2+y^3=10x^3-y^2=1其中x,y为方程组的未知量在Matlab的命名窗口中输进:syms x y(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/237^(1/2)/2 - 1/237^(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的符号向量。对以上五种输进输出都完全一样-5520对于不可符号求解的函数零点/方程解,solve抛出警告并返回一个数值解:exp1 = atan(x) - x - 1; % 不可符号求零点的表达式solve(exp1 == 0, x) % 命令行输进% 命令行输出:警告: Cannot solve symbolically. Returning a numeric approximation instead.ans =-2.132267725272885131625420696936
非线性方程组可以直接求解可以用fsolve函数求解,数值解法可以用牛顿法迭代求解。[1]中有很多程序例子可以直接调用里面的函数。对于计算机求解方程组,微分方程组的求解构思可以了解些算法的知识[2]。参考:[1] 龚纯, 王正林. (2008). MATLAB 语言常用算法程序集. 电子工业出版社. [2] 李信真,车刚明, 欧阳洁,等.(2010).计算 *** .西北工业大学出版社.