Matlab学习牛顿迭代法求解方程
本篇在课程学习中编程实现的迭代法计算非线性方程或者超越方程近似根的算法。
开始接触Matlab
第一次使用 Matlab ,虽然以前一直有听说Matlab,直到现在才发现,Matlab确实十分方便,因为刚开始使用这个,遇到了一些困难,于是便用博客的记录下来,方便后续复习。
变量创建
在MATLAB中提供了sym和syms函数用来定义符号变量,这两种定义符号变量的方式是不同的
>> x=sym('x');
>> y=sym('y') %此时就成功创建了xy变量
>> sym x;
>> sym y
ans =
y
%就只创建了最后那一个变量,所以使用sym要注意。
使用syms创建多个变量会比sym方便许多。
>> syms x y
%创建x,y变量。
多项式创建
构造带数值多项式最简单的方法就是直接带入向量。使用poly2sym来实现,调用格式为poly2sym(p) 其中p 是多项式的系数向量。
>> p=[1 2 3 4 5]%创建系数向量。
p =
1 2 3 4 5
>> poly2sym(p)
%创建多项式;
%也可以使用poly2sytm(p,y)创建以y为变量的函数,如果不指定就默认使用x作为变量;
%使用poly2sym(p,y)前提是要使用sym或者syms创建y的变量,否则会报错。
ans =
x^4 + 2*x^3 + 3*x^2 + 4*x + 5
多项式的导数
可以使用函数diff来实现,其调用格式为diff(fx),其中fx是待求的函数。
>> p=[1 2 3 4 5]
p =
1 2 3 4 5
>> fx=poly2sym(p)%生成待求函数
fx =
x^4 + 2*x^3 + 3*x^2 + 4*x + 5
>> fx0=diff(fx)
fx0 =
4*x^3 + 6*x^2 + 6*x + 4
>>
多项式的导数也可以使用函数polyder来实现,其调用格式为polyder(p) 其中p是多项式的系数向量。
>> p=[1 2 3 4 5]
p =
1 2 3 4 5
>> poly2sym(p)
ans =
x^4 + 2*x^3 + 3*x^2 + 4*x + 5
>> polyder(p)%得到一阶导数系数向量。
ans =
4 6 6 4
>> poly2sym(ans)%使用poly2sym创建一阶导数多项式。
ans =
4*x^3 + 6*x^2 + 6*x + 4
牛顿迭代法
首先创建一个m文件,可以在命令行中输入edit,或者在编辑器中点击新建脚本;
再使用cd命令更改当前文件夹,例如:
>> cd 'D:\matlab_pro\project'
牛顿迭代法是使用切线去靠近待求的函数的解,是不断的使用切线去靠近待求的点。
先选定一个初始点,然后求这个点的切线与x轴的交点,再用这个交点做垂线与曲线的交点作下一个切线,再连续迭代就会越来越靠近函数的解。
、
将求得的值与之前的值做差取绝对值,可以通过范数来实现即调用norm函数
clear
syms a
%定义待求方程 x^2-1
x=0:0.1:5
y=x.^2-1
b=solve(a^2-1,a)%使用solve函数直接求解该方程,与迭代法得到的值做比较
x0=input(" 输入估计值:") %输入初始的预估的值
for i=1:5
y0=x0.^2-1
k=2*x0
x0=x0-y0/k
error(i,1)=norm(x0-b) %存储每一次迭代之后的误差的绝对值
end
vpa(x0,3)
vpa(b,3)
plot(error)%打印error曲线