第六章 计算方法  非线性方程求解 多项式插值与曲线拟合 数值微分与数值积分 求常微分方程数值解命令.

Slides:



Advertisements
Similar presentations
数学建模计算 1 基于 MATLAB 的数学建模竞赛计算 计算在建模竞赛中的作用 数学建模竞赛中的数学软件 MATLAB 数学建模工具箱 数学建模 MATLAB 命令及建模应用.
Advertisements

数值分析 第五节 数值微分 在实际问题中,往往会遇到某函数 f(x) 是用表格 表示的, 用通常的导数定义无法求导, 因此要寻求其他 方法近似求导。常用的数值微分方法有 : 一. 运用差商求数值微分 二.运用插值函数求数值微分 三. 运用样条插值函数求数值微分 四. 运用数值积分求数值微分.
一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
第五节 全微分方程 一、全微分方程及其求法 二、积分因子法 三、一阶微分方程小结. 例如 所以是全微分方程. 定义 : 则 若有全微分形式 一、全微分方程及其求法.
第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
第九章 常微分方程数值解法 §1 、引言. 微分方程的数值解:设方程问题的解 y(x) 的存在区间是 [a,b] ,令 a= x 0 < x 1
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
高等数学一 主讲 杨俊 演示文稿制作 杨俊. 高等数学一 第 3 章 一元函数微分学的应用 第 4 章 一元函数 积分学及应用 第 1 章 函数、极限与连续 第 2 章 导数与微分.
第 4 章 数值微积分. 4.1 内插求积 Newton-Cotes 公式 第 4 章 数值微积分 4.1 内插求积 Newton-Cotes 公式.
一、会求多元复合函数一阶偏导数 多元复合函数的求导公式 学习要求: 二、了解全微分形式的不变性.
计算机数学基础(下) --数值分析 教师:孙继荣 电话: 028 -
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
1 、牛顿 - 莱布尼兹公式 另外若给出的函数 f(x) 是数据表,也不好求函数的积分。 计算定积分的方法: 但是求函数 f(x) 的原函数 F(x) 不一定比计算积分容易, 例如函数 找不到用初等函数表示的原函数。 一、数值求积的基本思想 实验 4 数值积分与微分 主讲人:魏志强.
第二章 导数与微分. 二、 微分的几何意义 三、微分在近似计算中的应用 一、 微分的定义 2.3 微 分.
全微分 教学目的:全微分的有关概念和意义 教学重点:全微分的计算和应用 教学难点:全微分应用于近似计算.
第五篇 社會變遷與社會問題 單元 16 人口變遷. 本課綱要 人口變遷的意涵 人口變遷的理論 世界人口變遷 世界人口問題及對策 社會變遷的適應 臺灣人口問題與政策 考題大觀園.
第6章 利用MATLAB语言 求解科学运算问题
一、能线性化的多元非线性回归 二、多元多项式回归(线性化)
一、二阶行列式的引入 用消元法解二元线性方程组. 一、二阶行列式的引入 用消元法解二元线性方程组.
第三章 函数逼近 — 最佳平方逼近.
例题 教学目的: 微积分基本公式 教学重点: 牛顿----莱布尼兹公式 教学难点: 变上限积分的性质与应用.
第3章 MATLAB数值计算 2017/9/9.
数 学 分 析 第九章 定积分 第二节 微积分学基本公式 主讲:师建国.
第四章 函数的积分学 第六节 微积分的基本公式 一、变上限定积分 二、微积分的基本公式.
§5.3 定积分的换元法 和分部积分法 一、 定积分的换元法 二、 定积分的分部积分法 三、 小结、作业.
第四章 一元函数的积分 §4.1 不定积分的概念与性质 §4.2 换元积分法 §4.3 分部积分法 §4.4 有理函数的积分
第四章 数值积分与数值微分 — 基本概念 — Newton-Cotes 公式.
定积分习题课.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
第三章 导数与微分 习 题 课 主要内容 典型例题.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
体育选项课件 健美操理论课 任课教师:黄明礼 湄洲湾职业技术学院.
第5章 §5.3 定积分的积分法 换元积分法 不定积分 分部积分法 换元积分法 定积分 分部积分法.
西南科技大学网络教育系列课程 数学软件 数学软件 第5讲 MATLAB数值计算二 主讲教师: 鲜大权 副教授 西南科技大学理学院数学系.
计算机数学基础 主讲老师: 邓辉文.
数学模型实验课(三) 插值与三维图形.
MATLAB及其应用 第三讲 数据处理 授课人:鲍文 在此幻灯片插入公司的徽标 从“插入”菜单 选择图片 找到徽标文件 单击“确定”
Matlab 选讲 二 上海交通大学数学系 刘小军
第三单元 第4课 Matlab数据插值 1.一维插值 2.二维插值 3.对非网格数据进行插值.
一.多项式构造及其运算 1、多项式构造 poly2str(p,’x’) 将表示多项式系数的行向量p转换为变量是x的多项式形式。
插值与拟合 一、插值的基本原理 二、拟合的基本原理 三、插值与拟合的关系 四、插值的MATLAB实现 五、拟合的Matlab实现.
导数的应用 ——函数的单调性与极值.
第二十二章 曲面积分 §1 第一型曲面积分 §2 第二型曲面积分 §3 高斯公式与斯托克斯公式.
第三单元 第3课 实验 多元函数的积分 实验目的:掌握matlab计算二重积分与三重积分的方法,提高应用重积分解决有关应用问题的能力。
概 率 统 计 主讲教师 叶宏 山东大学数学院.
实验一 计算复变函数极限、微分、积分、 留数、泰勒级数展开式 (一) 实验类型:验证性 (二) 实验类别:基础实验
作业 P158 习题 2 1(2)(4) (5). 2(1). 预习 P156— /5/2.
作业 P152 习题 复习:P 预习:P /5/2.
第三单元 第2课 实验 一元函数的积分 实验目的:掌握matlab求解有关不定积分和定积分的问题,深入理解定积分的概念和几何意义。
正切函数的图象和性质 周期函数定义: 一般地,对于函数 (x),如果存在一个非零常数T,使得当x取定义域内的每一个值时,都有
3.1 变化率与导数   3.1.1 变化率问题 3.1.2 导数的概念.
概 率 统 计 主讲教师 叶宏 山东大学数学院.
第三章 函数的微分学 第二节 导数的四则运算法则 一、导数的四则运算 二、偏导数的求法.
第 四 章 迴歸分析應注意之事項.
一 般 的 代 数 方 程 函数solve用于求解一般代数方程的根,假定S为符号表达式,命令solve (S)求解表达式等于0的根,也可以再输入一个参数指定未知数。例: syms a b c x S=a*x^2+b*x+c; solve(S) ans = [ 1/2/a*(-b+(b^2-4*a*c)^(1/2))]
建模常见问题MATLAB求解  .
第二章 函 数 插 值 — 三次样条插值.
正弦、余弦函数的性质 华容一中 伍立华 2017年2月24日.
第四节 第七章 一阶线性微分方程 一、一阶线性微分方程 *二、伯努利方程.
§2 方阵的特征值与特征向量.
教学大纲(甲型,54学时 ) 教学大纲(乙型, 36学时 )
Volterra-Lotka方程 1925年, A. Lotka(美)和V. Volterra(意)给出了第一个两物种间的捕食模型。
第三章 函数 逼近 — 曲线拟合的最小二乘法.
《偏微分方程》第一章 绪论 第一章 绪论 1.1.
百雞問題 製作者:張美玲 資料來源:數學誕生的故事—凡異出版社.
数学模型实验课(二) 最小二乘法与直线拟合.
西南科技大学网络教育系列课程 数学软件 数学软件 第7讲 MATLAB符号计算二 主讲教师: 鲜大权 副教授 西南科技大学理学院数学系.
§4.5 最大公因式的矩阵求法( Ⅱ ).
实验二 定积分的近似计算.
Matlab插值与拟合 插值 拟合.
Presentation transcript:

