Рабочая шпора
.doc1. Метод Ньютона (касательных)
Метод Хорд
Метод хорд и касательных
f(x)*f’’(x)>0 – хорды
f(x)*f’’(x)<0 – касательные
Метод половинного деления
a b f(a) f(b) c f(c)
Метод простой итерации для уравнения
Xn+1=x=ф(xn)
Если || ф’(x)||>1 то xn+1=xn-JI*f(xn)
JI=1/max{f’(a),f’(b)}
2. Метод сеток
i\j |
|
|
0 |
1 |
2 |
3 |
4 |
5 |
6 |
|
|
y\x |
|
a |
a+hx |
2 |
2,5 |
3 |
3,5 |
4 |
|
|
|
C |
|
|
|
|
|
|
|
D |
4 |
1 |
|
(Uac+Ucd)/2 |
Ucd |
Ucd |
Ucd |
Ucd |
Ucd |
2,125 |
|
3 |
0,75 |
|
Uac |
u |
1,993 |
2,622 |
3,250 |
3,878 |
Ubd |
|
2 |
0,5 |
|
Uac |
v+(Ucd-v)* (i-1)/(N-1)=u |
(u+(Ubd-u)*j)/(M-1) |
(u+(Ubd-u)*j)/(M-1) |
(u+(Ubd-u)*j)/(M-1 |
(u+(Ubd-u)*j)/(M-1 |
Ubd |
|
1 |
c+hy |
|
Uac |
Uac+(Ubd- Uac)*j/M=v |
Uac+(Ubd- Uac)*j/M |
Uac+(Ubd- Uac)*j/M |
Uac+(Ubd- Uac)*j/M |
Uac+(Ubd- Uac)*j/M |
Ubd |
|
0 |
c |
|
-0,5 |
Uab |
Uab |
Uab |
Uab |
Uab |
0,625 |
|
|
|
A |
|
|
|
|
|
|
|
B |
Альфа=hy2/(2*(hx2+hy2))
Бета=hx2/(2*(hx2+hy2))
1 |
2 |
3 |
4 |
5 |
7 |
8 |
9 |
10 |
6 |
11 |
12 |
13 |
14 |
15 |
1-6: Альфа*(верх+низ)+Бета*(лево+право)
7-15: Альфа*(верх_этой+низ)+Бета*(лево_этой+право)
3. Вычисление интегралов
J(h)=h*(f(x0)/2+CУММА f(x)+f(xk))
Метод Эйлера
xn+1=xn+h, yn+1=yn+h*f(x,y)
Модифицированный метод Эйлера
xn+1=xn+h, yn+1=yn+h*f1, f1=f(x1,y1)
x1= xi+h/2, y1=y+h/2*f(xi,yi)
Усовершенствованный метод Эйлера
xn+1=xn+h, yn+1=yn+h/2*f1, f1= f(xi,yi)+f(x1,y1)
x1= xi+h, y1=yn+h*f(xi,yi)
Метод простой итерации для нелинейных систем
x0,y0 – точка пересечения на графике
Ф(x,y)
Метод Ньютона для нелинейных систем
x0,y0 – точка пересечения на графике
f(x),f(x0),F(x),F(x0),определитель,F-1=1/опр*F(x0)
xn+1=xn-F-1(xn)*f(xn)
Метод Ньютона для линейных систем
Сделать нули по диагонали в матрице
x0=g, xn+1=x1*xn1+…+xk*xnk
Метод Зейделя
xn+1=x1*xn1+x2*xn+1…+xk*xn+1 k
Метод скорейшего спуска
f(x),f(x0),W(x),W(x0),WT(x0), WT(x0)*f(x0), W(x0)* WT(x0)*f(x0)
M0=(f(x0);W(x0)* WT(x0)*f(x0))/(W(x0)* WT(x0)*f(x0);W(x0)* WT(x0)*f(x0))
xn+1=xn-Mn(WT(xn),f(xn))
Метод Данилевского
M3j=x4j/x43
xiji+1=xiji+M3j/xi3
Метод наименьших квадратов
y = bx + a
b=(n*S(xiyi)-S(xi)*S(yi))/(n*S(xi2)-S2(xi))
a=(S(yi)-b*S(xi))/n
Метод Рунге-Кутта
(i) |
(x) |
y |
R=f(x,y) |
Δy |
1 |
xi |
yi |
R1=f(xi,yi)*H |
R1 |
2 |
xi+H/2 |
yi+R1/2 |
R2=f(x(i),y(i))*H |
2*R2 |
3 |
xi+H/2 |
yi+R2/2 |
R3=f(x(i),y(i))*H |
2*R3 |
4 |
xi+H |
yi+R3 |
R4=f(x(i),y(i))*H |
R4 |
|
|
Δyi=1/6*(R1+R2+R3+R4) |
||
i+1 |
yi+Δyi |
|
|
|
Метод прогонки
yi+1=miyi+niyi-1=f^ih2; mi=-(2-qih2)/(1+pih/2); ni=(1-pih/2)/(1+pih/2)
f^i=fi/(1+pih/2); c1=1/m1; d1=f^1h2-n1A; ci=1/(mi-nici-1);
di=f^ih2-nici-1di-1