Matlab 线性规划

最新推荐文章于 2024-07-30 23:24:40 发布

W | Z | H

最新推荐文章于 2024-07-30 23:24:40 发布

阅读量33

点赞数

文章标签： matlab 算法

本文链接：https://blog.csdn.net/wzh941009/article/details/135937199

版权

本文详细介绍了如何在Matlab中通过`linprog`函数解决线性规划问题，包括数据预处理、构造AUG矩阵、求解线性优化问题，并对结果进行平滑处理。

摘要由CSDN通过智能技术生成

function phi_res = linear_program(phi,m)
% 线性规划
% phi:数据
% m:阶数

% A_AUG
n = length(phi);
I_n = eye(n);
Z = zeros(n-m+1,n);
M = zeros(n-m+1,n);
for i = 1:size(M,1)
M(i,i:i+m-1) = CSG(m);
end
A_AUG = [I_n,-I_n;I_n,I_n;Z,M];
% phi_AUG
phi_AUG = [-phi,phi,zeros(1,n-m+1)]';

C = [ones(1,n),zeros(1,n)];
phi_c=linprog(C,-A_AUG,-phi_AUG);

phi_res = phi_c(n+1:end);

k = (m-1)/2;
sm = smth(k)';
for i = 1+k:size(phi_res,1)-k
phi_res(i) = sum(phi_res(i-k:i+k) .*sm);
end
end
%%
function res = CSG(m)
res = [];
for i = 1:m
res = [res,6*(2*i-m-1)/(m*(m+1)*(m-1))]; %#ok<*AGROW>
end
end

function sm = smth(k)
tmp = 0;
for i = 1:k
tmp = tmp + i*i;
end
tmp = tmp * 2;
dm = -k/tmp : 1/tmp : k/tmp;

for i = 1:k+1
sm(i) = -sum(dm(1:i-1)) - dm(i)/2;
end
for i = k+2:2*k+1
sm(i) = sum(dm(i+1:2*k+1)) + dm(i)/2;
end
end