dt时间一直到k个dt时间后的分析,可以得到:=-;=-=-;=-所以接下来我们可以进行MATLAB的编程。MATLAB程序如下:clear;clf;A=[0,0];B=[0,1];C=[1,1];D=[1,0];k=0;dt=0.004;v=1;grid;holdon;axis([0101]);whilek<10000k=k+1;plot(A(1),A(2),'g.');plot(B(1),B(2),'r.');plot(C(1),C(2),'y.');plot(D(1),D(2),'b.');e1=(B-A)/norm(B-A);A=A+v*dt*e1;e2=(C-B)/norm(C-B);B=B+v*dt*e2;e3=(D-C)/norm(D-C);C=C+v*dt*e3;e4=A-v*dt*e1-D;d4=norm(e4);e4=e4/d4;D=D+v*dt*e4;d1=norm(B-A);d2=norm(C-B);d3=norm(D-C);d4=norm(A-D);fprintf('k=%.0fd1=%.4fd2=%.4fd3=%.4fd4=%.4f\n',k,d1,d2,d3,d4);fprintf('A(%.2f,%.2f)B(%.2f,%.2f)C(%.2f,%.2f)D(%.2f,%.2f)\n',A(1),A(2),B(1),B(2),C(1),C(2),D(1),D(2));ifd1<=0.005breakendifd2<=0.005breakendifd3<=0.005breakendifd4<=0.005breakendpause(0.01) end 最终获得的图像如下图所示:七、结论与分析:该实验建立在我们事先做好的假设的基础上进行,结论也和我们最初的判断相吻合。同时通过这个实验,我们还可以看出,数学无处不存在着美,那是一种均衡、自然、对称的美。
全文预览
