数学建模算法之插值算法

例题

方案

以三次埃尔米特插值和三次样条插值法为例

代码

clear;clc
load data.mat

x = data(1,:);%第一列

[n,m] = size(data);%n为Z的行数，m为Z的列数

% 注意Matlab的数组中不能保存字符串，
% 如果要生成字符串数组，就需要使用元胞数组，其用大括号{}定义和引用
ylab={'周数','轮虫','溶氧','COD','水温','PH值','盐度','透明度','总碱度','氯离子','透明度','生物量'};  
% 要画的图形的标签

figure(1);
P = zeros(11,15);%结果矩阵赋初值
for i = 2:n
    y = data(i,:);
    new_x=1:15;
    p1 = pchip(x,y,new_x);
    subplot(4,3,i-1);%将所有图依次变现在4*3的一幅大图上
    plot(x,y,'ro',new_x,p1,'-');%画出每次循环处理后的图像
    axis([0 15,-inf,inf])  %设置坐标轴的范围，这里设置横坐标轴0-15，纵坐标不变化
    ylabel(ylab{i})%y轴标题  这里是直接引用元胞数组中的字符串
    P(i-1,:)=p1;%将每次插值之后的结果保存在P矩阵中
end
legend('原始数据','三次埃尔米特插值数据','Location','SouthEast')%加上标注
P = [1:15; P]  %把P的第一行加上周数

figure(2);
P = zeros(11,15);%结果矩阵赋初值
for i = 2:n
    y = data(i,:);
    new_x=1:15;
    p1 = spline(x,y,new_x);
    subplot(4,3,i-1);%将所有图依次变现在4*3的一幅大图上
    plot(x,y,'ro',new_x,p1,'-');%画出每次循环处理后的图像
    axis([0 15,-inf,inf])  %设置坐标轴的范围，这里设置横坐标轴0-15，纵坐标不变化
    ylabel(ylab{i})%y轴标题  这里是直接引用元胞数组中的字符串
    P(i-1,:)=p1;%将每次插值之后的结果保存在P矩阵中
end
legend('原始数据','三次样条插值数据','Location','SouthEast')%加上标注
P = [1:15; P]  %把P的第一行加上周数

结果