已给sin0.32=0.314567,sin0.34=0.333487,sin0.36=0.352274,用线性插值及抛物插值计算sin0.3367的值并估计截断误差。

1. 线性插值

clc;clear;
y=sin_L(0.32,0.314567,0.34,0.333487,0.3367);
function y=sin_L(x0,y0,x1,y1,x)
% sin_L输出sin(x)使用线性插值计算得到的函数值
% 例如: y=sin_L(0.32,0.314567,0.34,0.333487,0.3367)
% y =
% 0.3304
% R =
% 9.1892e-06 % 以下为判断输入值是否合法的代码
if nargin~= %nargin表示number of function input arguments,判断函数输入参数的个数
error('请输入一次插值的2个插值节点和待插值点')
end
if ~( isnumeric(x0)&&isnumeric(y0)&&isnumeric(x1)&&isnumeric(y1)&&isnumeric(x) ) %isnumeric()检测字符串是否只由数字组成,如果字符串中只包括数字,就返回Ture,否则返回False。
error('输入参数必须是数')
end % 核心计算的代码
%y=y0+(y1-y0)*(x-x0)/(x1-x0);%已知一次多项式,即直线上两点,可以用点斜式写出直线
y=y0*(x-x1)/(x0-x1)+y1*(x-x0)/(x1-x0); %两点式
% 以下为求解截断误差的代码
syms M2; % 因为sin(x)的二阶导数是本身,所以只需要挑出最大的y值,即可得到M2 %syms函数用于创建符号对象
if y0>y1
M2=y0;
else
M2=y1;
end
R=M2*abs((x-x0)*(x-x1))/;
disp(['插值结果是:y=',num2str(y)]);
disp(['最大误差是:R=',num2str(R)]);
end

2. 抛物插值

clc;clear;
%y1=sin_L(0.32,0.314567,0.34,0.333487,0.3367);%一次性插值
y2=sin_T(0.32,0.314567,0.34,0.333487,0.36,0.352274,0.3367)%二次抛物插值
function y=sin_L(x0,y0,x1,y1,x)
% sin_L输出sin(x)使用线性插值计算得到的函数值
% 例如: y=sin_L(0.32,0.314567,0.34,0.333487,0.3367)
% y =
% 0.3304
% R =
% 9.1892e-06 % 以下为判断输入值是否合法的代码
if nargin~= %nargin表示number of function input arguments,判断函数输入参数的个数
error('请输入一次插值的2个插值节点和待插值点')
end
if ~( isnumeric(x0)&&isnumeric(y0)&&isnumeric(x1)&&isnumeric(y1)&&isnumeric(x) ) %isnumeric()检测字符串是否只由数字组成,如果字符串中只包括数字,就返回Ture,否则返回False。
error('输入参数必须是数')
end % 核心计算的代码
%y=y0+(y1-y0)*(x-x0)/(x1-x0);%已知一次多项式,即直线上两点,可以用点斜式写出直线
y=y0*(x-x1)/(x0-x1)+y1*(x-x0)/(x1-x0); %两点式
% 以下为求解截断误差的代码
syms M2; % 因为sin(x)的二阶导数是本身,所以只需要挑出最大的y值,即可得到M2 %syms函数用于创建符号对象
if y0>y1
M2=y0;
else
M2=y1;
end
R=M2*abs((x-x0)*(x-x1))/;
disp(['插值结果是:y=',num2str(y)]);
disp(['最大误差是:R=',num2str(R)]);
end function y=sin_T(x0,y0,x1,y1,x2,y2,x)
% sin_T输出sin(x)使用抛物插值计算得到的函数值
% 例如: y=sin_T(0.32,0.314567,0.34,0.333487,0.36,0.352274,0.3367)
% y =
% 0.3304
% R =
% 2.0315e-07 % 以下为判断输入值是否合法的代码
if nargin~=
error('请输入线性插值的插值节点和插值点')
end
if ~( isnumeric(x0)&&isnumeric(y0)&&isnumeric(x1)&&isnumeric(y1)&&isnumeric(x2)&&isnumeric(y2)&&isnumeric(x) )
error('输入参数必须是数')
end % 核心计算的代码
y=y0*(x-x1)*(x-x2)/((x0-x1)*(x0-x2))+y1*(x-x0)*(x-x2)/((x1-x0)*(x1-x2))+y2*(x-x0)*(x-x1)/((x2-x0)*(x2-x1)); %{
% 以下为求解截断误差的代码
y_0=cos(x0); % 因为sin(x)的三阶导数是cos(x),那么只要求出x0,x1,x2的cos值,然后去最大即可得到M3
y_1=cos(x1);
y_2=cos(x2);
syms M3;
if y_0>y_1
M3=y_0;
else
M3=y_1;
end
if y_2>M3
M3=y_2;
end
R=M3*abs((x-x0)*(x-x1)*(x-x2))/
end
%} % 以下为求解截断误差的代码
%sysm M3;
M3=cos(x0);%因为在第一象限,cosx是减函数,所以节点处的最大值在第一个节点那里取得
R=M3*abs((x-x0)*(x-x1)*(x-x2))/
end

参考来源:https://blog.csdn.net/m0_37395228/article/details/80874393

