VS2010调用Matlab的m文件所生产的dll实现二次规划算法.docxVIP

VS2010调用Matlab的m文件所生产的dll实现二次规划算法

引言

二次规划算法可在SQP算法中成为通用的二次规划子问题程序。二次规划算法设计到众多的矩阵运算，用C++代码编程需要较大工作量，而矩阵运算在Matlab下那么相当方便。但SQP多用于非线性系统的优化，如航空发动机的优化，航空发动机模型通常在C代码下更为方便。所以本文采用的策略是在Matlab下编写实现二次规划算法的m文件，使用Matlab编译器将m文件编译成dll动态链接库，然后在VS下调用dll来实现二次规划算法。

1系统环境说明

操作系统：Window8X64.

Matlab版本：Matlab2012b(Matlab8.0).

VS版本：VS2010.

添加系统环境变量:

Path变量中增加：C:\ProgramFiles\MATLAB\R2012b\bin\win64;

重启电脑。

2二次规划算法的Matlab实现

2.1编程准备工作

首先为本文建立一个单独的文件夹\QPTEST，再在该文件夹下建立子文件夹\QPMAT。然后启动Matlab并将工作目录指定到\QPMAT下。

2.2m文件编程

新建m文件并保存文件名为qpsubp.m，编辑内容如下，代码分析参考文献[1]程序12.1：

function[d,mu,lam]=qpsubp(dfk,Bk,Ae,hk,Ai,gk)

n=length(dfk);l=length(hk);m=length(gk);

gamma=0.05;epsilon=1.0e-6;rho=0.5;sigma=0.2;

ep0=0.05;mu0=0.05*zeros(l,1);lam0=0.05*zeros(m,1);

d0=ones(n,1);u0=[ep0;zeros(n+l+m,1)];

z0=[ep0;d0;mu0;lam0];

k=0;

z=z0;ep=ep0;d=d0;mu=mu0;lam=lam0;

while(k=150)

dh=dah(ep,d,mu,lam,dfk,Bk,Ae,hk,Ai,gk);

if(norm(dh)epsilon)

break;

end

A=JacobiH(ep,d,mu,lam,dfk,Bk,Ae,hk,Ai,gk);

b=beta(ep,d,mu,lam,dfk,Bk,Ae,hk,Ai,gk,gamma)*u0-dh;

dz=A\b;

if(l0m0)

de=dz(1);dd=dz(2:n+1);du=dz(n+2:n+l+1);dl=dz(n+l+2:n+l+m+1);

end

if(l==0)

de=dz(1);dd=dz(2:n+1);dl=dz(n+2:n+m+1);

end

if(m==0)

de=dz(1);dd=dz(2:n+1);du=dz(n+2:n+l+1);

end

i=0;

while(i=20)

if(l0m0)

dh1=dah(ep+rho^i*de,d+rho^i*dd,mu+rho^i*du,lam+rho^i*dl,...

dfk,Bk,Ae,hk,Ai,gk);

end

if(l==0)

dh1=dah(ep+rho^i*de,d+rho^i*dd,mu,lam+rho^i*dl,...

dfk,Bk,Ae,hk,Ai,gk);

end

if(m==0)

dh1=dah(ep+rho^i*de,d+rho^i*dd,mu+rho^i*du,lam,...

dfk,Bk,Ae,hk,Ai,gk);

end

if(norm(dh1)=(1-sigma*(1-gamma*ep0)*rho^i)*norm(dh))

mk=i;break;

end

i=i+1;

if(i==20)

mk=10;

end

end

alpha=rho^mk;

if(l0m0)

ep=ep+