MATLAB在常微分方程上的應用 楊惠如 老師:王天楷教授 2005/8/30.

Slides:



Advertisements
Similar presentations
版 画 制 作版 画 制 作 版 画 种 类版 画 种 类 版 画 作 品版 画 作 品 刘承川.
Advertisements

Differentiation 微分 之二 以公式法求函數的微分. Type 函數形式 Function f (x) Derivative d f (x) /d x c=constant 常數 c0 Power of x xaxa a x a-1 Trigonometric 三角函數 sin x cos.
第五节 全微分方程 一、全微分方程及其求法 二、积分因子法 三、一阶微分方程小结. 例如 所以是全微分方程. 定义 : 则 若有全微分形式 一、全微分方程及其求法.
Chap 3 微分的應用. 第三章 3.1 區間上的極值 3.2 Rolle 定理和均值定理 3.3 函數的遞增遞減以及一階導數的判定 3.4 凹面性和二階導數判定 3.5 無限遠處的極限 3.6 曲線繪圖概要 3.7 最佳化的問題 3.8 牛頓法 3.9 微分.
第九章 常微分方程数值解法 §1 、引言. 微分方程的数值解:设方程问题的解 y(x) 的存在区间是 [a,b] ,令 a= x 0 < x 1
§ 2.2 线性微分方程与常数变易法 /Linear ODE and variation of constants Method/
663 Chapter 14 Integral Transform Method Integral transform 可以表示成如下的積分式的 transform  kernel Laplace transform is one of the integral transform 本章討論的 integral.
大家好,我是 NTOU/MSV 工數教學的 York 老師,今天跟大家講解工程數學的第 一課,一階常微分方程。 一般我們講一階 ODE , first order ordinary differential equation 。這邊的話我 們建議大家要先修過微積分,才能懂得我們要講甚麼。另外也有 PDE.
计算机数学基础(下) --数值分析 教师:孙继荣 电话: 028 -
CH2: 微分學 切 The definition of derivatives CH2: 微分學 Step1 :
工職數學 第四冊 第一章 導 數 1 - 1 函數的極限與連續 1 - 2 導數及其基本性質 1 - 3 微分公式 1 - 4 高階導函數.
实验目的 实验内容 MATLAB 2 、学会用 Matlab 求微分方程的数值解. 实验软件 1 、学会用 Matlab 求简单微分方程的解析解. 1 、求简单微分方程的解析解. 2 、求微分方程的数值解.
不定積分 不定積分的概念 不定積分的定義 16 不定積分的概念 16.1 不定積分的概念 以下是一些常用的積分公式。
大綱 1. 三角函數的導函數. 2. 反三角函數的導函數. 3. 對數函數的導函數. 4. 指數函數的導函數.
2-1 極限的概念 2-2 無窮等比級數 2-3 多項式函數的導數導函數 2-4 微分公式 2-5 微分的應用 2-6 積分的概念與反導函數 信樺文化.
第一單元 建立java 程式.
第六章 6.1 常微分方程及其求解概述 6.2 初值问题解法 6.3 边值问题解法
第九章 常微分方程的数值解法 主 要 内 容 §1、引言 §2、初值问题的数值解法--单步法 §3、龙格-库塔方法 §4、收敛性与稳定性
坚持唯物辩证法 反对形而上学 轻松学习.
色彩性格 -----发挥你的最佳本色 我相信生命是最富有情趣的旅程。 他可以成为你梦想不到、更美妙的旅程;
岡山區103年第12次 登革熱聯繫會報會議 岡山區公所 103年12月30日 1.
遞迴關係-爬樓梯.
第五章 法律价值 第一节 法律价值 一、价值和法律价值概念
运用Matlab GUI辅助大学物理实验 蒋志洁 中山大学 物理学院
增值评价 2014级 初中起点报告 解读培训 辽宁省基础教育质量监测与评价中心.
Ordinary Differential Equations
计算机数学基础(下) 第5编 数值分析 第14章 常微分方程的数值解法.
数值计算方法 第八章 常微分方程初值问题数值解法  重庆邮电大学.
Differentiation 微分 之一 微分的基本原理.
本章大綱 9.1 Sequence數列 9.2 Infinite Series無窮級數
MPLAB IDE 8 建立專案.
均勻降雨時坡地之水理分析 國立中興大學水土保持學系碩士班 Hydraulic analysis of flow down a slope under a uniform rainfall 指導教授:謝平城 教授 學 生:楊佳錫 日 期:
Differentiation 微分 之一 微分的基本原理.
Methods of Integration 積分的方法
数学软件 Matlab —— Matlab 符号运算.
内容: 1. 库模块简介 2.基本建模方法 3.模型举例 4.子系统与模块封装技术 5.函数的编写与应用
Simulink模擬基礎 主要內容 Simulink簡介 Simulink模組庫 Simulink的基本操作 S-函數.
線性一階微分方程與尤拉法 線性一階微分方程式求解 (Linear First-Order Differential Equations)
Methods 靜宜大學資工系 蔡奇偉副教授 ©2011.
CH5、SIMULINK仿真基础 在工程实际中,控制系统的结构往往很复杂,如果不借助专用的系统建模软件,则很难准确地把一个控制系统的复杂模型输入计算机,对其进行进一步的分析与仿真。 1990年,Math Works软件公司为MATLAB提供了新的控制系统模型图输入与仿真工具,并命名为SIMULAB,该工具很快就在控制工程界获得了广泛的认可,使得仿真软件进入了模型化图形组态阶段。但因其名字与当时比较著名的软件SIMULA类似,所以1992年正式将该软件更名为SIMULINK。
檔案與磁碟的基本介紹.
第二章 SPSS的使用 2.1 啟動SPSS系統 2.2 結束SPSS系統 2.3 資料分析之相關檔案 2.4 如何使用SPSS軟體.
第4章 非线性规划 基本概念 2011年11月.
MATLAB及其应用 第三讲 数据处理 授课人:鲍文 在此幻灯片插入公司的徽标 从“插入”菜单 选择图片 找到徽标文件 单击“确定”
用相频曲线测阻尼系数的探索 指导教师 陈乾 吉新程.
第四节 函数展开成幂级数 本节内容: 一、泰勒 ( Taylor ) 级数 二、函数展开成幂级数 第十二章 两类问题: 在收敛域内 求 和
虎克定律與簡諧運動 教師:鄒春旺 日期:2007/10/8
第一單元 建立java 程式.
Matlab及其应用 讲座之五 动态系统仿真——Simulink 主讲人:鲍文 在此幻灯片插入公司的徽标 从“插入”菜单 选择图片
第14章 總體經濟政策之爭論:法則與權衡性.
人工智慧應用蟲體辨識 亞洲大學 資訊工程學系 學生:殷聖展、廖哲毅、吳京育 指導教授:蔡志仁 教授
引用檔案.
貝氏刷牙法 (Bass Method) 外埔國小.
均勻降雨時坡地之 水流動力分析 指導教授:謝平城 教授 授課教授:林俐玲 教授 學 生:楊佳錫.
Video 影像 (VideoPlayer 影像播放器、Camcorder 錄影機) 靜宜大學資管系 楊子青
物理化學輔助學習工具 2018/12/04.
工程數學 Chapter 15 Power Series , Taylor Series 楊學成 老師.
 隐式欧拉法 /* implicit Euler method */
