How to Input a function equation in matlab

12 次查看（过去 30 天）
Ivy Shen
Ivy Shen 2018-10-3
编辑： Stephen23 2018-10-3
Hi, I have an equation as following:
F - hm*(a-b)*log(1+F/(hm*(a-b))) = kt
I am trying to use fzero to solve for F in matlab, but I don't know how to input this equation in matlab first?
  4 个评论
显示 2更早的评论
Ivy Shen
Ivy Shen 2018-10-3
I am trying to write a function as following:
function F = GreenAmp(Ks,hm,t,a,b)
F - hm*(a-b)*log(1 + F/(hm*(a-b))) = Ks*t;
And I don't know what's the next step. Is there any easier way to do it? Can you specify which documentation I should look at?
Thank you.
Stephen23
Stephen23 2018-10-3
编辑：Stephen23 2018-10-3
function F = GreenAmp(Ks,hm,t,a,b)
F - hm*(a-b)*log(1 + F/(hm*(a-b))) = Ks*t;
That will not work:
  • variable F is used but undefined. In order to solve for F, F must be one of the input arguments to the function.
  • the LHS of the = sign is not valid for assignment: the LHS should be one variable, which is the value of the function when it is evaluated.
"Can you specify which documentation I should look at?"
https://www.mathworks.com/help/matlab/getting-started-with-matlab.html
https://www.mathworks.com/help/matlab/ref/function.html
https://www.mathworks.com/help/matlab/math/parameterizing-functions.html
https://www.mathworks.com/help/matlab/ref/fzero.html

请先登录，再进行评论。

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

回答（1 个）

Stephen23
Stephen23 2018-10-3
编辑：Stephen23 2018-10-3
>> b = 1;
>> a = 2;
>> hm = 0.5;
>> kt = 0.9;
>> fun = @(F) F - hm*(a-b)*log(1+F/(hm*(a-b)));
>> y = fzero(@(x)fun(x)-kt,0)
y = 1.6230
>> fun(y)
ans = 0.90000
Note that this curve might have multiple solutions, depending on the parameter values: fzero only finds one solution, so you will have to adjust its starting value to find any other solutions. You should always plot the curve and check if the solution is reasonable:
X = -10:0.1:10;
Y = arrayfun(fun,X);
plot(X,Y)
hold on
plot(y,0.9,'*')

类别

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

标签

产品


版本

R2016b

Community Treasure Hunt

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

Start Hunting!

Translated by