王维颖应用数学班.docVIP

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
王维颖应用数学班.doc

《》 课程报告 选题名称:《Newton插值法》 学生姓名:王维颖 学 号:201318050122 班 级:应用数学131801班 完成时间: 2015 年 7 月17日 目 录 一法的发展历史 二的基本思想 的应用情况.数值实验 参考文献一的发展历史 插值法利用函数f (x)在某区间中若干点的函数值,作出适当的特定函数,在这些点上取已知值,在区间的其他点上用这特定函数的值作为函数f (x)的近似值。如果这特定函数是多项式,就称它为插值多项式。利用插值基函数很容易得到拉格朗日插值多项式,公式结构紧凑,在理论分析中甚为方便,但当插值节点增减时全部插值基函数均要随之变化,整个公式也将发生变化, 这在实际计算中是很不方便的,为了克服这一缺点,提出了牛顿插值。 牛顿插值通过求各阶差商,递推得到的一个公式: 牛顿插值与拉格朗日插值具有唯一性。 二的基本思想 将待求的n次插值多项式改写为具有承袭性的形式,然后利用插值条确定的待定系数,以求出所要的插值函数。。构造牛顿插值多项式。输入要计算的函数点,x并计算的值,利用牛顿插值公式,当增加一个节点时,只需在后面多计算一项,而前面的计算仍有用;另一方面的各项系数恰好又是各阶差商,而各阶差商可用差商公式来计算。 为的一阶差商, 为的n阶差商。 差商表: 一阶差商 二阶差商 ... n阶差商 ... ... ... ... ... ... ... ... 牛顿插值法计算步骤: 1.输入n值及;要计算的函数点x。 2.对给定的x,由 计算的值。 3. 输出。 四、算法的应用情况与凸轮轮廓接触,并传递动力和实现预定的运动规律的构件,一般做往复直线运动或摆动,称为从动件凸轮机构在应用中的基本特点在于能使从动件获得较复杂的运动规律。因为从动件的运动规律取决于凸轮轮廓曲线,所以在应用时,只要根据从动件的运动规律来设计凸轮的轮廓曲线就可以了凸轮机构广泛应用于各种自动机械、仪器和操纵控制装置。凸轮机构之所以得到如此广泛的应用,主要是由于凸轮机构可以实现各种复杂的运动要求,而且结构简单、紧凑。 盘形凸轮、圆柱凸轮和移动凸轮等,其中圆柱凸轮的凹槽曲线是空间曲线,因而属于空间凸轮。从动件与凸轮作点接触或线接触,有滚子从动件、平底从动件和尖端从动件等。尖端从动件能与任意复杂的凸轮轮廓保持接触,可实现任意运动,但尖端容易磨损,适用于传力较小的低速机构中。为了使从动件与凸轮始终保持接触,可采用弹簧或施加重力。具有凹槽的凸轮可使从动件传递确定的运动,为确动凸轮的一种。一般情况下凸轮是主动的,但也有从动或固定的凸轮。多数凸轮是单自由度的,但也有双自由度的劈锥凸轮。凸轮机构结构紧凑,最适用于要求从动件作间歇运动的场合。它与液压和气动的类似机构比较,运动可靠,因此在自动机床、内燃机、印刷机和纺织机中得到广泛应用。但凸轮机构易磨损,有噪声,高速凸轮的设计比较复杂,制造要求较高 .数值实验 n=length(X);A=zeros(n,n);A(:,1)=Y; s=0.0;p=1.0;q=1.0;c1=1.0; for j=2:n for i=j:n A(i,j)=(A(i,j-1)-A(i-1,j-1))/(X(i)-X(i-j+1)); end b=poly(X(j-1));q1=conv(q,b);c1=c1*j;q=q1; end C=A(n,n);b=poly(X(n));q1=conv(q1,b); for k=(n-1):-1:1 C=conv(C,poly(X(k)));d=length(C);C(d)=C(d)+A(k,k); end L(k,:)=poly2sym(C);Q=poly2sym(q1); syms M wcgs=M*Q/c1;Cw=q1/c1 在matlab工作窗口输入 X=[-2.15 -1.00 0.01 1.02 2.03 3.25]; Y=[17.03 7.24 1.05 2.03 17.06 23.05]; [A,C,L,wcgs,Cw]=newpoly1(X,Y) 运行后输出 Cw = 0.0014 -0.0044 -0.0079 0.0243 0.0061 -0.0202 0.0002 A = 17.0300 0 0 0 0 0 7.2400 -8.5130 0 0 0 0 1.0500 -6.1287 1.1039 0 0

文档评论(0)

整理王 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档