可以根据t的范围搜索
clc
h=4.3532;
k=93.5;
i=4*pi/180;
j=3*pi/180;
r=74.25;
a=(h-(r-sqrt(r*r-h*h))*tan(j)-k*tan(i))/tan(j);
b=k-(h-(r-sqrt(r*r-h*h))*tan(j))/tan(i);
c=h-(r-sqrt(r*r-h*h))*tan(j)-k*tan(i);
w=4000;
r1=69.75;
format short
i0=0;
for t=-10:0.5:10
i0=i0+1;
pp(i0)=(b-(b/a)*(r1+(r*cos((180*asin(h/r)/pi-w*t)*pi/180)-r-r1)*cos(3*w*t*pi/180))-((b/c)*r*sin((180*asin(h/r)/pi-w*t)*pi/180)+(r1-r*cos((180*asin(h/r)/pi-w*t)*pi/180)-r)*sin(3*w*t*pi/180)))/(cos(3*w*t*pi/180)+(b/a)*sin(3*w*t*pi/180))-93.5;
if pp(i0)<0.0001
t
end
end
本文如未解决您的问题请添加抖音号:51dongshi(抖音搜索懂视),直接咨询即可。