题目

有没有高手帮我用matlab解下这个方程：
syms m g k k1 t t1 v0 s;
m=1000 ;
g=9.80665;
k=2.95;
k1=0.64*s;
t=sqrt(2*3000*1000/k*g);
t1=7;
方程是：((m/(k+k1)*g)*tan(sqrt((k+k1)/m*g)*(t-t1)+atan(sqrt((k+k1)/k)*tan(sqrt(k*g/m)*t1))^2+1/(k^2*t^2/m^2+1/v0)^2+2*k*t/(m*v0)))=25
求 s

提问时间：2021-04-01

答案

syms s v0
m=1000;
g=9.80665;
k=2.95;
k1=0.64*s;
t=sqrt(2*3000*1000/k*g);
t1=7;
f=((m/(k+k1)*g)*tan(sqrt((k+k1)/m*g)*(t-t1)+atan(sqrt((k+k1)/k)*tan(sqrt(k*g/m)*t1))^2+1/(k^2*t^2/m^2+1/v0)^2+2*k*t/(m*v0)))-25;
>> s=solve(f,s)
Warning: Explicit solution could not be found.
> In C:Matlab6p5 oolboxsymbolicsolve.m at line 136
In C:Matlab6p5 oolboxsymbolic@symsolve.m at line 49

s =

[ empty sym ]
v0没有数值
此外方程没有解析解,可以取v0的值后,用数值解