这个题用matlab穷举法怎么实现阿？高手帮帮忙阿

CHOOT · 发表于 2005-8-11 00:57:44

<

align=center>1997年B题截断切割
<

> 
<

>B题 截断切割
 某些工业部门（如贵重石材加工等）采用截断切割的加工方式。这里“截断切割”是指将物体沿某个切割平面分成两部分。从一个长方体中加工出一个已知尺寸、位置预定的长方体（这两个长主体的对应表面是平行的）通常要经过6次截断切割。
 设水平切割单位面积的费用是垂直切割单位面积费用的r倍，且当先后两次垂直切割的平面（不管它们之间是否穿插水平切割）不平行时，因调整刀具需额外费用e.
 试为这些部门设计一种安排各面加工次序（称“切割方式”）的方法，使加工费用最少。（由工艺要求，与水平工作台接触的长方体底面是事先指定的）详细要求如下：
1、需考虑的不同切割方式的总数
2、给出上述问题的数学模型和求解方法。
1、试对某部门用的如下准则作出评价：每次选择一个加工费用最少的待切割面进行切割。
2、对于e=0的情形有无简明的优化准则。
3、用以下实例验证你的方法：待加工长方体和成品长方体的长、宽、高分别为10、14.5、19和3、2、4，二者左侧面、正面、底面之间的距离分别为6、7、9（单位均为厘米）。垂直切割费用为每平方厘米1元，r和e的数据有以下4组：
a.r=1 e=0 ;b.r=1.5 e=0 ; c.r=8 ,e=0 ; d.r=1.5;2≤e≤15
对最后一组数据应给出所有最优解，并进行讨论。

CHOOT · 发表于 2005-8-11 01:33:53

大家帮帮忙啊！今晚要交了阿 [em06]

jlw686 · 发表于 2005-8-11 22:11:11

matlab递归效率太低，故我是用c实现的。 
/*设水平切割单位面积的费用是垂直切割单位面积费用的r倍，且当先后两次垂直切割的平面（不管它们之间 
是否穿插水平切割）不平行时，因调整刀具需额外费用e. 
用以下实例验证你的方法：待加工长方体和成品长方体的长、宽、高分别为10、14.5、19和3、2、4， 
二者左侧面、正面、底面之间的距离分别为6、7、9（单位均为厘米）。 
垂直切割费用为每平方厘米1元，r和e的数据有以下4组： 
 
a.r=1 e=0 ;b.r=1.5 e=0 ; c.r=8 ,e=0 ; d.r=1.5;2≤e≤15 
*/ 
#include <stdio.h> 
#define OK 1 
#define ERROR 0 
typedef int Status; 
//min_cost为最小花费，cost为总花费。 
//min_trace[6]为最小花费时的步骤。 
//trace[6]用来记录步骤，n表示已走完的步数。 
//在步骤中，0--5分别代表前后左右上下。 
//vs1为当前的前垂直面积，vs2为当前的侧垂直面积，ls为当前的水平面积。 
//vnum为已垂直切割的次数，lnum为已水平切割的次数。 
//p为两次垂直切割平面不平行的概率。 
//x,y,z为长方体的位置(二者左侧面、正面、底面之间的距离)。 
 
Status Find(int trace[],int n,int k){ 
for(int i=0;i<n;i++) 
 if(trace==k)break; 
if(i==n)return OK; 
return ERROR; 
}//Find 
 
Status Print(FILE* fp,int trace[],float cost){ 
printf("cost=%f\n",cost); 
fprintf(fp,"cost=%f\n",cost); 
printf("该路径为："); 
for(int i=0;i<6;i++){ 
 switch(trace){ 
 case 0: 
 printf("前--"); 
 fprintf(fp,"前--"); 
 break; 
 case 1: 
 printf("后--"); 
 fprintf(fp,"后--"); 
 break; 
 case 2: 
 printf("左--"); 
 fprintf(fp,"左--"); 
 break; 
 case 3: 
 printf("右--"); 
 fprintf(fp,"右--"); 
 break; 
 case 4: 
 printf("上--"); 
 fprintf(fp,"上--"); 
 break; 
 case 5: 
 printf("下--"); 
 fprintf(fp,"下--"); 
 break; 
 }//switch 
}//for 
printf("END\n"); 
fprintf(fp,"END\n"); 
return OK; 
}//Print 
 
