數(shù)值分析實(shí)驗(yàn)題(-華科)-資料下載頁(yè)
2025-03-25 02:50本頁(yè)面
【正文】 9。,14)。grid onfunction testCharpt5_1%對(duì)數(shù)值分析實(shí)驗(yàn)題第3章第1題進(jìn)行分析%輸入?yún)?shù)：參數(shù)a，步長(zhǎng)h%精確解和數(shù)值解圖形對(duì)比%%第1問(wèn)輸入a=[2 1 1 2]% 輸入a的取值h=[ ]%輸入h的取值%%第2問(wèn)輸入% a=[20 20]% 輸入a的取值% h=[ ]%輸入h的取值%%func=inline(39。1+(yx).*a39。)。%定義函數(shù)for i=1:length(a) x=0:h(i):1。%求解區(qū)間 y=x。 N=length(x)。 y(1)=1。 for n=1:N1 k1=func(a(i),x(n),y(n))。 k2=func(a(i),x(n)+h(i)/2,y(n)+k1*h(i)/2)。 k3=func(a(i),x(n)+h(i)/2,y(n)+k2*h(i)/2)。 k4=func(a(i),x(n)+h(i),y(n)+k3*h(i)) 。 y(n+1)=y(n)+h(i)*(k1+2*k2+2*k3+k4)/6。%數(shù)值解 end y0=exp(a(i)*x)+x。%精確解% figure()%如果疊繪圖去掉此句命令 plot(x,y0) hold on plot(x,y,39。*39。) legend(39。精確解39。,39。數(shù)值解39。) xlabel(39。x39。)。 ylabel(39。y39。)。 title([39。微分方程數(shù)值解,a=39。,num2str(a(i)),39。,h=39。,num2str(h(i))],39。FontSize39。,14)。 grid onend
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1