how to use a piecewise function with fsolve?

Question

Ekin Ceyda Cetin 2017-1-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/322262-how-to-use-a-piecewise-function-with-fsolve

评论： Walter Roberson 2018-11-10

I want to use piecewise objective function with fsolve. If i define the objective function using an if statement I get 'Conversion to logical from sym is not possible.'. When I use piecewise built-in function it returns an error when I try convert the objective function to a matlabFunction.

Actually, the objective function consists of 3 functions, for a region i need fsolve to consider only two of the functions. I was thinking of doing this by defining one of the functions as NaN in a certain region. Is this possible?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2017-1-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/322262-how-to-use-a-piecewise-function-with-fsolve#answer_252341

It is not possible to use piecewise with fsolve. It is possible to use piecewise with vpasolve.

Also I discovered earlier today that if you have a symbolic expression that contains piecewise, that if you use matlabFunction() with the 'file' option then the code written will use 'if' to decide the case and do the appropriate computation. Unfortunately the generated code is not vectorized (which would require logical indexing rather than if), but provided you used the matlabFunction 'vars' option to put all of variables into a single argument, you would be able to use it with fsolve.

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

Ekin Ceyda Cetin 2017-1-28

在 MATLAB Online 中打开

Thanks for your answer.

Do I need to do logical indexing if I use 'file' and 'vars' options?

However I don't get why it returns the error about if statement. 'Illegal use of reserved keyword "if".'

syms param1 param2 param3;
param=[param1,param2,param3];
n=numel(peakdata);
p1=piecewise(param(1)~=0 & param(1)<1, param(2)+param(3)/param(1)*(gamma(1-param(1))-1),...
            param(1)==0,param(2)+param(3)*(-psi(1)),...
            param(1)>1,NaN);
p2=piecewise(param(1)~=0 & param(1)<0.5, (param(3)^2)*((gamma(1-2*param(1)))-(gamma(1-2*param(1)))^2)/(param(1))^2,...
            param(1)==0, (param(3)^2)*(pi^2)/6,...
            param(1)>0.5,NaN);
p3=piecewise(param(1)==0, 12*(6^0.5)*zeta(3)/pi^3,...
            sign(param(1))*(gamma(1-3*param(1))-3*gamma(1-param(1))*gamma(1-2*param(1))+2*(gamma(1-param(1)))^3)/((gamma(1-2*param(1))-(gamma(1-param(1)))^2)^(3/2)));
mean=(1/n)*sum(peakdata(1:end));
f1=piecewise(param(1)<1,p1-mean,param(1)>1,NaN);
a2=0;
for l=1:n
a2=a2+(peakdata(l)-((1/n)*sum(peakdata(1:end))))^2;
end
variance=(1/(n-1))*a2;
f2=piecewise(param(1)<0.5,p2-variance,param(1)>0.5,NaN);
a3=0;   
for l=1:n
    a3=a3+((peakdata(l)-((1/n)*sum(peakdata(1:end))))/(((1/(n-1))*a2)^0.5))^3;
end
skewness=1/n*a3;
f3=p3-skewness;
f=matlabFunction([f1;f2;f3],'vars',{param},'File','mom_gev');
param0=[0.2,2,1];
options=optimoptions('fsolve','MaxFunctionEvaluations',500,'Display','iter-detailed')
[xfinal,fval,exitflag,output] = fsolve(f,param0,options);

Cagdas Ozgenc 2018-11-10

Has this bug been fixed? Thanks.

Walter Roberson 2018-11-10

在 MATLAB Online 中打开

Unfortunately no, if you use matlabFunction with 'file' and piecewise() then although it generates appropriate if to handle the piecewise, it still does so assuming the input is a scalar. For example

F = [a, a^2, piecewise(a<0, -(-a)^(1/3), a^(1/3))]

is currently still turned into code that has

if (a < 0.0)
    t0 = -(-a).^(1.0./3.0);
else
    t0 = a.^(1.0./3.0);
end
F = [a,a.^2,t0];

instead of generating logical indexing or a loop.

请先登录，再进行评论。

how to use a piecewise function with fsolve?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

how to use a piecewise function with fsolve?

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

WeChat

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

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