题目与上一篇博客相同,求解第二问
非线性规划具体思路看上一篇
11行开始 使用蒙特卡罗法进行预估近似解 循环10的六次方次使其具有统计意义 生成0到n-1(参考题目设定大致范围)的随机数作为随机生成的解,由于后续需要进行等式约束条件判断(较难实现)为了提高效率前十二个变量随机生成整数,后四个变量根据题目随机生成0到9的数
在循环过程中先判断不等式约束是否成立,再判断等式约束是否成立(在判断语句中嵌套相关值的计算可以提高效率,减少无意义的计算)最后记录每一次循环的更优值,最终求得近似理想解,将近似理想解作为fmincon的初始值有助于减小非线性规划求解的误差,可以看出最终解明显是优于上一篇中的求解的
以下为所有变量