山野莽夫

  • 归档
    • 随笔
    • 建站资源
    • 分享
    • 代码
  • 地球物理学
    • 专业课
    • 概念解释
  • 计算机
  • 互联网
  • 教程
  • 规划
  • 实验室
    • 珍藏的软件
    • 贴吧云签到
    • A1账号自助申请
山野莽夫
小学生的挣扎的点点滴滴
  1. 首页
  2. 地球物理学
  3. 正文

计算方法实习----牛顿前插公式

2016年6月6日 3563点热度 0人点赞 0条评论
  1. 实验三 牛顿前插公式
    1. 实验要求(含有算法说明和程序说明)

    牛顿插值公式相比拉格朗日插值法,有一定的优越性,简单计算方便。我们使用牛顿插值公式的时候,如果遇到等距节点,即节点为xi=x0+ih(i=0,1,2……);h为步长这个时候可以简化计算公式,同时避免了除法运算,因为在计算机中做一次除法运算必做一次加法运算运算量大得多,精度也较难控制。因此我们提出了差分的概念

    ∆mfi=∆m-1fi+1-∆m-1fi为m阶差分。根据差分和差商的关系于是得到牛顿前插公式.根据牛顿前插公式,通过对公式的分解,然后把公式进行组合既可以得到这个程序。

    1. 源代码

    function ntqc3(x,y,xhat)

    x=[20 21 22 23 24];

    y=[1.30103 1.32222 1.34242 1.36173 1.38021];

    n=length(y);

    A=zeros(n,n);

    xhat=21.4;

    for i=1:n

    A(i,1)=y(i);

    end

    for j=2:n

    fprintf('\n ');


    for k=j:5

    A(k,j)=A(k,j-1)-A(k-1,j-1);


    end

    end

    t=(xhat-x(1))/(x(2)-x(1));

    N=A(1,1);

    for m=2:n-1

    M=1;


    for n=1:m-1

    M = M*(t-n+1);


    end

    P=1;


    for q=1:m-1

    P = P*q;


    end

    N = N+A(1,m)/P*M;

    end

    A

    ('当x=%f时,y的值是 %f\n',xhat,N);

    function y=f(x) %f Ϊ±»»ýº¯Êý±í´ïʽ

    % y =sin (x);

    y =1./sqrt(1+x.*x.*x);

    end

3、实验结果


A =

1.3010 0 0 0 0

1.3222 0.0212 0 0 0

1.3424 0.0202 -0.0010 0 0

1.3617 0.0193 -0.0009 0.0001 0

1.3802 0.0185 -0.0008 0.0001 -0.0000

当x=21.400000时,y的值是 1.330413

>>

标签: 牛顿前插公式 计算方法 试验
最后更新:2018年4月6日

小菜菜

菜鸟

打赏 点赞
< 上一篇
下一篇 >

文章评论

取消回复

此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据。

标签聚合
wordpress onedrive 模板 宝塔面板 c语言 地震学程序 虚拟机 ppt
最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
直立长方体重力异常及其一阶偏导数 本站提供免费office365教育版A1子号(5Tonedrive网盘),可以自助申请 wordpress添加水印---dx-watermark插件 认识测井软件--CifLog 微软新谷歌chromium内核Edge浏览器下载与简单使用体验 Phpstudy出品免费linux小皮面板简单安装和使用评测-让天下没有难配的环境

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2