参考北京师范大学的《计算物理基础》
第四章:求方程的零点(方程求根)
计算物理基础_中国大学MOOC(慕课)www.icourse163.org1.求方程零点的算法
- 二分法
- 切线法(比二分法快)
- 弦割法
- 迭代法
2.求方程零点的指令
(1)fzero:求单变量实函数
的零点(通过函数句柄调用函数)
x=fzero(fun,x0)
x=fzero(fun,x0,options,p1,p2,...)
[x,fval]=fzero(...)
x,fval:零点及其对应的函数值,没有零点则返回NaN
fun:求零点的单变量实函数
x0:搜寻零点的出发点或区间[x0(1),x0(2)]
options:控制算法的优化选项,[ ]是默认值
两种方法:
(2)roots(c) 求多项式的全部零点
将多项式系数按降幂排列,如求多项式
的零点。
>> roots([5,3,-4,1,7,9])
ans =
0.9530 + 0.7622i
0.9530 - 0.7622i
-1.3127 + 0.0000i
-0.5966 + 0.7515i
-0.5966 - 0.7515i
(3)fsolve求多变量非线性方程组的复数根
x=fsolve(fun,x0)
[x,fval]=fsolve(fun,x0,...)
x,fval:零点及其对应的函数值,没有零点则返回NaN
fun:求解的非线性方程组
x0:猜测的初始解
2020.12.5