本帖最后由 飞人周杰伦 于 2016-3-21 21:39 编辑

求目标函数

2.png (1.85 KB, 下载次数: 0)

2016-3-21 21:00 上传

,其中

11.png (29.86 KB, 下载次数: 0)

2016-3-21 21:37 上传

，目标函数中n=23为决策单元数量，t=2014为年份，i=4为投入种类，j=2为期望产出种类。x,y,b分别为投入，期望产出，非期望产出向量。x=(total asset; the number of employee; operating expenses total; cost of goods sale)y=(net sale; gross income)  b=(CO2)

22.png (21.15 KB, 下载次数: 0)

2016-3-21 21:39 上传

x和y都是已知的23×4维向量和23x2维向量，b是23×1维已知的向量。

5.png (444 Bytes, 下载次数: 0)

2016-3-21 21:04 上传

和

6.png (431 Bytes, 下载次数: 0)

2016-3-21 21:04 上传

之类的系数是未知量，x和y里面的元素见world里面的数据。

本人由于初学最优化编程，了解一些皮毛，写了一个只能意会但很有问题的代码，如下：

本人想问的问题如下：

1.二重求和能否用我代码中的symsum进行嵌套来表示？

2.约束条件中的偏导数如何用代码表示？

3.主函数写的很乱，自己都不知道对不对。

请高人能指点我，程序到处都是问题，题目意思易懂，但是表达起来就感觉有问题。

文件一：%函数表达式

function f=my_fun(x)

f=symsum((symsum(x(2)(j)*(Y(n,j))^t,1,2)

+x(6)*(b(n))^t

+1/2*symsum(symsum(x(1)(i,i')*(X(n,i))^t*(X(n,i'))^t, i, 1, 4),i',1,4)

+1/2*symsum(symsum(x(2)(j,j')*(Y(n,j))^t*(Y(n,j'))^t, j, 1, 2),j',1,2)

+1/2*x(7)*b(n)^(2t)

+symsum(x(4)(j)*((Y(n,j))^t)*b(n)^t,1,2)

+symsum(x(5)(i)*((X(n,i))^t)*b(n)^t,1,4)

+symsum(symsum(x(3)(i,j)*(X(n,i))^t*(Y(n,j))^t, i, 1, 4),j,1,2),1,23);

文件二：%约束条件

function [c,ceq]=xianzhi(x,i,i',j,j')

c(1)=-f;

c(2)=diff(f(X(n,i),Y(n,j),b),Y(n,j));

c(3)=-diff(f(X(n,i),Y(n,j),b),b);

c(4)=-diff(f(X(n,i),Y(n,j),b),X(n,i));

ceq(1)=symsum(x(2)(j),1,2)-x(6)+1;

ceq(2)=symsum(symsum(x(2)(j,j'),1,2),1,2)-symsum(x(4)(j),1,2);

ceq(3)=symsum(x(4)(j),1,2)-x(7);

ceq(4)=symsum(symsum(x(3)(i,j),1,4),1,2)-symsum(x(5)(j),1,4);

ceq(5)=symsum(symsum(x(2)(j,j'),1,2),1,2)+x(7)-symsum(x(4)(j),1,2);

ceq(6)=x(1)(i,i')-x(1)(i',i);

ceq(7)=x(2)(j,j')-x(2)(j',j);

文件三：%主程序

clc;clear;    % syms k symsum(k^2, 0, 10) 级数求和,symsum(x^k/sym('k!'), k, 0, Inf)

syms X Y b;  %建立符号变量

syms alfa  beta  fi  ;%  x(1,2,3)

syms  exci   yita    gama1   gama2;%x(4,5,6,7)

syms  n  t;

% syms      i  i'  j  j';

i  i'  j  j';

N=23; %决策单元数量

t=2014;%年份

alfa=[];

beta=[];

fi=[];

exci=[];

yita=[];

D=[];

x0=[1;1;1;1;1;1;1];

lb=[];

ub=[];

A=[];

B=[];

Aeq=[];

Beq=[];

[x,fval]=fmincon(@my_fun,x0,A,B,Aeq,Beq,lb,ub,@xianzhi)