, ,
. , Z, .. . Z . , , . :
.
:
r .
.
, .. , , Z , ( ), , r . , , , , , Z . , . , , , , . r , , .
: , .
, ∈ , , , . k=1 .
. k- .
. :
,
- , ; - .
:
1)
2)
.
, , .
, . .
. (k+1)- .
1) R=1/g(i);
, 2.
clc
clear
|
|
[x11 x22]=meshgrid([-15:0.1:15],[-15:0.1:15]);
syms x1 x2;
eps=10^(-5);
f=7*x1^2+2*x1*x2+5*x2^2-2*x1-10*x2;
z=7.*x11.^2+2.*x11.*x22+5.*x22.^2-2.*x11-10.*x22;
xx0=[5 10];
m=10000;
beta=0.5;
g1=x1+x2-6;
n=1;
k=0;
g=[g1];
plot(xx0(1),xx0(2),'o');
hold on
while(1)
plot(xx0(1),xx0(2),'.');
hold on
a=0;
for i=1:n
s=subs(g(i),[x1 x2],xx0);
if(abs(s)>=2*eps)
a=a+1/g(i);
end
end
if (a==0)
RES=xx0;
break;
end
F=f+m*a;
RES=method_NS(xx0,F,z,eps,g,n);
s=1/subs(g(i),[x1 x2],RES);
if (abs(m*s)<eps)
break;
end
m=m*beta;
k=k+1;
xx0=RES;
end
xx0
k
ff=subs(f,[x1 x2],xx0)
plot(RES(1),RES(2),'o');
:
Eps=10^(-5)
F=x1^2+x2^2
X0=[5; 10]
X = [2; 2]
F(x)=8
m | x | F(x) | Time (seconds) | k |
2.000004060440634 2.000007075540867 | 8.000044543992551 | 3.575780 | ||
2.000004717017259 2.000004719098306 | 8.000037744506781 | 5.27514 | ||
2.000004752794619 2.000004752790486 | 8.000038022385597 | 6.380273 | ||
2.000006127663761 2.000006127663761 | 8.000049021385184 | 9.462643 | ||
2.000006433557238 2.000006433557238 | 8.000051468540686 | 13.098821 |