MATLAB 曲线拟合

matlab曲线拟合

clc;
clear;
I= imread('E:\rc\09-27\xuan16-30result\16-30 (8).tiff');
figure,imshow(I,[]);
hold on
x=70;
y  = (4.24e+08)*x^(-3.589)+47.51;

CurtPtx=x;CurtPty=y;%%当前点的坐标
x_left = [];
y_left = [];
% for x=70:180
%     y  = (4.24e+08)*x^(-3.589)+47.51;
%     
%     NewPtx=x;
%     NewPty=y;%%新点的坐标
%     plot([CurtPtx NewPtx],[CurtPty  NewPty]);
%     CurtPtx=NewPtx;
%     CurtPty=NewPty;%%用新点的坐标代替当前点的坐标向前递推
% end
% hold off
% 
% hold on
% x=220;
% y  = (7.631e-13)*x^5.685+32.83;
% CurtPtx=x;CurtPty=y;%%当前点的坐标
% for x=220:310
%     y  = (7.631e-13)*x^5.685+32.83;
%     NewPtx=x;
%     NewPty=y;%%新点的坐标
%     plot([CurtPtx NewPtx],[CurtPty  NewPty]);
%     CurtPtx=NewPtx;
%     CurtPty=NewPty;%%用新点的坐标代替当前点的坐标向前递推
% end
% hold off
% figure,imshow(I,[]);

x_left = 72:180;
y_left = (4.24e+08)*x_left .^ (-3.589