Size of Q is undefined

2 次查看（过去 30 天）
feeroz babu
feeroz babu 2021-11-23
评论： feeroz babu 2021-11-23
function [j] = armijo_grad_direction_dnm_f1(P1,VG,sigma)
j = 0;
signal = 0;
nn = size(Q);
n = nn(1);
I = eye(n);
while( signal == 0 && j < MAX_j)
Q = P1^(1/2)*expm(P1^(-1/2)*2^(-j)*VG*P1^(-1/2))*P1^(1/2); % Geodesic with direction 2^(-j)*VG starting at P1.
PHIQ = 1/2*trace(I-2*Q^(-1)+(Q^(-1))^(2)); % Merit Function at Q.
PHIP1 = 1/2*trace(I-2*P1^(-1)+(P1^(-1))^(2));
ARM = PHIQ - PHIP1 - 2^(-j)*sigma*trace(VG*P1^(-1)*VG*P^(-1)); %Armijo Gradient Direction
if(ARM > 0)
j = j+ 1;
else
signal = 1;
end
end
if ( j >= MAX_j)
fprintf('\nArmijo: j not found with %d iterations.\n',MAX_j);
return;
end
end
  2 个评论
KSSV
KSSV 2021-11-23
How did you run the function? What inputs you have used?
feeroz babu
feeroz babu 2021-11-23
I am running this in octave then I got
error: 'P1' undefined near line 4, column 4
error: called from
armijo_newton_direction_dnm_f1 at line 4 column 4

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

KSSV
KSSV 2021-11-23
You need to define the input variables and then call the function. It seems you are striaght away hitting the f5/ run button so you are getting error.
P1 = define your value ; % enter value
VG = define your value ; % enter input
sigma = define your value ;
% now call the function
j = armijo_grad_direction_dnm_f1(P1,VG,sigma)
  7 个评论
显示 5更早的评论
KSSV
KSSV 2021-11-23
Your function name is _armijo_grad_direction_dnm_f1 and you are trying to call armijo_newton_direction_dnm_f2 it seems,
feeroz babu
feeroz babu 2021-11-23
Yes I have two functions. Ok, done. I selected wrong directory. Thank you for your kind help. My code runs.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Chebyshev 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by