牛顿迭代法 (应用:西北农林科技大学120 求解方程)
牛顿迭代法
编辑本段产生背景
牛顿迭代法(Newton's method)又称为牛顿-拉夫逊方法(Newton-Raphson method),它是牛顿在17世纪提出的一种在实数域和复数域上*似求解方程的方法。多数方程不存在求根公式,因此求精确根非常困难,甚至不可能,从而寻找方程的*似根就显得特别重要。方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程f(x) = 0的单根附*具有*方收敛,而且该法还可以用来求方程的重根、复根,此时线性收敛,但是可通过一些方法变成超线性收敛。另外该方法广泛用于计算机编程中。编辑本段牛顿迭代公式
设r是f(x) = 0的根,选取x0作为r初始*似值,过点(x0,f(x0))做曲线y = f(x)的切线L,L的方程为y = f(x0)+f'(x0)(x-x0),求出L与x轴交点的横坐标 x1 = x0-f(x0)/f'(x0),称x1为r的一次*似值。过点(x1,f(x1))做曲线y = f(x)的切线,并求该切线与x轴交点的横坐标 x2 = x1-f(x1)/f'(x1),称x2为r的二次*似值。重复以上过程,得r的*似值序列,其中x(n+1)=x(n)-f(x(n))/f'(x(n)),称为r的n+1次*似值,上式称为牛顿迭代公式。View Code
1 编辑本段C语言代码
2
3 double func(double x) //函数
4
5 {
6
7 return x*x*x*x-3*x*x*x+1.5*x*x-4.0;
8
9 }
10
11 double func1(double x) //导函数
12
13 {
14
15 return 4*x*x*x-9*x*x+3*x;
16
17 }
18
19 int Newton(double *x,double precision,int maxcyc) //迭代次数
20
21 {
22
23 double x1,x0;
24
25 int k;
26
27 x0=*x;
28
29 for(k=0;k<maxcyc;k++)
30
31 {
32
33 if(func1(x0)==0.0)//若通过初值,函数返回值为0
34
35 {
36
37 printf("迭代过程中导数为0!\n");
38
39 return 0;
40
41 }
42
43 x1=x0-func(x0)/func1(x0);//进行牛顿迭代计算
44
45 if(fabs(x1-x0)<precision || fabs(func(x1))<precision) //达到结束条件
46
47 {
48
49 *x=x1; //返回结果
50
51 return 1;
52
53 }
54
55 else //未达到结束条件
56
57 x0=x1; //准备下一次迭代
58
59 }
60
61 printf("迭代次数超过预期!\n"); //迭代次数达到,仍没有达到精度
62
63 return 0;
64
65 }
66
67 int main()
68
69 {
70
71 double x,precision;
72
73 int maxcyc;
74
75 printf("输入初始迭代值x0:");
76
77 scanf("%lf",&x);
78
79 printf("输入最大迭代次数:");
80
81 scanf("%d",&maxcyc);
82
83 printf("迭代要求的精度:");
84
85 scanf("%lf",&precision);
86
87 if(Newton(&x,precision,maxcyc)==1) //若函数返回值为1
88
89 printf("该值附*的根为:%lf\n",x);
90
91 else //若函数返回值为0
92
93 printf("迭代失败!\n");
94
95 getch();
96
97 return 0;
98
99 }
100
101 编辑本段C++代码
102
103 #include<iostream>
104
105
106 #include<cmath>
107
108 using namespace std;
109
110 int main()
111
112 {
113
114 double diedai(int a,int b,int c,int d,double x);
115
116 int a,b,c,d;
117
118 double x=10000.0;
119
120 cout<<"请依次输入方程四个系数:";
121
122 cin>>a>>b>>c>>d;
123
124 x=diedai(a,b,c,d,x);
125
126 cout<<x;
127
128 return 0;
129
130 }
131
132 double diedai(int a,int b,int c,int d,double x)
133
134 {
135
136 x=x-(a*pow(x,3.0)+b*pow(x,2.0)+c*pow(x,1.0)+d)/(3*a*pow(x,2.0)+2*b*pow(x,1.0)+c);
137
138 if(abs(a*x*x*x+b*x*x+c*x+d)<=0.000001)
139
140 return x;
141
142 else
143
144 return diedai(a,b,c,d,x);
145
146 }
147
148 编辑本段matlab代码
149
150 1.定义函数
151
152 function y=f(x)
153
154 y=f(x);%函数f(x)的表达式
155
156 function y=z(x)
157
158 y=z(x);%函数z(x)的表达式
159 2.主程序
160
161 x=X;%迭代初值
162
163 i=0;%迭代次数计算
164
165 while i<= I;%迭代次数
166
167 y=x-y(x)/z(x);%牛顿迭代格式
168
169 if abs(y-x)>ε;%收敛判断
170
171 x=y;
172
173 else break
174
175 end
176
177 i=i+1;
178
179 end
180
181 fprintf('\n%s%.4f\t%s%d','x=',x,'i=',i) %输出结果