close all
clc
clear
x=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30];
y=[64 111 156 155 135 113 78 104 85 77 42 83 56 63 85 81 92 99 80 40 56 43 48 63 53 82 92 52 33 66];
z=[19 18 15 26 27 29 31 31 32 25 24 26 27 32 28 23 21 24 22 25 28 33 35 35 34 28 26 22 29 32];
x=x';%日期
z=y';%PM2.5数据
y=z';%最高温度
steps=0:31;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%数据的拟合
[X,Y]=meshgrid(steps);
Z1=griddata(x,y,z,X,Y);
Z2=griddata(x,y,z,X,Y,'cubic');
Z3=griddata(x,y,z,X,Y,'nearest');
Z4=griddata(x,y,z,X,Y,'v4');
subplot(2,2,1);
mesh(X,Y,Z1);
hold on
plot3(x,y,z,'o');
hold off
subplot(2,2,2);
mesh(X,Y,Z2);
hold on
plot3(x,y,z,'o');
hold off
subplot(2,2,3);
mesh(X,Y,Z3);
hold on
plot3(x,y,z,'o');
hold off
subplot(2,2,4);
mesh(X,Y,Z4);
hold on
plot3(x,y,z,'o');
hold off