【MATLAB第75期】#源码分享 | 基于MATLAB的不规则数据插值实现时间序列数据扩充

【MATLAB第75期】#源码分享 | 基于MATLAB的不规则数据插值实现时间序列数据扩充

如时间数据以单位1为间隔排序， 可插间隔为0.5的数据 。

一、实现效果

1.规则间隔数据

2.非规则间隔数据

二、主程序代码

1.插值测试效果

%%  清空环境变量
warning off             % 关闭报警信息
close all               % 关闭开启的图窗
clear                   % 清空变量
clc                     % 清空命令行

%%  导入数据（时间序列的单列数据）
y= xlsread('数据集.xlsx');
N=size(y,1);
x=1:1:N;   % x也可以为不规则间隔插值  ， 如x=sort(rand(N,1))
d=0.5;%间隔设置  ，根据x变量的间隔设置
p=0.01;%插值密度大于0，p值越小精度越高。

%% 测试插值效果
xt1=(x)';%测试样本 插值扩充后的x轴
xt2=xt1+p;%测试样本  插值偏移位置为 xt1和xt2坐标的中点
yti=cuts(x,y,xt2,xt1);%测试样本 插值结果
error = sqrt(sum((yti -y ).^2) ./ N);
  %绘图 
  figure()
  plot(x,y,'b-x',reshape([xt1 xt2]',[],1),reshape([yti yti]',[],1),'ro',mean([xt1 xt2],2),yti,'ro', 'LineWidth', 2)
  legend('实际值','插值')
xlabel('x')
ylabel('y')
string = {
    
    '测试集插值误差';['RMSE=' num2str(error)]};
title(string)
grid

2.正式插值

%% 正式插值
x1=(0.5:1:N-0.5)';%测试样本 插值扩充后的x轴
x2=x1+p;%测试样本  插值偏移位置为 xt1和xt2坐标的中点
yi=cuts(x,y,x2,x1);%测试样本 插值结果
  %绘图 
  figure()
  plot(x,y,'b-x',reshape([x1 x2]',[],1),reshape([yi yi]',[],1),'ro',mean([x1 x2],2),yi,'ro', 'LineWidth', 2)
  legend('实际值','插值')
xlabel('x')
ylabel('y')
string2 = {
    
    '插值效果'};
title(string2)
grid

3.非规则数据插值

自定义x区间范围，比如同样是100个样本，x轴编号通过随机在1:500中抽取100个数字作为x轴，y不变。也可以包含小数 ，或者自己随机定义。

xx = randperm(500); %
x=sort(xx(1:N));   % x也可以为不规则间隔插值

正式插值时x1改为：

x1=(1:10:500)';%测试样本 插值扩充后的x轴

三、代码获取

CSDN后台私信回复“75期”即可。