python编写分段Hermite插值多项式

embedded/2025/1/11 14:20:38/
python">x = [0, 1, 2];   % x 坐标
y = [1, 2, 0];   % y 坐标
dy = [1, -1, 2]; % 导数值
X = linspace(0, 2, 100);  % 生成 100 个插值点n = length(x);  % 数据点的个数m = length(X);  % 插值点的个数H = zeros(1, m); % 用于存储结果for k = 1:mxi = X(k);  % 当前插值点% 计算每个分段的插值sum_val = 0;for i = 1:n-1% 构造 Hermite 插值的每一段if xi >= x(i) && xi <= x(i+1)% 计算 Hermite 插值多项式t = (xi - x(i)) / (x(i+1) - x(i)); % 归一化变量% Hermite插值的基函数h00 = 2*t^3 - 3*t^2 + 1;h10 = t^3 - 2*t^2 + t;h01 = -2*t^3 + 3*t^2;h11 = t^3 - t^2;% 计算当前段的插值sum_val = y(i)*h00 + dy(i)*(x(i+1)-x(i))*h10 + y(i+1)*h01 + dy(i+1)*(x(i+1)-x(i))*h11;endendH(k) = sum_val;end% 绘制插值结果
plot(X, H, 'r', x, y, 'bo');
legend('Hermite Interpolation', 'Data Points');
title('Piecewise Hermite Interpolation');
xlabel('x');
ylabel('y');

在这里插入图片描述
注:

在 Python 中,** 是一个 运算符,它有两种主要用途:

1. 幂运算(指数运算)

  • ** 用于计算一个数的幂,即对一个数字进行指数运算。例如:
python">result = 2 ** 3  # 计算 2 的 3 次方
print(result)    # 输出 8

在这个例子中,2 ** 3 计算的是 2 的 3 次方,结果为 8。


http://www.ppmy.cn/embedded/153031.html

相关文章

Jenkins git SSH获取code报错:git@github.com: Permission denied (publickey).

这个错误信息表明在尝试通过 SSH 连接到 GitHub 时&#xff0c;出现了权限被拒绝的问题&#xff0c;通常是由于 SSH 公钥未正确配置或未被 GitHub 识别。以下是解决此问题的步骤&#xff1a; 1. 确保 SSH 密钥已生成 首先&#xff0c;检查你是否已经生成了 SSH 密钥。如果没有…

移远BC28_opencpu方案_pin脚分配

先上图&#xff0c;BC28模块的pin脚如图所示&#xff1a; 下面看看GPIO的复用管脚 然后我自己整理了一份完整的pin功能列表

Docker能跑些什么呢?第三期

Docker 是一个开源的应用容器引擎&#xff0c;让开发者可以打包他们的应用以及依赖包到一个可移植的镜像中&#xff0c;然后发布到任何流行的 Linux或Windows操作系统的机器上&#xff0c;也可以实现虚拟化。容器是完全使用沙箱机制&#xff0c;相互之间不会有任何接口。 下面…

数据挖掘实训:基于CEEMDAN与多种机器学习模型股票预测与时间序列建模

股票市场的预测是金融领域中的一个重要研究方向&#xff0c;尤其是如何通过历史数据预测未来股价。这篇博客将带你走进如何通过多种方法&#xff0c;如CEEMDAN&#xff08;Complete Ensemble Empirical Mode Decomposition with Adaptive Noise&#xff09;分解技术与机器学习模…

腾讯云AI代码助手编程挑战赛-知识百科AI

作品简介 知识百科AI这一编程主要用于对于小朋友的探索力的开发&#xff0c;让小朋友在一开始就对学习具有探索精神。在信息化时代下&#xff0c;会主动去学习自己认知以外的知识&#xff0c;同时丰富了眼界&#xff0c;开拓了新的知识。同时催生了在大数据时代下的信息共享化…

Ubuntu 22.04 桥接配置

1. 网卡配置文件备份 sudo cp /etc/netplan/01-network-manager-all.yaml /etc/netplan/01-network-manager-all.yaml.bak 1.1 静态ip配置 1.1.3 解释 network: version: 2 ethernets: # 网卡名称 使用ifconfig查看 ens33: dhcp4: false dhcp6: false …

语音实战(一)中文语音识别

一、下载文件 from datasets import load_datasetname mozilla-foundation/common_voice_16_0 load_dataset(name, zh-CN, splittrain).save_to_disk(dataset/ name) from transformers import Wav2Vec2BertForCTCname lansinuote/Chinese_Speech_to_Text_CTC Wav2Vec2BertF…

毕业设计-模拟闯红灯的抓拍系统

1 前言 随着城市化进程的加速&#xff0c;交通安全问题日益凸显。闯红灯作为一种常见的交通违法行为&#xff0c;不仅严重威胁道路交通安全&#xff0c;还容易引发交通事故。为了有效遏制闯红灯行为&#xff0c;提高公众的交通安全意识&#xff0c;本文提出了一种模拟闯红灯抓…