27,580
社区成员
发帖
与我相关
我的任务
分享请问各位大神，哪错了，一直运行不出来结果？
clear
syms v0 o1 g m y b x L R Lm o0
m=800;v0=31.5;g=9.8;b=0.01;
while 1
Lm=-1;
for o1=30:0.01:70
R=m*v0*cos(1/180*o0*pi)/b*(1-exp(-(m*v0*sin(1/180*o0*pi)*b+m^2*g+m^2*g*lambertw(-1/m*(m*g+v0*sin(1/180*o0*pi)*b)/g*exp((-m*v0*sin(1/180*o0*pi)*b-m^2*g)/m^2/g)))/g/m^2));
f=sqrt (1+((m^2*g/b^2+m*v0*sin(1/180*o0*pi)/b)*b/m/v0/cos(1/180*o0*pi)-m^2*g/b/ (-x*b+m*v0*cos(1/180* o0*pi)))^2);
L=int(f,0,R);
if L<Lm
o0=o1-0.1
break
end
Lm=L
end
b=b+0.01
if b>2
break
end
end
...全文
请发表友善的回复…
发表回复
yolyry 2019-05-02
- 打赏
- 举报
用SQL来编程，兄弟，何必呢
