clear all for i=1:10 x1(i)=sqrt(25^2-(2.5*i)^2) end y1=2.5:2.5:25; %plot(x1,y1,'.') %pause %z1=zeros(1,10) h=1; z1=ones(1,10)*h; %plot3(x1,y1,z1,'.'); %圆上的点 for i=1:10 x2(i)=sqrt(30^2-h^2); x3(i)=sqrt(60^2-h^2); end y2=2.5:2.5:25; y3=2.5:2.5:25; z2=ones(1,10)*1/2*h; z3=zeros(1,10); plot3(x1,y1,z1,'.'); hold on plot3(x2,y2,z2,'.'); hold on plot3(x3,y3,z3,'.'); hold on for i=1:10 d(i)=60-x1(i)-sqrt((x2(i)-x1(i))^2+(y2(i)-y1(i))^2);%钢筋下方距离 lamda(i)=sqrt((x2(i)-x1(i))^2+(z2(i)-z1(i))^2)/(60-x1(i)); %z3(i)=(lamda(i)*h+z2(i)-h)/lamda(i); z3(i)=((lamda(i)-1)*z1(i)+z2(i))/lamda(i); x3(i)=(x2(i)+x1(i)*(lamda(i)-1))/lamda(i); end pause for i=1:10 x4=[x1(i) x2(i) x3(i)]; y4=[y1(i) y2(i) y3(i)]; z4=[z1(i) z2(i) z3(i)]; plot3(x4,y4,z4) hold on end