Finding x for y = f(x) in simulink
9 views (last 30 days)
Show older comments
I have an equation which is uel = 3 + 3 * iel + 3 * log (3 * iel +1) a.k.a uel = f(iel) and I want to find iel ( uel is already known), so I have used fsolve to solve the problem in matlab and it works. Here is the code
---------------------------------------------------------------------
function iel = test(uel)
f = @(x) uel-(3 + 3 * x + 5 * log(3 * x + 1));
iel = fsolve(f,uel);
end
----------------------------------------------------------------------
But when i try to add this code into matlab function block in simulink, the model just don't run. I want my input to be uel and output to be iel. Can anyone please help me on this problem. Thank you so so much
2 Comments
Answers (2)
Paul
on 19 Jan 2022
Edited: Paul
on 19 Jan 2022
According to the Symbolic Math Toolbox, the solution can be expressed as
%iel = wrightOmega(uel/3 - log(3) - 2/3) - 1/3
assuming that iel and uel are real. If they can be complex, the expression for iel is the same, but there are some other conditions that have to be satisfied.
Check
%doc wrightOmega
for more info.
Check a couple of test cases (note that I changed the 5 to a 3 in the function test() to match the equation in the question)
uel = 0.3
iel = test(uel)
iel = wrightOmega(uel/3 - log(3) - 2/3) - 1/3
uel = 5.4
iel = test(uel)
iel = wrightOmega(uel/3 - log(3) - 2/3) - 1/3
function iel = test(uel)
% f = @(x) uel-(3 + 3 * x + 5 * log(3 * x + 1)); % doesn't match equation in Question
f = @(x) uel-(3 + 3 * x + 3 * log(3 * x + 1));
iel = fsolve(f,uel);
end
0 Comments
Dat Trương Lê Trí
on 6 Jul 2022
Hi, I am facing the same problem, did you have any posible solution for this case?
0 Comments
See Also
Categories
Find more on Spectral Measurements in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!