Unable to find explicit solution for "solve"

3 次查看（过去 30 天）
Nate Weger
Nate Weger 2019-3-17
编辑： Stephan 2019-3-17
Code summary: Trying to calculate an array of 'Inet' values, one for every 'V' value. All of the other values are individul numbers. The code is below. The reason I believe it is not working is that there are 'Inet' variables on each side of the equation.
V = 0:.01:3;
syms Inet;
I = solve(Inet == Isc - Is.*(exp(q.*(V+Inet*Rs)./(n*kB*Tc))-1)-(V+Inet*Rs)./Rsh, Inet);
Here is the output warning:
Warning: Unable to find explicit solution. For options, see help.
After running the code, 'I' has no value. I would like it to be an array of 301 numbers, just like 'V'. Can anybody see what I'm doing wrong? Any help is appreciated

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

采纳的回答

Stephan
Stephan 2019-3-17
Hi,
solve numeric:
% Values for constants - I took some fantasy values...
Isc = 0.1;
Is = 1;
q = 2;
Rs = 0.0135;
Rsh = -2;
n = 4;
kB = 0.4;
Tc = 10;
% Definition of V as column vector
V = (0:.01:3)';
% Solve for Inet
x0 = ones(size(V,1),1);
Inet = fsolve(@(x)solve_Inet(x,Isc,Is,q,V,Rs,Rsh,n,kB,Tc),x0);
% Inet Function
function F = solve_Inet(x,Isc,Is,q,V,Rs,Rsh,n,kB,Tc)
Inet = x;
F = Isc - Is.*(exp(q.*(V+Inet*Rs)./(n*kB*Tc))-1)-(V+Inet*Rs)./Rsh - Inet;
end
Best regards
Stephan
  5 个评论
显示 3更早的评论
Nate Weger
Nate Weger 2019-3-17
Yes, that's exactly what I was looking for, thank you! Could you please explain what is going on with optimoptions and why you needed a function?
Stephan
Stephan 2019-3-17
编辑：Stephan 2019-3-17
the standard settings use 100 * number of variables. 30100 in your case. The message you recieved told us, that this limit was reached and there is no solution found. I used optimoptions to set it to 500* number of variables, which appears to be enough for a solution.
I choosed a function - also a function handle would work - no special reason. Fsolve tries to find x that statisfies F=0.
If this was useful please accept my answer.

请先登录，再进行评论。

更多回答（0 个）

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by