gradient-and-newton.zip_newton

% function [J,k1,k2,ee,time] = gradient(d,a,t) function [J,k1] = gradient(d,a,t) syms x y z; a1=a; e = 0.0015; k1 = 0; J=[]; t1=1; while(t1 > t) J1 = 0; J2 = []; k1=k1+1; for i = 1:20 J2 = (x*d(1,i)+y*d(2,i)+z*d(3,i)-1)^2; J1 = J1+J2; end J3 = [diff(J1,x); diff(J1,y); diff(J1,z)]; J4 = subs(J1,{x,y,z},{a1(1),a1(2),a1(3)}); J5 = subs(J3,{x,y,z},{a1(1),a1(2),a1(3)}); J = [J,J4]; a1 = a1 - e*J5; t1 = norm(J5)*e; if(k1>30) break; end end % k2 = []; % ee=0.001:0.0002:0.005; % % time=zeros(1,length(ee)); % for i=1:length(ee) % a0=a; % time1=clock; % e1 = ee(i); % k = 0; % t2=1; % while(t2>t) % k = k + 1; % J1 = 0; % J2 = []; % for j = 1:20 % J2 = (x*d(1,j)+y*d(2,j)+z*d(3,j)-1)^2; % J1 = J1+J2; % end % J3 = [diff(J1,x); % diff(J1,y); % diff(J1,z)]; % % J5 = subs(J3,{x,y,z},{a0(1),a0(2),a0(3)}); % % a0 = a0 - e1*J5; % t2 = norm(J5)*e1; % if(k>30) % break; % end % end % k2 = [k2,k]; % time(i) = etime(clock,time1); % end