数学建模学习-线性规划

最新推荐文章于 2024-09-16 13:36:19 发布

Xuuuuuuuuuuuuuuuuuuu

最新推荐文章于 2024-09-16 13:36:19 发布

阅读量128

点赞数

文章标签：数学建模学习

本文链接：https://blog.csdn.net/qq_30796747/article/details/132756025

版权

预备知识

在这里插入图片描述

1.待解决的问题

要解决的问题
表1.1投资的相关数据

s_i	r_i(%)	q_i(%)	p_i(%)	u_i(元)
s₁	28	2.5	1	103
s₂	21	1.5	2	198
s₁	23	5.5	4.5	52
s₁	25	2.6	6.5	40

2.符号规定

在这里插入图片描述

3.基本假设

在这里插入图片描述
好的假设可以节省很多力气

4.模型分析与建立

在这里插入图片描述
目标函数为
$\begin{cases} max \sum\limits_{i=0}^{n}(r_i -p_i)x_i \ 投资收益 \\ min \ max\{q_ix_i\} \ 最大风险 \end{cases}$
约束条件为
$\begin{cases} \sum\limits_{i=0}^{n}(1+p_i)x_i = M \ 资金总和 \\ x_i>=0, i=0,1,...,n \ 每一个项目的投资金额 \end{cases}$

在这里插入图片描述

以这种方案构建模型

min f =(-0.05,-0.27,-0.19,-0.185,-0.185)x ₀

$\begin{cases} x_0 + 1.01x_1 + 1.02x_2 + 1.045x_3 + 1.065x_4 = 1 \\ 0.025x_1<=a \\ 0.015x_2<=a \\ 0.055x_3<=a \\ 0.026x_4<=a \\ x_i>=0(i=0,1,...,4) \end{cases}$
a是任意给的风险度所以从0开始步长为Δa=0.001进行循环搜索

代码如下

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')