数据样本缺失前16天的偶数天数据,使用分段插值或埃尔米特插值进行数据样本恢复
实现代码:
clc;clear;
data = xlsread('Z.xlsx')
x = data(1,:);
xNew = 1:16;
dataNew = [xNew]
figure('Name','生成样本表');
for i=2:8
subplot(4,2,i-1)
y = data(i,:);
p0 = spline(x,y,xNew); % 分段样条
p1 = pchip(x,y,xNew); % 分段埃尔米特
dataNew = [dataNew;p0]; % 扩充数据
plot(x,y,'o',xNew,p0,'r-',xNew,p1,'b-')
legend('样本点','三次样条特插值','三次埃尔米特插值','Location','SouthEast')
end
dataNew
不同插值方式生成对应数据样本图示:
恢复数据样本:
数学建模纯小白,望指点