close all;
for i=1:5
for j=1:5
poly=['poly',num2str(i),num2str(j)];
fitobject1 = fit([x1,y1],z1,poly);
whole1nh=[fitobject1(x,y)];
goodnessfit(i,j) = goodnessOfFit(whole1nh(:),Wholeimage(:),'MSE');
mesh(Wholeimage)
hold on
mesh(whole1nh)
saveas(gcf,['G:\poly结果评估\',num2str(i),num2str(j),'.jpg']);
end
hold off
end
for i=1:5
for j=1:5
poly=['poly',num2str(i),num2str(j)];
fitobject1 = fit([x1,y1],z1,poly);
whole1nh=[fitobject1(x,y)];
goodnessfit(i,j) = goodnessOfFit(whole1nh(:),Wholeimage(:),'MSE');
mesh(Wholeimage)
hold on
mesh(whole1nh)
saveas(gcf,['G:\poly结果评估\',num2str(i),num2str(j),'.jpg']);
end
hold off
end