数学建模数值分析模型.pptxVIP

数学建模数值分析模型会计学 黑龙江科技学院 数 学 建 模第六章数值分析模型插值法非线性方程求根第2页/共61页迭代法弦截法和抛物线法建模举例重点:插值法和非线性方程求根难点:利用数值分析方法建立数学模型 理学院 黑龙江科技学院 数 学 建 模第3页/共61页 数值分析（numerical analysis）是研究用计算机求解各种数学计算问题的数值计算方法及其理论与软件实现的学科。数值分析就是介绍如何用计算机来解决数学问题，以各种各样的程序语言来设计出数值计算程序，然后依靠计算机的强大计算能力来求解这些数学问题，数值分析对数学理论与程序设计并重。 运用数值分析解决问题的过程可分为如下几步：实际问题→数学模型→数值计算方法→程序设计→上机计算求出结果。 数值分析这门学科有如下特点：（1）面向计算机（2）有可靠的理论分析（3）要有好的计算复杂性（4）要有数值实验（5）要对算法进行误差分析 理学院 黑龙江科技学院 数 学 建 模第4页/共61页函数逼近问题设y = f(x)，若对以函数y = f(x)来说 ① 其值是通过实验或观测得到，不知其解析表达式；② 解析表达式很复杂，不便分析。问题：能否构造一个较为简单的函数P(x)近似地表示f(x)。这就是函数逼近问题。上述函数f(x)称为被逼近函数，P(x)称为逼近函数。逼近方式有两种：插值和拟合。 理学院----被插值函数-----插值条件----插值函数----插值节点第5页/共61页 6.1 插值法在生产和科学研究中,经常出现这样的问题:由实验或测量得到的某一函数 在一系列点 处的值,需要构造一个简单函数 作为函数 的近似表达式: ,使得 这类问题称为插值问题.[min｛xi｝,max{xi} ]= [a,b]----插值区间插值函数：有各种类型，如代数多项式，三角函数，有理函数等。当插值函数为多项式时，称为（代数）插值多项式。 黑龙江科技学院 数 学 建 模第6页/共61页从几何上看，插值法就是要求一条曲线 它通过已知的n+1个点(xi,yi)(i=0,1, …,n)，并用 近似表示 f(x).（下图）y?y＝f(x)??yi??y0xx0xio 理学院 黑龙江科技学院 数 学 建 模第7页/共61页 一、插值基函数与Lagrange插值1. 简单情形 (1) n = 1时. 设 yi = f(xi) i = 0，1.作直线方程： 令：称 为两点式插值或线性插值。 理学院 黑龙江科技学院 数 学 建 模第8页/共61页 (2) n = 2时. 设yi = f(xi)i = 0，1，2. 令：称 为三点式插值或抛物插值。 理学院 黑龙江科技学院 数 学 建 模第9页/共61页2. 推广 n = 1时，记则 n=2时，记则 理学院 黑龙江科技学院 数 学 建 模第10页/共61页一般地令 则lj(x) (j = 0，1，2，…，n)为n次多项式称为Lagrange插值基函数， 为Lagrange插值多项式。 理学院 黑龙江科技学院 数 学 建 模7576777879902.7682.8332.9032.9793.0623.153（1）作一分段线性插值函数（2）用上述插值函数计算和的函数值。第11页/共61页例 给定数组 理学院 黑龙江科技学院 数 学 建 模解 由插值基函数的表达式，在75到80的6个点间有5个线性插值函数，以区间为例，此时则在区间上有第12页/共61页. 理学院 黑龙江科技学院 数 学 建 模第13页/共61页Matlab代码如下:function [Y,Phi]=FenDuanXianXingChaZhi(xx)clcx1=75:80;y=[2.768,2.833,2.903,2.979,3.062,3.153]; n=size(x1,2); syms x positivefor i=1:(n-1)Phi(i)=y(i)*(x-x1(i+1))/(x1(i)-x1(i+1))+y(i+1)*(x-x1(i))/(x1(i+1)-x1(i));endPhi=Phi; l=find(x1xx);Y=subs(Phi(l(1)-1),xx); end 理学院 黑龙江科技学院 数 学 建 模第14页/共61页函数的调用格式为xx=75.5[Y,Phi]=FenDuanXianXingChaZhi(xx)得到的结果为：Y =2.8005Phi =(13*x)/200 - 2107/1000 (7*x)/100 - 2487/1000 (19*x)/250 - 2949/1000 (83