lsqcurvefit-torque.m

Lsqcurvefit-Torque

for i=0:1:10
thvec=[0;10;20;30;40;50;60;70;80;90;100]
nevec=[950 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000]%(RPM)
emap=[0   -7.9498  -12.2427  -16.2176  -19.0795  -23.6904  -26.7113  -29.0962  -31.1632  -33.2301 -35.6151;
      8.2201   11.2728   12.2586   13.7690   15.9950   -4.2134   -8.2042  -11.3046  -14.0711  -16.8536  -20.0494;
      21.2100   22.8636   26.5841   27.8561   33.9138   15.2636   10.3029    6.4870    3.0209   -0.4770   -4.4837;
      40.5598   41.9272   42.6586   43.5330   51.8326   34.7406   28.8100   24.2787   20.1130   15.8996   11.0820;
      58.3197   56.8569   58.4151   59.2100   69.7356   54.2176   47.3171   42.0703   37.2050   32.2761   26.6477;
      72.8996   73.0586   76.0795   78.0669   87.6544   73.6946   65.8243   59.8619   54.2971   48.6527   42.2134;
      87.4795   89.2602   87.0661   96.9238  105.5732   93.1715   84.3314   77.6535   71.3891   65.0293   57.7791;
      102.0594  105.4619  111.5673  115.7807  123.4761  112.6485  102.8385   95.4452   88.4812   81.4058   73.3448;
      116.6393  121.6636  129.7088  134.6376  141.3949  132.1255  121.3456  113.2368  105.5732   97.7824   88.9104;
      131.2192  137.8652  147.2142  153.4945  159.3138  151.6025  139.8527  131.0284  122.6652  114.1590  104.4761;
      145.7991  154.0669  164.5606  172.3514  177.1213  171.0795  158.3598  148.8200  139.7573  130.5355  120.0418]%(N-m)

nevec=[950 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000]/100000;
emap(1)=[145.7991  154.0669  164.5606  172.3514  177.1213  171.0795  158.3598  148.8200  139.7573  130.5355  120.0418];
x0 = [1,1,1,1,1,1,1] ;
f=optimset;f.TolFun=1e-10;
[x,resnorm] = lsqcurvefit(@Torque,x0,nevec,emap(1),[],[],f)