Status Print_f2(FILE* fp,int trace[],float cost){ 
fprintf(fp,"cost=%f\n",cost); 
fprintf(fp,"该路径为："); 
for(int i=0;i<6;i++){ 
 switch(trace){ 
 case 0: 
 fprintf(fp,"前--"); 
 break; 
 case 1: 
 fprintf(fp,"后--"); 
 break; 
 case 2: 
 fprintf(fp,"左--"); 
 break; 
 case 3: 
 fprintf(fp,"右--"); 
 break; 
 case 4: 
 fprintf(fp,"上--"); 
 break; 
 case 5: 
 fprintf(fp,"下--"); 
 break; 
 }//switch 
}//for 
fprintf(fp,"END\n"); 
return OK; 
}//Print 
 
 
Status cutting(FILE*fp1,FILE*fp2,float& min_cost,int min_trace[],float cost,int trace[],int n,float vs1,float vs2,float ls, 

float unit_cost,int p,int e,int r,float x,float y,float z,float l,float
w,float h, 
 float L,float W,float H){ 
if(n>=6){ 
 Print(fp1,trace,cost); 
 if(min_cost>cost){ 
 fclose(fp2); 
 fp2=fopen("最优方案.txt","w"); 
 min_cost=cost; 
 for(int i=0;i<6;i++){ 
 min_trace=trace; 
 }//for 
 }//if 
 if(min_cost==cost) 
 Print_f2(fp2,trace,cost); 
 return OK; 
}//if 
float cost1; 
if(Find(trace,n,0)){ 
 trace[n]=0; 
 cost1=cost+unit_cost*vs1+p*e; 

cutting(fp1,fp2,min_cost,min_trace,cost1,trace,n+1,vs1,(W-y)*H,(W-y)*L,unit_cost,p,e,r,x,0,z,l,w,h,L,W-y,H); 
}//if 
if(Find(trace,n,1)){ 
 trace[n]=1; 
 cost1=cost+unit_cost*vs1+p*e; 

cutting(fp1,fp2,min_cost,min_trace,cost1,trace,n+1,vs1,(w+y)*H,(w+y)*L,unit_cost,p,e,r,x,y,z,l,w,h,L,w+y,H); 
}//if 
if(Find(trace,n,2)){ 
 trace[n]=2; 
 cost1=cost+unit_cost*vs2+p*e; 

cutting(fp1,fp2,min_cost,min_trace,cost1,trace,n+1,H*(L-x),vs2,W*(L-x),unit_cost,p,e,r,0,y,z,l,w,h,L-x,W,H); 
}//if 
if(Find(trace,n,3)){ 
 trace[n]=3; 
 cost1=cost+unit_cost*vs2+p*e; 

cutting(fp1,fp2,min_cost,min_trace,cost1,trace,n+1,H*(l+x),vs2,W*(l+x),unit_cost,p,e,r,x,y,z,l,w,h,l+x,W,H); 
}//if 
if(Find(trace,n,4)){ 
 trace[n]=4; 
 cost1=cost+r*unit_cost*ls+p*e; 

cutting(fp1,fp2,min_cost,min_trace,cost1,trace,n+1,L*(h+z),W*(h+z),ls,unit_cost,p,e,r,x,y,z,l,w,h,L,W,h+z); 
}//if 
if(Find(trace,n,5)){ 
 trace[n]=5; 
 cost1=cost+r*unit_cost*ls+p*e; 

cutting(fp1,fp2,min_cost,min_trace,cost1,trace,n+1,L*(H-z),W*(H-z),ls,unit_cost,p,e,r,x,y,0,l,w,h,L,W,H-z); 
}//if 
return OK; 
}//cutting 
 
void main(){ 
int min_trace[6],trace[6],n(0),p(0),e,r; 
float min_cost(10000),cost(0),unit_cost(1),x(6),y(7),z(9); 
float l(3),w(2),h(4),L(10),W(14.5),H(19),vs1,vs2,ls; 
FILE* fp1,*fp2; 
fp1=fopen("切割方案.txt","w"); 
fp2=fopen("最优方案.txt","w"); 
vs1=L*H; 
vs2=W*H; 
ls=L*W; 
printf("e="); 
scanf("%d",&e); 
printf("r="); 
scanf("%d",&r); 
cutting(fp1,fp2,min_cost,min_trace,cost,trace,n,vs1,vs2,ls,unit_cost,p,e,r,x,y,z,l,w,h,L,W,H); 
printf("最少花费为：%f\n",min_cost); 
fprintf(fp1,"最少花费为：%f\n",min_cost); 
printf("是请参看“切割方案.txt,最优方案.txt”\n"); 
fclose(fp1); 
fclose(fp2); 
}//main

nihaofeng · 发表于 2005-8-13 17:49:05

有些看不懂啊

mybliao · 发表于 2005-8-15 07:55:02

function qiege price=inf; x1=6; x2=10-3-6; x3=7; x4=14.5-2-7; x5=9; x6=19-4-9; n=zeros(1,6); num=0; best_n=[]; for n1=1:6, for n2=1:6, for n3=1:6, for n4=1:6, for n5=1:6, for n6=1:6, if n1~=n2&n1~=n3&n1~=n4&n1~=n5&n1~=n6&n2~=n3&n2~=n4&n2~=n5&n2~=n6... &n3~=n4&n3~=n5&n3~=n6&n4~=n5&n4~=n6&n5~=n6, n=[n1 n2 n3 n4 n5 n6]; num=num+1; a=10;b=14.5;c=19; cur_price=0; for i=1:6, switch n(i) case 1 a=a-x1; cur_price=b*c+cur_price; case 2 a=a-x2; cur_price=b*c+cur_price; case 3 b=b-x3; cur_price=a*c+cur_price; case 4 b=b-x4; cur_price=a*c+cur_price; case 5 c=c-x5; cur_price=a*b+cur_price; case 6 c=c-x6; cur_price=a*b+cur_price; end end if cur_price<price, price=cur_price; best_n=n'; elseif cur_price==price best_n=[best_n,n']; end end end end end end end end best_n price num

		自动登录	找回密码
密码			注-册-帐-号

这个题用matlab穷举法怎么实现阿？高手帮帮忙阿

浏览过的版块