最小二乘法辨识

%
clc;
clear all

close all

cu=xlsread（‘asxishu‘‘sheet1‘‘A2:A66‘）;  X1=cu;
xishu=xlsread（‘asxishu‘‘sheet1‘‘B2:B66‘）;  Y=xishu;
%diye=xlsread（‘cusolid‘‘sheet1‘‘C2:C66‘）;  Y=diye;


x=[X1 ]; 
y=Y;

P0=[0.02 0.03 ]; % P0的第一个的ORP的系数，第二个是入口钴（取对数后）的系数   1.0.*x（:3）
%  f=@（px）   x（:2）.*p（1）.*exp（-p（2）.*（x（:1）*（430/150））） 原，430是反应器体积，150是流量
f=@（px）  p（1）*log（x（:1）.*1）+p（2）%（x（:1）.^2）.*p（1）+x（:1）.*p（2）+p（3）