题主的问题属于极值最优化问题,可以用fmincon()函数求解。求解步骤:
1、创建自定义目标函数,fmincon_fun(k)
y=int((x-800)/200,0,q);
eq1=10000*(q-y)+500*y;
eq2=1000*q*(1+2.25/100)-3000*q*(1+2.25/100);
eq3=6000*q+30*(5*q-a*q-400+5*q-6*q-400)+1000*a^2+1000*b^2;
maxPI=-(eq1-eq2-eq3);
2、确定k0初值,k0=[0,0,0]
3、确定q、a、b的上下限,lb=[800,0,0];ub=[1000,5,5];
4、使用fmincon函数求解,q、a、b的值
5、运行代码,得到结果