例题
方案
以三次埃尔米特插值和三次样条插值法为例
代码
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的第一行加上周数