第六章 计算方法  非线性方程求解 多项式插值与曲线拟合 数值微分与数值积分 求常微分方程数值解命令

已知函数 f(x),求 x 使 f(x)=0 —— 非线性方程 例 x = tan x x y (k = 1,2,·······) f=inline('x-tan(x)') x=fzero(f,3.5) x = 4.4934 求一元函数零点 fzero(fun, x0) X = fzero(fun, x0)寻找fun在x0附近的零点

r x 例 x 3 – 30 x 2 + 2552 =0 p=[1 -30 0 2552]; roots(p) 半径r =10 cm的球体,密度 =0.638.浸入水深度 x=? ans = 26.3146 11.8615 -8.1761 15/18

求多项式零点: roots(p); 例2 求 x3 –8x2 +6x– 30=0的解 P=[1 -8 6 -30]; r=roots(P) r = 7.7260 0.1370 + 1.9658i 0.1370 - 1.9658i 例1.求方程 x3 – 4x + 5 = 0 的解 P=[1 0 -4 5]; roots(P) ans = -2.4567 1.2283 + 0.7256i 1.2283 - 0.7256i

P10(t) f(t) f(x) 多项式插值原理 已知 f(x) 在点 xi 上的函数值 yi=f(xi), (i=0,1,2,···,n) 求多项式 P(x)=a0 + a1x +···+ anxn 满足: P(xk)= yk (k = 0,1,…,n) P10(t) f(t) f(x)

zi=interp2(x, y, z, xi, yi, ‘method’) 一维插值: yi = interp1(x, y, xi, ‘method ’) x = 0:10; y = sin(x); xi = 0:.25:10; yi = interp1(x,y,xi); plot(x,y,'o',xi,yi) method nearest 最近点插值 linear 线性插值 spline 样条插值 cubic 立方插值 二维插值 zi=interp2(x, y, z, xi, yi, ‘method’) 三维插值 vi = interp3(x,y,z,v,xi,yi,zi, ‘method’) 17/18

例1.正弦曲线数据插值试验

离散数据的多项式拟合 x x1 x2 ·········· xm f(x) y1 y2 ·········· ym 使得 已知数据表 求拟合函数:  (x) = a0 + a1x + ······+ anxn 使得 达到最小

