Re: [繪圖] 聯立不等式繪圖

看板MATLAB作者 (外派越南工作者)時間11年前 (2013/12/09 22:42), 編輯推噓0(000)
留言0則, 0人參與, 最新討論串2/2 (看更多)
※ 引述《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
文章代碼(AID): #1IfTRIXz (MATLAB)
討論串 (同標題文章)
文章代碼(AID): #1IfTRIXz (MATLAB)