现代设计方法

机械设计方法实验报告

姓名：

学号：

成绩：

指导教师：

进退试算法实验报告

一、实验目的

1.加深对进退试算法的基本理论和算法步骤的理解。

2.培养独立编制、调试计算机程序的能力。

3.掌握常用优化程序的使用方法。

4.培养灵活运用优化设计方法解决工程实际问题的能力。

二、实验要求

1.明确进退试算法基本原理及程序框图。

2.编制进退试算法程序。

三．实验内容

计算实例：用进退试算法求函数())2

t

f的搜索区间。

(+

=t t

①.进退试算法基本原理简述

进退试算法的基本思想是：按照一定的规律给出若干试算点，一次比较各试算点的函数值的大小，直到找出相邻的三点的函数值按“高——低——高”变化的单峰区间为止。

②、程序的流程图

③.编制进退试算法程序

#include

#include

#define f(t) (t*(t+2))

void sb(double *a,double *b)

{

double t0,t1,t,h,alpha,f0,f1;

int k=0;

printf("请输入初始点t0=");

scanf("%lf",&t0);

printf("\n请输入初始步长h=");

scanf("%lf",&h);

printf("\n请输入加步系数alpha(需大于1)=");

scanf("%lf",&alpha);

f0=f(t0);

t1=t0+h;

f1=f(t1);

while(1)

{

printf("\nf1=%lf,f2=%lf,t0=%lf,t=%lf,h=%lf,k=%d",f0,f1,t0,t1,h,k);

if(f1

{

h=alpha*h;

t=t0;

t0=t1;

f0=f1;

k++;

}

else

{

if(k==0)

{h=-h;t=t1;}

else

{

*a=t

*b=t>t1?t:t1;

break;

}

}

t1=t0+h;

f1=f(t1);

}

}

main()

{

double a=0,b=0;

double *c,*d;

c=&a,d=&b;

sb(c,d);

printf("\na=%lf,b=%lf",a,b); }

④．程序运行结果

鲍威尔共轭方向法实验报告

一、实验目的

5. 加深对鲍威尔法的基本理论和算法步骤的理解。

6. 培养独立编制、调试计算机程序的能力。

7. 掌握常用优化程序的使用方法。

8. 培养灵活运用优化设计方法解决工程实际问题的能力。

二、实验要求

3. 明确鲍威尔法基本原理及程序框图。

4. 编制鲍威尔法程序。

三．实验内容

计算实例：用鲍威尔法求函数()212

221212141060,x x x x x x x x f -++--=的极

小值。

①.鲍威尔法基本原理简述

任选一初始点X 0，再选两个线性无关的向量。从X 0出发，顺

次沿e 1、e 2作一维搜索得01X 、02X ，两点连线得一新方向d 1，用d

1代替e 1形成两个线性无关向量e 2、d 1，作为下一轮搜索方向。再从02

X 出发，沿d 1作一维搜索得点01X ，作为下一轮迭代的初始点。从X 1出发，顺次沿e 2、d 1作一维搜索，得到点11X 、12X ，两点的连线得一新方向d 2。10X 、12X 两点是从不同点X 0、11X 出发，分别沿d 1方向进

行一维搜索而得到的极小点。再从12X 出发，沿d 2作一维搜索得点X 2，即是二维问题的极小点X *。

③.编制鲍威尔法程序

#include "stdio.h"

#include "stdlib.h"

#include "math.h"

double objf(double x[])

{double ff;

ff=60-10*x[0]-4*x[1]+x[0]*x[0]+x[1]*x[1]-x[0]*x[1];

return(ff);

}

void jtf(double x0[],double h0,double s[],int n,double a[],double b[]) {int i;

double *x[3],h,f1,f2,f3;

for(i=0;i<3;i++)

x[i]=(double *)malloc(n*sizeof(double));

h=h0;

for(i=0;i

*(x[0]+i)=x0[i];

f1=objf(x[0]);

for(i=0;i

*(x[1]+i)=*(x[0]+i)+h*s[i];

f2=objf(x[1]);

if(f2>=f1)

{ h=-h0;

for(i=0;i

*(x[2]+i)=*(x[0]+i);

f3=f1;

for(i=0;i

{*(x[0]+i)=*(x[1]+i);

*(x[1]+i)=*(x[2]+i);

}

f1=f2;

f2=f3;

}

for(;;)

{h=2*h;

for(i=0;i

*(x[2]+i)=*(x[1]+i)+h*s[i];

f3=objf(x[2]);

if(f2

else

{ for(i=0;i

{*(x[0]+i)=*(x[1]+i);

*(x[1]+i)=*(x[2]+i);