最簡單求解一個微分方程數(shù)值解得方法：Euler法

function [x,y]=Euler_method(dufun,span,h,x0,y0)

%EuLer格式,

%求解方程y'=dufun(x,y);其中x \in[a,b];y0為初始值；n為自變量的離散個數(shù)；y為求解結(jié)果

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

n=length(x);

y=zeros(1,n);%存放數(shù)值的解

x(1)=x0;

y(1)=y0;

for i=1:n-1

y(i+1)=y(i)+h.*feval(dufun,x(i),y(i));%Euler格式

end

２．改進的Euler格式

function[x,y]=Gaijin_Euler(func,span,h,x0,y0)

%改進的Euler格式；

%span為區(qū)間

%y(x0)=y0;

%n表示區(qū)間等份；

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

n=length(x);

y=zeros(1,n);

x(1)=x0;

y(1)=y0;

for i=1:n

y1=y(i)+h*feval(func,x(i),y(i));

y2=y(i)+h*feval(func,x(i+1),y1);

y(i)=(y1+y2)/2;

end

方法二：常見的Rungerkutta3與Rungerkutta4格式

１.Rungerkutta3

function [x,y]=Rungekutta3(f1,span,h,x0,y0)

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

n=length(x);

y=zeros(1,n);

x(1)=x0;

y(1)=y0;

for j=1:n-1

K1=feval(f1,x(j),y(j));

K2=feval(f1,x(j)+h/2,y(j)+h/2*K1);

K3=feval(f1,x(j)+h,y(j)-h*K1+h*2*K2);

y(j+1)=y(j)+(h/6)*(K1+4*K2+K3);

end

１.Rungerkutta4

function [x,y]=Rungekutta4(f1,span,h,x0,y0)

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

n=length(x);

y=zeros(1,n);

x(1)=x0;

y(1)=y0;

for j=1:n-1

K1=feval(f1,x(j),y(j));

K2=feval(f1,x(j)+h/2,y(j)+h/2*K1);

K3=feval(f1,x(j)+h/2,y(j)+h/2*K2);

K4=feval(f1,x(j)+h,y(j)+h*K3);

y(j+1)=y(j)+(h/6)*(K1+2*K2+2*K3+K4);

end

調(diào)試和運行matlab代碼(源程序)的技巧和教程

轉(zhuǎn)載請標明出處:專注matlab代碼下載的網(wǎng)站http://www.downma.com/ 本文主要給大家分享使用matlab編寫代碼,完成課程設計.畢業(yè)設計或者研究項目時,matlab調(diào)試程序的技巧 ...

如何加速MATLAB代碼運行

學習筆記 V1.0 2015/4/17 如何加速MATLAB代碼運行 概述 本文源于LDPCC的MATLAB代碼,即.由于代碼的問題,在信息位長度很長 ...

多分類問題中，實現(xiàn)不同分類區(qū)域顏色填充的MATLAB代碼(demo：Random Forest)

之前建立了一個SVM-based Ordinal regression模型,一種特殊的多分類模型,就想通過可視化的方式展示模型分類的效果,對各個分類區(qū)域用不同顏色表示.可是,也看了很多代碼,但基本都是 ...

卷積相關公式的matlab代碼

取半徑=3 用matlab代碼實現(xiàn)上式公式: length=3;for Ki = 1:length for Kj = 1:length for Kk = 1:length Ksigma(Ki,Kj,K ...

JAVA調(diào)用matlab代碼

做實驗一直用的matlab代碼,需要嵌入到java項目中,matlab代碼拼拼湊湊不是很了解,投機取巧采用java調(diào)用matlab的方式解決. 1. ? ?matlab版本:matlabR2014a ...

直方圖均衡化與Matlab代碼實現(xiàn)

昨天說了,今天要好好的來解釋說明一下直方圖均衡化.并且通過不調(diào)用histeq函數(shù)來實現(xiàn)直方圖的均衡化. 一.直方圖均衡化概述 直方圖均衡化(Histogram Equalization)?又稱直方圖平 ...

將labelme 生成的.json文件進行可視化的代碼+label.png 對比度處理的matlab代碼

labelme_to_dataset 指令的代碼實現(xiàn): show.py文件 #!E:\Anaconda3\python.exe import argparse import json import o ...

關于本科畢業(yè)論文《Laguerre小波在數(shù)值積分與微分方程數(shù)值解中的應用》存在的問題與小結(jié)

本科的畢業(yè)設計是通過Laguerre小波函數(shù)來近似表達某個需要求積分或解微分方程的函數(shù),將原函數(shù)很難求得函數(shù)用小波函數(shù)表達出來,這 ...

SVM實例及Matlab代碼

******************************************************** ***數(shù)據(jù)集下載地址 :http://pan.baidu.com/s/1geb8CQf ...

隨機推薦

Enterprise Solution 企業(yè)資源計劃管理軟件 C/S架構，支持64位系統(tǒng)，企業(yè)全面應用集成，制造業(yè)信息化

Enterprise Solution是一套完整的企業(yè)資源計劃系統(tǒng),功能符合眾多制造業(yè)客戶要求.系統(tǒng)以.NET Framework技術作為開發(fā)架構,完善的功能可有效地幫助企業(yè)進行運營策劃,減低成本,如 ...

關于xml的使用。

使用的常用類: XmlSerializer ParaMapping StreamReader DirectionaryInfo FileInfo using as object 例子: public ...

radio被選中，但是重復點擊后事件不觸發(fā)

網(wǎng)上找了好多帖子,都沒用,在百度知道發(fā)現(xiàn)以下答案 知識點:使用 jq的prop才能設置 html $('.ss').click(function(){ $(this).find("input ...

Codeforces Round #249 (Div. 2)

A.水題. #include #include #include #include

分享一些Comet開發(fā)經(jīng)驗

前言 本comet技術主要用于數(shù)據(jù)庫持久層的 穿越防火墻 遠程訪問.只要有一臺中繼網(wǎng)站,任意地點的數(shù)據(jù)庫都能被訪問. Comet概念介紹 WebIM.網(wǎng)頁的客服.meebo等大家聽說過了.最近還有個兄 ...

字體圖標 icon font

Icon font icon font 指的是用字體文件代替圖片文件,來展示圖標.特殊字體等元素的方法. 應用場景: iconfont的優(yōu)缺點 大小能夠自由地變化 顏色能夠自由地改動 加入陰影效果 * ...

Mysql的Error 1364

在安裝mysql 5.5.8的最后一步,應用安全配置時, 彈出錯誤對話框,?Error Nr.1364 authentication_string 解決辦法: 1.暫時不要關閉該窗口,cmd命令,執(zhí)行 ...

杭電ACM水仙花數(shù)

水仙花數(shù) Time Limit: 2000/1000 MS (Java/Others)????Memory Limit: 65536/32768 K (Java/Others) Total Submi ...

Java后臺實現(xiàn)方法

Java后臺實現(xiàn)方法 首先后臺結(jié)構分為四個部分(以表schedule為例) entity>mapper>service>controller 1. 在entity里面寫好實體,新建目 ...

Python學習第四天

一.數(shù)字 int 二.字符串 str #以下均為補充內(nèi)容 #對于空字符串是假 #數(shù)字0是假 #數(shù)字和字符串可以相互轉(zhuǎn)換 # a="123" # b=int(a) # a=123 # ...

總結(jié)

以上是生活随笔為你收集整理的微分方程数值解 matlab,本学期微分方程数值解课程总结（matlab代码）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。