- 1、本文档共65页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
9_基于Mathematica的数值计算
第九讲 数值计算 6.1 求近似函数 在生产和实验中, 人们经常遇到需要通过某个未知的函数f(x)在有限个给定点的函数值:{xi, yi}, i=1,2,…., n, 这里 f(xi) = yi 去获得函数f(x)的近似函数?(x), 求近似函数?(x)的方法主要有拟合方法和插值方法。 6.1.1 曲线拟合 曲线拟合主要用来求一元近似函数, 它是根据最小二乘原理的意义下获得近似函数的, 此近似函数具有在数据点处的误差平方和最小的特点。记函数集合: M=Span[?0, ?1, ?2,…, ?m]= {?(x)|?(x)= a0 ?0(x)+a1?1(x)+…+am?m(x), ai ?R} 称集合M为函数?0, ?1, ?2,…, ?m张成的空间,m+1个函数?0(x), ?1(x), ?2(x),…, ?m(x)称为拟合基函数集合, 它们都是已知的函数。 Mathematica曲线拟合的一般形式为: Fit[{数据点集合}, {拟合基函数集合}, 自变量名] 具体的拟合命令有: 命令形式1:Fit[{{x1,y1},{x2,y2},...,{xn,yn}},{ ?0, ?1, ?2,…, ?m },x] 功能:根据数据点集{{x1,y1},{x2,y2},...,{xn,yn}}求出具有拟合函数为 ?(x)= a 0 ?0(x)+a1?1(x)+…+a m?m(x) 形式的近似函数?(x) 命令形式2:Fit[{y1,y2,...,yn},{ ?0, ?1, ?2,…, ?m },x] 功能:根据数据点集{{1,y1},{2,y2},...,{n,yn}}求出具有拟合函数为 ?(x)= a 0 ?0(x)+a1?1(x)+…+a m?m(x) 形式的近似函数?(x) 命令形式3:Fit[{{x1,y1},{x2,y2},...,{xn, yn}}, Table[x^i,{i,0,m}] ,x] 功能:根据数据点集{{x1,y1},{x2,y2},...,{xn,yn}}求出拟合函数为m次多项式的近似函数?(x) =a 0 +a1x+ a2x2 +…+a mx m 多项式拟合算法 输入n+1个拟合点: (xi, yi),i=0,1,…,n 根据散点图确定拟合多项式的次数m 计算相应正规线性方程组的系数和右端项 解正规正规线性方程组,得解:a0*,a1*,…,a m* 写出拟合多项式?*(x)= a0*+ a1*x+ a2*x2+ …+ am*xm 求m次多项式拟合程序 Clear[xi,xx,yi]; xi=Input[xi=] yi=Input[yi=] n=Length[xi]; h=ListPlot[Table[{xi[[i]],yi[[i]]},{i,1,n}],PlotStyle-PointSize[0.04]] m=Input[多项式次数m=] s=Table[Sum[xi[[k]]^i,{k,1,n}],{i,0,2m}]; a=Table[s[[i+j-1]],{i,1,m+1},{j,1,m+1}]; Print[a=,MatrixForm[a]]; b=Table[Sum[xi[[k]]^i*yi[[k]],{k,1,n}],{i,0,m}]; Print[b=,MatrixForm[b]]; xx=Table[x[i],{i,1,m+1}]; g=Solve[a.xx==b,xx]; fa=Sum[x[i]*t^(i-1),{i,1,m+1}]/.g[[1]]; p=fa//N p1=Plot[p,{t,xi[[1]],xi[[n]]},DisplayFunction-Identity]; Show[{p1,h},DisplayFunction-$DisplayFunction]; 说明:本程序用于求m次多项式拟合。程序执行后,按要求通过键盘输入拟合基点xi:{x0 , x1, ... , xn }、对应函数值yi:{ y0 , y1 , … , yn }后,计算机给出散点图和请求输入拟合多项式次数的窗口,操作者可以根据散点图确定拟合多项式的次数通过键盘输入,程序即可给出对应的正规方程组系数矩阵a、常数项b、m次拟合多项式和由拟合函数图形和散点图画在一起的图形。 程序中变量说明 xi:存放拟合基点{x0 , x1, ... , xn } yi: 存放对应函数值{y0 , y1 , … , yn} m: 存放拟合多项式次数 a: 存放正规方程组系数矩阵 b: 存放正规方程组常数项 p: 存放m次拟合多项式 h: 存放散点图 p1:存放拟合函数图形 xx:定义正规方程组变量,存放m次拟合多项式的系数 注:语句s=Ta
您可能关注的文档
- 施工现场灭火应急救援预案.doc
- 建筑施工现场三宝四口五临边如何使用防范好.doc
- 时代家园应急预案.doc
- 机电事故应急预案.doc
- 暂存间管理规程.doc
- 机械设备事故调查分析处理规定.doc
- 景区防汛预案.doc
- 查勘定损初级试题.doc
- 案例分析题3.doc
- 校园安全教育系列.ppt
- 2023年江苏省镇江市润州区中考生物二模试卷+答案解析.pdf
- 2023年江苏省徐州市邳州市运河中学中考生物二模试卷+答案解析.pdf
- 2023年江苏省苏州市吴中区中考冲刺数学模拟预测卷+答案解析.pdf
- 2023年江苏省南通市崇川区田家炳中学中考数学四模试卷+答案解析.pdf
- 2023年江西省吉安市中考物理模拟试卷(一)+答案解析.pdf
- 2023年江苏省泰州市海陵区九年级(下)中考三模数学试卷+答案解析.pdf
- 2023年江苏省苏州市高新二中中考数学二模试卷+答案解析.pdf
- 2023年江苏省南通市九年级数学中考复习模拟卷+答案解析.pdf
- 2023年江苏省南通市海安市九年级数学模拟卷+答案解析.pdf
- 2023年江苏省泰州市靖江外国语学校中考数学一调试卷+答案解析.pdf
最近下载
- 篮球进校园方案.pdf VIP
- CESA--2023-105《算力服务可信性评估规范》.pdf
- 部编版四年级语文上册《精卫填海》教材分析及相关素材.docx VIP
- 2024年共青团知识自测题库及答案.pdf VIP
- 环环相扣PPT立体图表.ppt
- 2024-2025学年八年级英语上学期期中考试(上海卷)(含解析)(牛津上海版).pdf VIP
- 药学专业知识一.pdf
- 湖北省武汉市洪山高级中学2024-2025学年高一上学期9月考试化学试卷.docx VIP
- 2024年公共卫生与预防医学继续教育平台“大学习”活动线上培训栏目考试题库(1392题).docx
- 地方标准审评工作细则.pdf VIP
文档评论(0)