班級:博碩子一甲 授課老師:鐘國家 助教:陳國政
第八章 服務部門成本分攤.
第一章 直角坐標系 1-3 函數及其圖形.
補充 數值方法 數值方法.
幂函数.
What is “this”? 在物件導向程式設計中,類別的定義就是在說明如果創建了“這個物件”的話,它會具有那些屬性與功能,以及這些功能是如何實現的。 而所謂的“這個物件”就以 this 來表示。 當我們在JavaScript與jQuery中寫 script 程式(函式)時,“誰”呼叫這個函式,這個“誰”就是該函式中所謂的.
教学大纲(甲型,54学时 ) 教学大纲(乙型, 36学时 )
教学大纲(甲型,54学时 ) 教学大纲(乙型, 36学时 )
Chapter 4 Multi-Threads (多執行緒).
9.5 函数的幂级数展开式 通过上节的学习知道:任何一个幂级数在其收敛区间 内,均可表示成一个函数(即和函数).但在实际中为了便于
張智星 台大資工系 MIR實驗室 第23章 程式碼的重複使用 張智星 台大資工系 MIR實驗室.
物理化學輔助學習工具 2018/12/04.
Presentation transcript:

MATLAB在常微分方程上的應用 楊惠如 老師:王天楷教授 2005/8/30

數值分析法 Taylor series Runge-Kutta Method

y後=y前+(x後-x前)y’前+ y’’前 Taylor’s series 一階 二階 Runge-Kutta Method A first-order Runge-Kutta method: uses a first-order Taylor series expansion 一階 y後=y前+(x後-x前)y’前 二階 y後=y前+(x後-x前)y’前+ y’’前

MATLAB在常微分方程上的應用 函數解法的特點 M檔案 ODE範例

函數解法的特點 ode45 ode23 ode113 ode23t ode23s ode15s ode23tb 4階/5階Rung-Kutta(4,5)法。採單步法,速度快,作為解問題的初試法 ode23 2階/3階Rung-Kutta(2,3)法。對於誤差容許範圍較寬且稍帶剛性問題比ode45效果好 ode113 採可變階的Adams-Bashforth-Moulton PECE法,屬於多步解法,比ode45更適合於誤差容許範圍要求較嚴格時 ode23t 採自由內插法實現的梯形規則,適用剛性並要求無數值衰減的問題 ode23s 採用改進的2階Rosenbrock法,採單步法,比ode15s更適於誤差範圍容許較寬的情況。 ode15s 可變階的數值微分公式法(NDFs),屬於多步解法,如果用在ode45效果很差或失敗,可考慮用此方式。可解剛性問題。 ode23tb 採用TR-BDF2方法,和ode23s一樣,對於誤差容許範圍較寬的情況,比ode15s效果好,可解剛性問題 資料來源:蔡宏睿

M檔案 以MATLAB程式碼所撰寫的檔案通常以 “ m ” 為副檔名,所以這些檔案又稱為M檔案 M檔案又可分為兩類—”函式” 和”主程式” 資料來源:蔡宏睿

ODE範例 y’=3x2 函式 function dy=ff(x,y) dy=3*x.^2 執行 ode45(‘ff’,[2 7],0.5) 檔名 y的初值 X的範圍

ODE範例 函式 function yprime=cc(x,y) yprime=[y(2);-sin(y(1))] 執行 ode45(‘cc’,[0 10],[1;1.5])

ODE範例 y1(0)=1 y2(0)=1 y3(0)=1 函式 function yprime=bb(t,y,a,b,c) yprime=[-y(2)-y(3);y(1)+a*y(2);b+y(3)*(y(1)-c)] options=odeset(‘AbsTol’,1e-7,’Reltol’,1e-4) 絕對誤差 相對誤差 1·10-4 執行 ode45(‘bb’,[0 100],[1;1;1],options,a,b,c)

bb.m

報告完畢 謝 謝