I am getting this preallocation warning

Vânia Gonçalves

2022 11 19

2 个回答

2 次查看（30 天）

0 个投票

5 个评论
显示 3更早的评论隐藏 3更早的评论

Vânia Gonçalves 2022-11-19

在 MATLAB Online 中打开

clear; clc; close all
% Function Definition (Enter your Function here):
syms X Y;
f = 100*(Y-(X^2))^2 + (1-X)^2;
% Initial Guess (Choose Initial Guesses):
x(1) = -1.1;
y(1) = 0.5;
e = 10^(-3); % Convergence Criteria
i = 1; % Iteration Counter
% Gradient and Hessian Computation:
df_dx = diff(f, X);
df_dy = diff(f, Y);
J = [subs(df_dx,[X,Y], [x(1),y(1)]) subs(df_dy, [X,Y], [x(1),y(1)])]; % Gradient
ddf_ddx = diff(df_dx,X);
ddf_ddy = diff(df_dy,Y);
ddf_dxdy = diff(df_dx,Y);
ddf_ddx_1 = subs(ddf_ddx, [X,Y], [x(1),y(1)]);
ddf_ddy_1 = subs(ddf_ddy, [X,Y], [x(1),y(1)]);
ddf_dxdy_1 = subs(ddf_dxdy, [X,Y], [x(1),y(1)]);
H = [ddf_ddx_1, ddf_dxdy_1; ddf_dxdy_1, ddf_ddy_1]; % Hessian
S = inv(H); % Search Direction
% Optimization Condition:
while norm(J) > e
    I = [x(i),y(i)]';
    lam = J*J'./(J*H*J');
    x(i+1) = I(1)-lam*S(1,:)*J';
    y(i+1) = I(2)-lam*S(2,:)*J';
    i = i+1;
    J = [subs(df_dx,[X,Y], [x(i),y(i)]) subs(df_dy, [X,Y], [x(i),y(i)])]; % Updated Jacobian
    ddf_ddx_1 = subs(ddf_ddx, [X,Y], [x(i),y(i)]);
    ddf_ddy_1 = subs(ddf_ddy, [X,Y], [x(i),y(i)]);
    ddf_dxdy_1 = subs(ddf_dxdy, [X,Y], [x(i),y(i)]);
    H = [ddf_ddx_1, ddf_dxdy_1; ddf_dxdy_1, ddf_ddy_1]; % Updated Hessian
    S = inv(H); % New Search Direction
end

Jan 2022-11-19

@Vânia Gonçalves: A general hint is to search in the net at first. Asking an internet search engine for "Matlab preallocate" shows useful explanations already.

Vânia Gonçalves 2022-11-19

i already searched in the internet, but i didn´t find any usefull explanation, so i hoped find help here.

but thanks anyawy

请先登录，再进行评论。

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

Follow Question

回答（2 个）

Steven Lord 2022-11-19

1 个投票

This is a Code Analyzer warning not a run-time warning.

If you're not sure what preallocation is, why it can be important, or how to do it see the first example on this documentation page.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Vânia Gonçalves 2022-11-19

Thank you very much !!

请先登录，再进行评论。

Jan 2022-11-19

移动：Jan 2022-11-19

在 MATLAB Online 中打开

0 个投票

Seriously? If I ask Google, I get e.g.:

For your code:

...
maxIter = 1000;
x = zeros(maxIter, 1);
y = zeros(maxIter, 1);
while norm(J) > e && i < maxIter
    ...
end
x = x(1:i);  % Crop unused elements
y = y(1:i);