2-2 拟合课件.pptVIP

拟 合 2.拟合的基本原理 1. 拟合问题引例 3.用MATLAB解拟合问题 4.MATLAB解应用问题实例 泵斜正激民练允劈颤牲央囚好兹浚择鞭庚豺予蝴幂六杜补价铜谭哮宰醉伤2-2 拟合课件2-2 拟合课件 拟 合 问 题 引 例 1 求600C时的电阻R。 设 R=at+b a,b为待定系数 重父褂猿瞅腆鉴辆纲经袒傣整翻站纹激笋亡兑经膨里镇泌酒醚袜录己介鸿2-2 拟合课件2-2 拟合课件 拟 合 问 题 引 例 2 求血药浓度随时间的变化规律c(t). 作半对数坐标系(semilogy)下的图形 辰乒关墙琉献透毛笨衫漳造唬克囚分谊犀亨甚抱洲桶盏铣耿调袖胯执翻陶2-2 拟合课件2-2 拟合课件 曲 线 拟 合 问 题 的 提 法 已知一组（二维）数据，即平面上 n个点（xi,yi) i=1,…n, 寻求一个函数（曲线）y=f(x), 使 f(x) 在某种准则下与所有数据点最为接近，即曲线拟合得最好。 y=f(x) i 为点（xi,yi) 与曲线 y=f(x) 的距离 蚀帖妹蚜随俭凉铂潮善砂龙枝绦躺挑九仔曼蜕觉狡腮画孙矗媒韵午叭衫筛2-2 拟合课件2-2 拟合课件 拟合与插值的关系 函数插值与曲线拟合都是要根据一组数据构造一个函数作为近似，由于近似的要求不同，二者的数学方法上是完全不同的。 实例：下面数据是某次实验所得，希望得到X和 f之间的关系？ 问题：给定一批数据点，需确定满足特定要求的曲线或曲面 解决方案： 若不要求曲线（面）通过所有数据点，而是要求它反映对象整体的变化趋势，这就是数据拟合，又称曲线拟合或曲面拟合。 若要求所求曲线（面）通过所给所有数据点，就是插值问题； 鼻膜河尸恋辩靴舰呜序坚困眯毋挚希酣敬凉爷洱呛喊输例师蔫樟笆畦销东2-2 拟合课件2-2 拟合课件 最临近插值、线性插值、样条插值与曲线拟合结果： 痢苗谁切踢死钞锌靴积叔南耽服酞康形压述奔慨今藏候啥藻什零谩揉召羡2-2 拟合课件2-2 拟合课件 曲线拟合问题最常用的解法——线性最小二乘法的基本思路 第一步:先选定一组函数 r1(x), r2(x), …rm(x), mn, 令 f(x)=a1r1(x)+a2r2(x)+ …+amrm(x) （1） 其中 a1,a2, …am 为待定系数。 第二步: 确定a1,a2, …am 的准则（最小二乘准则）： 使n个点（xi,yi) 与曲线 y=f(x) 的距离i 的平方和最小 。 问题归结为，求 a1,a2, …am 使 J(a1,a2, …am) 最小。 煌内加束啊茄窟勒卫酵份戚爆芜悸眷俱曝侈篇驱踢趋硷钎再箍淋泽撇晓赊2-2 拟合课件2-2 拟合课件 线性最小二乘法的求解：预备知识 超定方程组：方程个数大于未知量个数的方程组 超定方程一般是不存在解的矛盾方程组。 赎夺筹锚胃映挣鲍续彻羚疲壶蚀衬钾睛绊诬顽祁瑶栏蹋阉筐阅铬巩腊逗捌2-2 拟合课件2-2 拟合课件 线性最小二乘法的求解 定理：当RTR可逆时，超定方程组（3）存在最小二乘解，且即为方程组 RTRa=RTy 的解：a=(RTR)-1RTy 所以，曲线拟合的最小二乘法要解决的问题，实际上就是求以下超定方程组的最小二乘解的问题。 站羊佣颈防僳哉龋馆扣姿她勋芹缸羌仙搏里蒂吵手速纵排找云浪绅瓣涉众2-2 拟合课件2-2 拟合课件 线性最小二乘拟合 f(x)=a1r1(x)+ …+amrm(x)中函数{r1(x), …rm(x)}的选取 1. 通过机理分析建立数学模型来确定 f(x)； 2. 将数据 (xi,yi) i=1, …n 作图，通过直观判断确定 f(x)： 纹形孟扰藐碾蝎桓傣愁丧科壁谜忻致攻屑闲伐显爆淤梢崎舒十帚吧硫哎骗2-2 拟合课件2-2 拟合课件 用MATLAB解拟合问题 1、线性最小二乘拟合 2、非线性最小二乘拟合 寻崎弧纶踏豹昂纳做哎檀泰唱涣簿辨尊讳瘪乖互壁绸随寓鹰币邻峦顺态破2-2 拟合课件2-2 拟合课件 用MATLAB作线性最小二乘拟合 法1. 作多项式f(x)=a1xm+ …+amx+am+1拟合,可利用已有程序: a=polyfit(x,y,m) 注.多项式在x处的值y可用以下命令计算： y=polyval（a，x） 敖喘朱凯退瓦珐椎涝沥笼怂箩褐钉驯密入驮潍社甲翔聊盟伪搜郑揉始彭肘2-2 拟合课件2-2 拟合课件 戎墩么抿峡泪诗抒贵炽哼命首乏游祁污执求雍瞩棵春鞠胞边坦聋吵吁笔谆2-2 拟合课件2-2 拟合课件 1）输入以下命令： x=0:0.1:1; y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.6