五,优点和缺点

  拉格朗日插值法的公式结构整齐紧凑,在理论分析中十分方便,然而在计算中,当插值点增加或减少一个时,所对应的基本多项式就需要全部重新计算,于是整个公式都会变化,非常繁琐。这时可以用重心拉格朗日插值法或牛顿插值法来代替。此外,当插值点比较多的时候,拉格朗日插值多项式的次数可能会很高,因此具有数值不稳定的特点,也就是说尽管在已知的几个点取到给定的数值,但在附近却会和“实际上”的值之间有很大的偏差。这类现象也被称为龙格现象,解决的办法是分段用较低次数的插值多项式。(https://www.cnblogs.com/fengzhiyuan/p/8645246.html)

拉格朗日插值matlab实现的更多相关文章

  1. 多项式函数插值:全域多项式插值(一)单项式基插值、拉格朗日插值、牛顿插值 [MATLAB]

    全域多项式插值指的是在整个插值区域内形成一个多项式函数作为插值函数.关于多项式插值的基本知识,见“计算基本理论”. 在单项式基插值和牛顿插值形成的表达式中,求该表达式在某一点处的值使用的Horner嵌 ...

  2. 【matlab】 拉格朗日插值

    第一个函数  "lagrange1.m" 输入:X Y 与点x0 输出:插值函数对应函数值 y0 function y = lagrange1(X,Y,x0) n = length ...

  3. Educational Codeforces Round 7 F - The Sum of the k-th Powers 拉格朗日插值

    The Sum of the k-th Powers There are well-known formulas: , , . Also mathematicians found similar fo ...

  4. 常系数齐次线性递推 & 拉格朗日插值

    常系数齐次线性递推 具体记在笔记本上了,以后可能补照片,这里稍微写一下,主要贴代码. 概述 形式: \[ h_n = a_1 h_{n-1}+a_2h_{n-2}+...+a_kh_{n-k} \] ...

  5. 快速排序 and 拉格朗日插值查找

    private static void QuictSort(int[] zu, int left, int right) { if (left < right) { ; ; ]; while ( ...

  6. BZOJ3601 一个人的数论 莫比乌斯反演、高斯消元/拉格朗日插值

    传送门 题面图片真是大到离谱-- 题目要求的是 \(\begin{align*}\sum\limits_{i=1}^N i^d[gcd(i,n) == 1] &= \sum\limits_{i ...

  7. 【XSY1537】五颜六色的幻想乡 数学 生成树计数 拉格朗日插值

    题目大意 ​ 有一个\(n\)个点\(m\)条边的图,每条边有一种颜色\(c_i\in\{1,2,3\}\),求所有的包括\(i\)条颜色为\(1\)的边,\(j\)条颜色为\(2\)的边,\(k\) ...

  8. 【BZOJ2655】calc DP 数学 拉格朗日插值

    题目大意 ​ 一个序列\(a_1,\ldots,a_n\)是合法的,当且仅当: ​ 长度为给定的\(n\). ​ \(a_1,\ldots,a_n\)都是\([1,m]\)中的整数. ​ \(a_1, ...

  9. 【Luogu4781】【模板】拉格朗日插值

    [Luogu4781][模板]拉格朗日插值 题面 洛谷 题解 套个公式就好 #include<cstdio> #define ll long long #define MOD 998244 ...

随机推荐

  1. C语言单链表简单实现(简单程序复杂化)

    PS: goto还是很好玩的. #include <stdio.h> #include <stdlib.h> typedef struct _node{ int value; ...

  2. java笔记——java加载与执行

    java的加载与执行过程: 编译阶段: 首先,我们先创建一个java源文件xxx.java 2. 经过javac.exe编译(如果源文件没有语法错误),生成一个或多个.class文件,注意的是此时的字 ...

  3. python 之 网络编程(基于TCP协议的套接字通信操作)

    第八章网络编程 8.1 基于TCP协议的套接字通信 服务端套接字函数 s.bind() 绑定(主机,端口号)到套接字 s.listen() 开始TCP监听 s.accept() 被动接受TCP客户的连 ...

  4. JVM运行参数优化详细教程

    获取设置的参数str的值:  常用的-X参数有以下这些: 手动调用GC执行垃圾回收操作:(-XX:+DisableExplicitGC 手动调用将会失效) 查看tomcat的进程ID: 或者:

  5. vscode 连接远程服务器 sftp

    1.在vscode 应用商店搜索 sftp 下载安装 2.ctrl+shift+p 搜索sftp:config 生成sftp.json 3.配置你的服务器信息{ "name": & ...

  6. Java 并发框架Disruptor(七)

    Disruptor VS BlockingQueue的压测对比: import java.util.concurrent.ArrayBlockingQueue; public class ArrayB ...

  7. linux BufferedImage.createGraphics()卡住不动

    项目应用服务器tomcat7,在开发(windows).测试环境(linux 64bit)均正常.在生产环境(linux 64bit)一直启动不起来,也没有报错. 最终定位问题:执行到buffered ...

  8. docker的学习总结

    一 docker的8个使用场景1.简化配置虚拟机的最大好处是能在你的硬件设施上运行各种配置不一样的平台(软件, 系统), Docker在降低额外开销的情况下提供了同样的功能. 它能让你将运行环境和配置 ...

  9. easyui-datagrid 假分页

    假分页就是将数据一下全查出来,利用前端来把所有数据进行分页

  10. Java文件字符流

    1.字符编码(Character encoding)和编码集(Character set) 字符编码(Character encoding)是将字符转为字节或字节数组的过程. 字符集(Characte ...