1、姓名徐競(jìng)17評(píng)分實(shí)驗(yàn)報(bào)告課程名稱： 動(dòng)態(tài)規(guī)劃 實(shí)驗(yàn)名稱： 二維分配問(wèn)題 專(zhuān) 業(yè)： 信息與計(jì)算科學(xué) 指導(dǎo)教師： 滕宇 完成日期： 2010年 04月 07日 綜合實(shí)驗(yàn)報(bào)告一、實(shí)驗(yàn)名稱二維分配問(wèn)題實(shí)驗(yàn)報(bào)告 二、 實(shí)驗(yàn)?zāi)康模?） 理解多維動(dòng)態(tài)規(guī)劃的含義及應(yīng)用；（2） 掌握一維、二維以及多維火力分配問(wèn)題的計(jì)算方法與含義；（3） 上機(jī)編寫(xiě)、調(diào)試二維火力分配問(wèn)題的算法程序；三、 實(shí)驗(yàn)要求（1） 理解一維及二維火力分配問(wèn)題的計(jì)算方法；（2） 寫(xiě)出二維火力分配問(wèn)題的模型；（3） 理解二維火力分配問(wèn)題中拉格朗日乘數(shù)法，逐次逼近法，粗格子點(diǎn)法（疏密法）三種計(jì)算方法；（4） 選擇C+計(jì)算機(jī)語(yǔ)言編寫(xiě)一種算法的程序來(lái)實(shí)

2、現(xiàn)該實(shí)驗(yàn).四、 實(shí)驗(yàn)原理（1） 火力分配問(wèn)題，更一般提法，它是屬于資源分配問(wèn)題一類(lèi)的。所謂資源分配問(wèn)題，就是供應(yīng)量有限的一種或若干種資源(例如，原材料，資金，機(jī)器設(shè)備，勞力，食品，武器，兵力等)，分配給若干個(gè)使用者(或打擊目標(biāo))，而使最終結(jié)果達(dá)到最優(yōu)。（2） 二維火力分配問(wèn)題的非線性規(guī)劃的數(shù)學(xué)模型，有動(dòng)態(tài)規(guī)劃基本方程： 五、 實(shí)驗(yàn)題目設(shè)有3架敵機(jī), 價(jià)值分別為=3, =2, =5， 用二種導(dǎo)彈射擊，第一種導(dǎo)彈有6枚(m1=6)，第二種導(dǎo)彈有10枚(m2=10), 其單發(fā)擊毀概率如下表所示，求最優(yōu)分配方案(xk, yk) ，k=1，2，3,。 目標(biāo)兵器12310.40.10.520.20.40.

3、2六、 實(shí)驗(yàn)步驟（1） 由題意建立模型，設(shè)對(duì)空中敵機(jī)射擊總效果為V，則有： ； （2） 根據(jù)拉格朗日乘數(shù)法降維方法，在計(jì)算機(jī)上編寫(xiě)C程序求解。七、 程序及其運(yùn)行結(jié)果#include<iostream>#include<math.h>using namespace std;float maxc;int y,xf;float v17,v27,v37,f27,f17;int y117,y127,y137,xfi7,xfj7; /定義數(shù)據(jù)float factorial(float a,int x) if (x=0) return 1.0; else if (x=1) retur

4、n a; else return a*factorial(a,x-1);void comparemax(int a,float b,float c,int x,float para)int i;float ku11; for(i=0;i<=10;i+) kui=a*(1-factorial(b,x)*factorial(c,i)-para*i; /對(duì)第二種導(dǎo)彈（10顆）對(duì)目標(biāo)的射擊情況進(jìn)行列舉 maxc=ku0; for(i=0;i<=10;i+) if(kui>=maxc) maxc=kui; y=i; /找出各種情況的最大解 void funck2(float para)

5、 for(int i=0;i<=6;i+) comparemax(2,0.9,0.6,i,para); v2i=maxc; y12i=y; comparemax(5,0.5,0.8,i,para); v3i=maxc; y13i=y; float f7;int j=0; for(int x=0;x<=6;x+) for(int x2=0;x2<=x;x2+) fx2=v2x2+v3x-x2; maxc=f0; for(int i=0;i<=x;i+) if(fi>=maxc) maxc=fi; xf=i; f2j=maxc;xfij=xf;+j; void fun

6、ck1(float para) funck2(para); for(int i=0;i<=6;i+) comparemax(3,0.6,0.8,i,para); v1i=maxc; y11i=y; float f7;int j=0; for(int x=0;x<=6;x+) for(int x1=0;x1<=x;x1+) fx1=v1x1+f2x-x1; maxc=f0; for(int i=0;i<=x;i+) if(fi>=maxc) maxc=fi; xf=i; f1j=maxc;xfjj=xf;+j; void liebiao(float para) in

7、t i; for(i=0;i<=6;i+) comparemax(3,0.6,0.8,i,para); v1i=maxc; y11i=y; comparemax(2,0.9,0.6,i,para); v2i=maxc; y12i=y; comparemax(5,0.5,0.8,i,para); v3i=maxc; y13i=y; int excellre(float para) int x1,x2,x3,x,xd,j; int X7=0,1,2,3,4,5,6; liebiao(para); cout<<"X"<<'t'<

8、<"V1(x)"<<'t'<<"Y1"<<'t'<<"V2(x)"<<'t'<<"Y2"<<'t'<<"V3(x)"<<'t'<<"Y3"<<endl; for(j=0;j<=6;j+) cout<<Xj<<'t'&

9、lt;<v1j<<'t'<<y11j<<'t'<<v2j<<'t'<<y12j<<'t'<<v3j<<'t'<<y13j<<endl; cout<<"<<<<<<<<<<<<.>>>>>>>>>>"<<en

10、dl; funck1(para); cout<<"X"<<'t'<<"f3(x)"<<'t'<<"x3"<<'t'<<"f2(x)"<<'t'<<"x2"<<'t'<<"f1(x)"<<'t'<<"x1"&l

11、t;<endl; for(j=0;j<=6;j+) cout<<Xj<<'t'<<v3j<<'t'<<Xj<<'t'<<f2j<<'t'<<xfij<<'t'<<f1j<<'t'<<xfjj<<endl; x=6; x1=xfjx; xd=x-x1; x=xd; x2=xfix; xd=x-x2; x=xd; x3=Xxd

12、; cout<<"得出最后結(jié)果！"<<endl; cout<<x1<<'t'<<y11x1<<endl; cout<<x2<<'t'<<y12x2<<endl; cout<<x3<<'t'<<y13x3<<endl; if(y11x1+y12x2+y13x3=10) cout<<"para = "<<para<&

13、lt;"滿足條件"<<endl; cout<<"得出最優(yōu)解！"<<endl; cout<<x1<<'t'<<y11x1<<endl; cout<<x2<<'t'<<y12x2<<endl; cout<<x3<<'t'<<y13x3<<endl; return 1; else cout<<"para 不符合要求

14、"<<endl; return 0; int main() float para;int k; cout<<"請(qǐng)輸入?yún)?shù)para（0<=para<=1）= " cin>>para;cout<<endl; do k=excellre(para); if(k=0) cout<<"請(qǐng)輸入?yún)?shù)para（0<=para<=1）= " cin>>para;cout<<endl; else break; if(para>1|para<0) brea