|
发表于 2009-8-31 15:49:46
|
显示全部楼层
引用第272楼speary于2009-08-31 15:46发表的 :
结果是等于边长。
但是用计算机算好像不符合题目规定的吧,
而且这个活动不是要结果,主要是需要推理过程。 你这道题目用微分一起的话太烦了,而且本来就是计算题,为啥不可以?又不是闹经济转弯类的!!计算机条件下算起来很方便~~一下是我的程序:
Ax=zeros(1,1000);Ay=Ax;% 存放四个点的轨迹坐标
Bx=zeros(1,1000);By=Bx;
Cx=zeros(1,1000);Cy=Cx;
Dx=zeros(1,1000);Dy=Dx;
Ax(1)=0;Ay(1)=0; %初始坐标
Bx(1)=1;By(1)=0;
Cx(1)=1;Cy(1)=1;
Dx(1)=0;Dy(1)=1;
for t=1:999
Ax(t+1)=Ax(t)+0.001*(Bx(t)-Ax(t))/sqrt((Bx(t)-Ax(t))^2+(By(t)-Ay(t))^2);%A向B前进了千分之一
Ay(t+1)=Ay(t)+0.001*(By(t)-Ay(t))/sqrt((Bx(t)-Ax(t))^2+(By(t)-Ay(t))^2);
Bx(t+1)=Bx(t)+0.001*(Cx(t)-Bx(t))/sqrt((Cx(t)-Bx(t))^2+(Cy(t)-By(t))^2);
By(t+1)=By(t)+0.001*(Cy(t)-By(t))/sqrt((Cx(t)-Bx(t))^2+(Cy(t)-By(t))^2);
Cx(t+1)=Cx(t)+0.001*(Dx(t)-Cx(t))/sqrt((Dx(t)-Cx(t))^2+(Dy(t)-Cy(t))^2);
Cy(t+1)=Cy(t)+0.001*(Dy(t)-Cy(t))/sqrt((Dx(t)-Cx(t))^2+(Dy(t)-Cy(t))^2);
Dx(t+1)=Dx(t)+0.001*(Ax(t)-Dx(t))/sqrt((Ax(t)-Dx(t))^2+(Ay(t)-Dy(t))^2);
Dy(t+1)=Dy(t)+0.001*(Ay(t)-Dy(t))/sqrt((Ax(t)-Dx(t))^2+(Ay(t)-Dy(t))^2);
n(t+1)=sqrt((Ay(t+1)-Ay(t))^2+(Ax(t+1)-Ax(t))^2);
end
figure;
plot(Ax,Ay,Bx,By,Cx,Cy,Dx,Dy);
m=sum(n);
结果m=0.9990 |
|