数学建模与算法分析第一次打卡

Hello大家好，经过两天的学习，我弄完了第一章的部分，下面说说我的心得以及源代码。

---

心得体会：

第一章讲的是线性规划，总的来说就是给你一个目标和一些条件来求解最优解的过程。
线性规划分为单目标规划和多目标规划，一般来说，都是讲多目标规划转化为单目标规划来求解。
这里介绍三种方法：
①最常见的一种是使用系数s，将多个目标整合为一个目标，然后对整合后的目标进行线性规划求解。这有点类似于高数中的拉格朗日求解过程。这一种相对比较容易理解，但是系数s的选取是关键，决定了规划效果的好坏。笔者觉得，可以先使用循环将多个效果绘制在图上，然后来选取最好的那一个。
②二三两种方法都具有一定的技巧性，需要根据题目来确定，这里不过多展开，可以参考书上的解法，多加体会。

我觉得这一章还是比较容易理解的，不是非常难，主要就是需要注意目标函数要转化为min的形式，而且主要矩阵和列向量，在编写程序的时候不要搞混了就可以了。

内容不是很多，多加思考，是可以掌握的！

---

源代码：

f=[-2;-3;5];%价值向量，列向量
a=[-2,5,-1;1,3,1];%线性不等式约束 ， 矩阵
b=[-10;12];%线性不等式约束 ， 矩阵
aeq=[1,1,1];%线性等式约束  列向量
beq=7;%线性等式约束  列向量
lb=[0;0;0];%决策向量的下界向量
%[x,y]=linprog(f,a,b,aeq,beq,zeros(3,1));
[x,y]=linprog(f,a,b,aeq,beq,lb);
x,y=-y;

f=[2;3;1];
a=[-1,-4,-2;-3,-2,0];
b=[-8;-6];
lb=[0;0;0]
[x,y]=linprog(f,a,b,[],[],lb)

clc,clear
c=1:4;
c=[c,c]';%构造价值列向量
a=[1 -1 -1 1;1 -1 1 -3;1 -1 -2 3];
a=[a,-a];%构造变换后新的系数矩阵
b=[-2 -1 -0.5]';
[y,z]=linprog(c,a,b,[],[],zeros(8,1));%这里没有等式的约束，对应的矩阵为空矩阵
x=y(1:4)-y(5:end)%变换到原问题的解

clc,clear
a=0;
hold on
while a<0.05
    c=[-0.05,-0.27,-0.19,-0.185,-0.185];
    A=[zeros(4,1),diag([0.025,0.015,0.055,0.026])];
    b=a*ones(4,1);
    Aeq=[1,1.01,1.02,1.045,1.065];
    beq=1;
    LB=zeros(5,1);
    [x,Q]=linprog(c,A,b,Aeq,beq,LB);
    Q=-Q;
    plot(a,Q,'*k');
    a=a+0.001;
end
xlabel('a'),ylabel("Q")

---

本次讲解到此结束，谢谢大家，记得点赞哟！