【MATLAB编程实例练习】-(17)绘制含不等式的三维图形

需绘制的三维图形方程如下:
在这里插入图片描述
实现代码:

clear
clc
[x,y] = meshgrid(0:0.01:2);
z = -x.^2-y.^2+2*y+2*x.*y-x;
z(x+y>3) = NaN;
z(x>y) = NaN;
mesh(x,y,z);

结果:
在这里插入图片描述

2参数一次不等式组围成的区域的图表示-patch_NoEqual.m 如标题所示,这里要解决的是给出一个2元一次不等式组,画出其可行性区域,也就是类似与LMI中的可行性解的区域。        在File Exchange上也有实现类似功能的函数,但是我这个人比较懒,懒得看别人写的东西,因为可读性实在是太差,而且一般情况下鲁棒性也太差,只能适用于他们自己的应用。举个例子,如果不等式组中存在这么两个不等式:0.001x 0.002y 10<0;  3e16 x 9e15 y 1e16>0;我试了一下他们写的函数,会出现大数吃小数的现象,根本就画不出要求的区域。 函数说明: function patch_NoEqual %************************************************************************** % 功能:二元线性规划的可行性区域。将代数不等式组所围成的区域用patch表示出来,仅适用于线性不等式 % vtr_NonEuations: sym类型的不等式左边的表达式,右边都为0,为列向量 % flag:为一个表示vtr_NonEuations是大于零或者小于零的行向量,1表示大于零,-1表示小于零 % varargin{1}:Facecolor值 % varargin{2}:Edgecolor值 % varargin{3}:第三维的坐标值,我的应用中只需考虑为一个常数的情况;不为常数的情况下,只需保证它的维数与最后可行点的维数相同即可。 如何使用: 1) 不等式组的表示:所有不等式均采用符号表达式输入,vtr_NonEquations中存储不等式的左边,flag中标记不等号。 2)区域颜色设置:varargin{1,2}中设定区域颜色,Facecolor是区域的颜色;Edgecolor是边界的颜色。 3) varargin{3}:设置为0即可。 例子:%求解 x-y>0 x y<1  x 2y>1这3个的可行区域            syms x y real;          vtr_NonEquations=[x-y;x y-1;x 2*y-1];          flag = [1 -1 1]; patch_NoEqual; 优点:1)完全采用符号运算,不用考虑精度问题造成的大树吃小数              2)提前对平行边界、相同边界和多条边界交于一点这3种情况进行处理,减少了计算量和错误。 缺点:1) 由于我的应用中不涉及半开平面的情况,所以没有对其考虑。              2)只考虑了<和>两种情况,没有考虑不等式符号<=和>=,也就是没有考虑是否包边界的问题。 解决方案:1)可以设定一下xlim,ylim,然后加在原不等式组的基础上追加4个不等式xlim的下限<x<xlim的上限,ylim的下限<y<ylim的上限。                      2)对与每一个不等式可以追加一个flag,让它来决定是否包边界,然后再画图的时候采用不同的linestyle
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值