牛顿插值实验报告C语言,牛顿插值C语言.doc
《数值分析》实验报告
学号: 20120921101 姓名: 孙腾 班级: 计算1201 日期:2012/10/20
题目:Newton插值法的函数表如下
x
0.0 0.1 0.2 0.3 0.4
0.5000 0.5398 0.5793 0.6179 0.7554
(1)编程实现求Newton插值(2)用Newton算法求的近似值,统计计算时间,并与lagrange算法进行比较。
1、算法:step 1 输入[x0,x1,…xn],[y1,y2,…yn];
step 2 对j=0,1,2…n;
di
step 3 对k=0,1,2…n-1;
dj=(dj-d(j-i))/(xi-x(j-k-1));
step 4 输出[d0,d1,…dn]
程序:(Matlab)
①建立自定义函数Chashang.m
function f=Chashang(x,y,X)
syms t; %定义符号变量t,进行公式的化简和计算;
n=length(x); %测量向量x的长度,赋给n;
m=length(y); %测量向量y的长度,赋给m;
if m~=n %判断m和n是否相等,就是判断x与y是否一一对应;
error('样本数据中的x与y的对应个数不匹配');
end
A=zeros(n,n); %定义一个n行n列的零矩阵;
A(:,1)=y'; %把向量y转置,赋给零矩阵的第一列;
for j=2:n %第一个循环,变量为j,用来表示第几行;
for i=1:(n-j+1) %第二个循环,变量为i,用来表示第几列;
A(i,j)=(A(i+1,j-1)-A(i,j-1))/(x(i+j-1)-x(i)); %差商公式,A(i,j)表示零矩阵的第j行,第i列;
end
end
A %得到差商矩阵A;
②编写函数调用。
x=[]; %对应x值;
y=[]; %已知的函数值;
X=[]; %插值节点;
f=Chashang(x,y,X) %调用函数得到差商值;
2、㈠利用1中的求差商的算法,把求出的差商矩阵的提取对角线元素,利用秦九韶算法,得到牛顿插值多项式。
①首先建立自定义函数Newton.m
function f=Newton(x,y,X)
syms t; %定义符号变量t,进行公式的化简和计算;
n=length(x); %测量向量x的长度,赋给n;
m=length(y); %测量向量y的长度,赋给m;
if m~=n %判断m和n是否相等,就是判断x与y是否一一对应;
error('样本数据中的x与y的对应个数不匹配');
end
A=zeros(n,n); %定义一个n行n列的零矩阵;
A(:,1)=y'; %把向量y转置,赋给零矩阵的第一列;
for j=2:n %第一个循环,变量为
牛顿插值实验报告C语言,牛顿插值C语言.doc相关推荐
- mysql的触发器实验报告_数据库原理实验报告s11-数据库触发器的创建.doc
数据库原理实验报告s11-数据库触发器的创建.doc 数据库管理系统SQLSERVER实验报告第1页2011年5月4日实验11数据库触发器的建立实验日期和时间20141128实验室软件工程室班级12计 ...
- 计算机图形学直线线型实验报告,计算机图形学实验报告-直线中点bresenham算法的实现资料.doc...
计算机图形学实验报告-直线中点bresenham算法的实现资料.doc (10页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 29.90 积分 计算机图形 ...
- 上海应用技术学院c语言实验报告9,上海工程技术大学C语言实验报告
<上海工程技术大学C语言实验报告>由会员分享,可在线阅读,更多相关<上海工程技术大学C语言实验报告(7页珍藏版)>请在人人文库网上搜索. 1.实 验 报 告高级语言程序设计20 ...
- c语言实验分支程序设计二,C语言程序实验报告分支结构的程序设计(0页).doc
C语言程序实验报告分支结构的程序设计(0页) 数学与软件科学学院 实验报告 学期:11至12___ 第_1 学期 2011年10 月 17 日 课程名称:程序设计基础教程-C语言 专业:2010级5_ ...
- c语言设计实验报告答案,武汉理工大学《C语言程序设计》实验报告答案
武汉理工大学<C语言程序设计>实验报告答案 注:在Visual C++ 6.0编译环境中亲自调试通过,但不保证在Turbo C中通过. 实验二 选择结构的程序设计 (题目当初没抄下来,这是 ...
- 2015华北水院c语言实验报告5,华北水利水电大学C语言实验5.doc
华北水利水电大学C语言实验5 C语言程序设计实验报告 实验5 使用数组的程序设计 班级 2013156 学号 201315613 姓名 吴浩 一.实验目的二.实验内容及步骤#include int m ...
- 离散数学c语言实验报告,离散数学数理逻辑C++或C语言实验报告
离散数学实验报告 专业班级:12级计算机本部一班 姓名:鲍佳珍 学号: 1016 实验成绩: 1.[实验题目] 命题逻辑实验一 2.[实验目的] 熟悉掌握命题逻辑中的联接词,实现二元合取.析取.蕴涵 ...
- c语言第9章实验报告,《经典规范C语言程序设计》实验报告(第九章)文件程序设计.doc...
<经典规范C语言程序设计>实验报告(第九章)文件程序设计 <高级语言程序设计>实验报告 班级: 学号: 姓名: 成绩: 实验日期: 年 月 日 实验题目: 实验九 文件程序设计 ...
- c语言综合实验报告包括5个子系统,C语言综合实验课程设计
实验7 综合实验 一.实验目的 ? ? ? ? 掌握系统功能划分的方法 掌握模块设计的方法 熟练掌握在Code::Blocks中的跟踪调试与断点调试方法 培养良好的编程风格 二.实验环境 ? 个人计算 ...
最新文章
- php获得帮助类数据_PHP实现的一个时间帮助类
- EWSD命令输入格式解释
- SSRS动态设置文本框属性
- mysql中的参数如何调试_mysql 查询优化 ~ 查询参数调节
- 民科微服务小程序怎么注册_民科微服务小程序(医师护士电子化注册系统入口)公众号...
- 程序员被纳入新生代农民工;“腾讯视频崩了”上热搜;英特尔发布全新独立显卡品牌 Arc|极客头条...
- jfianl中 解决jetty Form too large
- Spring实战4:面向切面编程
- 【统计分析】2 地理统计
- java中如何进入代码块_Java中的匿名代码块
- java面向对象程序设计实验指导答案,Java面向对象程序设计实验指导与习题解答...
- 盘点:2012中国互联网大会十大亮点
- java开发一年后学习计划
- python的神经网络模块接法图解_关于pytorch中全连接神经网络搭建两种模式详解...
- 趋势跟踪系统的形成历程
- 1. 系统符号与正则符号
- win7网络上计算机进不去,Win7电脑开机进不了系统怎么办?
- uni-app 商城 的sku算法(vue)
- a标签跳转的四种方法
- 五子棋游戏图形化实现
热门文章
- stm32看门狗_STM32单片机:独立看门狗、窗口看门狗的配置
- Shell脚本实现每个工作日定时执行指定程序
- HDU - 1525
- Ubuntu 16.04安装Bless十六进制编辑器
- 使用 HTML5 File API 实现client log
- 数据结构与算法之排序(归纳总结三)
- ios进度条Demo一个
- java.lang.ClassNotFoundException: net.sf.ezmorph.Morpher
- U811.1接口EAI系列之二-BOM构成-委外BOM构成--VB语言
- rvm的安装, 使用rvm, 安装ruby, 以及gem的使用 (转)