FZERO error with matrix dimensions

Question

Girija Bahety 2021-4-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/802541-fzero-error-with-matrix-dimensions

编辑： Matt J 2021-4-19

I am trying to evaluate a function with fzero and getting the following error message - I am not sure why its throwing an error:

Error using fzero (line 246)

FZERO cannot continue because user-supplied function_handle ==> xsld failed with the error below.

Incorrect dimensions for raising a matrix to a power. Check that the matrix is square and the power is a scalar. To perform elementwise matrix

powers, use '.^'.

Error in moments (line 80)

[w,fv]= fzero(@xsld,w0,options);

The codes are below:

tol=10^-5; %tolerance 
w0=[0.1,200]; %wage interval
options=optimset('TolX',tol,'Display','iter');
options=optimset(options,'PlotFcns',{@optimplotx,@optimplotfval});
try    
[w,fv]= fzero(@xsld,w0,options);
    catch  
    cr = 50;
    return
end

And the function @xsld has the following codes, saved as a function file (I am only showing the codes that have the argument w):

function[xld] = xsld(w)
some codes...
RW = repmat(w.*(s.^gamma),1,ns); %wage income 
A = (1-((1-nu)*(1-alpha))).*(((1-nu)*(1-alpha))/w)^(((1-nu)*(1-alpha))/(1-((1-nu)*(1-alpha))));
lu = (((theta_r.^nu).*(1-nu).*(1-alpha).*(k.^(alpha*(1-nu))))/w).^(1/(nu + alpha*(1-nu)));
Ld = sum(mean(lu,2)) + ((B*psi)/w)^(1/(1-psi));
xld = (sum(mean((((theta_r.^nu).*((1-nu)*(1-alpha)).*(k.^(alpha*(1-nu))))/w).^(1/(nu + alpha*(1-nu))),2)) + ((B*psi)/w)^(1/(1-psi))) - Ls; 