clear;
y=xlsread('E:\works\项目\MC100\电机台架数据\鲨湾96v电机145外特性.xlsx','MAP1','A2:A9'); % torque
x=xlsread('E:\works\项目\MC100\电机台架数据\鲨湾96v电机145外特性.xlsx','MAP1','W2:W9'); % D
z=xlsread('E:\works\项目\MC100\电机台架数据\鲨湾96v电机145外特性.xlsx','MAP1','X2:X9'); % Q
y=abs(y);
% Te = 1.5*4*(a*iq + (ld - lq)*idiq);
ft=fittype('6*(a*iq-(l)*id*iq)', 'independent',{'id','iq'}, 'coefficients', {'a', 'l'});
%ld=0.000070;
%lq=0.000140;
%ft=fittype('6*(a*iq+(ld-lq)*id*iq)', 'independent',{'id','iq'}, 'coefficients', {'a', 'ld', 'lq'});

Y=fit([x,z],y,ft,'StartPoint',[0.001, 0.00070])