资源简介
Boyd凸优化资料(内含中英文课本教材、习题、PPT、Matlab习题)建议大家配合官网上的视频食用,效果更佳哟~
代码片段和文件信息
% Generates figures 10.1 10.2 10.3 10.4 of Boyd & Vandenberghe
% Convex Optimization
%
% Infeasible start Newton method for solving
%
% minimize -sum_i log x_i
% subject to Ax = b
% Generate a feasible problem.
rand(‘state‘3);
randn(‘state‘4);
n = 100;
m = 50;
A = randn(mn);
x0 = rand(n1);
b = A*x0;
% Apply infeasible Newton method to KKT conditions
%
% -1./x + A‘*w = 0
% A*x = b
alpha = 0.011;
beta = 0.5;
tol = 1e-12;
maxiters = 20;
x = ones(n1);
w = zeros(m1);
resp = zeros(1 maxiters);
resd = zeros(1 maxiters);
steps = zeros(1 maxiters);
for k=1:maxiters
rp = A*x-b; rd = -1./x + A‘*w;
resp(k) = norm(rp); resd(k) = norm(rd);
normr = norm([rd; rp]);
if normr < 1e-12 break; end;
dw = - (A * diag(x.^2) * A‘) \ (A*((x.^2).*rd) - rp);
dx = - (x.^2) .* (rd + A‘*dw);
s = 1;
while min(x+s*dx) <= 0 s = s*beta; end;
while norm([ -1./(x+s*dx) + A‘*(w+s*dw); A*(x+s*dx)-b ]) > ...
(1-alpha*s) * normr
s = s*beta;
end;
steps(k) = s;
x = x + s*dx;
w = w + s*dw;
end;
% Figure 10.1
figure(1)
semilogy([0:k-1] resp(1:k) ‘o‘ [0:k-1] resp(1:k) ‘-‘ ...
[0:k-1] resd(1:k) ‘o‘ [0:k-1] resd(1:k) ‘--‘);
xlabel(‘iteration number‘)
ylabel(‘residuals‘)
% Figure 10.2
figure(2);
plot([1:k-1] steps(1:k-1) ‘o‘ [1:k-1] steps(1:k-1) ‘-‘);
xlabel(‘iteration number‘)
ylabel(‘step length‘)
axis([0 12 0.4 1.1]);
% Generate an infeasible problem.
rand(‘state‘3);
randn(‘state‘4);
A = randn(mn);
x0 = randn(n1);
b = A*x0;
% Apply infeasible Newton method to KKT conditions
%
% -1./x + A‘*w = 0
% A*x = b
alpha = 0.011;
beta = 0.5;
tol = 1e-12;
maxiters = 21;
x = ones(n1);
w = zeros(m1);
resp = zeros(1 maxiters);
resd = zeros(1 maxiters);
steps = zeros(1 maxiters);
for k=1:maxiters
rp = A*x-b; rd = -1./x + A‘*w;
resp(k) = norm(rp); resd(k) = norm(rd);
normr = norm([rd; rp]);
if normr < 1e-12 break; end;
dw = - (A * diag(x.^2) * A‘) \ (A*((x.^2).*rd) - rp);
dx = - (x.^2) .* (rd + A‘*dw);
s = 1;
while min(x+s*dx) <= 0 s = s*beta; end;
while norm([ -1./(x+s*dx) + A‘*(w+s*dw); A*(x+s*dx)-b ]) > ...
(1-alpha*s) * normr
s = s*beta;
end;
steps(k) = s;
x = x + s*dx;
w = w + s*dw;
end;
% Figure 10.3
figure(3)
semilogy([0:k-1] resp(1:k) ‘o‘ [0:k-1] resp(1:k) ‘-‘ ...
[0:k-1] resd(1:k) ‘o‘ [0:k-1] resd(1:k) ‘--‘);
xlabel(‘iteration number‘)
ylabel(‘residuals‘)
% Figure 10.4
figure(4);
plot([1:k-1] steps(1:k-1) ‘o‘ [1:k-1] steps(1:k-1) ‘-‘);
xlabel(‘iteration number‘)
ylabel(‘step length‘)
axis([0 20 0 0.3]);
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 0 2018-12-05 17:31 Boyd鍑镐紭鍖栬祫鏂?
文件 1737448 2018-11-29 16:25 Boyd鍑镐紭鍖栬祫鏂?slides-Convex Optimization.pdf
目录 0 2018-12-05 17:31 __MACOSX\
文件 0 2018-12-05 17:31 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?
文件 558 2018-11-29 16:25 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?._slides-Convex Optimization.pdf
文件 1121365 2018-11-29 16:24 Boyd鍑镐紭鍖栬祫鏂?bv_cvxbook_extra_exercises.pdf
文件 530 2018-11-29 16:24 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?._bv_cvxbook_extra_exercises.pdf
目录 0 2018-11-29 11:39 Boyd鍑镐紭鍖栬祫鏂?chapter 9\
文件 4525 2018-11-29 11:39 Boyd鍑镐紭鍖栬祫鏂?chapter 9\figure_9_25.m
目录 0 2018-12-05 17:31 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 9\
文件 628 2018-11-29 11:39 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 9\._figure_9_25.m
文件 2488 2018-11-29 11:39 Boyd鍑镐紭鍖栬祫鏂?chapter 9\figure_9_7.m
文件 627 2018-11-29 11:39 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 9\._figure_9_7.m
文件 9320 2018-11-29 11:39 Boyd鍑镐紭鍖栬祫鏂?chapter 9\figure_9_3.m
文件 627 2018-11-29 11:39 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 9\._figure_9_3.m
文件 5007 2018-11-29 11:39 Boyd鍑镐紭鍖栬祫鏂?chapter 9\figure_9_6.m
文件 627 2018-11-29 11:39 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 9\._figure_9_6.m
文件 553 2018-11-29 11:39 Boyd鍑镐紭鍖栬祫鏂?chapter 9\figure_9_2.m
文件 627 2018-11-29 11:39 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 9\._figure_9_2.m
文件 1473 2018-11-29 11:39 Boyd鍑镐紭鍖栬祫鏂?chapter 9\figure_9_23.m
文件 628 2018-11-29 11:39 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 9\._figure_9_23.m
文件 176 2018-11-29 11:39 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?._chapter 9
目录 0 2018-11-29 11:40 Boyd鍑镐紭鍖栬祫鏂?chapter 10\
文件 3802 2018-11-29 11:40 Boyd鍑镐紭鍖栬祫鏂?chapter 10\figure_10_6.m
目录 0 2018-12-05 17:31 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 10\
文件 629 2018-11-29 11:40 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 10\._figure_10_6.m
文件 1570 2018-11-29 11:40 Boyd鍑镐紭鍖栬祫鏂?chapter 10\figure_10_5.m
文件 629 2018-11-29 11:40 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 10\._figure_10_5.m
文件 2702 2018-11-29 11:39 Boyd鍑镐紭鍖栬祫鏂?chapter 10\figure_10_1.m
文件 629 2018-11-29 11:39 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?chapter 10\._figure_10_1.m
文件 176 2018-11-29 11:40 __MACOSX\Boyd鍑镐紭鍖栬祫鏂?._chapter 10
............此处省略23个文件信息
相关资源
- 数字信号处理-理论、算法与实现 第
- 冈萨雷斯—数字图像处理中文与英文
- 《MATLAB数值分析与应用》各章实验和
- cvx优化问题
- 凸优化公开课convex optimization,Boyd等著
- MATLAB与编程英文Chapman之习题解答
- 障碍法和原对偶内点法解决含有不等
- 数字信号处理理论算法与实现课后习
- 数学实验习题答案完整版
- matlab图像处理中英文翻译文献
- 纠错编码原理及MATLAB实现刘爱莲习题
-
simuli
nk练习题个人编写 - MATLAB程序设计与应用(刘卫国-第3版)
- 天线理论与技术第二版_钟顺时_第一章
- 椭球法求解凸优化
-
simuli
nk建模与仿真 与习题 - l1-Magic工具箱
- 上海电力大学-数值计算方法上机习题
- MATLAB程序设计与应用(刘卫国-第3版)
- OCR MATLAB OCR中英文段落的提取、分割
- 《卡尔曼滤波理论与实践matlab版的笔
- 最优阵列处理随书程序 (optimum arra
- Machine-Learning-exercises_finished andrew NG上
- cvx CVX工具包
- cvx 国外著名的凸优化程序包
- SSR_DOA 稀疏表示doa估计
- 数字信号处理matlab习题答案
- MATLAB仿真及电子信息应用第二版课后
- 《MATLAB教程及实训第2版》曹弋习题答
评论
共有 条评论