Question

使用matlab画出二元一次不等式组图：3x1+x2<=60000,2x1+x2<=50000,x1<=15000,x1、x2>=0.求代码。

如题

Answer 1

用MATLAB画满足不等式约束的区域，我至今没发现什么好办法。

对于题主现在的特殊情况（约束均为线性），应该稍微好一些。

我编写了一段代码，实现思路不算太好，贴出来抛砖引玉吧。

照理说应该根据各条边交点的位置自动判断坐标范围，但有点嫌麻烦，先凑合着用吧。

 

参考代码：

% 把线性不等式约束统一写成 Ax<=b 的形式
A = [3 1; 2 1; 1 0; -1 0; 0 -1];
b = [60000; 50000; 15000; 0; 0];
n1 = 600;
n2 = 400;
x1 = linspace(0,20000,n1);
x2 = linspace(0,60000,n2);
[x1,x2] = meshgrid(x1,x2);
X = ones(n2,n1);
L = any( A*[x1(:)'; x2(:)']>repmat(b,1,n1*n2) );
X( L ) = NaN;
pcolor(x1,x2,X)
shading interp

 

===============

PS: 晕。。。提交后发现已经有人回答并被采纳了，早知道就不献丑了

Answer 2

思路如下：

先画出两条直线，然后用fill命令对可行域进行填充。

代码如下：

clear all;
close all;
clc;
h1=ezplot('60000-3.*x-y',[0,25000,0,60000]);
hold on
h2=ezplot('50000-2.*x-y',[0,25000,0,60000]);
set(h1,'Color','blue')
set(h2,'Color','green')
title('不等式组图')
x=0:100:10000;
y1=zeros (1,size(x,2));
y2=-2.*x+50000;
X=[x,fliplr(x)];
Y=[y1,fliplr(y2)];
fill(X,Y,'r');
x=10000:100:15000;
y1=zeros (1,size(x,2));
y2=-3.*x+60000;
X=[x,fliplr(x)];
Y=[y1,fliplr(y2)];
fill(X,Y,'r');
xlabel('x1')
ylabel('x2')

生成图像如下：

追问

谢谢，能不能在这基础上再画一条H=250x1+100x2。还有注明所有线段的名称。

追答

最后这条直线（H=250x1+100x2）也影响约束区域么？

它有没有不等号方向 还是就单纯的显示

Answer 3

室友你要加油