Solve equation with multiple values for V

Question

Alexander Madsen 2021-10-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1570683-solve-equation-with-multiple-values-for-v

评论： Star Strider 2021-10-25

Hi, i'm trying so solve an equation with 19 different values for V and Pv. It's possible when i do it one at a time, but not when I'm trying with a loop.

Vr=12;
Vi=3;
Pr=max(P_test_data4);
Pv=P_test_data4;
V=WS_test_data4;
syms n
for i=0:length(WS_test_data4);
    
    solve((Pr*((V.^n)-(Vi.^n)))/((Vr.^n)-(Vi.^n))==Pv,n)
end
V=3.53;
Pv=73;
solve((Pr*((V^n)-(Vi^n)))/((Vr^n)-(Vi^n))==Pv,n)

The abbove code does'nt give an answer just " Empty sym: 0-by-1 "

The correct for the first one should be:

V=3.53;
Pv=73;
solve((Pr*((V^n)-(Vi^n)))/((Vr^n)-(Vi^n))==Pv,n)
ans = 2.6088704281837088027315795939332

Thanks in advance

CODE:

Vr = 12;
Vi = 3;
Pr = 5000;
syms n
V = [3;3.5;4;4.51;5;5.5;6;6.5;7;7.5;8.01;8.5;9;9.5;10;10.49;10.99;11.49;12]
Pv = [13;73;201;385;619;864;1174;1508;1894;2271;2710;3150;3523;3918;4244;4496;4711;4834;5000]
for i = 1:length(V)
    
    n = vpasolve((Pr*((V.^n)-(Vi.^n)))/((Vr.^n)-(Vi.^n))==Pv,n)
    
end
n = mean(n)

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2021-10-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1570683-solve-equation-with-multiple-values-for-v#answer_815843

Another approach —

Vr = 12;
Vi = 3;
Pr = 5000;
syms n
V = [3;3.5;4;4.51;5;5.5;6;6.5;7;7.5;8.01;8.5;9;9.5;10;10.49;10.99;11.49;12];
Pv = [13;73;201;385;619;864;1174;1508;1894;2271;2710;3150;3523;3918;4244;4496;4711;4834;5000];
for i = 1:length(V)
    
    nc{i} = vpasolve((Pr*((V(i).^n)-(Vi.^n)))/((Vr.^n)-(Vi.^n)) == Pv(i),n);
    
end
idx = cellfun(@(x)~isempty(x), nc)
idx = 1×19 logical array
   0   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1
nm = sprintf('%23.15E',mean(double([nc{idx}]),'omitnan'))
nm = '  1.173348363730592E+00'