求解二次规划问题
c++ osqp例子
c++ osqp例子
以下是一个使用C++和OSQP库求解二次规划问题的例子：
1. 安装OSQP库并在C++代码中包含osqp头文件：#include "osqp.h"
2. 创建二次规划问题的数据。
例如，考虑以下二次规划问题：
min 0.5*x^2 + 2*y
s.t. x + y <= 3
x, y >= 0
在这个例子中，目标函数是一个二次函数，约束条件是线性的。
你可以根据实际问题修改目标函数和约束条件。
要使用OSQP库求解这个问题，你需要创建一个OSQP问题对象，并将目标函数和约束条件传递给它。
然后，调用OSQP的求解函数求解问题。
最后，你可以检查求解结果，以确定问题是否有解以及解的质量如何。
请注意，这只是一个简单的例子，实际问题可能会更加复杂。
你可能需要根据问题的具体情况进行适当的修改和调整。
有关OSQP库的更详细信息和其他示例，请参考OSQP的官方文档或其他相关资源。
qp算法的c语言实现
qp算法的c语言实现简介qp算法（Quadratic Programming Algorithm）是一种用于求解二次规划问题的优化算法。
二次规划问题是指目标函数为二次型，约束条件为线性等式或不等式的优化问题。
qp算法通过迭代寻找目标函数的极小值点，以及满足约束条件的最优解。
在本文中，我们将介绍如何使用C语言来实现qp算法，并提供一个简单的示例来说明其使用方法。
算法原理qp算法通过将二次规划问题转化为一系列线性规划子问题来求解。
具体而言，它将原始问题转化为一个等价的对偶问题，并利用对偶问题中的特性来进行求解。
qp算法基于拉格朗日乘子法和KKT条件（Karush-Kuhn-Tucker conditions）。
拉格朗日乘子法通过引入拉格朗日乘子来构建一个增广拉格朗日函数，将原始问题转化为一个无约束优化问题。
KKT条件则是指在最优解处，原始问题和对偶问题满足一系列约束条件。
qp算法主要分为两个阶段：预处理阶段和迭代阶段。
预处理阶段主要用于计算一些必要的矩阵和向量，并进行一些初始化操作。
迭代阶段则是通过不断更新变量来逐步接近最优解。
C语言实现下面我们将介绍如何使用C语言来实现qp算法。
首先，我们需要定义一些必要的数据结构和变量。
#include <stdio.h>#include <stdlib.h>typedef struct {double *H; // Hessian矩阵double *f; // 线性项double *A; // 不等式约束系数矩阵double *b; // 不等式约束常数向量int n; // 变量个数int m; // 不等式约束个数} QPProblem;typedef struct {double *x; // 最优解向量double obj; // 目标函数值} QPSolution;在定义了必要的数据结构后，我们可以开始实现qp算法的各个函数。
求解二次规划问题的拉格朗日及有效集方法
求解二次规划问题的拉格朗日及有效集方法——最优化方法课程实验报告学院：数学与统计学院班级：硕2041班姓名：王彭学号：3112054028指导教师：阮小娥同组人：钱东东求解二次规划问题的拉格朗日及有效集方法求解二次规划问题的拉格朗日及有效集方法摘要二次规划师非线性优化中的一种特殊情形，它的目标函数是二次实函数，约束函数都是线性函数。
由于二次规划比较简单，便于求解（仅次于线性规划），并且一些非线性优化问题可以转化为求解一些列的二次规划问题，因此二次规划的求解方法较早引起人们的重视，称为求解非线性优化的一个重要途径。
二次规划的算法较多，本文仅介绍求解等式约束凸二尺规划的拉格朗日方法以及求解一般约束凸二次规划的有效集方法。
关键字：二次规划，拉格朗日方法，有效集方法。
- 1 -《最优化方法》课程实验报告- 2 - 【目录】摘要........................................................................................................................... - 1 -1 等式约束凸二次规划的解法............................................................................... - 3 -1.1 问题描述.................................................................................................... - 3 -1.2 拉格朗日方法求解等式约束二次规划问题............................................ - 3 -1.2.1 拉格朗日方法的推导...................................................................... - 3 -1.2.2 拉格朗日方法的应用...................................................................... - 4 -2 一般凸二次规划问题的解法............................................................................... - 5 -2.1 问题描述.................................................................................................... - 5 -2.2 有效集法求解一般凸二次规划问题........................................................ - 6 -2.2.1 有效集方法的理论推导.................................................................. - 6 -2.2.2 有效集方法的算法步骤.................................................................. - 9 -2.2.3 有效集方法的应用........................................................................ - 10 -3 总结与体会......................................................................................................... - 11 -4 附录..................................................................................................................... - 11 -4.1 拉格朗日方法的matlab程序................................................................. - 11 -4.2 有效集方法的Matlab程序 .................................................................... - 11 -求解二次规划问题的拉格朗日及有效集方法- 3 -1 等式约束凸二次规划的解法1.1 问题描述我们考虑如下的二次规划问题⎪⎩⎪⎨⎧=+b Ax t s x c Hx x T T ..,21min (1.1) 其中n n R H ⨯∈对称正定，n m R A ⨯∈行满秩，n R x c,∈，m R b ∈。
裘布依公式范文
裘布依公式范文裘布依公式（also known as Juebu Yi formula）是一个用来解决二次多项式最优化问题的公式。
它最早由英国数学家艾伦·喬治·裘布依（Alan George Joubert）在1911年提出，得到了广泛的应用和研究。
裘布依公式在工程、经济学和运筹学等领域都有重要的应用。
x = \frac{-b ± \sqrt{b^2-4ac}}{2a}\]其中，a、b和c都是常数，且$a \neq 0$。
这是一个二次方程的解析解，可以得到方程的根x。
裘布依公式基于这个二次方程的特点，对于给定的a、b和c，可以通过计算出根的值来解决相应的最优化问题。
\underset{x}{\text{minimize}} \quad \frac{1}{2} x^TAx + bx \]其中，A是一个对称矩阵，x是一个向量，b是一个常数。
约束条件可以表示为：\begin{align*}Ax &\geq b \\Cx&=d\\Fx &\leq g\end{align*}\]其中，C和F是矩阵，d和g是向量。
裘布依公式可以用来求解这个二次规划问题的最优解。
具体来说，裘布依公式可以通过将二次规划问题转化为一个等价的线性规划问题来解决。
线性规划问题是一个目标函数为线性函数，约束条件为线性函数的最优化问题。
通过裘布依公式，可以将二次规划问题转化为一个等价的线性规划问题，然后利用线性规划的常规方法来求解最优解。
裘布依公式的应用不仅限于二次规划问题。
在经济学中，裘布依公式可以用来求解最大化利润的方程，最小化成本的方程等。
在工程学中，裘布依公式可以应用于优化问题的设计和优化，如最优化控制问题、最优化参数估计问题等。
总之，裘布依公式是一个重要的数学工具，可以用来解决二次多项式最优化问题。
它的应用广泛，并且不仅限于二次规划问题，还可以推广到更高次的多项式。
裘布依公式为解决最优化问题提供了一种有效的方法和工具。
二次规划基本介绍
二次规划基本介绍二次规划（Quadratic Programming，简称QP）是数学规划的一种特殊形式，它的目标函数是二次函数，约束条件是线性函数。
在实际应用中，二次规划被广泛应用于经济学、运筹学、工程学等领域，具有重要的理论和实际意义。
二次规划的一般形式可以表示为：$$\begin{aligned}\min_{x} \quad & \frac{1}{2} x^T Q x + c^T x \\\text{s.t.} \quad & Ax \ge b \\&Cx=d\end{aligned}$$其中，$x$是优化变量，$Q$是一个对称正定的矩阵，$c$、$b$、$d$都是列向量，$A$、$C$是约束矩阵。
在约束条件中，$Ax \ge b$表示一组不等式约束，$Cx = d$表示一组等式约束。
二次规划的优化目标是寻找满足约束条件的$x$，使得目标函数最小。
目标函数由两部分组成，一部分是二次项，一部分是线性项，其中$Q$是二次项的系数矩阵，$c$是线性项的系数向量。
由于$Q$是一个对称正定矩阵，所以二次项是凸的，使得问题具有良好的性质。
二次规划的解法有多种方法，以下介绍其中两种常用的方法：内点法和激活集方法。
内点法是一种迭代求解二次规划问题的方法。
它通过将二次规划问题转化为一系列等价的线性规划问题来求解。
在每一次迭代中，内点法通过将问题的方向限制在可行域的内部，逐渐逼近最优解。
使用内点法求解二次规划问题的一个优点是，可以在多项式时间内找到最优解，尤其适用于大规模的问题。
激活集方法是一种基于约束的求解方法。
它通过不断修改约束条件，从而求解二次规划问题。
在每一次迭代中，激活集方法选取一个子集，称为“激活集”，包含满足等式约束、不等式约束等的约束条件。
然后通过解析方法或数值方法求解这个子问题，得到对应的最优解。
该方法的优点是，可以很好地处理不等式约束和等式约束，并且收敛性良好。
除了内点法和激活集方法，还有其他的求解方法，如：序列二次规划、信赖域算法、光滑方法等。
MATLABquadprog函数求解二次规划问题
MATLABquadprog函数求解⼆次规划问题【例】求如下⼆次规划问题。
【分析】⾸先应该把⽬标函数表⽰成如下矩阵形式：这⾥要细说⼀下如何写成矩阵形式。
⾸先，向量x是很容易写出的，因为f(x)包含两个变量x1和x2，因此其次，向量f只与两个变量x1和x2的⼀次项有关，所以f T x=-2x1-6x2，因此最后，矩阵H只与两个变量x1和x2的⼆次项有关，所以，这⾥要注意的是不同于⼆次型，这⾥有个系数1/2，所以矩阵H的元素是⼆次型中的矩阵元素⼤⼩的两倍。
给出⼀个规律：设矩阵H第i⾏第j列的元素⼤⼩为H(i,j)，⼆次项x i x j的系数为a(i,j)，则本例中，，这是由于x1的平⽅项（即x1x1）系数为1/2，所以第1⾏第1列的元素为1=2*(1/2)，x2的平⽅项（即x2x2）系数为1，所以第2⾏第2列的元素为2=2*1，x1x2项（即x2x1）的系数为-1，所以第1⾏第2列和第2⾏第1列的元素均为-1。
⽬标函数搞定之后，下⾯来看约束条件部分，约束条件应该写成如下形式：本例中约束条件只有不等式约束，因此A eq和b eq为空，对于A和b很容易就可以得出来：⽽约束条件中对变量x1和x2只给出下限，没有给上限，因此ub为空，得到了所有的参数，将参数输⼊MATLAB，编程如下：（代码是直接在Command Window中⼀⾏⼀⾏录⼊的，所以每⾏前⾯有符号“>>”）>> H = [1 -1; -12];>> f = [-2; -6];>> A = [11; -12; 21];>> b = [2; 2; 3];>> lb = [0; 0];>> [x,fval,exitflag,output,lambda] = quadprog(H,f,A,b,[],[],lb)输出以下结果：Warning: Large-scale algorithm does not currently solve this problem formulation,using medium-scale algorithm instead.> In quadprog at 291Optimization terminated.x =0.66671.3333fval =-8.2222exitflag =1output =iterations: 3constrviolation: 0algorithm: 'medium-scale: active-set'firstorderopt: []cgiterations: []message: 'Optimization terminated.'lambda =lower: [2x1 double]upper: [2x1 double]eqlin: [0x1double]ineqlin: [3x1 double]参考⽂献：【1】孙⽂瑜, 徐成贤,朱德通.最优化⽅法(第⼆版)[M]. 北京:⾼等教育出版社, 2010.【2】龚纯，王正林. 精通MATLAB最优化计算[M].北京: 电⼦⼯业出版社,2009.【3】lnsunqingshen, 464518439., 百度知道，2011-06-20.【4】李明强.⼏类特殊凸⼆次规划问题的求解算法研究[D].⼭东科技⼤学,2013 .【5】于绍慧.边界约束凸⼆次规划的求解[D].南京航空航天⼤学,2005.。
quadprog函数的介绍和应用
二次规划问题1.二次规划及其基本思想二次规划问题是最简单的非线性规划问题，它是指约束为线性，目标函数为二次函数的优化问题，这类优化问题在非线性规划中研究得最早，也研究得最成熟。
二次规划迭代法的基本思想是把一般的非线性规划问题转化为一系列二次规划问题进行求解，并使得迭代点能逐渐向最优点逼近，最后得到最优解。
如果某非线性规划的目标函数为自变量的二次函数，约束条件全是线性函数，就称这样规划为二次规划。
其数学模型为：2.二次规划问题的数学模型⎪⎩⎪⎨⎧≤≤=≤+ubx lb beq x Aeq bAx t s x f Hx x TT x ·..21min ，式中，H ，A 和Aeq 为矩阵；f,b, beq, lb, ub, 和x 为向量。
3.QuadProg 函数quadprog 函数可以求解二次规划问题。
quadprog 函数的几种调用格式：•x=quadprog(H,f,A,b)这个函数的功能是：用来解最简单，最常用的模型：x f Hx x T T +21Subject tobAx ≤•x=quadprog(H,f,A,b,Aeq,beq)仍然求解上面的问题，但添加了等式约束条件Aeq*x=beq。
•x=quadprog(H,f,A,b,lb,ub,)定义设计变量的下届Ib和上界ub,使得lb<=x<=ub。
•x=quadprog(H,f,A,b,lb,ub,x0)功能同上，并设置初值x0。
•x=quadprog(H,f,A,b,lb,ub,x0,options)根据options参数指定的优化参数进行最小化。
以上通用格式为x=quadprog(problem)√problem 结构字段说明字段说明字段说明H二次规划中的二次项矩阵lb自变量下界约束f二次规划中的一次项向量ub自变量上界约束Aineq线性不等约束的系数矩阵x0初始点bineq线性不等约束的右端向量solver求解器，为“quadprog”Aeq线性等式约束的系数矩阵options options结构beq线性等式约束的右端向量•[x,fvaI]=quadprog(H,f,A,b)这个函数的功能是，返回解x 处的目标函数值，即二次规划的极值fval=•[x,fvaI,exitfIag]=quadprog(H,f,A,b)返回exitfIag 参数，描述计算的退出条件。
二次规划资料
向。
内点法的改进
• 修正内点法：引入正则化项，提高内点法的稳定性和收敛性。
• 梯度投影法：利用梯度的投影性质，简化内点法的计算。
• 并行内点法：利用多核处理器并行计算，提高计算速度。
修正牛顿法
修正牛顿法原理
• 基本思想：引入正则化项，使得海森矩阵具有更好的条件数。
• 更新公式：^(k+1) = ^k - ^(-1)^k - ^(-1)，其中为步长因子。
• 敏感性分析图：绘制模型结构与最优解的关系图，直观
的可行域，从而影响最优解的值和位置。
展示模型结构变化对最优解的影响。
06
二次规划问题的拓展与推广
多目标二次规划问题
多目标二次规划问题
• 定义：多目标二次规划问题是一类求解多个目标函数的二次规划问题，目标函数
之间可能存在冲突或权衡。
• 决策变量：多目标二次规划问题需要求解一组决策变量的最优值。
非线性二次规划问题
• 定义：非线性二次规划问题是一类目标函数或约束条件为非线性函数的二次规划
问题。
• 决策变量：非线性二次规划问题需要求解一组决策变量的最优值。
• 目标函数：非线性二次规划问题的目标函数是一个非线性二次多项式函数，通常
表示为最小化形式。
非线性二次规划问题的求解方法
• 转化为线性问题：通过变量替换或线性化方法，将非线性二次规划问题转化为线性
参数变化对最优解的影响
敏感性分析的方法
• 目标函数系数变化：目标函数系数的变化会影响最优解
• 参数扫描：遍历参数取值范围，观察最优解的变化情况。
的值和位置。
• 敏感性分析图：绘制参数与最优解的关系图，直观展示
• 约束条件变化：约束条件的变化会影响最优解的可行域，
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验2求解二次规划问题LINDO 可以求解二次规划（QP ）问题。
例如：⎪⎩⎪⎨⎧<=+>++-+=7.011.19.02.1..4.03min 22y y x y x t s yxy y x f由LAGRANGE 乘子法，得()()()7.011.19.02.14.0322-+-++-+-+-+y C y x B y x A y xy y x ，分别对x 、y 求偏导，得到两个约束条件：4.09.0202.16->++-->+--C B A x y B A y x在LINDO 中输入下列命令： MIN X+Y+A+B+C ST6X-Y-1.2A+B>02Y-X-0.9A+B+C>-0.4 1.2X+0.9Y>1.1 X+Y=1 Y<0.7 END QCP 4注释：MIN X+Y+A+B+C 一句只代表变量的出场顺序；QCP 4 一句代表前4行不是原问题真正的约束，原问题真正的约束从第5行开始。
LINDO 运行后输出以下结果：STATUS OPTIMALQP OPTIMUM FOUND AT STEP 7OBJECTIVE FUNCTION V ALUE1) 1.355556V ARIABLE V ALUE REDUCED COST X 0.666667 0.000000 Y 0.333333 0.000000A 10.888889 0.000000B 9.400000 0.000000C 0.000000 0.366667ROW SLACK OR SURPLUS DUAL PRICES2) 0.000000 -0.6666673) 0.000000 -0.3333334) 0.000000 -10.8888895) 0.000000 9.4000006) 0.366667 0.000000NO. ITERATIONS= 7这个结果说明：LINDO求解此二次规划问题（QP）共用7步迭代得到最优解fmin = 1.355556，X = 0.666667，Y = 0.333333。
第5个松弛变量取值0.366667，其它松弛变量都取0值，即，这个最优解使得前4个约束条件都取等号；其对偶问题的最优解（影子价格）DUAL PRICES为Y1 = -0.666667，Y2 = -0.333333，Y3 = -10.888889，Y4 = 9.4，Y5 = 0。
农户生产的优化模型本文内容取自生产实践，豫东一个普通农户，该农户所在地区的农业生产条件、气候状况属于中等。
下列各变量的假设均建立在农村一般农业生产条件、气候状况之上。
假设（面积单位：亩）：X1 = 用于完成上缴国家任务的小麦一年总种植面积X2 = 用于生产、生活的小麦一年总种植面积X3 =用于生产、生活的油菜一年总种植面积X4 =用于生产、生活的红薯一年总种植面积X5 =用于完成上缴国家任务的棉花一年总种植面积X6 =用于生产、生活的棉花一年总种植面积X7 =用于完成上缴国家任务的玉米一年总种植面积X8 =用于生产、生活的玉米一年总种植面积X9 =用于生产、生活的芝麻一年总种植面积X10 =用于生产、生活的花生一年总种植面积X11 =用于生产、生活的大豆一年总种植面积X12 =用于生产、生活的西瓜一年总种植面积X13 =用于生产、生活的番茄一年总种植面积X14 =用于生产、生活的白菜一年总种植面积X15 =用于生产、生活的辣椒一年总种植面积X16 =用于生产、生活的茄子一年总种植面积X17 =用于生产、生活的豆角一年总种植面积X18 =用于生产、生活的萝卜一年总种植面积F ( X1，X2，。
，X18 ) = 该农户一年总的纯收入根据上表情况，得到：小麦：600X1≥780，600X2≥4380，4月底收获后再种其它作物，有X1 + X2+X3 =利润 0.35（600X 1-780）+ 0.32*780 - 75X 1 = 135X 1-30； 0.35（600X 2-4380）- 75X 2 = 135X 2-1533；油菜：200X 3 ≥30，利润 200X 3 –30 –50X 3；= 150 X 1-30；红薯：2000 X 4 ≥600+1500，利润 0.19（2000X 4 –2100）– 80X 4 = 300 X 4-399； 棉花：100X 5≥300，10%*100X 6≥12，利润 2.95*100X 5-80X 5-80X 6 = 215 X 5-80 X 6 ； 玉米：700X 7≥600，700X 8≥2600，利润 0.27（700X 7-600）+0.22*600-36X 7 = 153X 7 –30, 0.27（700X 8-600）-36X 8 = 153X 8 –702,芝麻：40%*120X 9≥6，利润 2（120X 9-6/0.4）- 20X 9 = 220 X 9-30； 花生：利润 0.90*300X 10 - 20X 10 = 250X 10 ； 大豆：利润 0.90*150X 11 - 10X 11 = 125X 11 ； 西瓜：利润 0.12*3500X 12 - 100X 12 = 320X 12 ； 番茄：利润 0.08*10000X 13 - 120X 13 = 680X 13 ； 白菜：利润 0.05*10000X 14 - 120X 14 = 380X 14 ； 辣椒：利润 5*200X 15 - 90X 15 = 910X 12 ；茄子：利润 0.06*10000X 16 - 100X 16 = 500X 12 ； 豆角：利润 0.20*6250X 17 - 120X 17 = 1130X 17 ； 萝卜：利润 0.05*10000X 18 - 50X 18 = 450X 18 ；得到模型如下：⎪⎪⎪⎩⎪⎪⎪⎨⎧≥==++≥==≥≥≥≥≥≥-++++++++++++-++++=∑=0,12,12321,819,7268,767,12610,35,21420,3320,438260,78160..4.2015184501711301650015910143801368012320111251025092208153715368052154300315021351135max 184xj xj x x x x x x x x x x x x t s x x x x x x x x x x x x x x x x x x f j在LINDO 中送入以下命令：MAX 135X1+135X2+150X3+300X4+215X5-80X6+153X7+153X8+220X9+250X10+125X11+320X12+680X13+380X14+910X15+500X16+1130X17+450X18 ST60X1>78 60X2>438 20X3>3 20X4>21 X5>37X7>67X8>268X9<1X1+X2+X3=12X4+X5+X6+X7+X8+X9+X10+X11+X12+X13+X14+X15+X16+X17+X18=12 END得到结果：LP OPTIMUM FOUND AT STEP 8OBJECTIVE FUNCTION V ALUE1) 7148.214V ARIABLE V ALUE REDUCED COSTX1 1.300000 0.000000X2 7.300000 0.000000X3 3.400000 0.000000X4 1.050000 0.000000X5 3.000000 0.000000X6 0.000000 1210.000000X7 0.857143 0.000000X8 3.714286 0.000000X9 0.000000 910.000000X10 0.000000 880.000000X11 0.000000 1005.000000X12 0.000000 810.000000X13 0.000000 450.000000X14 0.000000 750.000000X15 0.000000 220.000000X16 0.000000 630.000000X17 3.378572 0.000000X18 0.000000 680.000000ROW SLACK OR SURPLUS DUAL PRICES2) 0.000000 -0.2500003) 0.000000 -0.2500004) 65.000000 0.0000005) 0.000000 -41.5000006) 0.000000 -915.0000007) 1.200000 0.0000008) 0.000000 -139.5714269) 0.000000 -139.57142610) 1.000000 0.00000011) 0.000000 150.00000012) 0.000000 1130.000000NO. ITERATIONS= 8RANGES IN WHICH THE BASIS IS UNCHANGED:OBJ COEFFICIENT RANGESV ARIABLE CURRENT ALLOWABLE ALLOWABLECOEF INCREASE DECREASE X1 135.000000 14.999999 INFINITYX2 135.000000 14.999999 INFINITYX3 150.000000 INFINITY 14.999999X4 300.000000 830.000000 INFINITYX5 215.000000 915.000000 INFINITYX6 -80.000000 1210.000000 INFINITYX7 153.000000 976.999939 INFINITYX8 153.000000 976.999939 INFINITYX9 220.000000 910.000000 INFINITYX10 250.000000 880.000000 INFINITYX11 125.000000 1005.000000 INFINITYX12 320.000000 810.000000 INFINITYX13 680.000000 450.000000 INFINITYX14 380.000000 750.000000 INFINITYX15 910.000000 220.000000 INFINITYX16 500.000000 630.000000 INFINITYX17 1130.000000 INFINITY 220.000000X18 450.000000 680.000000 INFINITYRIGHTHAND SIDE RANGESROW CURRENT ALLOWABLE ALLOWABLERHS INCREASE DECREASE2 78.000000 195.000000 77.9999923 438.000000 195.000000 438.0000004 3.000000 65.000000 INFINITY5 21.000000 67.571426 20.9999986 3.000000 3.378572 3.0000007 1.200000 INFINITY 1.2000008 6.000000 23.650000 6.0000009 26.000000 23.650000 25.99999810 1.000000 INFINITY 1.00000011 12.000000 INFINITY 3.25000012 12.000000 INFINITY 3.378572 最后得到F MAX = 7148.214 – 2015.4 = 5132.814（元）。