How to solve this equation Numerically [VPASOLVE]

Question

Meddour Aissam riad 2020-10-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/603277-how-to-solve-this-equation-numerically-vpasolve

回答： Alan Stevens 2020-10-1

采纳的回答： Alan Stevens

在 MATLAB Online 中打开

Hi there,

When i try to numerically solve my equation i got the following error :

Error using mupadengine/feval_internal (line 172)

More equations than variables is only supported for polynomial systems.

I simplified the equation to give you an example, and here is the equation example (the matlab code):

syms Idn
%%data
    Ism=200;
    Wb=1500;
    speedend=6000;
    Vsm=100;
    Ws=Wb:15:speedend;
    
    Eq2=Idn==sqrt(Idn)+(((2.4*10.^(-5))*(Idn).^2-0.013*(Idn+3.3)*10^(-3)))*Ws/Vsm %%The equation
    Idn=vpasolve(Eq2,Idn)

Can anyone help me to solve this kind of equation?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Star Strider 2020-10-1

在 MATLAB Online 中打开

Probably not.

When I ran this:

syms Idn
%%data
Ism=200;
Wb=1500;
speedend=6000;
Vsm=100;
Ws=Wb:15:speedend;
Eq2=Idn==sqrt(Idn)+(((2.4*10.^(-5))*(Idn).^2-0.013*(Idn+3.3)*10^(-3)))*Ws/Vsm %%The equation
Idns=solve(Eq2,Idn)
Idn = simplify(vpa(Idns), 'Steps',150)

the result was:

Idns =
Empty sym: 0-by-1
 
Idn =
Empty sym: 0-by-1

.

请先登录，再进行评论。

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

Answer 1

Alan Stevens 2020-10-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/603277-how-to-solve-this-equation-numerically-vpasolve#answer_503725

在 MATLAB Online 中打开

You can solve it numerically using fzero.

Ism=200;
Wb=1500;
speedend=6000;
Ws=Wb:15:speedend;
Idn = zeros(size(Ws));
Idn0 = 1;
for i = 1:numel(Ws)
Idn(i) = fzero(@f,Idn0,[],Ws(i));
end
plot(Ws,Idn),grid
xlabel('Ws'),ylabel('Idn')
function  F = f(Idn,Ws)
    Vsm=100;
    F = sqrt(Idn)+(((2.4*10.^(-5))*(Idn).^2-0.013*(Idn+3.3)*10^(-3)))*Ws/Vsm-Idn; %%The equation
 end
    