數(shù)值分析之帶初值的常微分方程數(shù)值解法(二)

【實(shí)驗(yàn)原理】

一迫吐、改進(jìn)的歐拉法(預(yù)報(bào)校正)

1)預(yù)報(bào)校正法基本思路

改進(jìn)歐拉法先用歐拉法求出預(yù)報(bào)值,再利用梯形公式求出校正值溯祸,局部截?cái)嗾`差比歐拉法低了一階叹俏,可以較大程度地提高計(jì)算精度。

2)預(yù)報(bào)校正法算法描述

二嚷往、四階規(guī)范龍格-庫塔公式

【實(shí)驗(yàn)內(nèi)容】

一葛账、回答下面的問題

1.什么是常微分方程的解析解和數(shù)值解?

2.歐拉左矩形公式和右矩形公式是什么,有什么區(qū)別皮仁。

3.解函數(shù)y(x)的光滑性表現(xiàn)在什么地方籍琳,解函數(shù)y(x)“變化劇烈”如何體現(xiàn),請舉例說明贷祈。

4.寫出一個二階龍格-庫塔公式趋急、四階標(biāo)準(zhǔn)龍格-庫塔公式

二、編程計(jì)算課后習(xí)題中規(guī)定題目势誊,交回實(shí)驗(yàn)報(bào)告與計(jì)算

實(shí)驗(yàn)代碼

%定義函數(shù)
function z=f1(x,y)
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
z=1/(x^2)-y/x;
end
%改進(jìn)歐拉法
clear;clc;
h=0.05; 
x=1:h:2; %定義函數(shù)范圍
y(1)=1; %初值

for i=1:20 
    yp=y(i)+h*f1(x(i),y(i));%改進(jìn)(迭代)公式
    yc=y(i)+h*f1(x(i+1),yp);
    y(i+1)=(yp+yc)/2;%校正值
end 

%-----------輸出--------------------
disp('i=0...10的值為:'); 
for i=1:2:21 
    fprintf('%f\n',y(i)) 
end
%-----------輸出--------------------
%四階RK 
clear all
a=1;b=2;h=0.1;%邊界及步長
x0=a;y0=1;%初值
 n=(b-a)/h;
N=[0:1:n]';
x=zeros(n+1,1);
y=zeros(n,1);%用矩陣存儲
x(1)=x0;y(1)=y0;
 for i=1:n
    x(i+1)=x(i)+h;
end
%----------------四階RK--------------------
for i=1:n
    k1=h*f1(x(i),y(i));
    k2=h*f1(x(i)+h/2,y(i)+k1/2);
    k3=h*f1(x(i)+h/2,y(i)+k2/2);
    k4=h*f1(x(i)+h,y(i)+k3);
    y(i+1)=y(i)+(k1+2*k2+2*k3+k4)/6;
end
%----------------四階RK--------------------
[N,x,y]%輸出

實(shí)驗(yàn)結(jié)果

運(yùn)行結(jié)果

實(shí)驗(yàn)代碼(僅改動部分)

%--------------------函數(shù)定義-------------------
function z=f1(x,y)
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
z=-50*y+50*x^2+2*x;
end
%--------------------函數(shù)定義-------------------

a=0;b=1;h=0.1; %邊界及步長
x0=a;y0=1/3; %初值
T=[N,x,y] %輸出解
dy=y-(1/3*exp(-50*x)+x.^2) %輸出與精確解的誤差

實(shí)驗(yàn)結(jié)果過多呜达,此處略去


相關(guān)文章:
MATLAB數(shù)值分析之?dāng)?shù)值積分(一)

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市粟耻,隨后出現(xiàn)的幾起案子查近,更是在濱河造成了極大的恐慌,老刑警劉巖勋颖,帶你破解...
    沈念sama閱讀 211,743評論 6 492
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件嗦嗡,死亡現(xiàn)場離奇詭異,居然都是意外死亡饭玲,警方通過查閱死者的電腦和手機(jī)侥祭,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 90,296評論 3 385
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來茄厘,“玉大人矮冬,你說我怎么就攤上這事〈喂” “怎么了胎署?”我有些...
    開封第一講書人閱讀 157,285評論 0 348
  • 文/不壞的土叔 我叫張陵,是天一觀的道長窑滞。 經(jīng)常有香客問我琼牧,道長,這世上最難降的妖魔是什么哀卫? 我笑而不...
    開封第一講書人閱讀 56,485評論 1 283
  • 正文 為了忘掉前任巨坊,我火速辦了婚禮,結(jié)果婚禮上此改,老公的妹妹穿的比我還像新娘趾撵。我一直安慰自己,他們只是感情好共啃,可當(dāng)我...
    茶點(diǎn)故事閱讀 65,581評論 6 386
  • 文/花漫 我一把揭開白布占调。 她就那樣靜靜地躺著暂题,像睡著了一般。 火紅的嫁衣襯著肌膚如雪究珊。 梳的紋絲不亂的頭發(fā)上薪者,一...
    開封第一講書人閱讀 49,821評論 1 290
  • 那天,我揣著相機(jī)與錄音苦银,去河邊找鬼啸胧。 笑死,一個胖子當(dāng)著我的面吹牛幔虏,可吹牛的內(nèi)容都是我干的纺念。 我是一名探鬼主播,決...
    沈念sama閱讀 38,960評論 3 408
  • 文/蒼蘭香墨 我猛地睜開眼想括,長吁一口氣:“原來是場噩夢啊……” “哼陷谱!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起瑟蜈,我...
    開封第一講書人閱讀 37,719評論 0 266
  • 序言:老撾萬榮一對情侶失蹤烟逊,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后铺根,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體宪躯,經(jīng)...
    沈念sama閱讀 44,186評論 1 303
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,516評論 2 327
  • 正文 我和宋清朗相戀三年位迂,在試婚紗的時候發(fā)現(xiàn)自己被綠了访雪。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 38,650評論 1 340
  • 序言:一個原本活蹦亂跳的男人離奇死亡掂林,死狀恐怖臣缀,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情泻帮,我是刑警寧澤精置,帶...
    沈念sama閱讀 34,329評論 4 330
  • 正文 年R本政府宣布,位于F島的核電站锣杂,受9級特大地震影響脂倦,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜元莫,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,936評論 3 313
  • 文/蒙蒙 一赖阻、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧柒竞,春花似錦、人聲如沸播聪。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,757評論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽森书。三九已至茴迁,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間夹供,已是汗流浹背霎俩。 一陣腳步聲響...
    開封第一講書人閱讀 31,991評論 1 266
  • 我被黑心中介騙來泰國打工哀军, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人打却。 一個月前我還...
    沈念sama閱讀 46,370評論 2 360
  • 正文 我出身青樓杉适,卻偏偏與公主長得像,于是被迫代替她去往敵國和親柳击。 傳聞我的和親對象是個殘疾皇子猿推,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 43,527評論 2 349

推薦閱讀更多精彩內(nèi)容