山野莽夫

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

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

2016年6月6日 4275点热度 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日

小菜菜

菜鸟

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

文章评论

razz evil exclaim smile redface biggrin eek confused idea lol mad twisted rolleyes wink cool arrow neutral cry mrgreen drooling persevering
取消回复

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

联系方式

QQ群 | TG群 | 邮箱

最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
免费不限速大容量云盘服务-小麦云盘-520免费升级专业版账户活动还有免费临时存储-云链服务 宝塔面板配置onedrive索引工具OLAINDEX 数据结构(C语言版)第2版参考答案 申请google voice 方法-纪念我撸到的第一个gv号 水平多层介质反射波时距曲线地震检波器组合方向频率特性曲线 提供一个科学上网工具
标签聚合
wordpress ppt 宝塔面板 c语言 模板 地震学程序 onedrive 虚拟机
最近评论
小菜菜 发布于 7 个月前(11月24日) 这玩意已经废了,成收割工具了,不能再用了。
eamon 发布于 8 个月前(11月07日) 我一年不用了才发现这个休眠管理费每月15,一共扣了我135元,然后我消费还消费不了,我宁愿消费掉也不...
magic 发布于 12 个月前(07月03日) 请问账号不注销会有什么影响吗?
magic 发布于 12 个月前(07月01日) 我想问一下 如果不注销账号就留着会怎么样
qwp6601 发布于 1 年前(06月04日) 有没有方法改为bing

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2