東南大學(xué)計(jì)算方法與實(shí)習(xí)實(shí)驗(yàn)報(bào)告-資料下載頁(yè)
2025-08-14 21:43本頁(yè)面
【正文】 } coutsetw(15)f2setw(15)fabs(f2q)endl。 if(fabs(f2q))break。 q=f2。 n2++。 f2=0。 } coutsetprecision(7) 用 復(fù) 化 辛 卜 生 公 式 計(jì) 算 得 f1=f2, 迭 代 次 數(shù) 為 ： n=n2endl。 return 0。 }l 運(yùn)行窗口實(shí)習(xí)題七用四階龍格庫(kù)塔公式（取h=）解下列微分方程組初值問(wèn)題：（1） l 實(shí)驗(yàn)代碼C語(yǔ)言程序如下：includevoid Runge_Kutta(float(*f1)(float x,float y1,float y2),float(*f2)(float x,float y1,float y2),float a,float b,float y0,int N){ float x=a,K1,K2,K3,K4,L1,L2,L3,L4,y1=y0,y2=y0。 float h=(ba)/N。 int i。 printf(x[0]=%f\ty1[0]=%f\ty2[0]=%f\n,x,y1,y2)。 for(i=1。i=N。i++){ K1=(*f1)(x,y1,y2)。 L1=(*f2)(x,y1,y2)。 K2=(*f1)(x+h/2,y1+h*K1/2,y2+h*L1/2)。 L2=(*f2)(x+h/2,y1+h*K1/2,y2+h*L1/2)。 K3=(*f1)(x+h/2,y1+h*K2/2,y2+h*L2/2)。 L3=(*f2)(x+h/2,y1+h*K2/2,y2+h*L2/2)。 K4=(*f1)(x+h,y1+h*K3,y2+h*L3)。 L4=(*f2)(x+h,y1+h*K3,y2+h*L3)。 y1=y1+h*(K1+2*K2+2*K3+K4)/6。 y2=y2+h*(L1+2*L2+2*L3+L4)/6。 x=a+i*h。 printf(x[%d]=%f\ty1[%d]=%f\ty2[%d]=%f\n,i,x,i,y1,i,y2)。 }}float f1(float x,float y1,float y2){ return 1202*y1+2*y2。}float f2(float x,float y1,float y2){ return 2*y15*y2。}void main(){ float a=0,b=1,y0=0。 Runge_Kutta(f1,f2,a,b,y0,10)。}l 運(yùn)行窗口小結(jié)通過(guò)對(duì)計(jì)算方法一個(gè)學(xué)期的學(xué)習(xí)，我受益匪淺。理論和上機(jī)實(shí)踐的結(jié)合讓我對(duì)知識(shí)有了更進(jìn)一步的理解和掌握。并讓我能學(xué)會(huì)運(yùn)用C++的功能來(lái)代替有些筆算起來(lái)很復(fù)雜的計(jì)算。不僅方便了個(gè)人，而且還緊緊跟隨現(xiàn)代科技發(fā)展的腳步。12
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1