Re: [繪圖] 聯立不等式繪圖
※ 引述《JeremyKSKGA (GIANTS!!!)》之銘言:
: 從來沒做過聯立不等式繪圖...囧
: 譬如我想做以下的聯立不等式
: 2 x^2 + 4 y^2 - 4 x^2 * y^2 < 1
: 1 - x^2 - 2 y^2 > 0
: 請問我該如何求出x對y作圖
: 以及指令求出x和y的最大值??
利用fmincon函數求解
------------------------------------------------------
function pttex140
% 解法
[xy target] = fmincon(@obj,[0;0],[],[],[],[],[],[],@myfun);
x = xy(1)
y = xy(2)
maxfun = target*(-1)
% 驗算
%
a = 2.*x.^2 + 4.*y.^2 - 4.*(x.^2).*(y.^2)-1
b = x.^2 + 2.*y.^2 -1
function [f,feq] = myfun(x)
f(1) = 2.*x(1).^2 + 4.*x(2).^2 - 4.*(x(1).^2).*(x(2).^2)-1;
f(2) = x(1).^2 + 2.*x(2).^2 -1;
feq = [];
function f = obj(x)
f = -x(1)-x(2);
-------------------------------------------------------------
Active inequalities (to within options.TolCon = 1e-06):
lower upper ineqlin ineqnonlin
1
x =
0.5880
y =
0.3434
maxfun =
0.9313
a =
4.2277e-10
b =
-0.4185
先看看. 有問題再討論.
另外繪圖的部分就自己疊代後繪圖
這邊不著墨太多
--
1.MATLAB programming 2.ASPEN process simulation package
3.FORTRN programming 4.Advance Engineering Mathematics
5.Process Control Theory
6.Chemical Engineering Basic Theory(Kinetic.thermodynamics.transport)
7.Numerical Method and Analysis
8.MATLAB Toolbox.Simulink system basic design
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 113.186.203.142
討論串 (同標題文章)
本文引述了以下文章的的內容:
完整討論串 (本文為第 2 之 2 篇):
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章