亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        四階龍格-庫塔方法的程序設(shè)計(jì)與應(yīng)用

        2019-10-31 04:14:07羅麗珍吳慶軍
        科教導(dǎo)刊·電子版 2019年24期
        關(guān)鍵詞:龍格常微分方程

        羅麗珍 吳慶軍

        摘 要 本文通過介紹四階龍格-庫塔方法,通過預(yù)報(bào)斜率和泰勒展開式推導(dǎo)出龍格—庫塔格式。了解它的基本思想與算法步驟、MATLAB語言編寫的程序。列舉一些例子,運(yùn)用四階龍格-庫塔方法的MATLAB程序在軟件中運(yùn)行,求解出常微分方程的數(shù)值解,同時(shí)將求解出的數(shù)值解與精確解進(jìn)行比較。

        關(guān)鍵詞 龍格-庫塔方法 常微分方程 數(shù)值解

        中圖分類號(hào):TP337文獻(xiàn)標(biāo)識(shí)碼:A

        0引言

        從17世紀(jì)以來國內(nèi)外數(shù)學(xué)家對(duì)常微分方程的研究取得了很多的成果.歐拉在研究中指出常微分方程存在唯一解和無數(shù)解,他用近似值求解微分方程,發(fā)現(xiàn)用積分因子求解微積分方程的特殊算法。拉格朗日建立了一階微分方程理論,他將參數(shù)變法應(yīng)用到四階非齊次方程的求解。

        我們生活中許多問題的解決都運(yùn)用到常微分方程,常微分方程的數(shù)值解法中經(jīng)常使用的方法是四階龍格-庫塔方法。各個(gè)領(lǐng)域和工程問題中的原理和演變規(guī)律都是用常微分方程來描述的,如在物理方面的電路中電流變化的規(guī)律、航天航空方面衛(wèi)星運(yùn)轉(zhuǎn)問題、經(jīng)濟(jì)方面物品供給以及需求與物價(jià)的之間的關(guān)系、軍事方面研究深水炸彈在水下的運(yùn)動(dòng)等。對(duì)這些事物、現(xiàn)象變化規(guī)律的描述、認(rèn)知和分析,需要運(yùn)用常微分方程來解決。人們使用常微分方程數(shù)值解法的四階龍格-庫塔方法去研究這些問題很實(shí)用,而且具有很重要的應(yīng)用價(jià)值。

        目前,常微分方程在解決我們生活中的問題很實(shí)用,許多問題都運(yùn)用常微分方程來求解。中國科學(xué)技術(shù)大學(xué)學(xué)者倪興在常微分方程的研究中寫了關(guān)于歐拉法、方法等幾種方法,他運(yùn)用常微分計(jì)算衛(wèi)星運(yùn)動(dòng)的初軌,把方法運(yùn)用到衛(wèi)星軌道改進(jìn)的例子中;揚(yáng)州大學(xué)學(xué)者馮建強(qiáng)和孫詩一研究四階方法的推導(dǎo),他寫出了如何推導(dǎo)的過程。在高校數(shù)值分析、數(shù)值計(jì)算方法與實(shí)驗(yàn)等教材中,許多作者都出版關(guān)于常微分方程初值問題數(shù)值解法的教材書,歐拉方法、改進(jìn)歐拉法和方法等,同時(shí)在教材書中寫入各種實(shí)際問題的例子,運(yùn)用這些方法去解決常微分方程的初值問題。本文主要介紹常微分方程數(shù)值解法的四階方法,使用四階方法求解常微分方程。介紹四階方法的基本思想、算法步驟和MATLAB程序,同時(shí)使用四階方法在案例中解決問題。

        1四階方法

        1.1四階方法的思想與算法

        對(duì)于一階常微分方程初值問題,根據(jù)拉格朗日微分中值定理可知存在使

        記為在區(qū)間上的平均斜率。

        在給定的區(qū)間里預(yù)測(cè)個(gè)點(diǎn)處的斜率,再把它們的加權(quán)平均作為平均斜率的近似值。這里取在區(qū)間上若干個(gè)點(diǎn)的斜率值或者預(yù)報(bào)斜率值的加權(quán)平均值,作為平均斜率的近似值,令區(qū)間上若干個(gè)點(diǎn)的斜率值或預(yù)報(bào)斜率值為,以及權(quán)系數(shù)為,使差分格式

        當(dāng)上式為階時(shí),把它叫做階龍格—庫塔格式.通過上式能夠推導(dǎo)得精確度比較高的公式來求解常微分方程。

        當(dāng)時(shí),龍格-庫塔的格式為

        在上式中,在點(diǎn)的斜率是在點(diǎn)的預(yù)報(bào)斜率為是,參數(shù)都是待定參數(shù)。當(dāng)待定參數(shù)使上式是二階格式時(shí),稱上式是二階龍格-庫塔格式。

        函數(shù)在點(diǎn)的泰勒展開式為

        令是準(zhǔn)確的,有,根據(jù)二元函數(shù)泰勒展開有

        要使上面的截距誤差為,把與泰勒展開式比較可得。通過以上的推導(dǎo)可以得出待定參數(shù)的值。

        由可知二階龍格-庫塔格式是一個(gè)系列的差分格式。當(dāng)取時(shí)得

        這是改進(jìn)的歐拉公式。取,可得

        其中

        這是常用的二階龍格-庫塔方法,稱其為中心格式。

        三階龍格-庫塔格式:

        三階龍格-庫塔格式是比二階格式高一階的格式,它是在二階龍格-庫塔格式的基礎(chǔ)上進(jìn)一步構(gòu)造得到的格式。當(dāng)時(shí),龍格-庫塔的格式如下

        在上式里,在點(diǎn)的斜率是在點(diǎn)和點(diǎn)的預(yù)報(bào)斜率分別是、,當(dāng)待定參數(shù)和使上式為三階格式時(shí),則稱其為三階龍格—庫塔格式。

        類似二階龍格-庫塔格式的推導(dǎo),得到三階龍格-庫塔格式為:

        稱其為三階龍格-庫塔格式。

        四階龍格-庫塔格式:

        同樣是根據(jù)上述的推導(dǎo),得到四階龍格-庫塔格式為:

        稱上式為四階龍格-庫塔格式。在解決許多問題中常常使用的龍格-庫塔格式是四階龍格-庫塔格式。

        四階龍格-庫塔方法的算法步驟如下:

        第一步輸入?yún)^(qū)間等分?jǐn)?shù),初值;

        第二步輸出在的個(gè)點(diǎn)處的近似值;

        第三步令;

        第四步計(jì)算

        令,輸出;

        第五步如果,令,轉(zhuǎn)向第四步;否則停機(jī)。

        1.2四階方法的MATLAB程序

        四階龍格-庫塔法MATLAB函數(shù)文件nark4.m如下:

        function [x,y]=nark4(dyfun,xspan,y0,h)

        x=xspan(1):h:xspan(2);

        y(1)=y0;

        for n=1:length(x)-1

        k1=dyfun(x(n),y(n));

        K2=dyfun(x(n)+h/2,y(n)+h/2*k1);

        K3=dyfun(x(n)+h/2,y(n)+h/2*k2);

        k4=dyfun(x(n+1),y(n)+h*k3);

        y(n+1)=y(n)+h*(k1+2*k2+2*k3+k4)/6

        end

        x=x';y=y';

        1.3四階方法的應(yīng)用

        例2.1:用四階方法求解初值問題的數(shù)值解,取步長(zhǎng)。其精確解為,并將解得的數(shù)值解與精確解進(jìn)行比較。

        猜你喜歡
        龍格常微分方程
        磁懸浮式雙自由度軌道車輛軸箱振動(dòng)能量采集器研究
        懸索橋吊索尾流致振非定常理論分析
        基于Labview的龍格—庫塔法估算Pu內(nèi)照射攝入量
        非結(jié)構(gòu)網(wǎng)格的有限體積法研究
        淺談常微分方程與大學(xué)生創(chuàng)新培養(yǎng)的聯(lián)系
        東方教育(2016年18期)2017-01-16 11:17:03
        基于常微分方程的數(shù)學(xué)建模問題的求解
        常微分方程的教學(xué)探討
        考試周刊(2016年82期)2016-11-01 12:54:50
        數(shù)學(xué)建模在高數(shù)教學(xué)中的應(yīng)用
        常微分方程新的可解類型及在經(jīng)濟(jì)學(xué)中的應(yīng)用
        商(2016年22期)2016-07-08 15:04:37
        常微分方程的思想方法及對(duì)經(jīng)濟(jì)管理的影響
        商(2016年22期)2016-07-08 14:57:18
        久久午夜无码鲁丝片直播午夜精品| 人人爽人人澡人人人妻| 色婷婷五月综合亚洲小说| 真实国产网爆门事件在线观看| 成人亚洲av网站在线看| 制服丝袜一区二区三区| 国产一卡2卡3卡四卡国色天香| 一本久道久久综合五月丁香| 日韩精品一级在线视频| 亚洲乱码av中文一区二区| 最新亚洲人成网站在线观看| 亚洲中文无码久久精品1| 亚洲中文字幕不卡一区二区三区 | 中文字幕丰满伦子无码| 日韩精品电影在线观看| 在线视频日韩精品三区| 女优av一区二区三区| 日本熟妇色xxxxx欧美老妇| 国产一区二区精品久久凹凸| 中文字幕文字幕视频在线| 精品av熟女一区二区偷窥海滩| 狠狠躁夜夜躁人人爽天天不卡软件| 男人天堂AV在线麻豆| 免费国产一区二区视频| 国产精品久久久国产盗摄| 国产小视频网址| 亚洲一区二区三区av无| 妺妺窝人体色www婷婷| 亚洲成a人v欧美综合天堂麻豆| 国产极品美女到高潮视频| 91九色最新国产在线观看| 无码乱人伦一区二区亚洲一| 国产91吞精一区二区三区 | 亚洲九九九| 亚洲丰满熟女乱一区二区三区| 国产精品久线在线观看| 国产啪精品视频网给免丝袜| 日本人妻三级在线观看| 伊人大杳焦在线| 熟妇人妻AV中文字幕老熟妇| 国产熟女乱综合一区二区三区 |