定义函数:

function yhat = myfun(x,y,xhat)
n = length(y);
h = x(2)-x(1);%等距节点的步长
t = (xhat-x(1))/h;%预测点与初值点相距的步长
c = y(:);
z = zeros(n);
z(:,1) = x';z(:,2) = y';%初始化差分表
for j = 2:nfor i = n:-1:jc(i) = c(i)-c(i-1);z(i-j+1,j+1) = c(i);end
end
y1 = ones(1,n);
r = ones(n,n);
for i = 0:n-1y1(i+1) = c(i+1)/factorial(i);
end
for i =1:n-1r(i+1,1:i) = (t-i+1):t ;b = prod(r,2);
end
yhat = y1*b;
disp("差分表:");disp(z);
fprintf("N_%d(%f) = %f\n",n-1,xhat,yhat);

调用:

x = [20 21 22 23];
y = [1.30103 1.32222 1.34242 1.36173];
xhat = 21.4;
yhat= myfun(x,y,xhat);

牛顿前插公式计算程序例题与matlab代码(计算方法)相关推荐

  1. 牛顿前插公式与后插公式简单例子 BIT数值分析5.5

    题目 思路 公式: 牛顿前插公式 牛顿后插公式 尽量方法一般化,而不是针对某个题目 应该有更简单的方法,但是我没有查资料,直接就上手了 (就这写了好久,脑壳疼,不说了) 代码 # 牛顿前插后插公式 i ...

  2. python:实现newton_forward_interpolation牛顿前插算法(附完整源码)

    python:实现newton_forward_interpolation牛顿前插算法 from __future__ import annotationsimport math# for calcu ...

  3. 计算方法 matlab代码,计算方法及其MATLAB实现

    1.课程简介 计算方法是一门研究求解数学问题数值近似解的专业基础课.作为一门数学课程,计算方法与其它基础数学课程有着本质上的区别,它不仅研究自身的理论,而且更多地与实际问题相结合,提供具有应用价值的理 ...

  4. matlab数理基础之 牛顿-莱布尼兹公式具体应用有哪些,用matlab如何实现

    一.什么是牛顿莱布尼兹公式 牛顿-莱布尼兹公式也被称为牛莱公式,它是微积分学中一个非常重要的公式,用于计算一个区间内函数的定积分. 具体而言,牛顿-莱布尼兹公式描述了定积分与不定积分之间的关系,即一个 ...

  5. 数学建模之微分方程(符实现例题和MATLAB源码)

    微分方程的基本概念 微分方程:一般的,凡表示未知函数.未知函数的导数与自变量之间的关系的方程,叫做微分方程,有时也简称方程. 微分方程的阶:微分方程中所出现的未知函数的最高阶导数的阶数,叫做微分方程的 ...

  6. 艾特肯法方程解matlab程序,牛顿迭代法matlab代码

    牛顿法 迭代公式: x(k1) xk [2 f (x(k) )]1f (x(k) ) Matlab 代码: function [x1,k] =newton(x1,eps) hs=inline('(x ...

  7. matlab用牛顿差值计算三次差值多项式,计算方法用Newton插值多项式求函数的近似值.docx...

    计算方法用Newton插值多项式求函数的近似值 计算方法课程设计题 目: 用Newton插值多项式 处理磁化曲线学 院: 理学院班 级: 学 生 姓 名: 学 生 学 号: 指 导 教 师: 2017 ...

  8. 数学建模微分方程导弹问题matlab求解,数学建模之微分方程(符实现例题和MATLAB源码)...

    微分方程的基本概念 微分方程:一般的,凡表示未知函数.未知函数的导数与自变量之间的关系的方程,叫做微分方程,有时也简称方程. 微分方程的阶:微分方程中所出现的未知函数的最高阶导数的阶数,叫做微分方程的 ...

  9. 电网络节点电压matlab,电力网潮流电压计算例题与matlab程序技术总结.docx

    电力网潮流电压计算例题与MATLAB程序编写佘名寰本文介绍了两个电力网潮流计算例题.一例为5个节点5条支路:另一例有6个节点7条支路,有PQ负荷节点也有PV发电机节点,变压器支路标么变比不为1.本文给 ...

最新文章

  1. PL/SQL 的一些用法
  2. 一分钟教你学会python_十分钟教你学会python编写小游戏
  3. drop 很慢 物化视图_终于解决了物化视图复制的问题
  4. H5实例 移动端页面练习
  5. post url 后面跟参数_都2019年了,还问GET和POST的区别
  6. java面试准备题目
  7. 数据分析的3大阶段6大步骤
  8. 关于蓝天准系统换京东方屏后不显bios的解决方案
  9. 坯子库无法一键安装插件没用_坯子插件库的下载地址、安装以及运用办法
  10. Ubuntu16.04安装网易云音乐后打不开
  11. 陈彤离职,新浪在门户竞争中将继续被边缘化
  12. 如何编辑制作并发送手机报?
  13. 大数据相关书籍(包含Java, Scala, R, Linux, Spark, Hadoop, Hive, Hbase, Sqoop, Flume, Strom)
  14. 读书印记 - 《历代经济变革得失》
  15. Jaeger知识点补充
  16. 积累20180604
  17. 字符编码:ANSI、ASCII、Unicode、UTF-8、UTF-16、UTF-32概念和格式转换
  18. 【郑轻邀请赛 G】密室逃脱
  19. 深入理解JUC——AQS 及 ReentrantLock
  20. 上下Div固定,中间div高度自适应(随浏览器高度变化而变化)

热门文章

  1. python catia harness_介绍一些CATIA汽车线束设计模块的常用命令
  2. Hadoop 基本常用命令
  3. Java搭建订单状态机模型
  4. 2015年二本投档线
  5. RHCSA操作第四次作业
  6. cdn缓存服务器有网站图片,CDN缓存服务器图片存储一致性hash算法的理解
  7. SQL SERVER练习题及答案2
  8. 计算机网络基础常见名词缩写汇总
  9. 态℃ |一条毛巾在拼多多的新消费奇遇记
  10. Python文本变量与函数的解析执行,增强自动化测试数据驱动