微分方程數(shù)值解-展示頁(yè)
2025-07-02 00:43本頁(yè)面
【正文】 寫“關(guān)鍵詞”三字，提煉35個(gè)關(guān)鍵詞，用分號(hào)隔開，小4號(hào)黑體)； 正文部分采用三級(jí)標(biāo)題；第1章 (小二號(hào)黑體居中，) 小三號(hào)黑體（段前、）（段前、）參考文獻(xiàn)（黑體小二號(hào)居中，），參考文獻(xiàn)用五號(hào)宋體，參照《參考文獻(xiàn)著錄規(guī)則（GB/T 7714－2005）》。=2*t*u2u0=12 實(shí)驗(yàn)數(shù)據(jù)與實(shí)驗(yàn)結(jié)果1）實(shí)驗(yàn)結(jié)果歐拉法的值與真實(shí)的值改進(jìn)歐拉法的值與真實(shí)的值4階Rungekutta法的值與真實(shí)的值歐拉法的值與真實(shí)的值改進(jìn)歐拉法的值與真實(shí)的值4階Rungekutta法的值與真實(shí)的值從上圖可以看出，Euler法計(jì)算出的結(jié)果與其他方法所得結(jié)果相對(duì)來(lái)說精度太低。歐拉法function[t,u]=euler(f,t0,u0,tf,h)n=(tft0)/h。t(1)=t0。 u(i+1)=u(i)+h*feval(f,t(i),u(i))。u(1)=u0。for i=1:n t(i+1)=t0+i*h。 uc=u(i)+h*feval(f,t(i+1),up)。end四階Rungekutta法function [T,Y]=rk4(f,a,b,ya,N)h=(ba)/N。Y(1)=ya。 k2=h*feval(f,T(j)+h/2,Y(j)+k1/2)。 k4=h*feval(f,T(j)+h,Y(j)+k3)。end實(shí)驗(yàn)二 求解二階常微分初值問題（二）1 實(shí)驗(yàn)內(nèi)容分別用差分法、Jacobi迭代法求解初值問題dudx2=sinx 0x1u0=u1=02 實(shí)驗(yàn)數(shù)據(jù)與實(shí)驗(yàn)結(jié)果（可用文字描述或貼圖的方式進(jìn)行說明） 1）實(shí)驗(yàn)結(jié)果3 程序代碼清單（可直接將可運(yùn)行源代碼粘貼在下面的方框中）Matrix:function [A,b] = matrix(x0,xf,h)n = (xfx0) / h。b = zeros(n1,1)。 A(i, i+1) = 1。endb(n1) = h^2*((n1)*h)。xold=x0。E=eye(n)。d=D\b。while i=max xnew=C*xold+d。 disp(39。252。213。178。)。 else xold=xnew。])。enddisp(39。252。168。187。213。178。)。238。243。252。206。253。243。196。225。170。186。)。Gauss:function [x, k] = Gauss(A, b, x0, to1, max)D = diag(diag(A))。U = triu(A, 1)。f = (D+L) \ b。k = 1。 x = T*x0 + f。 if(k = max) disp(39。252。206。253。171。224。172。201。220。187。213。178。)。 end % disp([k, x39。 207。202。195。178。189。185。endSor:function [x, k] = sor(A, b, x0, w, to1, max)if(w = 0 || w =2) erro。endD = diag(diag(A))。U = t
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1