Question

matlab 的solve函数求解等式中的未知数，未知数形式为N^-小数，如下，求解多个小时，怎么办？

clc
clear
d=2.1;
w=11;
R=0
Smax=570;
Smin=0;
Sr=Smax-Smin
Sm=Smax*(1+R)/2
Kt=1.96;
r1=(w/2)^2/d
Kf=1+(Kt-1)/(1+(0.46904415*r1)^0.5)
syms N
e = (Kf*Sr)^2/(107000*Sr*2)
eqn=((1286-Sm)/107000)*(2*N)^-0.0875 + 0.15*(2*N)^-0.6594 ==e
Ni=solve(eqn,N)
Ni

Answer 1

题主给出的方程，用solve函数求解要花费多个小时，为了解决该问题，你可以考虑其他的解方程函数（如vpasolve,二分法，牛顿法等）。虽然solve和vpasolve是属于同一类函数，但vpasolve函数擅长求解方程的数值解，而solve函数擅长求解方程的解析解和数值解，所以vpasolve求解速度要高于solve函数。

对于本例，只要vpasolve替代solve函数，就可以顺利得到其数值解。

。。。。。。

syms N

eqn=((1286-Sm)/107000)*(2*N)^-0.0875 + 0.15*(2*N)^-0.6594 ==e

Ni=vpasolve(eqn,N)