《机械优化设计》课程实践报告.doc

提供各专业全套设计
(课程实践报告封面模版)
合肥工业大学
《机械优化设计》课程实践
研究报告
班级:
学号:
姓名:
授课老师:
日期:
目录
一主要内容
1、一维搜索程序作业
= ..........................................1
..........................2
2、单位矩阵程序作业............................................4
3、其他工程优化问题..................................9
4连杆机构问题.....................................12
二实践心得体会...............................15
一: 主要内容
1. 一维搜索程序作业:
λ= (y2 > y1)
证明:、α2 的位置相对于区间[a,b] 两端点具有对称性,即
已知 a1=a2 , 要求α11=α22
由于α1=b-λ(b-a)
α2=a+λ(b-a)
若使α11=α22
则有:b1-λ(b1-a1)=a2+λ(b2-a2)= a1+λ2(b1-a1)
因此: b1- a1=(λ2+λ)( b1- a1)
( b1- a1)( λ2+λ-1)=0
因为: b1= a1
所以: λ2+λ-1=0
则有: 取方程正数解得
若保留下来的区间为[α1,b] ,根据插入点的对称性,也能推得同样的λ的值。
:


例:(1)a=0 ,b=2π,f(x)=cox(x)
(2)a=0 ,b=10, f(x)=(x-2)2+3
(1)
#include<>
#include<>
void main(void)
{
int i;
float a1,a2,aa,y1,y2,ymin,e;
float a=0,b=2*,n=;
a1=b-n*(b-a);
a2=a+n*(b-a);
print(“输入精度:”);
scanf(“%f”,&e);
for(i=0;i=10000;i=i++)
{
y1=cos(a1);
y2=cos(a2);
if(y1<y2)
{
a=a1;
a1=a2;
a2=a+n*(b-a);
}
If(y1<y2)
b=a2;
a2=a1;
a1=b-n*(b-a);
}
if(fabs(b-a)/b<e)&&(fabs(y1-y2)/y2<e))
{
aa=(a+b)/2;
ymin=cos(aa);
printf(“x=%\tf(x)=%\n”),aa,ymin);
break;
}
}
}
运行结果:
(2)
#include<>
#include<>
void main(void)
{
int i;
float a1,a2,aa,y1,y2,ymin,e;
float a=0,b=10,n=;
a1=b-n*(b-a);
a2=a+n*(b-a);
print(“输入精度:”);
scanf(“%f”,&e);
for(i=0;i=10000;i=i++)
{
y1=(a1-2)*(a1-2)+3;
y2=(a2-2)*(a2-2)+3;
if(y1>=y2)
{
a=a1;
a1=a2;
a2=a+n*(b-a);
}
If(y1<y2)
b=a2;
a2=a1;
a1=b-n*(b-a);
}
if(fabs(b-a)/b<e)&&(fabs(y1-y2)/y2<e))
{
aa=(a+b)/2;
ymin=(aa-2)*(aa-2)+3;
printf(“x=%\tf(x)=%\n”),aa,ymin);
break;
}
}
}
运行结果:

编写生成单位矩阵的程序
程序文本
#include <>
void main(void)
{
int a[100][100];
int N,i,j;
printf("请输入所要输出矩阵的阶数(最多100阶):");
scanf("%d",&N);
printf("输出的矩阵阶数