春數(shù)學(xué)實(shí)驗(yàn)基礎(chǔ)實(shí)驗(yàn)報(bào)告(1)常微分方程-展示頁
2025-01-27 22:28本頁面
【正文】 1)。y39。x39。2*x+y^239。t39。Dy(0)=139。y(0)=039。*(Dy)^2+2*y=sin(t)39。)。plot(t1,y1(:,1),39。% ode45[t1,y1]=ode45(ex1_2fun,[0,5],[0。f(2)=*y(2).^22*y(1)+sin(t)。Euler法只有一階精度，所以實(shí)際應(yīng)用效率比較差，而ode45的效果比較好，很接近真實(shí)值。解析解39。euler法39。ode4539。pausehold off。)。,39。% 解析解y0=dsolve(39。)。plot(x2,y2,39。168。pausehold on。ko39。[x1,y1]=ode45(odefun,[0,3],1)。y39。x39。x+y39。y=y39。for i=1:length(t)1y(i+1)=y(i)+h*feval(odefun,t(i),y(i))。年級(jí)、專業(yè) 姓名 學(xué)號(hào) 名單序號(hào) 實(shí)驗(yàn)時(shí)間 2013年 3月 日 使用設(shè)備、軟件 PC, MATLAB 注： 實(shí)驗(yàn)報(bào)告的最后一部分是實(shí)驗(yàn)小結(jié)與收獲 實(shí)驗(yàn)一 常微分方程1. 分別用Euler法和ode45解下列常微分方程并與解析解比較： (1) 編寫Euler法的matlab函數(shù)，function [t,y]=euler(odefun,tspan,y0,h)t=tspan(1):h:tspan(2)。y(1)=y0。endt=t39。下面比較三者的差別：% ode45odefun=inline(39。,39。,39。)。plot(x1,y1,39。)。% Euler[x2,y2]=euler(odefun,[0,3],1,)。r+39。pausehol
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1