二阶龙格库塔法迭代公式用Matlab怎么编程
四阶龙格库塔法精度为4，属于单步递推法，单步递推法的基本思想是从（x(i),y(i)）点出发，以某一斜率沿直线达到（x(i+1),y(i+1)）点，从上述定义可以看出，龙格库塔实质上是求一阶微分方程，对其进行排列后就可以进入Matlab进行编程。
但是如果将一阶导看作变量，则二阶导也不过是这个变量的一阶导而已，对于下述二阶方程：基本思想如下令位移为q的一阶导，即位移的一阶导（速度）为q的二阶导。
是用一个简单的迭代式来画分叉混沌图的，从这个例子里，不难学到如何如何画分叉图，其原理不难，两个循环来在2维平面作图，一个if筛选迭代的结果。
还是类似上面那个例子，这里给出不一样的代码实现方法。没看懂的可以再看一遍，里面主要的还是两个for循环，下面这个代码是用第三个for循环来实现对迭代结果的筛选的，这里不必纠结细节，功能实现用if或者for都可以。
依次类推，如果在区间[xi,xi+1]内多预估几个点上的斜率值K1、K2、……Km，并用他们的加权平均数作为平均斜率K*的近似值，显然能构造出具有很高精度的高阶计算公式。经数学推导、求解，可以得出四阶龙格－库塔公式，也就是在工程中应用广泛的经典龙格－库塔算法这样就完成了matlab的编程。
绛旓細鍥闃堕緳鏍煎簱濉旀硶绮惧害涓4锛屽睘浜庡崟姝ラ掓帹娉曪紝鍗曟閫掓帹娉曠殑鍩烘湰鎬濇兂鏄粠锛坸(i),y(i)锛夌偣鍑哄彂锛屼互鏌愪竴鏂滅巼娌跨洿绾胯揪鍒帮紙x(i+1),y(i+1)锛夌偣锛屼粠涓婅堪瀹氫箟鍙互鐪嬪嚭锛岄緳鏍煎簱濉斿疄璐ㄤ笂鏄眰涓闃跺井鍒嗘柟绋嬶紝瀵瑰叾杩涜鎺掑垪鍚庡氨鍙互杩涘叆Matlab杩涜缂栫▼銆備絾鏄鏋滃皢涓闃跺鐪嬩綔鍙橀噺锛屽垯浜岄樁瀵间篃涓嶈繃鏄繖涓彉閲忕殑...
绛旓細杩欓噷K锛漟(Xn+th,Y(Xn+th))绉颁负骞冲潎鏂滅巼锛岄緳鏍煎簱濉旀柟娉曞氨鏄眰寰桲鐨勪竴绉嶇畻娉曘傚埄鐢ㄨ繖鏍风殑鍘熺悊锛岀粡杩囧鏉傜殑鏁板鎺ㄥ锛堣繃浜庣箒鐞愮渷鐣ワ級锛屽彲浠ュ緱鍑烘埅鏂宸负O(h^5)鐨勫洓闃堕緳鏍煎簱濉斿叕寮锛欿1锛漟(Xn,Yn);K2=f(Xn+h/2,Yn+(h/2)*K1);K3=f(Xn+h/2,Yn+(h/2)*K2);K4=f(Xn+h,Yn+h*K3...
绛旓細鍥闃堕緳鏍-搴撳鏄眰瑙ｅ父寰垎鏂圭▼锛堝父寰垎鏂圭▼缁勶級绮惧害鏈楂樼殑涓绉嶆暟鍊兼柟娉曘傚洓闃堕緳鏍-搴撳杩唬鍏紡涓 鏍规嵁鍥涢樁榫欐牸-搴撳杩唬鍏紡锛屽彲鐢╩atlab杞欢缂栧啓寰楀埌濡備笅璁＄畻浠ｇ爜锛氳皟鐢ㄦ柟娉曪細clc;close all;y0=1;h=0.1;a=0;b=1;[x,y] = runge_kutta(@func,y0,h,a,b);disp(' x y')[x',y']...
绛旓細鏁板煎垎鏋愪腑锛榫欐牸锛搴撳娉锛圧unge-Kutta methods锛夋槸鐢ㄤ簬闈炵嚎鎬у父寰垎鏂圭▼鐨勮В鐨勯噸瑕佺殑涓绫婚殣寮忔垨鏄惧紡杩唬娉銆傝繖浜涙妧鏈敱鏁板瀹跺崱灏斅烽緳鏍煎拰椹竵路濞佸皵娴峰路搴撳浜1900骞村乏鍙冲彂鏄庛傞緳鏍-搴撳(Runge-Kutta)鏂规硶鏄竴绉嶅湪宸ョ▼涓婂簲鐢ㄥ箍娉涚殑楂樼簿搴﹀崟姝ョ畻娉曪紝鍏朵腑鍖呮嫭钁楀悕鐨勬鎷夋硶锛岀敤浜庢暟鍊兼眰瑙ｅ井鍒嗘柟绋嬨傜敱浜庢...
