Октавный индекс за пределами

Question

Программа представляет собой оптимизированный спуск градиента.

Вот код:

clear all 
    close all 
    [x,y] = meshgrid(-2:0.1:2);            
    z = x.^2 + 100*y.^2;           
    n = 1; 
    k(n)=0.01; 
    arret = 0.0001;               
    mesh(x,y,z);                
    [x1(n),y1(n)] = ginput(1);            
    diff_x(n) = 2*x1(n);             
    diff_y(n) = 200*y1(n);             
    while sqrt(diff_x(n)^2 + diff_y(n)^2) > arret && n < 30 
    k(n) = sqrt(diff_x(n)^2 + diff_y(n)^2)/(8*x1(n).^2+2*10.^6*y1(n).^2);   
    x1(n+1) = x1(n) - k(n)*diff_x(n);          
    y1(n+1) = y1(n) - k(n)*diff_y(n);  
    n = n+1; 
    diff_x(n) = 2*x1(n); 
    diff_y(n) = 200*y1(n); 
    z1(n) = x1(n).^2 + 100*y1(n).^2;          
    plot3(x1(n),y1(n),z1(n));            
    end 

    x1(n) 
    y1(n) 
    n 

Так что я получил это, и я не понимаю, почему.

error: GradientPasOptFinal2: A(I): index out of bounds; value 2 out of bound 1

error: called from:

error: error: C:\Octave\octave-3.8.2\GradientPasOptFinal2.m at line 13, column 8

решаемые: к (п) между y1 (п) и п была причина, я не знаю, почему, но теперь программа работает, спасибо!

Answer 1

Вы определили только k (1). k (2) не определено. Переместить строку:

k(n) = sqrt(diff_x(n)^2 + diff_y(n)^2)/(8*x1(n).^2+2*10.^6*y1(n).^2); 

Впереди этих линий:

x1(n+1) = x1(n) - k(n)*diff_x(n);          
y1(n+1) = y1(n) - k(n)*diff_y(n); 

И он должен работать.