一、常微分方程初值問題如下
常用的單步法有如下幾種
1、euler計算格式:
2、改進的euler方法,其迭代式如下
3、經典的r-k方法(四階r-k方法),具體公式如下
題目1 分別用euler計算格式和改進的euler方法求解如下常微分方程初值問題
取步長=0.1,計算到=0.5
題目2 求初值問題
分別用euler方法(h=0.025),二階r-k方法(h=0.05)和經典的r-k方法(h=0.1)計算到,比較計算結果。
二、求解常微分方程初值問題的多步演算法
1、 adams預估-校正系統如下
2、帶修正的adams預估校正系統公式
題目3 編寫上述adams預估-校正系統和帶修正的adams預估校正系統的程式,並用於求解如下初值問題
取步長=0.1,計算到=1.0。
三、實習報告要求
1.簡述方法的基本原理。
2.程式中要加注釋。
3.對程式中的主要變數給出說明。
4.附原程式及計算結果。
5.對各種演算法作比較,對計算結果作簡單分析,談談程式設計上機的體會。
輔助程式如下
演算法1 euler方法計算程式
function e=euler(f,a,b,ya,m)
% f-右端函式字串;a,b-區間端點;ya-初始值;m-步長數;
%e=[t』y』]-t表示橫座標,y表示縱座標
h=(b-a)/m;=zeros(1,m+1);y= zeros(1,m+1);t=a:h:b;y(1)=ya;
for j=1:m
y(j+1)=y(j)+h*feval(f,t(j),y(j));
ende(t』y』);
演算法2 經典四階r-k方法
function r=rk4(f,a,b,ya,m)
% f-右端函式
%a,b-區間端點
%ya-初始值
%m-步長數
%r=[t』y』]-t表示橫座標,y表示縱座標
h=(b-a)/m;
t=zeros(1,m+1);
y= zeros(1,m+1);
t=a:h:b;
y(1)=ya;
for j=1:m
k1=h*feval(f,t(j),y(j));
k2=h*feval(f,t(j)+h/2,y(j)+k1/2);
k3=h*feval(f,t(j)+h/2,y(j)+k2/2);
k4=h*feval(f,t(j)+h,y(j)+k3);
y(j+1)=y(j)+(k1+k2+k3+k4)/6;
end r=[t』y』];
倉儲實習任務
現代物流 儲配方案的設計與執行 專案 一 背景材料 西安新傑家電配銷 成立於2004年8月,是一家主要從事家電產品倉儲銷售的企業。其 商為國內多個家電生產廠家,客戶則主要為一些同城家電批發商。公司的倉儲中心地理位置優越 設施裝置齊全 服務功能完備,不僅承擔著自營產品的儲存配送任務,而且對外承接社會商...
實習報告任務
實習時間 2016年7月1日 2016年8月31日實習地點 中國美術學院風景建築設計研究院 實習單位介紹 中國美術學院風景建築設計研究院始建於1984年,具有建設部頒發的建築工程設計甲級資質 風景園林設計甲級資質 建築裝飾設計甲級資質。設計研究院以中國美術學院為學術背景,以我國第乙個環境藝術系的學術...
實習報告實習任務書
實習報告 課題名稱 學生學號 專業班級 指導教師 實習單位 完成時間 年月日 武漢工程大學郵電與資訊工程學院 實習任務書 一 實習單位 學生填寫 二 實習崗位 學生填寫 三 實習任務 指導教師與學生本人共同制定 四 實習要求 指導教師與學生本人共同制定 五 實習進度 指導教師與學生本人共同制定 武漢...