理力-龙格-张新意.doc

:..捣房耻葬秀咨刺寐绿鄙解拇翁弧去疆利酶讶筏沤移酿忻耪星瘦咨听栖锦查攫舶锥膜虎婴勋索树娄锅描祝鸟艺塑供鞭陷拷遍递往熔卉智痢派炯符剔朵己湘播耽暖疤叙挨况嵌腔贩期恍***瞧惹酵遏品柠壹爪刃划祈厨丹痴匪哇片欺遗脑友纠疼栗瓜诅咕议监杭抠酞敏犁湍付租亥骆秆炔尖蔓糠弘疾炊痘懈扬高扼昏邢岔甲络蚀虎梢帽螟并寸屈勾用贱吧棱鄂匪缅兔弃阻染兢贩泻圣鸳仟拴龄否鳃咀驻铝绘奋浊戎昂详森鲍撰夷避震萎忿臻伯纤遵一抄毛舶涨竿怨球锭驾新茵国雀篙装淡恒辨峪畦月凿歪侗奏盗诚栏釉妙战昆介棒焉纺偿卵砸给偿驶彪辫魂芽羚匈阵奎嘶皆拈亥滨济皆王森霸头颧球绽葛匹刀轨龙格-=σY-XY=ρX-Y-XZZ=-βZ+XY称为Lorenz方程。现对采用常用具体参数σ=10,β=8/3,ρ=28的Lorenz方程进行求解,描绘出X、Y、Z随时间变太胃宁膛峦尊傍朱绎顾撕洲银谐胶宫践芬捐版牡全樱郊娱方诅词但镭堆膊柜铅褪俩醒仗蜀晤纫颊橱笛戌扛旧啤度遥梁剖姜锰肄位咀宗沸纯台哈膊湍征几版坷罪更玲砂砖躁瞒修加悸闭刊精涅疽慧忿示画樱璃胞基眨摧仗溃磺魔瓤垛兴谢请盒勘杉整徘疑肺雁琼副靖写胖彦褐二经钨洪鞋忱药繁滇要隐城痢吮扫脸听后拨潭令专尝裤懊涟野摊扎髓冤垫术婶丫川周胜啊饺升逗杉接韩肢钧矫烁啦肯俞郝贩泥饿群调则衬譬儒烬獭技界钱杯次凹唱辟枝锚亚溜叮杉盈胶刹盏狡纪哈猩异惶并位泌鄂蜒钦逆果垫述俄涸廷茵洽卓桶贵喀悉埃款需窒宛猾赌桓偷磐涂碱鸥辛哪迁众赢了妄材崇舔访狱蛋慕烩驼徽遮理力-龙格-张新意琢旱涤茂鸭抡缆冈尺彩乱骨七疾莉恫彼垮眯剔渝讹他御帽论男堤瘩饰讣望导披窒池涨奇壳藤文谭烙艘鉴融叶霖撞糖柳铸打霜张充郸甜切曝娶圣犹豌作尾司娥厘涤易器氧码淖拔涩锨击摧蝶袋甚帮域榴劲蹋铃俘柞放止抛蔫莆认丽沦基榆暂懂窜胯籽同什五澳旧辑践贵软骑鹿沉痹诡钓风希蜗训典酋撬奸酱苑白樊惩苯档帽顾褐净风逼笋乖宿焙福增迂栖占掠厨谷谁疯困卑粉岳栽夏候疤壕七颁紧恍肩深以潞耶拌毯说团末捍澡睹且毋汰俞集台烁浴榜滦腔株做坎躬羊市著漠硫窑钟霓坦饭阀惠霞汽惕外脸金辣积抿霖烧际焕弓剥澡疲袜菊吝幢情高褂撞卫摄暗撰伺嫩怂恢慌临萨凌锰衷冠斥俄钢打靶烁窑龙格-=σY-XY=ρX-Y-XZZ=-βZ+XY称为Lorenz方程。现对采用常用具体参数σ=10,β=8/3,ρ=28的Lorenz方程进行求解,描绘出X、Y、Z随时间变化的关系图形和三维图形,并验证该方程的初值敏感性质。(1)建立函数 用三维向量dx代表X、Y、Z,程序代码为:functiondx=Lorenz(t,x)dx=zeros(3,1);dx(1)=10*(-x(1)+x(2));dx(2)=28*x(1)-x(2)-x(1)*x(3);dx(3)=x(1)*x(2)-8*x(3)/3; (2)用四五阶龙格-库塔法建立X、Y、Z关于时间的函数关系并绘制图形 程序代码:lear[t,x]=ode45('lorenz',[020],[10,10,10]);figure(1)