用MATLAB龙格库塔法解决二阶微分方程y''+ay=0，a为常数，可以随便设，初值为y0=0,程序怎么写啊？？ Matlab四阶龙格库塔法求解常微分方程
matlab \u600e\u6837\u7528\u9f99\u683c\u5e93\u5854\u6cd5\u6c42\u4e8c\u9636\u5fae\u5206\u65b9\u7a0b\u521d\u503c\u8981\u786e\u5b9a
\u5047\u8bbea=b=c=d=1;
\u7ed3\u679c\u4e3a\uff1a
\u7528Matlab\u56db\u9636\u9f99\u683c\u5e93\u5854\u6cd5\u6c42\u5e38\u5fae\u5206\u65b9\u7a0b\u53ef\u4ee5\u6309\u7167\u4ee5\u4e0b\u65b9\u6cd5\u53bb\u5b9e\u73b0\u3002
1\u3001\u9996\u5148\u5efa\u7acb\u81ea\u5b9a\u4e49\u5fae\u5206\u65b9\u7a0b\u51fd\u6570
function f = ode_fun(x,y)
f=y+2*x/y^2;
end
2\u3001\u7136\u540e\u7528\u56db\u9636\u9f99\u683c\u5e93\u5854\u6cd5\u6c42\u5176\u6570\u503c\u89e3
figure(2)
y0=[1]; %\u521d\u503cy\uff080\uff09=1
h=0.1;
a=0;
b=5;
[x,y] = runge_kutta(@(x,y)ode_fun(x,y),y0,h,a,b);
disp(' x y')
A=[x',y']
plot(x,y,'LineWidth',1.5),grid on
xlabel('x'),ylabel('y(x)');
3\u3001\u7ed8\u5236y-x\u7684\u66f2\u7ebf\u56fe
[T,Y] = ode45(@vdp1000,[0 3000],[2 0]); 这是龙格库塔4阶算法的示例，[2,0]为初值。
绛旓細鍐matlab鍛戒护绐楀彛杈撳叆 doc ode45 鍙互鏌ョ湅榫欐牸搴撳绠楁硶鐨勮缁嗚В閲婂拰鐢ㄦ硶 [T,Y] = ode45(@vdp1000,[0 3000],[2 0]); 杩欐槸榫欐牸搴撳4闃剁畻娉曠殑绀轰緥锛孾2,0]涓哄垵鍊笺
绛旓細璇ヤ簩闃跺井鍒嗘柟绋嬬敤榫欐牸搴撳娉鍙互杩欐牱鏉ユ眰瑙ｃ傜涓姝ワ紝鏍规嵁璇ヤ簩闃跺井鍒嗘柟绋嬶紝鑷畾涔夊井鍒嗘柟绋嬪嚱鏁帮紝func(t,y)绗簩姝ワ紝鏍规嵁鍒濆鏉′欢锛岀‘瀹歽鍜寉'鐨勫垵鍊硷紝鍗硑0=[0,0]绗笁姝ワ紝浣跨敤ode45鍑芥暟姹傝В銆恡,y銆戠殑鏁板艰В锛屽嵆 [t,y] = ode45(@func,[0 0.0005],y0);绗洓姝ワ紝鏍规嵁t銆亂銆亂'鍊硷紝缁樺埗t鈥攜...
绛旓細闇瑕佺敤榫欐牸-搴撳娉姹傝В,涔熷氨鏄痬atlab閲岀殑ode23/ode45鍑芥暟姹傝В澶氫釜浜岄樁寰垎鏂圭▼缁勪负绠鍖栭棶棰,姹傚姪澶х鎻愮ず濡備綍鐢╩atlab鐢ㄦ暟鍊间及璁＄殑鏂规硶姹傝В绫讳技鍥剧墖涓殑鏂圭▼缁勩... 闇瑕佺敤榫欐牸-搴撳娉曟眰瑙,涔熷氨鏄痬atlab閲岀殑ode23/ode45鍑芥暟姹傝В澶氫釜浜岄樁寰垎鏂圭▼缁勪负绠鍖栭棶棰,姹傚姪澶х鎻愮ず濡備綍鐢╩atlab鐢ㄦ暟鍊间及璁＄殑鏂规硶姹傝В绫讳技鍥剧墖涓殑...
绛旓細2銆佸缓绔榫欐牸搴撳绠楁硶鍑芥暟runge_kutta锛堬級璋冪敤鏍煎紡锛歔t,x] = runge_kutta(@(t,x)func(t,x),x0,h,a,b);3銆佺劧鍚庢牴鎹畑鍜寈'鏁版嵁锛岀粯鍒跺嚭x(t)銆亁鈥(t)鐨勫浘褰俻lot锛坸锛堬細锛1锛,x锛堬細锛2锛夛級
绛旓細缂栧啓M鏂囦欢,鐢榫欐牸搴撳娉缂栧埗绋嬪簭,鍗冲彲姹傝В.鎴栬呯敤娆ф媺娉,鎴栬呰绠鍗曠偣鐨勫涓:>> y=dsolve('D2y=70*t-80*t^2', 'y(0)=0,y(20)=0')y = -20/3*t^4+35/3*t^3+146000/3*t >> y=subs(y,'t','x')y = -20/3*(x)^4+35/3*(x)^3+146000/3*(x)...
绛旓細姹傝В浜岄樁寰垎鏂圭▼锛屽垵濮嬫潯浠惰繕闇瑕佺粰鍑簓1'(0)鍜寉2'(0)銆傝繖閲屾殏鏃舵寜鐓0澶勭悊銆俧unction zd530003514 a=0.1;b=0.1;Y0 = [b-1; 0; b; 0];瑙ｆ柟绋 [t,Y]= ode45(@ode,[0 10],Y0);y1=Y(:,1);y2=Y(:,3);缁樺浘 subplot 211 plot(t,y1);subplot 212 plot(t,y2);寰垎鏂圭▼瀹氫箟...
绛旓細鍥涢樁榫欐牸搴撳娉绮惧害涓4锛屽睘浜庡崟姝ラ掓帹娉曪紝鍗曟閫掓帹娉曠殑鍩烘湰鎬濇兂鏄粠锛坸(i),y(i)锛夌偣鍑哄彂锛屼互鏌愪竴鏂滅巼娌跨洿绾胯揪鍒帮紙x(i+1),y(i+1)锛夌偣锛屼粠涓婅堪瀹氫箟鍙互鐪嬪嚭锛岄緳鏍煎簱濉斿疄璐ㄤ笂鏄眰涓闃跺井鍒嗘柟绋嬶紝瀵瑰叾杩涜鎺掑垪鍚庡氨鍙互杩涘叆Matlab杩涜缂栫▼銆備絾鏄鏋滃皢涓闃跺鐪嬩綔鍙橀噺锛屽垯浜岄樁瀵间篃涓嶈繃鏄繖涓彉閲忕殑...
绛旓細榫欐牸搴撳娉姹備簩闃跺井鍒嗘柟绋嬫楠ゅ涓嬶細1銆侀鍏堥渶瑕侀夋嫨涓涓傚綋鐨勫垵濮嬪硷紝浠ュ強涓涓闀縣銆傚垵濮嬪肩殑閫夋嫨閫氬父鍩轰簬闂鐨勫垵濮嬫潯浠讹紝鑰屾闀縣鍒欏喅瀹氫簡杩唬鐨勭簿搴︺2銆佷娇鐢ㄩ夋嫨鐨勫垵濮嬪煎拰姝ラ暱h锛岃绠楀嚭绗竴涓偣鐨勫紋锛坔锛夛紝杩欏彲浠ラ氳繃鍏紡y锛0锛+h*f锛0锛夊緱鍑猴紝鍏朵腑f锛0锛夋槸鍒濆鏉′欢鍑芥暟銆3銆佽绠楃涓涓...
绛旓細鐢榫欐牸搴撳姹傝В浜岄樁寰垎鏂圭▼缁勭殑瀹炵幇鏂规硶:鈪0=[0.5 0];h=0.1;a=0;b=1;[t,x] = runge_kutta(@(t,x)func(t,x),x0,h,a,b);disp(' t x(t) x鈥(t)')A=[t',x'];disp(A)figure,plot(t,x,'.-'),grid on xlabel('t'),ylabel('x(t),x鈥(t)');title(...
绛旓細濡傛灉棰勫厛姹涓涓偣鐨勬枩鐜囧氨鏄簩闃榫欐牸搴撳娉锛屽鏋滈鍏堝彇鍥涗釜鐐瑰氨鏄洓闃堕緳鏍煎簱濉旀硶銆備竴闃跺父寰垎鏂圭▼鍙互鍐欎綔锛歽'=f(x,y),浣跨敤宸垎姒傚康銆(Yn+1-Yn)/h= f(Xn,Yn)鎺ㄥ嚭锛堣繎浼肩瓑浜庯紝鏋侀檺涓篩n'锛塝n+1=Yn+h*f(Xn,Yn)鍙﹀鏍规嵁寰垎涓煎畾鐞嗭紝瀛樺湪0<t<1,浣垮緱 Yn+1=Yn+h*f(Xn+th,Y(Xn+...
