close all
clc
clear all
L1=1000;
L2=1100;
W1=200;
W2=200;
d=1500;
l=(d-L1):round(min(sqrt(L2^2+W2^2/4),L1+d));
thet1=imag(acos((L1.^2+d.^2-l.^2)/2*L1.*d));
thet2u=imag(acos((L1.^2+l.^2-d.^2)/2*L1.*l)+asin((W2/2.*l)-3.14));
thet2l=imag(acos((L1.^2+l.^2-d.^2)/2*L1.*l)-asin((W2/2.*l)-3.14));
plot(thet1,thet2u);
hold on;
plot(thet1,thet2l);
可以编译成功,但不知道为什么和预想的曲线相差太大