Lagrange 乘子法求最优解
2017-04-18 16:24
120 查看
clc clear syms x y z r1 r2 w f=x^2+y^2+z^2+w^2; g1=2*x-y+z-w-1; g2=x+y-z+w-1; h=f-r1*g1 -r2*g2; hx=diff(h,x); hy=diff(h,y); hz=diff(h,z); hw=diff(h,w); hr1=diff(h,r1); hr2=diff(h,r2); r=solve([hx==0,hy==0,hz==0,hw==0,hr1==0,hr2==0],[x,y,z,w,r1,r2]); arr_x=double(r.x) arr_y=double(r.y) arr_z=double(r.z) arr_w=double(r.w) arr_fv=[]; for i=1 : length(arr_x) xv=arr_x(i); yv=arr_y(i); zv=arr_z(i); wv=arr_w(i); fv=subs(f,[x,y,z,w],[xv,yv,zv,wv]); arr_fv(i)= double(fv); g1v=subs(g1,[x,y,z,w],[xv,yv,zv,wv]) g2v=subs(g2,[x,y,z,w],[xv,yv,zv,wv]) end arr_fv disp('after sort:') sort(arr_fv)
View Code
相关文章推荐
- 波松分酒问题 C++求最优解.
- 广度搜索——分酒问题(最优解)
- Lagrange插值问题
- [0612]数字反转C++最优解
- 解空间树搜索 及 最优解
- 树上最优解+贪心+dfs+LA3902
- [物理学与PDEs]第4章第3节 一维反应流体力学方程组 3.2 一维反应流体力学方程组的 Lagrange 形式
- Hoj1412 求两个集合的和{A} + {B}——java实现最优解的辛酸路
- hdoj1003注意在寻找最优解对应的时候寻找的顺序
- Lagrange 四平方定理
- Project Euler 103:Special subset sums: optimum 特殊的子集和:最优解
- 使用贪心算法求解多机调度问题的近似最优解。
- 拉格朗日(Lagrange)插值曲线
- leetcode最优解整理2(Linked list/ Stack/Queue/Array/HashTable/ Tree/bfs/union find)
- HDU 2639 Bone Collector II(01背包变形【第K大最优解】)
- 求和小于k的最长子数组的最优解
- 黄金分割法推导及单变量函数近似最优解
- 防止SQL注入的方法和最优解
- 动态规划——最长公共子序列LIS(dp+最优解)
- 波松分酒问题 C++求最优解.