ieee14.c

一种求解目标函数最小化的MATLAB仿真程序.

      J[701] = 0.0;      J[702] = 0.0;      J[703] = 0.0;      t1282 = -J[189];      J[704] = t1087*t867+J[189]*t870+t1282*t1008;      J[705] = 0.0;      t1289 = t663+t662+t666+t665+t671-t670;      t1293 = -t293*t288-t295*t290;      t1295 = -t671+t670;      J[706] = t1017*t1000+t1108*t867+t1180*t870-2.0*t490*t1005+t1289*t1008+t1293*t1047+t1295*t1050; 
      t1297 = -J[192];      t1301 = t311*t288-t295*t308;      J[707] = t1297*t1008+t1301*t1047+J[192]*t1050;      J[708] = 0.0;      J[709] = 0.0;      J[710] = 0.0;      J[711] = 0.0;      J[712] = 0.0;      J[713] = 0.0;      J[714] = 0.0; 
      J[715] = 0.0;      J[716] = 0.0;      J[717] = t1007*t1000+t1255*t1005+t999*t1008;      J[718] = t1042*t867+t1137*t870+J[198]*t1005+J[221]*t1008;      J[719] = 0.0;      t1317 = -J[246];      t1319 = -J[269];      J[720] = t1013*t1000+t1064*t867+t1153*t870+t1265*t1005-J[223]*t1008+t1317*t1047+t1319*t1050;      t1323 = -t720+t719;      t1325 = -J[224]; 
      J[721] = J[201]*t1005+J[224]*t1008+t1323*t1047+t1325*t1050;      J[722] = 0.0;      J[723] = 0.0;      J[724] = 0.0;      J[725] = 0.0;      J[726] = 0.0;      J[727] = t1089*t867+J[212]*t870+t1282*t1005-J[212]*t1008;      J[728] = 0.0;      J[729] = J[76]*t1000+J[99]*t867+J[122]*t870+t1289*t1005+J[214]*t1008+J[237]*t1047+J[260]*t1050;      t1342 = -t698-t697; 
      J[730] = t1297*t1005-J[215]*t1008+t1342*t1047+J[215]*t1050;      J[731] = 0.0;      J[732] = 0.0;      J[733] = 0.0;      J[734] = 0.0;      J[735] = 0.0;      J[736] = 0.0;      J[737] = 0.0;      J[738] = 0.0;      J[739] = 0.0; 
      J[740] = 0.0;      J[741] = t1046*t867+t1139*t870+J[244]*t1050;      J[742] = 0.0;      J[743] = t1269*t1005+t1317*t1008+J[246]*t1050;      t1357 = t717-t716+t720-t719+t725+t724+t728+t727;      t1361 = -t323*t320+t328*t326;      t1362 = v[12];      t1364 = -t725-t724;      t1365 = v[13];      t1369 = -t336*t333+t341*t339; 
      t1370 = v[20];      t1372 = -t728-t727;      t1373 = v[21];      J[744] = t1069*t867+t1157*t870+t1276*t1005+t1323*t1008+2.0*t316*t1047+t1357*t1050+t1361*t1362+t1364*t1365+t1369*t1370+t1372*t1373;      t1378 = t346*t320+t328*t349;      t1380 = -J[248];      J[745] = J[248]*t1050+t1378*t1362+t1380*t1365;      J[746] = 0.0;      J[747] = 0.0;      J[748] = 0.0; 
      t1385 = t443*t333+t341*t446;      t1387 = -J[252];      J[749] = J[252]*t1050+t1385*t1370+t1387*t1373;      t1391 = -J[235];      J[750] = t1093*t867+J[235]*t870+t1391*t1050;      J[751] = 0.0;      t1395 = -J[237];      J[752] = t1293*t1005+J[237]*t1008+t1395*t1050;      t1403 = t695+t694+t698+t697+t703-t702+t706-t705;      t1407 = -t326*t320-t328*t323; 
      t1409 = -t703+t702;      t1413 = -t339*t333-t341*t336;      t1415 = -t706+t705;      J[753] = t1115*t867+t1186*t870+t1301*t1005+t1342*t1008-2.0*t499*t1047+t1403*t1050+t1407*t1362+t1409*t1365+t1413*t1370+t1415*t1373;      t1417 = -J[239];      t1421 = t349*t320-t328*t346;      J[754] = t1417*t1050+t1421*t1362+J[239]*t1365;      J[755] = 0.0;      J[756] = 0.0;      J[757] = 0.0; 
      t1424 = -J[243];      t1428 = t446*t333-t341*t443;      J[758] = t1424*t1050+t1428*t1370+J[243]*t1373;      J[759] = 0.0;      J[760] = 0.0;      J[761] = 0.0;      J[762] = 0.0;      J[763] = 0.0;      J[764] = t1049*t867+t1141*t870+J[244]*t1047+J[267]*t1050;      J[765] = 0.0; 
      J[766] = t1271*t1005+t1319*t1008+J[246]*t1047+J[269]*t1050;      t1446 = -J[293];      t1448 = -J[316];      t1450 = -J[477];      t1452 = -J[500];      J[767] = t1071*t867+t1159*t870+t1278*t1005+t1325*t1008+t1357*t1047-J[270]*t1050+t1446*t1362+t1448*t1365+t1450*t1370+t1452*t1373;      t1456 = -t749+t748;      t1458 = -J[271];      J[768] = J[248]*t1047+J[271]*t1050+t1456*t1362+t1458*t1365;      J[769] = 0.0; 
      J[770] = 0.0;      J[771] = 0.0;      t1462 = -t851+t850;      t1464 = -J[275];      J[772] = J[252]*t1047+J[275]*t1050+t1462*t1370+t1464*t1373;      J[773] = t1095*t867+J[258]*t870+t1391*t1047-J[258]*t1050;      J[774] = 0.0;      J[775] = t1295*t1005+J[260]*t1008+t1395*t1047-J[260]*t1050;      J[776] = J[100]*t867+J[123]*t870+J[192]*t1005+J[215]*t1008+t1403*t1047+J[261]*t1050+J[284]*t1362+J[307]*t1365+J[468]*t1370+J[491]*t1373;      t1489 = -t737-t736; 
      J[777] = t1417*t1047-J[262]*t1050+t1489*t1362+J[262]*t1365;      J[778] = 0.0;      J[779] = 0.0;      J[780] = 0.0;      t1495 = -t839-t838;      J[781] = t1424*t1047-J[266]*t1050+t1495*t1370+J[266]*t1373;      J[782] = 0.0;      J[783] = 0.0;      J[784] = 0.0;      J[785] = 0.0; 
      J[786] = 0.0;      J[787] = 0.0;      J[788] = 0.0;      J[789] = 0.0;      J[790] = t1361*t1047+t1446*t1050+J[293]*t1365;      t1505 = t749-t748+t754+t753;      t1509 = -t356*t359+t362*t364;      t1511 = -t754-t753;      J[791] = t1378*t1047+t1456*t1050+2.0*t354*t1362+t1505*t1365+t1509*t938+t1511*t941;      t1516 = t376*t359+t362*t381; 
      t1518 = -J[295];      J[792] = J[295]*t1365+t1516*t938+t1518*t941;      J[793] = 0.0;      J[794] = 0.0;      J[795] = 0.0;      J[796] = 0.0;      J[797] = 0.0;      J[798] = 0.0;      t1522 = -J[284];      J[799] = t1407*t1047+J[284]*t1050+t1522*t1365; 
      t1528 = t737+t736+t742-t741;      t1532 = -t364*t359-t362*t356;      t1534 = -t742+t741;      J[800] = t1421*t1047+t1489*t1050-2.0*t508*t1362+t1528*t1365+t1532*t938+t1534*t941;      t1536 = -J[286];      t1540 = t381*t359-t362*t376;      J[801] = t1536*t1365+t1540*t938+J[286]*t941;      J[802] = 0.0;      J[803] = 0.0;      J[804] = 0.0; 
      J[805] = 0.0;      J[806] = 0.0;      J[807] = 0.0;      J[808] = 0.0;      J[809] = 0.0;      J[810] = 0.0;      J[811] = 0.0;      J[812] = 0.0;      J[813] = t1364*t1047+t1448*t1050+J[293]*t1362+J[316]*t1365;      t1552 = -J[340]; 
      t1554 = -J[363];      J[814] = t1380*t1047+t1458*t1050+t1505*t1362-J[317]*t1365+t1552*t938+t1554*t941;      t1558 = -t776+t775;      t1560 = -J[318];      J[815] = J[295]*t1362+J[318]*t1365+t1558*t938+t1560*t941;      J[816] = 0.0;      J[817] = 0.0;      J[818] = 0.0;      J[819] = 0.0;      J[820] = 0.0; 
      J[821] = 0.0;      J[822] = t1409*t1047+J[307]*t1050+t1522*t1362-J[307]*t1365;      J[823] = J[239]*t1047+J[262]*t1050+t1528*t1362+J[308]*t1365+J[331]*t938+J[354]*t941;      t1576 = -t766-t765;      J[824] = t1536*t1362-J[309]*t1365+t1576*t938+J[309]*t941;      J[825] = 0.0;      J[826] = 0.0;      J[827] = 0.0;      J[828] = 0.0;      J[829] = 0.0; 
      J[830] = 0.0;      t1580 = -t937;      J[831] = t937*t927+t1580*t941;      J[832] = 0.0;      J[833] = 0.0;      J[834] = 0.0;      J[835] = 0.0;      J[836] = 0.0;      J[837] = t1509*t1362+t1552*t1365+J[340]*t941;      t1590 = t773+t772+t776-t775; 
      J[838] = t965*t927+t1516*t1362+t1558*t1365+2.0*t385*t938+t1590*t941;      J[839] = 0.0;      J[840] = 0.0;      J[841] = 0.0;      J[842] = 0.0;      J[843] = 0.0;      J[844] = 0.0;      J[845] = 0.0;      t1594 = -J[331];      J[846] = t1532*t1362+J[331]*t1365+t1594*t941; 
      t1601 = t763-t762+t766+t765;      J[847] = t986*t927+t1540*t1362+t1576*t1365-2.0*t517*t938+t1601*t941;      J[848] = 0.0;      J[849] = 0.0;      J[850] = 0.0;      J[851] = 0.0;      J[852] = 0.0;      J[853] = 0.0;      J[854] = t940*t927+t1580*t938-t940*t941;      J[855] = 0.0; 
      J[856] = 0.0;      J[857] = 0.0;      J[858] = 0.0;      J[859] = 0.0;      J[860] = t1511*t1362+t1554*t1365+J[340]*t938+J[363]*t941;      J[861] = t967*t927+t1518*t1362+t1560*t1365+t1590*t938-J[364]*t941;      J[862] = 0.0;      J[863] = 0.0;      J[864] = 0.0;      J[865] = 0.0; 
      J[866] = 0.0;      J[867] = 0.0;      J[868] = 0.0;      J[869] = t1534*t1362+J[354]*t1365+t1594*t938-J[354]*t941;      J[870] = J[56]*t927+J[286]*t1362+J[309]*t1365+t1601*t938+J[355]*t941;      J[871] = 0.0;      J[872] = 0.0;      J[873] = 0.0;      J[874] = 0.0;      J[875] = 0.0; 
      J[876] = 0.0;      t1628 = -t945;      J[877] = t945*t927+t1628*t949;      J[878] = 0.0;      J[879] = 0.0;      J[880] = 0.0;      J[881] = 0.0;      J[882] = 0.0;      J[883] = 0.0;      J[884] = 0.0; 
      t1633 = t793+t792+t798+t797;      t1637 = -t402*t400+t407*t405;      t1639 = -t798-t797;      J[885] = t970*t927+2.0*t397*t946+t1633*t949+t1637*t954+t1639*t957;      t1644 = t420*t400+t407*t423;      t1646 = -J[389];      J[886] = J[389]*t949+t1644*t954+t1646*t957;      J[887] = 0.0;      J[888] = 0.0;      J[889] = 0.0; 
      J[890] = 0.0;      J[891] = 0.0;      J[892] = 0.0;      J[893] = 0.0;      t1651 = t783-t782+t788-t787;      t1655 = -t405*t400-t407*t402;      t1657 = -t788+t787;      J[894] = t991*t927-2.0*t522*t946+t1651*t949+t1655*t954+t1657*t957;      t1659 = -J[380];      t1663 = t423*t400-t407*t420; 
      J[895] = t1659*t949+t1663*t954+J[380]*t957;      J[896] = 0.0;      J[897] = 0.0;      J[898] = 0.0;      J[899] = 0.0;      J[900] = t948*t927+t1628*t946-t948*t949;      J[901] = 0.0;      J[902] = 0.0;      J[903] = 0.0;      J[904] = 0.0; 
      J[905] = 0.0;      J[906] = 0.0;      J[907] = 0.0;      t1674 = -J[434];      t1676 = -J[457];      J[908] = t972*t927+t1633*t946-J[411]*t949+t1674*t954+t1676*t957;      t1680 = -t825+t824;      t1682 = -J[412];      J[909] = J[389]*t946+J[412]*t949+t1680*t954+t1682*t957;      J[910] = 0.0; 
      J[911] = 0.0;      J[912] = 0.0;      J[913] = 0.0;      J[914] = 0.0;      J[915] = 0.0;      J[916] = 0.0;      J[917] = J[57]*t927+t1651*t946+J[402]*t949+J[425]*t954+J[448]*t957;      t1692 = -t810-t809;      J[918] = t1659*t946-J[403]*t949+t1692*t954+J[403]*t957;      J[919] = 0.0; 
      J[920] = 0.0;      J[921] = 0.0;      J[922] = 0.0;      t1696 = -t953;      J[923] = t953*t927+t1696*t957;      J[924] = 0.0;      J[925] = 0.0;      J[926] = 0.0;      J[927] = 0.0;      J[928] = 0.0; 
      J[929] = 0.0;      J[930] = 0.0;      J[931] = t1637*t946+t1674*t949+J[434]*t957;      t1706 = t822+t821+t825-t824+t830+t829;      t1710 = -t433*t431+t438*t436;      t1712 = -t830-t829;      J[932] = t975*t927+t1644*t946+t1680*t949+2.0*t428*t954+t1706*t957+t1710*t1370+t1712*t1373;      t1717 = t451*t431+t438*t454;      t1719 = -J[436];      J[933] = J[436]*t957+t1717*t1370+t1719*t1373; 
      J[934] = 0.0;      J[935] = 0.0;      J[936] = 0.0;      J[937] = 0.0;      J[938] = 0.0;      J[939] = 0.0;      t1723 = -J[425];      J[940] = t1655*t946+J[425]*t949+t1723*t957;      t1730 = t807-t806+t810+t809+t815-t814;      t1734 = -t436*t431-t438*t433; 
      t1736 = -t815+t814;      J[941] = t996*t927+t1663*t946+t1692*t949-2.0*t531*t954+t1730*t957+t1734*t1370+t1736*t1373;      t1738 = -J[427];      t1742 = t454*t431-t438*t451;      J[942] = t1738*t957+t1742*t1370+J[427]*t1373;      J[943] = 0.0;      J[944] = 0.0;      J[945] = 0.0;      J[946] = t956*t927+t1696*t954-t956*t957;      J[947] = 0.0; 
      J[948] = 0.0;      J[949] = 0.0;      J[950] = 0.0;      J[951] = 0.0;      J[952] = 0.0;      J[953] = 0.0;      J[954] = t1639*t946+t1676*t949+J[434]*t954+J[457]*t957;      t1759 = -J[481];      t1761 = -J[504];      J[955] = t977*t927+t1646*t946+t1682*t949+t1706*t954-J[458]*t957+t1759*t1370+t1761*t1373; 
      t1765 = -t854+t853;      t1767 = -J[459];      J[956] = J[436]*t954+J[459]*t957+t1765*t1370+t1767*t1373;      J[957] = 0.0;      J[958] = 0.0;      J[959] = 0.0;      J[960] = 0.0;      J[961] = 0.0;      J[962] = 0.0;      J[963] = t1657*t946+J[448]*t949+t1723*t954-J[448]*t957; 
      J[964] = J[58]*t927+J[380]*t946+J[403]*t949+t1730*t954+J[449]*t957+J[472]*t1370+J[495]*t1373;      t1784 = -t842-t841;      J[965] = t1738*t954-J[450]*t957+t1784*t1370+J[450]*t1373;      J[966] = 0.0;      J[967] = 0.0;      J[968] = 0.0;      J[969] = 0.0;      J[970] = 0.0;      J[971] = 0.0;      J[972] = 0.0; 
      J[973] = 0.0;      J[974] = t1369*t1047+t1450*t1050+J[477]*t1373;      J[975] = 0.0;      J[976] = 0.0;      J[977] = 0.0;      J[978] = t1710*t954+t1759*t957+J[481]*t1373;      t1799 = t851-t850+t854-t853;      J[979] = t1385*t1047+t1462*t1050+t1717*t954+t1765*t957+2.0*t459*t1370+t1799*t1373;      J[980] = 0.0;      J[981] = 0.0; 
      J[982] = 0.0;      t1803 = -J[468];      J[983] = t1413*t1047+J[468]*t1050+t1803*t1373;      J[984] = 0.0;      J[985] = 0.0;      J[986] = 0.0;      t1807 = -J[472];      J[987] = t1734*t954+J[472]*t957+t1807*t1373;      t1815 = t839+t838+t842+t841;      J[988] = t1428*t1047+t1495*t1050+t1742*t954+t1784*t957-2.0*t540*t1370+t1815*t1373; 
      J[989] = 0.0;      J[990] = 0.0;      J[991] = 0.0;      J[992] = 0.0;      J[993] = 0.0;      J[994] = 0.0;      J[995] = 0.0;      J[996] = 0.0;      J[997] = t1372*t1047+t1452*t1050+J[477]*t1370+J[500]*t1373;      J[998] = 0.0; 
      J[999] = 0.0;      J[1000] = 0.0;      J[1001] = t1712*t954+t1761*t957+J[481]*t1370+J[504]*t1373;      J[1002] = t1387*t1047+t1464*t1050+t1719*t954+t1767*t957+t1799*t1370-J[505]*t1373;      J[1003] = 0.0;      J[1004] = 0.0;      J[1005] = 0.0;      J[1006] = t1415*t1047+J[491]*t1050+t1803*t1370-J[491]*t1373;      J[1007] = 0.0;      J[1008] = 0.0; 
      J[1009] = 0.0;      J[1010] = t1736*t954+J[495]*t957+t1807*t1370-J[495]*t1373;      J[1011] = J[243]*t1047+J[266]*t1050+J[427]*t954+J[450]*t957+t1815*t1370+J[496]*t1373; 

   return;
}

/*         Gateway Routine         */

void mexFunction(
   int             nlhs,
   mxArray  *plhs[],
   int             nrhs,
   const mxArray  *prhs[]
   )
{
   double          *f;
   double          *J;
   double          *data;
   double          *x;
   double          *param;
   double          *v;
   /* Create a matrix for the return argument */

   plhs[0] = mxCreateDoubleMatrix(  23, 1, mxREAL);

   plhs[1] = mxCreateDoubleMatrix(  23,  44, mxREAL);

   /* Assign pointers to the various parameters */

   f = mxGetPr(plhs[0]);

   J = mxGetPr(plhs[1]);

   data = mxGetPr(prhs[0]);
   x = mxGetPr(prhs[1]);
   param = mxGetPr(prhs[2]);
   v = mxGetPr(prhs[3]);
   /* Do the actual computations in a subroutine */

  IEEE14(f,J,data,x,param,v);
   return;
}