^사용하여 오류 하나의 인수는 정사각형 행렬이어야하고 다른 하나는 스칼라 여야합니다. POWER를 사용합니다. (^) elementwise 전력 내가이 함수에서 최소화 함수를 작성하기 위해 노력하고있어

Question

,

(4-2.1*x1^2+x1^4/3)*x1^2+x1*x2+(-4+4*x2^2)*x2^2 

, -10 < = X1 < = 10 -10 < = X2 < = 10 .

내가 쓰고있는이게 맞습니까? 그에서

min = -10; 
max = 10; 
x1 = min+max*rand(); 
x2 = min+max*rand(); 
T = 1.0; 
Tmin = 0.0000001; 
refresh = 0.9; 
finalMin = (4-2.1*x1^2+x1^4/3)*x1^2+x1*x2+(-4+4*x2^2)*x2^2; 
finalCoordinate=[x1,x2]; 
while(T>Tmin) 
    i=0; 
    while(i<=300) 
     x1= min+max*rand(1,1); 
     x2= min+max*rand(1,1); 
     newMin = (4-2.1*x1^2+x1^4/3)*x1^2+x1*x2+(-4+4*x2^2)*x2^2; 
     delta1 = finalMin-newMin; 
     if(newMin < finalMin) 
      finalMin=newMin; 
      finalCoordinate=[x1,x2]; 
     elseif(2.71828^(delta/T))>(0+(rand(1)*1)) 
      finalCoordinate=[x1,x2]; 
     end; 
     i=i+1; 
    end; 
    T = T*refresh; 
end; 

그러나 왜,

elseif(2.71828^(delta/T))>(0+(rand(1)*1)) 




Error using^One argument must be a square matrix and the other must be a scalar. Use POWER (.^) for elementwise power. 

어떻게에서 얻기 위해 100 %의 정확도는 좌표?

Answer 1

글쎄, 오류 메시지에 따르면 delta은 스칼라가 아닐 수도 있습니다. 여기에 코드에서 정의하지 않았습니다.

fun = @(x) (4-2.1*x(1)^2+x(1)^4/3)*x(1)^2+x(1)*x(2)+(-4+4*x(2)^2)*x(2)^2; 
x0 = rand(2,1)*20-10; 
[x,fval] = fmincon(fun, x0, [],[], [],[], [-10;-10], [10;10])