大连理工c语言第五次上机作业参考答案
大连理工c语言第五次上机
作业参考答案
-标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII
第五次上机作业：（请用指针实现）
1． 数据移位
有n （n<=10）个整数，使前面各数顺序向后移m （m<=n ）个位置，最后m 个数变成最前面的m 个数，如图所示。编写函数shift(a,n,m)实现上述功能，在主函数中输入n 个整数和输出调整后的n 个数。
SAMPLE OUTPUT Input n(n<=10)：10
Input 10 integers ：1 2 3 4 5 6 7 8 9 10 Output the original ：1 2 3 4 5 6 7 8 9 10
Input m(m<=10)：4
Output the shifted ：7 8 9 10 1 2 3 4 5 6
#include
#define N 10
void shift(int a[ ],int n,int m);
void main( )
{ int i,n,m;
int a[N];
printf("Input n(n<=10):");
scanf("%d",&n);
printf("Input %d integers:\n",n);
for(i=0;i scanf("%d",a+i); printf("Output the original:\n"); for(i=0;i printf("%d ",*(a+i)); printf("\n"); printf("Input m(m<=%d):",n); scanf("%d",&m); shift(a,n,m); printf("Output the shifted:\n"); for(i=0;i printf("%d ",*(a+i)); printf("\n"); } void shift(int a[ ],int n,int m) { int i,j,t; for(i=0;i { t=a[n-1]; for(j=n-2;j>=0;j--) a[j+1]=a[j]; a[0]=t; } } n 2．字符串排序 在主函数中输入10个不等长的字符串放入二维字符数组中，编写函数sort()利用指针数组对其排序，在主函数中输出排好序的字符串。 #include #include #define N 10 #define M 80 void sort(char *s[ ]); void main() { char s[N][M],*strp[N]; int i; printf("Input 10 strings:\n"); for(i=0;i gets(s[i]); for(i=0;i strp[i]=s[i]; sort(strp); printf("Output sorted strings:\n"); for(i=0;i puts(strp[i]); } void sort(char *s[ ]) { int i,j,k; char *p; for(i=0;i { k=i; for(j=i+1;j if(strcmp(s[k],s[j])>0) k=j; f(i!=j) { p=s[k]; s[k]=s[i]; s[i]=p; } } } 3．学生成绩统计 有一个班级30个学生5门课，编写函数实现下列功能： （1）函数average()求第一门课的平均分； （2）函数findfail()找出有2门以上（含2门）课程不及格的学生，输出其学号和全部课程成绩及平均成绩； （3）函数findgood()找出平均成绩在90分以上或全部课程成绩在85分以上的学生。 在主函数中输入数据并输出结果。 #include #define N 30 #define M 5 int average(int a[ ][M+3]); void findfail(int a[ ][M+3]); void findgood(int a[ ][M+3]); void main( ) { int s[N][M+3],i,j; int ave; printf("Input 5 scores for each student:\n"); for(i=0;i { s[i][0]=i+1; printf("%d:",i+1); for(j=1;j<=M;j++) scanf("%d",&s[i][j]); printf("\n"); } ave=average(s); printf("The average of the first course is %d\n",ave); findfail(s); printf("The students who failed in 2 or more courses:\n"); for(i=0;i { if(s[i][M+2]=='N') { for(j=0;j printf("%-5d",s[i][j]); printf("\n"); } } findgood(s); printf("The students whose average>90 or each course>85:\n"); for(i=0;i { if(s[i][M+2]=='Y') { for(j=0;j printf("%-5d",s[i][j]); printf("\n"); } } } int average(int a[ ][M+3]) { int i,sum=0; for(i=0;i sum+=a[i][1]; return(sum/N); } void findfail(int s[][M+3]) { int i,j,count,sum; for(i=0;i { count=0; sum=0; for(j=1;j<=M;j++) { sum+=s[i][j]; if(s[i][j]<60) count++; } s[i][M+1]=sum/M; s[i][M+2]=0; if(count>=2) s[i][M+2]='N'; } } void findgood(int s[ ][M+3]) { int i,j,count; for(i=0;i { count=0; for(j=1;j<=M;j++) if(s[i][j]>=85) count++; if(count==M || s[i][M+1]>=90) s[i][M+2]='Y'; } } 4．矩阵乘法 在主函数中输入m和n (0 C=A XB。并在主函数中输出矩阵C。 SAMPLE OUTPUT Input m and n (1~5)：4 3 Input A4 X 3： 1 2 3 4 5 6 7 8 9 10 11 12 Input B3 X 4： 1 2 3 4 4 5 6 8 7 8 9 10 C=AXB is： 30 36 42 50 66 81 96 116 102 126 150 182 138 171 204 248 #include #define N 5 void multiply(int a[][N],int b[][N],int c[][N],int m,int n); void main() { int a[N][N],b[N][N],c[N][N]; int i,j,m,n; printf("Input m and n(1～5):"); scanf("%d%d",&m,&n); printf("Input A %d*%d:\n",m,n); for(i=0;i for(j=0;j scanf("%d",&a[i][j]); printf("Input B %d*%d:\n",n,m); for(i=0;i for(j=0;j 1.输入十个整数，分别输出这10个数再计算并输出它们的和。 源程序： #include } printf("Input x:"); scanf("%d",&x); for(j=0;j<10;j++) { if(x==a[j]) { printf("%d\n",j); goto L1; } } printf("not found\n"); L1: ; } 运行结果： 3.输入一个整数n（1 工科数学分析上机作业 说明：以下两道题均是使用Matlab 语言，且在Matlab 7.0中运行通过。 1.（两个重要极限）计算下列函数的函数值并画出图形，观察两个重要极限值。 (1)y=f(x)=; (2)y=f(x)=. sin x x (1+x)1x 解：（1）求解过程如下： >> syms x >> y=limit(sin(x)/x) y = 1 >> ezplot(sin(x)/x,[-10*pi,10*pi]) >> ezplot(sin(x)/x,[-1*pi,1*pi]) 其图形如下： （2）求解过程如下：>> syms x >> y=(1+x)^(1/x) y = (1+x)^(1/x) >> y=limit((1+x)^(1/x)) y = exp(1) >> ezplot((1+x)^(1/x),[-1000,1000]) >> ezplot((1+x)^(1/x),[-10,10]) >> ezplot((1+x)^(1/x),[-1,1]) 其图像如下： 分析如下：（1）当x 取值为[-30,30]时，由该题的第一个图像可以看到，函数值在不断震荡，一会为正数，一会为负数。 而当x 取值为[-3，3]时，函数值始终大于0。当x 趋近于0时，由该题的第二个图像可以得到函数值为1。 另外，该结论也可以由夹逼法则证明，结果不变，当x 趋近于0时，函数值仍为1。 （2）由该题的三个图像可以知道，该函数在定义域内为单调递减函数。且由该题的第一和二个图像知道，当x 在 [0,10]区间内，函数递减趋势非常迅速。由该题的第三个图像知道，当x 趋于0 时，函数值为自然对数的底数 e ,即约为2.71828. 3.计算f(x)=, 12+1√2π ∫x 0e ?t 2/2dt 1?x ?3的函数值{f （0.1k ）;k=1,2,…,30}.计算结果取7位有效数字。 解：计算过程为： >> f1= @(t) exp(-(t).^2/2) f1 = @(t) exp(-(t).^2/2) >> for i=1:30 实验一： 1、已知华氏温度，根据转换公式计算相应的摄氏温度。 摄氏温度= 5*（华氏温度-32）/9 2．任意输入一个有五位数字的正整数，逆序输出每一数位上的数字。 如： 输入 12345 输出 5 4 3 2 1 #include"stdio.h" main() { double t,T; printf ("已知华氏温度，根据转换公式计算相应的摄氏温度"); scanf ("%lf",&t); T=(5*(t-32)/9); printf ("摄氏温度=%lf",T); } #include"stdio.h" main() { int value,a,b,c,d,e,m,n,h; printf ("任意输入一个有五位数字的正整数，逆序输出每一数位上的数字"); scanf ("%d",&value); a=value/10000; m=value%10000; b=m/1000; n=m%1000; c=n/100; h=n%100; d=h/10; e=h%10; printf ("%d%d%d%d%d",e,d,c,b,a); } 实验二： 1．一个成年人的身高与其父母的身高密切相关，根据以下人体标准身高预测公式计算自己或某个同学的身高，并跟实际身高比较，给出结论是否符合该公式（误 差在3厘米范围内即视为相符）。 男性身高=(父亲身高+母亲身高)×1.08÷2(厘米) 女性身高=(父亲身高×0.923+母亲身高)÷2(厘米） 2．从键盘上输入税前收入，根据以下公式及规则，计算个人所得税税额并输出。应纳个人所得税税额=(税前收入—起征点额度)*适用税率—速算扣除数 其中，起征点额度是3500元，其他如下： 应纳税所得额税率速算扣除数 不超过1500元的部分 3% 0 超过1500元至4500元的部分 10% 105 超过4500元至9000元的部分 20% 555 超过9000元至35000元的部分 25% 1005 超过35000元至55000元的部分 30% 2755 超过55000元至80000元的部分 35% 5505 超过80000的部分 45% 13505 如：某人月收入5000元，则应纳个人所得税税额为（5000-3500）* 10% - 105 = 45元。 #include 实验5 函数（1） 【实验题5-1】以下程序的功能：从键盘输入两个数，计算出他们的和，并将结果输出在屏幕上。调试运行程序，观察出现了什么错误，提出两种修改方案。 【实验题5-2】下程序的功能：输出3~1000之的间的素数。请选择合适的内容填入程序的空白处，并调试运行程序。 【实验题5-3】从键盘输入一个整数n,计算n^3，并将结果输出在屏幕上。要求计算n^3用函数实现。 【实验题5-4】从键盘输入一个整数n，计算2n+1与++n之和。要求计算两数之和用函数 实现。 注意写计算表达式的时候3*n之间的乘号不要丢掉，否则会提示：invalid suffix “n” on integer constant. 调用用户自定义函数是只写ifun=fun(n); 就可以了，不要在括号里加上int。Ifun=fun(int n) × 【实验题5-5】编写程序，计算Cmn=n!/(n-m)!m!。要求如下： 1）编写一个函数计算n！； 2）编写主函数，输入n和m, 调用第1）项中的函数完成计算； 3）输出并计算C58、C712的结果。 【实验题6-1】以下程序功能：用已知三角形三边，求三角形面积。但程序有错，请上机调试并改正错误。 【实验题6-2】以下程序运行的结果应该是a=14, b=15,temp=0。在程序空白处填写适当的编译预处理命令。 【实验题6-3】用for循环输出20个0~99之间的随机数。要求没输出五个随机数换一行。 【实验题6-4】从键盘输入一个字母，如果是大写字母就转换成小写字母输出；如果是小写字母就转换成大写字母输出。要求用库函数实现。 2016年理工大学优化方法上机大作业学院： 专业： 班级： 学号： ： 上机大作业1： 1.最速下降法： function f = fun(x) f = (1-x(1))^2 + 100*(x(2)-x(1)^2)^2; end function g = grad(x) g = zeros(2,1); g(1)=2*(x(1)-1)+400*x(1)*(x(1)^2-x(2)); g(2) = 200*(x(2)-x(1)^2); end function x_star = steepest(x0,eps) gk = grad(x0); res = norm(gk); k = 0; while res > eps && k<=1000 dk = -gk; ak =1; f0 = fun(x0); f1 = fun(x0+ak*dk); slope = dot(gk,dk); while f1 > f0 + 0.1*ak*slope ak = ak/4; xk = x0 + ak*dk; f1 = fun(xk); end k = k+1; x0 = xk; gk = grad(xk); res = norm(gk); fprintf('--The %d-th iter, the residual is %f\n',k,res); end x_star = xk; end >> clear >> x0=[0,0]'; >> eps=1e-4; >> x=steepest(x0,eps) 2.牛顿法： function f = fun(x) f = (1-x(1))^2 + 100*(x(2)-x(1)^2)^2; end function g = grad2(x) g = zeros(2,2); 第一次上机作业（2014）： 要求： ?独立完成，上机调试通过后上交 ?作业提交邮箱： ?邮件主题命名：班级学号姓名（第1次作业）， 例如：电1301班2王鹏（第1次作业） ?以附件形式将每次作业的所有程序源程序压缩打包后提交，压缩文件命名格式同上。 ?程序必须采用缩进格式 1．大写字母转换成小写字母 从键盘输入一个大写英文字母，输出相应的小写字母。 例：输入G 输出g 2．求平方根 输入一个实数x，计算并输出其平方根（保留1 位小数）。 例：输入17 输出The square root of 17.0 is 4.1 3．温度转换 设计一个程序将华氏温度转换成摄氏温度c = 5/9(f-32) a)输入华氏温度（实型） b)输出的摄氏温度结果保留两位小数 例：Please input Fahrenheit temperature: 76.8 The corresponding Celsius temperature is 24.89 4. 计算旅途时间 输入二个整数time1 和time2，表示火车的出发时间和到达时间，计算并输出旅途时间。 (输入的有效的时间范围是0000 到2359，不需要考虑出发时间晚于到达时间的情况。) 例：输入712 1411 （出发时间是7：12，到达时间是14：11） 输出The train journey time is 6 hours 59 minutes. 5. 数字加密 输入一个四位数，将其加密后输出。方法是将该数每一位上的数字加9，然后除以10 取余，做为该位上的新数字，最后将第1 位和第3 位上的数字互换，第2 位和第4 位上的数字互换，组成加密后的新数。 例：输入1257 输出The encrypted number is 4601 思考题：你能否编程找出谁做的好事？ 有四位同学中的一位做了好事，不留名，表扬信来了之后，校长问这四位是谁做的好事。 ?A说：不是我。 ?B说：是C。 ?C说：是D。 ?D说：他胡说。 已知三个人说的是真话，一个人说的是假话。现在要根据这些信息，找出做了好事的人。 第二次上机作业（2014）： 要求： ?独立完成，程序必须采用缩进格式，上机调试通过后上交 实验7 数组（1） 【实验题7-1】此程序的功能：给数组元素赋值并进行输出。但程序中有错误，请查找原因，并调试程序以实现程序功能。 【实验题7-2】以下程序的功能：求整型数组中元素的最大值及其下标。请选择合适的内容填入程序中空白处。 【实验题7-3】以下程序的功能：把整型数组里的元素按逆序进行排序并输出。请选择合适的内容填入程序中空白处。 【实验题7-4】编写一个程序，定义一个有五个元素的一维数组，从键盘对该数组的元素进行赋值后再显示在屏幕上。 【实验题7-5】编写程序对整型数组iA的各元素进行从大到小排序。 【实验题7-6】在数组iA和数组iB中分别存放了na和nb个由小到大的有序数。编写程序，完成把两个数组中的元素按从小到大的顺序归并存放在数组iC中。 【实验题7-7】假设数组iA中的元素按由小到大的顺序存放。编写程序把数组iA中相同的元素删除至一个，然后按照每行六个元素的形式输出数组iA中的各个元素值。 【实验题7-8】在二维数组中，若某一位置上的元素在该行中最大，而在该列中最小，则该元素为该二维数组的一个鞍点。编写程序要求从键盘输入一个二维数组，当鞍点存在时，把鞍点找出来。 #include printf("M",a[i][j]); printf("\n"); } for(i=0;i k=0; for(j=1;j if(a[i][j]>a[i][k]) k=j; for(j=0;j if(a[j][k] break; if(j==M) printf("%d,%d,%d\n",a[i][k],i,k); } return 0; } 上面这个题代码运行不了，可又找不到问题所在，头疼！！！ 实验8 数组（2） 【实验题8-1】此程序的功能：过滤字符数组中数字并进行输出。但程序中有错误，请查找原因，并调试程序以实现程序功能。 【实验题8-2】此程序的功能：统计字符数组中出现字符“s”的次数，并进行输出。但程序中有错误，请查找原因，并调试程序以实现程序功能。 【实验题8-3】以下程序的功能：连续输入一批“#”结束的字符串，遇“##”则全部结束，统计并输出每个字符串的有效长度。请选择适当的内容填入空白处。 【实验题8-4】 大连理工大学《工程抗震》大作业 题目1：底部剪力法。 钢筋混凝土5层框架经质量集中后计算简图如下图所示，各层高均为3m ， 集中于各楼层的重力荷载代表值分别为： 1500kN G =，2550kN G =，3580kN G =，4600kN G =，5450kN G =。结构阻尼比0.05ξ=，自振周期为10.55s T =，Ⅰ1类 场地类别，设计地震分组为第一组，抗震设防烈度为8度（设计基本地震加速度为0.30g ）。按底部剪力法计算结构在多遇地震时的水平地震作用及地震剪力。 3580kN =2550kN =1500kN =（a ）计算简图 4600kN =5450kN = 解：查《建筑设计抗震规范》表5.1.4知，8度多遇地震，αmax=设计地震分组为第一组， Ι类场地，取Tg= Tg=＜T1=＜5Tg= α1=（Tg/T1）r η2αmax =（）××=≈ 查《建筑设计抗震规范》表5.2.1知，T 1=＞=×= 取δn=T1+=×+= 总水平地震作用标准值： F EK =α1Geq=×（500+550+580+600+450）×85%= 各楼层水平地震作用标准值： Fi=G i H i F EK (1-δn)/∑G j H j (i=1,2,3n) ∑G j H j =500×3 +550×6+580×9+600×12+450×15=23970KN ·m F 1=[500×3××]/23970= F 2=[550×6××]/23970= F 3=[580×9××]/23970= F 4=[600×12××]/23970= F 5=[450×15××]/23970= 计算各楼层的层间地震剪力 V 1= F 1+ F 2+ F 3+ F 4+ F 5=++++= V 2= F 2+ F 3+ F 4+ F 5=+++=152KN V 3= F 3+ F 4+ F 5=++= V 4= F 4+ F 5=+= V 5=F 5= 题目3：怎样判断土的液化如何确定土的液化严重程度，并简述抗液化措施。 答：饱和松散的砂土或粉土（不含黄土），地震时易发生液化现象，使地基承载力丧失或减弱，甚至喷水冒砂，这种现象一般称为砂土液化或地基土液化。其产生的机理为：地下水位以下的饱和砂土和粉土颗粒在地震作用下，土颗粒之间有变密的趋势。因空隙水不能及时排出，土颗粒就处于悬浮状态，形成如同液体一样的现象，即所谓的土的液化现象。地基土液化判别过程可以分为初步判断和标准贯入试验判别两大步骤。下面分别予以介绍。 1、初步判断 饱和的砂土或粉土（不含黄土）当符合下列条件之一时，可初步判别为不液化或不考虑液化影响： （1）地质年代为第四纪晚更新世（Q3）及其以前时且处于烈度7度或者8度地区时可判为不液化土。 （2）粉土的粘粒（粒径<0.005mm ）含量百分率当烈度为7度时大于10％、当烈度为8度时大于13％、当烈度为9度时大于16％，可判为不液化土。 （3）浅埋天然地基，当地下水位深度和覆盖非液化土层厚度满足下式之一时，可不考虑液化影响。 03w b d d d >+- 02 u b d d d >+- 函数定义： % 目标函数 function f = fun(x) fm=0; for i=1:499 fmi = (1-x(i))^2 + 100*(x(i+1)-x(i)^2)^2; fm=fm+fmi; end f =fm; end % 梯度 function g = grad(x) g = zeros(500,1); g(1)=2*(x(1)-1)+400*x(1)*(x(1)^2-x(2)); for i=2:499 g(i)=2*(x(i)-1)+400*x(i)*(x(i)^2-x(i+1))+200*(x(i)-x(i-1)^2); end g(500) = 200*(x(500)-x(499)^2); end % 二阶梯度 function g = grad2(x) g = zeros(500,500); g(1,1)=2+400*(3*x(1)^2-x(2)); g(1,2)=-400*x(1); for i=3:500 g(1,i)=0; end for i=1:498 g(500,i)=0; end g(500,499)=-400*x(499); g(500,500)=200; for i=2:499 for j=1:500 if j==i-1 g(i,j)= -400*x(i-1); elseif j==i g(i,j)= 2+400*(3*x(i)^2-x(i+1))+200; elseif j==i+1 g(i,j)= -400*x(i); else g(i,j)=0; end end end end 1.最速下降法 function x_star = steepest(x0,eps) gk = grad(x0); res = norm(gk); k = 0; while res > eps && k<=50000 dk = -gk; 答案仅供参考 实验一 4．设计main函数实现从键盘上输入两个整型变量a、b的值，交换它们的值并输出。#include printf("x=%d,y=%d\n",x,y); getch(); } 实验二 1．编写程序，从键盘上输入一个整数（例如560）表示分钟，将其换算成用小时和分钟表示，然后输出至屏幕。 #include 1．将顺序表逆置，要求用最少的附加空间。 参考答案 #include 第一次上机作业 1.利用Matlab自带命令产生1000个均匀随机变量服从U(0,1)。 >>unifrnd(0,1,20,50) ans= Columns1through10 0.81470.65570.43870.75130.35170.16220.10670.85300.78030.5470 0.90580.03570.38160.25510.83080.79430.96190.62210.38970.2963 0.12700.84910.76550.50600.58530.31120.00460.35100.24170.7447 0.91340.93400.79520.69910.54970.52850.77490.51320.40390.1890 0.63240.67870.18690.89090.91720.16560.81730.40180.09650.6868 0.09750.75770.48980.95930.28580.60200.86870.07600.13200.1835 0.27850.74310.44560.54720.75720.26300.08440.23990.94210.3685 0.54690.39220.64630.13860.75370.65410.39980.12330.95610.6256 0.95750.65550.70940.14930.38040.68920.25990.18390.57520.7802 0.96490.17120.75470.25750.56780.74820.80010.24000.05980.0811 0.15760.70600.27600.84070.07590.45050.43140.41730.23480.9294 0.97060.03180.67970.25430.05400.08380.91060.04970.35320.7757 0.95720.27690.65510.81430.53080.22900.18180.90270.82120.4868 0.48540.04620.16260.24350.77920.91330.26380.94480.01540.4359 0.80030.09710.11900.92930.93400.15240.14550.49090.04300.4468 0.14190.82350.49840.35000.12990.82580.13610.48930.16900.3063 0.42180.69480.95970.19660.56880.53830.86930.33770.64910.5085 0.91570.31710.34040.25110.46940.99610.57970.90010.73170.5108 0.79220.95020.58530.61600.01190.07820.54990.36920.64770.8176 0.95950.03440.22380.47330.33710.44270.14500.11120.45090.7948 Columns11through20 0.64430.31110.08550.03770.03050.05960.17340.95160.03260.2518 0.37860.92340.26250.88520.74410.68200.39090.92030.56120.2904 0.81160.43020.80100.91330.50000.04240.83140.05270.88190.6171 0.53280.18480.02920.79620.47990.07140.80340.73790.66920.2653 0.35070.90490.92890.09870.90470.52160.06050.26910.19040.8244 0.93900.97970.73030.26190.60990.09670.39930.42280.36890.9827 0.87590.43890.48860.33540.61770.81810.52690.54790.46070.7302 C语言上机考试题—01 1. //#include 绪论部分： 2、简述环境问题的分类？（10分） 答：环境问题是多方面的，但大致可分为两类：原生环境问题和次生环境问题。由自然力引起的为原生环境问题，也称为第一环境问题。由于人类生产和生活引起生态系统破坏和环境污染，反过来又危及人类自身和生存和发展的现象，为次生环境问题，也叫第二环境问题。原生环境问题和次生环境问题很难截然分开，它们之间常常存在着某种程度的因果关系和相互作用。 4、什么是环境化学，学习环境化学有什么意义？（10分） 答：环境化学是一门研究有害化学物质在环境介质中的存在、化学特性、行为和效应及其控制的化学原理和方法的科学。 意义：用来掌握污染来源，消除和控制污染，确定环境保护决策，以及提供科学依据诸方面都起着重要的作用。 5、简述环境化学的分支学科。（10分） 答：主要包括6类。 ①环境分析化学：是研究化学品的形态、价态、结构、样品前处理和痕量分析的学科。 ②环境污染化学：大气、水体和土壤环境化学，元素循环的化学过程。 ③污染控制化学：主要研究与污染控制有关的化学机制及工艺技术中化学基础性问题。 ④污染生态化学：是研究化学污染物在生态系统中产生生态效应的化学过程的学科。 ⑤环境计算化学：主要利用有效的数学近似以及电脑程序计算分子的性质。 ⑥环境生物化学：是研究环境化学品对生命影响的学科。 第一章： 1、地球环境主要由哪些圈层构成？英文单词？各之间有什么联系？各有哪些性 质？（10分） 答：地球环境主要由大气圈（atmosphere）、水圈（hydrosphere）、土壤圈（pedosphere）、岩石圈（lithosphere）和生物圈（biosphere）构成。 联系：大气圈、水圈、土壤圈和生物圈共同组成了地球环境系统，每个圈层都离不开 2016年大连理工大学优化 方法上机大作业 -标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII 2016年大连理工大学优化方法上机大作业学院： 专业： 班级： 学号： 姓名： 上机大作业1： 1.最速下降法： function f = fun(x) f = (1-x(1))^2 + 100*(x(2)-x(1)^2)^2; end function g = grad(x) g = zeros(2,1); g(1)=2*(x(1)-1)+400*x(1)*(x(1)^2-x(2)); g(2) = 200*(x(2)-x(1)^2); end function x_star = steepest(x0,eps) gk = grad(x0); res = norm(gk); k = 0; while res > eps && k<=1000 dk = -gk; ak =1; f0 = fun(x0); f1 = fun(x0+ak*dk); slope = dot(gk,dk); while f1 > f0 + 0.1*ak*slope ak = ak/4; xk = x0 + ak*dk; f1 = fun(xk); end k = k+1; x0 = xk; gk = grad(xk); res = norm(gk); fprintf('--The %d-th iter, the residual is %f\n',k,res); end x_star = xk; end >> clear >> x0=[0,0]'; >> eps=1e-4; >> x=steepest(x0,eps) 第五次上机作业：（请用指针实现） 1．数据移位 有n（n<=10）个整数，使前面各数顺序向后移m（m<=n）个位置，最后m 个数变成最前面的m个数，如图所示。编写函数shift(a,n,m)实现上述功能，在主函数中输入n个整数和输出调整后的n个数。 SAMPLE OUTPUT Input n(n<=10)：10 Input 10 integers：1 2 3 4 5 6 7 8 9 10 Output the original：1 2 3 4 5 6 7 8 9 10 Input m(m<=10)：4 Output the shifted：7 8 9 10 1 2 3 4 5 6 #include 矩阵与数值分析上机作业 学校：大连理工大学 学院： 班级： 姓名： 学号： 授课老师： 注：编程语言Matlab 程序： Norm.m函数 function s=Norm(x,m) %求向量x的范数 %m取1,2,inf分别表示1,2，无穷范数 n=length(x); s=0; switch m case 1 %1-范数 for i=1:n s=s+abs(x(i)); end case 2 %2-范数 for i=1:n s=s+x(i)^2; end s=sqrt(s); case inf %无穷-范数 s=max(abs(x)); end 计算向量x，y的范数 Test1.m clear all; clc; n1=10;n2=100;n3=1000; x1=1./[1:n1]';x2=1./[1:n2]';x3=1./[1:n3]'; y1=[1:n1]';y2=[1:n2]';y3=[1:n3]'; disp('n=10时'); disp('x的1-范数:');disp(Norm(x1,1)); disp('x的2-范数:');disp(Norm(x1,2)); disp('x的无穷-范数:');disp(Norm(x1,inf)); disp('y的1-范数:');disp(Norm(y1,1)); disp('y的2-范数:');disp(Norm(y1,2)); disp('y的无穷-范数:');disp(Norm(y1,inf)); disp('n=100时'); disp('x的1-范数:');disp(Norm(x2,1)); disp('x的2-范数:');disp(Norm(x2,2)); disp('x的无穷-范数:');disp(Norm(x2,inf)); disp('y的1-范数:');disp(Norm(y2,1)); disp('y的2-范数:');disp(Norm(y2,2)); disp('y的无穷-范数:');disp(Norm(y2,inf)); disp('n=1000时'); disp('x的1-范数:');disp(Norm(x3,1)); disp('x的2-范数:');disp(Norm(x3,2)); disp('x的无穷-范数:');disp(Norm(x3,inf)); disp('y的1-范数:');disp(Norm(y3,1)); disp('y的2-范数:');disp(Norm(y3,2)); disp('y的无穷-范数:');disp(Norm(y3,inf)); 运行结果： n=10时 x的1-范数:2.9290；x的2-范数:1.2449； x的无穷-范数:1 y的1-范数:55； y的2-范数:19.6214； y的无穷-范数:10 n=100时 x的1-范数:5.1874；x的2-范数: 1.2787； x的无穷-范数:1 y的1-范数:5050； y的2-范数:581.6786； y的无穷-范数:100 n=1000时 x的1-范数:7.4855； x的2-范数:1.2822； x的无穷-范数:1 y的1-范数: 500500； y的2-范数:1.8271e+004；y的无穷-范数:1000 程序 Test2.m clear all; clc; n=100;%区间 h=2*10^(-15)/n;%步长 x=-10^(-15):h:10^(-15); %第一种原函数 上机考试习题集 1.编写程序，一次输入5个学生的3门课成绩，每输入一个学生的3门成绩后，立即统计并输出该学生的总分和平均分。 #include "stdio.h" void main() { int score[5][3], i,j; float sum[5],average[5]; printf("请分别输入三个学生的成绩：\n"); for(i=0;i<5;i++) { sum[i]=0; for(j=0;j<3;j++) {scanf("%d",&score[i][j]); average[i]=0; sum[i]+=score[i][j]; average[i]=sum[i]/5; } } for(i=0;i<5;i++) printf("第%d个学生的总分为：%f,平均成绩为：%f\n",i+1,sum[i],average[i]); } 2 、编写程序，提示用户从键盘输入一个double类型的数，使该数保留小数点后两位，对第三位小数进行四舍五入处理，然后在显示器上输出此数，以便验证处理是否正确。 #includeC语言上机作业
Removed_大连理工大学工科数学分析上机作业
C语言上机作业题及答案
5-6次C语言上机作业
2016年大连理工大学优化方法上机大作业
C语言上机作业试题5套含答案
7~8次C语言上机作业
大连理工大学(工程抗震)大作业
大连理工大学优化方法上机大作业程序
扬州大学C语言上机作业1-9整理
大连理工大学数据结构(一)上机作业答案——张老师
大连理工大学概率上机作业
C语言期末上机考试题及答案
大连理工大学-环境化学-所有作业答案
大连理工大学优化方法上机大作业
大连理工c语言第五次上机作业参考答案
大连理工大学矩阵与数值分析上机作业
c语言上机考试训练题库及答案
相关文档
- 大连理工大学优化方法上机作业
- 大连理工 程序设计 第5次上机作业
- 大连理工大学线代上机
- 2016年大连理工大学优化方法上机大作业
- 第一次上机作业参考答案
- 大连理工大学线性代数上机作业PPT
- 2020年大连理工大学概率与统计上机作业答案
- 大连理工大学矩阵与数值分析上机作业
- 大连理工大学矩阵与数值分析上机作业13388
- 大连理工大学概率上机作业
- 大连理工c语言第一次上机作业参考答案
- 2016年大连理工大学优化方法上机大作业
- 大连理工大学矩阵与数值分析上机作业
- Removed_大连理工大学工科数学分析上机作业
- 大连理工大学概率上机作业(含答案)
- 大连理工大学矩阵与数值分析上机作业
- 大连理工大学矩阵分析matlab上机作业
- 大连理工大学优化方法上机作业
- 大连理工大学线代上机作业三原文件
- 大连理工大学数据结构(一)上机作业答案——张老师