如果matlab有不等式工具箱,可以用linpro函数,具体用法百度一下
我没有这个工具箱,所以写了一段计算方法
k=1;
for k= 1:3600/9 %遍历
x(k) = k; % 甲种K件
tp = [(3600-9*k)/4, (2000-k*4)/5,(3000-k*3)/10];
y(k)= min(tp);% 乙种y件
z(k)= 70*k+120*y(k);%总利润
end
plot(x(1:400),z(1:400))
[mmax jia] = max(z)
%mmax = 42800
% jia = 200