Matlab非线性规划——求解周长固定四边形的最大面积

问题提出：

已知四边形的周长固定为4cm，四边形边长为a，b，c，d，对角线长度为x。求四边形面积最大时的每一条边长和面积。编写MATLAB程序进行求解。

问题分析：

四边形ABCD问题可将其划分为🔺ABC和🔺ACD两个三角形的问题

如上图所示四边形ABCD的面积公式如下所示：

根据三角形两边之和大于第三边以及周长固定为4cm建立约束条件：

fmincon函数用法参见： https://blog.csdn.net/qq_41301570/article/details/112428785

---

编写代码：

（1）编写约束函数的function，注意fmincon函数为求解最小值的函数，所以首先要在表示四边形面积的函数前加上负号。建立的第一个fun1.m如下：

function f=fun3(x)
%a表示为x(1),b表示为x(2);c表示为x(3);d表示为x(4);x表示为x(5);
s=(x(1)+x(2)+x(3)+x(4))/2;
f=-sqrt((s-x(1))*(s-x(2))*(s-x(3))*(s-x(4)));

（2）然后编写约束条件的function，注意不等式要写成小于等于0的形式，然后依次写出左侧表示的式子。建立的第二个fun2.m如下：

function [g,h]=fun4(x)
g=[-x(1)-x(2)+x(5)
   -x(1)+x(2)-x(5)
    x(1)-x(2)-x(5)
   -x(3)-x(4)+x(5)
   -x(3)+x(4)-x(5)
    x(3)-x(4)-x(5)];
h=[x(1)+x(2)+x(3)+x(4)-4];

（3）编写主函数

clear
clc
%fun1 fun2
[x,y]=fmincon('fun1',rand(5,1),[],[],[],[],zeros(5,1),[],'fun2');
x  %输出x矩阵，a，b，c，d，x的值
y=-y %fun1中目标函数加了负号，这里加回来。

---

运行结果：

最后：

不定期发布一些matlab设计内容,敬请期待。包括但不限于如下内容:信号处理、通信仿真、gui设计、matlab appdesigner，simulink仿真。有任何有关MATLAB的问题可加QQ：2802009708进行咨询。或扫码进行添加。