Finding a root with interval constraint

4 次查看（过去 30 天）

显示更早的评论

Chien-Chia Huang 2014-1-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/113633-finding-a-root-with-interval-constraint

评论： Walter Roberson 2014-1-25

采纳的回答： Walter Roberson

在 MATLAB Online 中打开

Hello there!

I am trying to find a point x within the time interval [t-1,t] (for some t, say t = 3) so that the function attains value zero. That is, I want to solve "Q_0 + integral(a+b*sin(c*t+d)-mu,t-1,x) = 0" for x in [t-1,t]. My code is the following

    y = fsolve(@(x) Q_0+(a-mu)*(x-t+1)-(b/c)*cos(c*x+d)+(b/c)*cos(c*(t-1)+d),0,optimset('Display','off'))

wherein (a,b,c,d) satisfy a + b*sin(c*t+d), and Q_0 and mu are constants. This code has no problem. However, the solution may sometimes be outside the time interval [t-1,t], which is not what I want.

So, my question is if there is a way to restrict the routine to find a solution that lies within [t-1,t] exactly?

Thanks!

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

请先登录，再进行评论。

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

采纳的回答

Walter Roberson 2014-1-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/113633-finding-a-root-with-interval-constraint#answer_122114

As your x0 is a scalar (0), your x are scalar, and that implies you can use fzero() instead of fsolve(). With fzero() you can pass the interval [t-1 t] as your x0.

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

Chien-Chia Huang 2014-1-25

编辑：Chien-Chia Huang 2014-1-25

在 MATLAB Online 中打开

Thanks, Walter. My code now goes like this (value of mu changes)

[FirstVanish,~,exitflag] = fzero(@(x) queuelength(counter_qln)+(a(j)-mu)*(x-t+1)-(b(j)/c(j))*cos(c(j)*x+d(j))+(b(j)/c(j))*cos(c(j)*(t-1)+d(j)),[t-1,t])

However, it showed the error msg

Error using fzero (line 274) The function values at the interval endpoints must differ in sign.

This is why I will need to know how to get things going without seeing the above.