非线性动力学学习报告.docVIP

非线性动力学学习报告 在课堂上老师以生动活泼的方式介绍了分形的相关知识，特别是展现了一些美丽的分形图案，我对此十分感兴趣，所以课后找了一些相关资料，学会了用仿射变换的循环迭代方法，在MATLAB平台下，实现了一些简单的飞行图案的绘制。具体内容见项目一。其中的数学原理由于我还不是特别清楚，所以在此进仅做一简要汇报，下面会具体叙述用MATLAB绘制分形图案的过程。 在项目二中，探讨了对于一根细长压杆，端部的压力大小与杆件变形之间的关系。这里的端部压力是较大的载荷（即大于临界力），那么经典的材料力学理论便束手无策，这里构建了一个压杆变形的微段迭代模型，把一个大变形非线性问题转化为有限个小变形的迭加，用MATLAB编程迭代计算的结果较好的吻合了铁木辛哥弹性稳定理论中有关压杆弹性屈曲中的一些成果。 项目一：用MATLAB绘制美丽的分形图案 上个世纪60年代，B.Mandelbrot对一个具有复杂几何性质但局部看起来仍然一样的几何对象提出了分形概念。在很多非线性动力学系统等血多领域都会看到分形的例子，随着电子计算机的发展，我们绘制出了很多分形图案。 在这个项目中，实现了用MATLAB来绘制蕨类植物枝叶和著名的Sierpinski三角形；另外还给出了一个通过编程绘制树枝的例子没有用到仿射变换，只是复杂的循环。 经过翻阅相关资料（考文献[1]），我了解到数学中的仿射变换的定义如下：设x是一个n维向量，A是n*n的矩阵，b是与x同维的向量，那么变换称作仿射变换，去不同的A，b就会得到不同的变换结果。如果打印前k次（k应该取较大的值）迭代过程中向量x在坐标系中所表示的所有点，那么就可以得到一幅漂亮的分形图案。其中矩阵A和向量b的取法涉及到很复杂的数学理论，在这里不做详细介绍。 基于前面的理论分析很容易得到MATLAB绘图程序代码及其运行结果。 1.、使用数学中的仿射变换理论，绘制蕨类植物枝叶 程序：%fenxing_juelei.m %蕨类植物模拟 x = [.5; .5]; %初值 h = plot(x(1),x(2),.); %绘制初值点 %设置用于后面随机数的判别向量 p = [ .85 .92 .99 1.00]; b1 = [0; 1.6]; b2 = [0; 1.6]; b3 = [0; .44]; %------仿射变换矩阵 A1 = [ .85 .04; -.04 .85]; A2 = [ .20 -.26; .23 .22]; A3 = [-.15 .28; .26 .24]; A4 = [ 0 0 ; 0 .16]; for i=1:20000 r = rand; %产生随机数 if r p(1) x = A1*x + b1; elseif r p(2) x = A2*x + b2; elseif r p(3) x = A3*x + b3; else x = A4*x; end plot(x(1),x(2),g),hold on %采用绿色绘制 end axis off %取消坐标轴 把该m文件放置到Matlab的当前工作目录下在命令行中输入fenxing_juelei，变得到了下面的运行结果。 在具体操作中，迭代了20000次，电脑计算了十多秒钟才计算绘制完毕，着看起来好像迭代次数，可是得到的图形靠下的部分还是不够细密。而实际上再多尺度上分型图形保持了其固有的形态，就是说当把分形图形的一部分无穷放大后仍然保留了其原来的形状特性。如果增加迭代次数显然可以得到更为细密的蕨类植物枝叶图形。 2.、绘制Sierpinski三角形分形图案 程序：fenxing_sie.m %Sierpinski三角形分形图案 %初值 x = [0; 0]; h = plot(x(1),x(2),.); %绘制初值 A=[1/2 0;0 1/2]; %仿射变换矩阵 b1=[0 0]; b2=[1/2 0]; b3=[1/4 sqrt(3)/4]; for i=1:10000 r = rand; %产生随机数 if r 0.66666666667 x = A*x + b1; elseif r 0.333333333333 x=A*x+b2; else x = A*x+b3; end plot(x(1),x(2),-),hold on %保持图形 end axis off %取消坐标轴 程序编制完成后将其放在MATLAB当前工作目录下，在命令行中输入fen