请求MATLAB解方程高手：MATLAB解超越方程,详情见问题补充.方程：f=4*((z1-hatx+xt1)^2+((hatx-xt1)*cotd(at))^2-C1)^0.5+((A2-xt1)^2-C2)^0.5-D=0其中：z1 = 20;z2 = 80;hatx = 1;at = 20;at1= [at-5:0.1:at+15];A2 = (z1+z2)*cosd(at)/cosd(at
网友回答
z1 = 20;
z2 = 80;
hatx = 1;
at = 20;
AT1= [at-5:0.1:at+15];
for i=1:length(AT1)
at1=AT1(i);
A2 = (z1+z2)*cosd(at)/cosd(at1)+2*hatx-z1;
C1 = (z1*cosd(at))^2;
C2 = (z2*cosd(at))^2;
D = (z1+z2)*cosd(at)*tand(at1);
[email protected](x)4*((z1-hatx+x)^2+((hatx-x)*cotd(at))^2-C1)^0.5+((A2-x)^2-C2)^0.5-D;
xx=fzero(f,0); %这个命令只能找出一个实数解,at1值较小时,好像没有实数解
x1(i)=xx;
endplot(AT1,x1)