


下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
多项式拟合1数表的拟合计算所谓拟合是指寻找一条平滑的曲线,最不失真地去表现测量数据。反过来说,对测量的实验数据,要对其进行公式化处理,也就是用一种计算方法,构造一个函数来近似表达数表的函数关系。由于函数构造方法的不同,有许多的逼近方法,机械设计中常用最小平方逼近(最小二乘法理论)来实现曲线的拟合。根据该理论可推导出计算公式,而MATLAB在此数学基础上用一个函数命令polyfit即可实现,命令格式为:P=polyfit(x,y,n)式中:x、y为已知数据,n为拟合多项式的阶次,p为返回所得多项式的系数向量,通常多项式拟合中阶数越大,拟合的精度就越高。例1:在工程技术中,通过实验获得一组y二f(x)实验数据如下表1:ii表1实验实测数据x0123456789y01.23.8&517.120.234.845.067.685.0下面程序分别设n=1,n=2进行一阶和二阶的拟合,结果如图1所示。%曲线拟合(Curvefitting)x=[0,1,2,3,4,5,6,7,8,9]y=[0.0,1.2,3.8,8.5,17.1,20.2,34.8,45.0,67.6,85.0]%绘实验节点数据plot(x,y,'*r')holdongrid%绘一阶拟合曲线p1=polyfit(x,y,1)py1=polyval(p1,x)plot(x,py1,'g')holdongrid%绘二阶拟合曲线py2=polyval(p2,x)plot(x,py2,'b')gridholdonp2=polyfit(x,y,2)py2=polyval(p2,x)plot(x,py2,'b')gridholdon图1曲线拟合根据所编制的程序,在MATLAB的工作空间可得n=1时,p=[3.3188-2.4545],1线性拟合得一条直线,即直线方程式为:y=3.3188x-2.45451n=2时,p=[0.16481.8358-0.4773],拟合得一条二次曲线,即曲线方程2式为:y二0.1648x2+1.8358x-0.477322数表的插值计算插值是在已知的节点数据基础上,平滑地估算出节点数据之间其它点对应的函数值,也就是通过实验或将一张已知图表处理成数表后,在实际使用时,自变量值不是数表中给定的节点时,须用插值计算方法求出相应的值。其数学基础是差分。根据数表维数不同可分为一维插值、二维插值、三维插值,根据构造插值函数的不同可分为线性插值、抛物线插值、样条插值等。同样在MATLAB中,可以直接用函数来实现。一维插值:Y=interp1(X,Y,X,'method')ii二维插值:Z=interp2(X,Y,乙X,Y,'method')iii三维插值:V=interp3(X,Y,Z,V,X,Y,Z,'method')iiii三次样条插值:Y二spline(X,Y,X)ii式中:对一维X为插值点,Y为插值点对应的函数值,二维、三维以此类推,等式右ii边为插值点,左边为对应的插值函数值。参数'method'表示用指定的插值方法进行插值,可取'nearest'(最近插值)、'linear'(线性插值)、'spline'(三次样条插值)、'cubic'(三次插值),默认方法是线性插值。例2:对一维插值,设有一组数据:x012345678y00.90.61.00.1-0.3-0.7-0.9-0.2要求在上述十个样本点间,产生对应的插值点以使数据进行更为平滑的效果处理。采用一维样条插值spline使输出结果平滑,上述问题的程序如下,结果如图2所示。%====一维样条插值====
y=[00.90.610.1-0.3-0.7-0.9-0.2]x=0:length(y)-1;x1=0:0.1:length(y)-1;y1=interp1(x,y,x1,spl
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 体育中心改造可行性研究报告
- 铁路车站客运管理课件
- 月幼儿园保健工作方案模板
- 2025-2030年中国站用屏项目投资可行性研究分析报告
- 2024-2030年中国泉州房地产行业市场全景分析及投资策略研究报告
- 2025年间苯二胺市场调研报告
- 4s店前台接待工作总结
- 中国桌面型资料架行业市场发展前景及发展趋势与投资战略研究报告(2024-2030)
- 秋季学期学校德育工作方案模板
- 中国微晶玻璃行业发展监测及市场发展潜力预测报告
- 用电办理授权委托书(新)
- 电气值班员技师实操题
- 手术物品清点课件
- 物理word版2023山东高考答题卡涂准考证号和条形码
- WS/T 512-2016医疗机构环境表面清洁与消毒管理规范
- GB/T 36089-2018丙烯腈-丁二烯橡胶(NBR)
- GB/T 35694-2017光伏发电站安全规程
- GB/T 26746-2011矿物棉喷涂绝热层
- GA/T 828-2009电子物证软件功能检验技术规范
- 中耳疾病(分泌性中耳炎、急慢性中耳炎)
- 20XX贵州省高职(专科)分类真习题卷
评论
0/150
提交评论