P(x)=a0xn + a1xn-1 + ······ + an-1x + an 多项式拟合命令 P=polyfit(X,Y,n) 求出(最小二乘意义下)n次拟合多项式 P(x)=a0xn + a1xn-1 + ······ + an-1x + an 计算结果为系数: P=[ a0, a1, ······ , an-1, an ] 多项式求值命令: y1=polyval(P, x) 其中,P是n次多项式的系数,x是自变量的值,y1是多项式在x处的值 17/18

定积分与积分和式 Sn 5.2908 5.1044 4.9835 ······ 4.8999 右矩形和 h 1 0.5 0.2 ······ 4/16

梯形数值积分公式 左矩形 梯形 右矩形 4.4429 4.8669 5.2908 4.6804 4.8924 5.1044 4.8139 4.8987 4.9835 4.8572 4.8996 4.9420 5/16

梯形数值积分命令 trapz() clear x=0:pi/100:pi; y=sin(x); trapz(x,y) ans = 1.9998 clear x=sort(rand(1,101)*pi); y=sin(x); trapz(x,y) ans = 1.9981 rand(1,101)产生101个均匀随机数,每个数都介于0-1之间

辛卜生求数值积分命令quad() clear fun=inline('1./(x.^3-2*x-5)') ezplot(fun,[0,2]) [q,n]=quad(fun,0,2) q = -0.4605 n = 53

二重积分dblquad()与三重积分 fun=inline('y*sin(x)+x*cos(y)') Q=dblquad(fun,pi,2*pi,0,pi) Q = -9.8698 [x,y]=meshgrid(pi:.1:2*pi,0:.1:pi); z=fun(x,y); mesh(x,y,z)

物理量在空间的分布情况其及随时间变化的规律。 研究具有实际应用背景的数学函数 物理量在空间的分布情况其及随时间变化的规律。 物理现象  物理定律  微分方程  数值解 建模型: 物理学家、化学家、生物学家、社会学家 2/18

 y(t)=y0 e r t 马尔萨斯(Malthus)人口模型 某时刻人口数为y0, 在以后某一时间 t 人口数y 的变化率正比于y Thomas Malthus (1766-1834) MATLAB符号求解方法 y=dsolve('Dy=r*y','t') y =C1*exp(r*t) 弹性力学问题 单摆问题 天体力学二体问题 电磁场计算············ 令 r = a, y0 = exp( b )  y(t) = exp( at + b) a = ?, b = ?  人口数据拟合实验 3/18

MATLAB求常微分方程数值解 简例 f=inline('y-x.*y.^2'); [x,y]=ode23(f,[0,2],1) 0 1.0000 0.0800 1.0796 0.2621 1.2485 0.4437 1.3755 0.6437 1.4402 0.8437 1.4209 1.0437 1.3374 1.2437 1.2196 1.4437 1.0923 1.6437 0.9709 1.8437 0.8620 2.0000 0.7870 f=inline('y-x.*y.^2'); [x,y]=ode23(f,[0,2],1) 4/18

ode23(),ode45() 求常微分方程初值问题数值解函数 解一阶常微分方程组初值问题命令格式 [T,y] = ode23(' F ',Tspan,y0) 其中, F是ODE文件, 表示 y' = F(t,y)中右端函数 Tspan = [t0 Tfinal] —— 求解区域; y0 —— 初始条件 注: 函数F(t,y) 必须返回列向量. 数值解 y 的每一行对应于列向量T中的每一时刻 5/18

捕食者和被捕食者模型 ( 0≤ t ≤ 15 ) y1(0)=20, y2(0)=20 数会急增,数学家将鱼分成捕食者和被捕食者,分别 生物学家研究发现,战争期间捕鱼量减少,食肉鱼总 数会急增,数学家将鱼分成捕食者和被捕食者,分别 以y2(t) , y1(t)表示它们在时刻t 的总数。假定被捕食者 (食用鱼)有充分的食物。 ( 0≤ t ≤ 15 ) y1(0)=20, y2(0)=20 6/18

function yp = lotka(t,y) % predator-prey model. MATLAB求解: (1)定义函数; (2)调用ode23 函数文件,文件名:lotka.m function yp = lotka(t,y) % predator-prey model. yp=diag([1 - .01*y(2), -1+.02*y(1)])*y; t0=0;tf=15; y0=[20 20]; [t,y]=ode23(‘lotka’,[t0,tf],y0); plot(t,y(:,1),t,y(:,2),'r') —— 捕食者y2 —— 被捕食者y1 7/18

洛伦兹模型与“蝴蝶效应” x(0)=0,y(0)=0,z(0)=0.01 function z=flo(t,y) Y-X平面 x(0)=0,y(0)=0,z(0)=0.01  = 8/3  = 10  =28 function z=flo(t,y) A=[-8/3 0 y(2); 0 -10 10; -y(2) 28 -1]; z=A*y; y0=[0;0;0.01]; [x,y]=ode23(‘flo’,[0, 80],y0); plot(y(:,2),y(:,1)) 8/18