大连理工离线作业《可编程控制器》

网络教育学院《可编程控制器》大作业

题目：十字路口交通灯控制设计

学习中心：]

层次：高起专

专业：电力系统自动化技术

年级： 1 6年秋季

学号：

学生姓名：

题目五：十字路口交通灯控制设计

起动后，南北红灯亮并维持30s。在南北红灯亮的同时，东西绿灯也亮，东西绿灯亮25s 后闪亮，3s后熄灭，东西黄灯亮，黄灯亮2s后，东西红灯亮，与此同时，南北红灯灭，南北绿灯亮。南北绿灯亮25s后闪亮，3s后熄灭，南北黄灯亮，黄灯亮2s后，南北红灯亮，东西红灯灭，东西绿灯亮。依次循环。

十字路口交通灯控制示意图及时序图如下图所示。

设计要求：（1）首先对可编程序控制器（PLC）的产生与发展、主要性能指标、分

类、特点、功能与应用领域等进行简要介绍；

（2）设计选用西门子S7－200 系列PLC，对其I/O口进行分配，列出

PLC控制程序（梯形图进行截图，语句表可直接拷贝）并对程序作

出解释；

（3）总结：需要说明的问题以及设计的心得体会。

十字路口交通灯控制设计

摘要

随着汽车进入家庭步伐的加快和城市汽车数量的增多，城市道路交通问题显得越来越重要。解决好十字路口交通信号灯控制问题是保障交通有序、安全、快速运行的重要环节。但现有的十字路口交通信号灯控制系统大都采用继电器或单片机实现，且都是单一的固有时序控制，不能够根据实际路面车流量情况进行调节控制，存在着功能少、可靠性差、维护量大等缺点。为了弥补原交通信号灯系统存在的种种缺点，本文设计了基于PLC控制的交通信号灯控制系统。该系统选用的可编程逻辑控制器是德国西门子公司的S7-200，具有一定的智能性，即可以根据路面车流量大小对十字路口的交通信号灯按高峰期、正常期和晚间几个时段进行分时控制。

1 设计背景

1.1 背景概述

随着汽车进入家庭步伐的加快和城市汽车数量的增多，城市道路交通问题显得越来越重要。马路上经常会看到这种现象:一旦整个路口的交通信号灯出现故障，若没有交警的及时疏导，该路口就会塞得一塌糊涂，甚至造成严重的交通事故。原交通信号控制大都采用继电器或单片机实现，存在着功能少、可靠性差、维护量大等缺点,越来越不能适应城市道路交通高速发展的要求。另外，根据人车流量的多少，可能随时增加路口的交通信号，比如增加转弯或人行道交通信号，原有系统的制约性就更加明显了。为了弥补原交通信号灯系统存在的以上缺点，我们引入了基于PLC控制的交通信号灯控制系统。

本文对十字路口交通信号灯控制系统，运用可编程逻辑器件PLC做了软件与硬件的设计，能基本达到控制要求。系统仅实现了小型PLC系统的一个雏形，在完善各项功能方面都还需要进一步的分析、研究和调试工作。如果进一步结合工业控制的要求，形成一个较为成型的产品，则需要作更多、更深入的研究。

1.2 可编程逻辑控制器简介

可编程逻辑控制器(Programmable Logic Controller,简称 PLC) 根据国际电工委员会(IEC)在1987年的可编程控制器国际标准第三稿中，对其作了如下定义：“可编程控制器是一种数字运算操作的电子系统，专为在工业环境应用而设计的。它采用可编程序的存储器，用来在其内部存储执行逻辑运算，顺序控制，定时，计数与算术运算等操作的指令，并通过数字式、模拟式的输入和输出，控

制各种类型的机械或生产过程。可编程控制器及其有关外部设备，都应按易于使工业控制系统形成一个整体，易于扩充其功能的原则设计。”可编程控制器作为目前工业自动化的重要基础设备，被称为“工业自动化三大支柱性产业之一”，在各工业生产领域发挥着愈来愈大的作用。

2 十字路口交通信号灯PLC控制系统简介

2.1 控制对象及要求

2.1.1 控制对象

本系统的控制对象有八个，分别是：

东西方向红灯（R—EW）两个；

南北方向红灯 (R—SN) 两个；

东西方向黄灯（Y—EW）两个；

南北方向黄灯 (Y—SN) 两个；

东西方向绿灯（G—EW）两个；

南北方向绿灯 (G—SN) 两个；

东西方向左转弯绿灯(L—EW)两个；

南北方向左转弯绿灯(L—SN)两个。

2.1.2 控制要求

1、系统工作受开关控制，起动开关 ON 则系统工作；起动开关 OFF 则系统停止工作；

2、交通信号灯按高峰时段、正常时段及晚上时段进行控制，这三个时段的的时序分配如图1所示；

3、在高峰时段，交通信号灯按图2所示时序控制；

4、在正常时段，交通信号灯按图3 所示时序控制；

5、晚上时段按提示警告方式运行，规律为：东、南、西、北四个黄灯全部闪亮，其余灯全部熄灭，黄灯闪亮按亮 0.4 秒，暗 0.6 秒的规律反复循环。

2.2 系统简介

本系统是一个十字路口交通灯的PLC控制系统，利用西门子公司的S7-200可编程逻辑控制器对十字路口的交通灯进行控制。本系统具有一定的智能性，即它可以对交通灯按高峰期、正常期及晚间几个时段进行分段控制。高峰期的控制方案为：

（1）南北方向左转弯灯和南北南北方向红灯同时亮10秒，同时东西方向红灯亮；

（2）南北方向绿灯亮35秒，东西方向红灯继续亮；

（3）南北方向黄灯闪烁5秒；东西方向红灯继续亮；

（4）东西方向左转弯绿灯和南北方向红灯同时亮10秒，东西方向红灯继续亮；

（5）东西方向绿灯亮25秒，南北方向红灯继续亮；

（6）东西方向黄灯闪烁5秒，南北方向红灯继续亮，然后跳至第（1）步依次循环。

正常期的控制方案为：

（1）南北方向左转弯灯和南北南北方向红灯同时亮10秒，同时东西方向红灯亮；

（2）南北方向绿灯亮30秒，东西方向红灯继续亮；

（3）南北方向黄灯闪烁5秒；东西方向红灯继续亮；

（4）东西方向左转弯绿灯和南北方向红灯同时亮10秒，东西方向红灯继续亮；

（5）东西方向绿灯亮30秒，南北方向红灯继续亮；

（6）东西方向黄灯闪烁5秒，南北方向红灯继续亮，然后跳至第（1）步依次循环。

晚间的控制方案为：

东、南、西、北四个黄灯全部闪亮，其余灯全部熄灭，黄灯闪亮按亮 0.4 秒，暗 0.6 秒的规律反复循环。

2.3 硬件选型

城市道路交通信号控制是典型的开关量顺序控制，采用PLC能充分利用它的优点。在这里我们采用德国西门子公司的S7-200可编程控制器，它是积木式结构，安装比较方便，中央处理单元和信号模板有多种类型，另外还具有如位控单元、PD调节等特殊功能模块。根据本系统输入点数及控制要求，中央处理单元可选用CPU224，该CPU板上本身具有10个数字量输入点，6个非隔离数字量输出点，最多能够带8个数字量信号模板。

电源模块将交流电源转换成供CPU，存储器等所有扩展模块使用的直流电源，是整个PLC系统的能源供给中心，它的好坏直接影响到PLC的稳定性和可靠。

S7-200属于小型PLC，电源模块与CPU模块封装在一起，通过连接总线为本机和扩展模块提供+5V（DC）电源。同时，还可通过端子向外输出一个+24V（DC）电源，供本机输入点和扩展模块继电器线圈使用。需注意的是，从资料中我们了解到，外部电源不可与S7-200的传感器电源并联使用。否则，交会导致两个电源的竟争而影响它们各自的输出，缩短其使用寿命，使得一个或两个电源同时失效，使PLC系统产生不正确的操作。正确的使用方法是S7-200的传感器电源和外部电源应该在不同的点上提供电源，而两者之间只能有一个会共连接点。

由于根据控制要求所确定的输入输出点分别人二个和九个，由于我们是以一个路口信号单独控制为例，考虑到够用为准。所以我们选择了CPU224这一具有较强控制功能的控制器。

另外，在硬件选型时，不要忘记完成现场测试及软件编程时所需的一些设备。

综上，得到系统硬件配置如表1所示：

表1 硬件配置表

名称数量

DC24V电源1

CPU224 1

PC/PPI编程电缆 1

STEP7编程软件 1

PC机 1

3 系统I/O分配

分析PLC的输入和输出信号,在满足控制要求的前提下,要尽量减少占用PLC 的I/O点。由系统控制要求可见，由控制开关输入的启、停信号是输入信号。由PLC的输出信号控制各指示灯的亮、灭。在交通灯布置图中，南北方向的三色灯共六盏，同颜色的灯在同一时间亮、灭；所以，可将同色灯两两并联，用一个输出信号控制。同理，东西方向的三色灯也依次设计。再加上东西方向左转的三色灯共九盏，所以其占9个输出点。由此可得系统I/O分配如表2所示：表2 系统I/O分配表

输入/输出设备/器件名称 I/O地址符号名数据类型

输入校正当前时钟I0.0 SB0 BOOL

程序启停按钮I0.1 SB1 BOOL

输出东西方向绿灯Q0.0 Q1 DINT

东西方向黄灯Q0.1 Q2 DINT

东西方向红灯Q0.2 Q3 DINT

南北方向绿灯Q0.3 Q4 INT

南北方向黄灯Q0.4 Q5 INT

南北方向红灯Q0.5 Q6 INT

东西方向左转弯灯Q0.6 Q7 INT

南北方向左转弯灯Q0.7 Q8 INT

4 PLC端子接线图

根据I/O表及PLC的配置图很容易就可以得到PLC端子接线图4如下所示：

5 软件设计

本控制系统的控制原理是：用一路数字量的不同输入状态来判定是否对时钟进行初始化，用一路数字量的不同输入状态分别用作程序的启动和停止控制，每一方向有红、黄、绿及转弯四种信号灯，分别对应四位数字量输出，两个方向共有8位数定量输出；在某一方向用两个延时脉冲定时器分别控制该方向黄灯闪烁的亮、灭时间，根据道路人车流量多少，分别设置各信号灯亮灭时间的长短，通过6个定时器依次交替工作，就可实现各方向交通信号灯的顺序工作。本文所设计的软件由一个主程序和四个子程序（时钟初始化子程序，晚间时段交通灯控制子程序，正常时段交通灯控制子程序和高峰时段交通灯控制子程序）组成。主程序主要任务包括：读取两个开关状态，根据开关的不同状态做出相应的处理，当开关SB0闭合时则对时钟进行初始化，反之则不对时钟进行初始化；当开关SB1闭合时，则读取时钟值，并做处理，根据处理后的时钟值的大小判定当前时间是属于哪个时间段，并调用相应的交通灯控制子程序，反之，则停止程序的运行主程序流程图如图5所示。晚间时段的控制规律为：两个方向的四个黄灯均按亮0.4秒灭0.6秒的规律闪烁，其余的交通灯全灭程序中将用到两个定时器T37和T38，各定时器的功能如表3所示。正常时段的控制方案结构图如图6所示，程序中将用到8个定时T37-T44，各定时器的功能如表4所示。高峰时段的控制方案结构图如图7所示，程序中将用到8个定时T37-T44，各定时器的功能如表5所示。

该程序实现了信号由东西左转、东西直行、南北直行依次循环变化。其优势思路简单，容易理解，对时钟的校正以及各时段的起始时间和终止时间的修改方便。如路口要求在晚上10：00以后实行各方向黄色信号灯闪烁功能，只需要将

实时采集PLC的时钟信号作为一个子程序的跳转条件，再增加一段闪光程序即可。如果需要将几个路口集中到一台PLC控制，根据实际需要的I/O点数，硬件上再增加相应的数字量输出模板即可。需要指出的是，用PLC实现城市道路关通信号控制，最好几个路口共用一套PLC，这样可以大大降低工程成本。

表3 晚间时段各定时器一个循环中的功能明细表

定时器t0 t1 T2

T37

定时0.4秒开始定时，黄灯亮定时到，输出ON且保持；黄灯灭开始下一次循环的定时

T38

定时1秒开始定时继续定时定时到，输出ON，随即复位开始下一次循环的定时，黄灯亮。

表4 正常时段各定时器一个循环中的功能明细表

定时器t0 t1 t2 t3 t4 t5 t6

T37

定时10秒开始定时,南北转弯灯、南北红灯、东西红灯亮。定时到，输出ON且保持；南北转弯灯灭，南北绿灯亮，东西红灯继续亮。 ON ON ON ON 开始下一个循环定时

T38

定时40秒开始定时继续定时定时到，输出ON且保持；南北绿灯灭，南北黄灯闪烁，东西红灯继续亮。ON ON ON 开始下一个循环定时T39

定时45秒开始定时继续定时继续定时定时到，输出ON且保持；南北黄灯灭，东西转弯灯、南北红灯亮，东西红灯继续亮。ON ON 开始下一

个循环定时

T40

定时55秒开始定时继续定时继续定时继续定时定时到，输出ON且保持；东西转弯、东西红灯灭，东西绿灯亮，南北红灯继续亮。ON 开始下一个循环定时

T41

时到，输出ON且保持；东西绿灯灭，东西黄灯闪烁，南北红灯继续亮。开始下一个循环定时

T42

定时90秒开始定时继续定时继续定时继续定时继续定时继续定时定时到，输出ON，随即复位开始下一次循环定时；东西黄灯、南北红灯灭，南北转弯灯、南北红灯、东西红灯亮。

表5 高峰时段各定时器一个循环中的功能明细表

定时器t0 t1 t2 t3 t4 t5 t6

T37

定时10秒开始定时,南北转弯灯、南北红灯、东西红灯亮。定时到，输出ON且保持；南北转弯灯灭，南北绿灯亮，东西红灯继续亮。 ON ON ON ON 开始下一个循环定时

T38

定时45秒开始定时继续定时定时到，输出ON且保持；南北绿灯灭，南北黄灯闪烁，东西红灯继续亮。ON ON ON 开始下一个循环定时T39

定时50秒开始定时继续定时继续定时定时到，输出ON且保持；南北黄灯灭，东西转弯灯、南北红灯亮，东西红灯继续亮。ON ON 开始下一

个循环定时

T40

定时60秒开始定时继续定时继续定时继续定时定时到，输出ON且保持；东西转弯、东西红灯灭，东西绿灯亮，南北红灯继续亮。ON 开始下一个循环定时

T41

定时85秒开始定时继续定时继续定时继续定时继续定时定时到，输出ON且保持；东西绿灯灭，东西黄灯闪烁，南北红灯继续亮。开始下一个循环定时

T42

续定时定时到，输出ON，随即复位开始下一次循环定时；东西黄灯、南北红灯灭，南北转弯灯、南北红灯、东西红灯亮。

心得体会：

本学期在老师的带领下，学习了一门新的课程《可编程控制器》也叫PLC，短暂的十几周课程里，经过老师的耐心讲解、指导和学习，我受益匪浅，学到了不少本专业的知识。

在没有接触这门课之前，也曾经在和自动化专业的同学交流时听过PLC这个词，并且从社会工作中已经接触过PLC控制是当今自动控制的主流，目前自动化机床大多采用PLC控制。但当时并未很清楚的了解它，以及怎样运用。

通过本学期学习，使我对PLC有了一定认识，为我今后的工作提供了很好的理论支持，相信在今后的工作中我将理论与实践进一步结合，会有更好的效果。最后感谢老师的辛勤付出。

大连理工大学优化方法上机大作业程序

函数定义： % 目标函数 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;

Removed_大连理工大学工科数学分析上机作业

工科数学分析上机作业 说明：以下两道题均是使用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

2016年大连理工大学优化方法上机大作业

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);

大连理工大学(工程抗震)大作业

大连理工大学《工程抗震》大作业

题目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 >+-

大工15春《建筑制图》在线作业1及答案

大工15春《建筑制图》在 线作业1及答案 -标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

建筑制图大工15春《建筑制图》在线作业1及答案 一,单选题 1. A. B. C. D. 若空间一点位于某一坐标轴上，则该点的三个坐标值（）。均不为零只有一个为零必有两个为零必有三个为零 正确答案：C 2.A. B. C. D. 某平面图形在 V 面、W 面上的投影均为直线，该平面为（）。水平面正平面侧平面铅垂面 正确答案：A 3.A. B. C. D. 下列选项中，不可能是圆锥体截交线的是（）。等腰三角形圆椭圆正方形 正确答案：D 4. A. B. C. D. 下列选项中，哪项不是投影三要素的内容。（）空间被投影物体投射线投影面投射中心 正确答案：D 5. A. B. C. D. 某点的坐标值为（10,20,30），则该点到 H 面的距离为（）。 10 20 25 30 正确答案：D 二,多选题 1. A. B. C. D. 下列选项中，可能是圆锥上的截交线的有（）。三角形圆形椭圆形矩形 正确答案：ABC 2. A. B. C. D. 平面截切圆球所得的截交线总是圆，但是截交线圆的投影可以是（）。直线圆椭圆三角形 正确答案：ABC 三,判断题

1. 轴测投影分为正轴测投影和斜轴测投影两类，正轴测投影属于平行投影，而斜轴测投影不属于平行投影。 A. 错误 B. 正确 正确答案：A 2. 投影轴上的点的投影中，有一个投影在投影轴上；另两个投影在原点处。 A. 错误 B. 正确 正确答案：A 3. 棱柱体的截交线可能是椭圆形。 A. 错误 B. 正确 正确答案：A 4. 空间平行的两直线段，其轴测投影可能不再平行了。 A. 错误 B. 正确 正确答案：A 5. 平面投影法包括斜投影、平行投影、中心投影三种。 A. 错误 B. 正确 正确答案：A 6. 三棱锥的截交线可能是三角形。 A. 错误 B. 正确 正确答案：B 7. 某两条直线在 H 面和 V 面上的投影均平行，则这两条直线在空间上一定互相平行。 A. 错误 B. 正确 正确答案：A 8. 一般位置直线的三面投影均小于其实长。 A. 错误 B. 正确 正确答案：B 9. 两个立体互相贯穿时称为互贯。 A. 错误 B. 正确 正确答案：B 10. 圆锥的母线是直线段。 A. 错误 B. 正确 正确答案：B 11. A 点不在直线 BC 上，则 A 点的投影不可能在直线 BC 的投影上。 A. 错误 B. 正确

大连理工大学数据结构(一)上机作业答案——张老师

1．将顺序表逆置，要求用最少的附加空间。 参考答案 #include #include #include #define OK 1 #define ERROR 0 #define INFEASIBLE -1 #define OVERFLOW -2 typedef int ElemType; typedef int Status; #define LIST_INIT_SIZE 100 #define LISTTINCREMENT 10 typedef struct{ ElemType *elem; int length; int listsize; }SqList; //创建空顺序表 Status InitList_Sq(SqList &L){ L.elem=(ElemType*)malloc(LIST_INIT_SIZE*sizeof(ElemType)); if(!L.elem)exit(OVERFLOW); L.length=0; L.listsize=LIST_INIT_SIZE; return OK; } //创建顺序表，插入元素 void ListInput_Sq(SqList &L){ int n,i; printf("input the length of Sqlist:"); scanf("%d",&n); L.length=n; for(i=0;i

大连理工大学概率上机作业

第一次上机作业 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

大连理工大学-环境化学-所有作业答案

绪论部分： 2、简述环境问题的分类？（10分） 答：环境问题是多方面的，但大致可分为两类：原生环境问题和次生环境问题。由自然力引起的为原生环境问题，也称为第一环境问题。由于人类生产和生活引起生态系统破坏和环境污染，反过来又危及人类自身和生存和发展的现象，为次生环境问题，也叫第二环境问题。原生环境问题和次生环境问题很难截然分开，它们之间常常存在着某种程度的因果关系和相互作用。 4、什么是环境化学，学习环境化学有什么意义？（10分） 答：环境化学是一门研究有害化学物质在环境介质中的存在、化学特性、行为和效应及其控制的化学原理和方法的科学。 意义：用来掌握污染来源，消除和控制污染，确定环境保护决策，以及提供科学依据诸方面都起着重要的作用。 5、简述环境化学的分支学科。（10分） 答：主要包括6类。 ①环境分析化学：是研究化学品的形态、价态、结构、样品前处理和痕量分析的学科。 ②环境污染化学：大气、水体和土壤环境化学，元素循环的化学过程。 ③污染控制化学：主要研究与污染控制有关的化学机制及工艺技术中化学基础性问题。 ④污染生态化学：是研究化学污染物在生态系统中产生生态效应的化学过程的学科。 ⑤环境计算化学：主要利用有效的数学近似以及电脑程序计算分子的性质。 ⑥环境生物化学：是研究环境化学品对生命影响的学科。 第一章： 1、地球环境主要由哪些圈层构成？英文单词？各之间有什么联系？各有哪些性 质？（10分） 答：地球环境主要由大气圈（atmosphere）、水圈（hydrosphere）、土壤圈（pedosphere）、岩石圈（lithosphere）和生物圈（biosphere）构成。 联系：大气圈、水圈、土壤圈和生物圈共同组成了地球环境系统，每个圈层都离不开

大工15春管理学在线作业1答案

大工15春《管理学》在线测试1 一、单选题（共10 道试题，共30 分。） 1. （）是整个管理系统的基础。 A. 高层管理者 B. 中层管理者 C. 基层管理者 D. 综合管理者 -----------------选择：C 2. （）是指利用非推理因素采取迂回探索的办法来激发专家创造力的特殊创新会议。 A. 头脑风暴法 B. 名义群体法 C. 德尔菲法 D. 畅谈会法 -----------------选择：B 3. （）是指综观全局、认清为什么要做某事的能力，也就是洞察企业与环境要素间相互影响和相互作用的复杂性的能力。 A. 概念技能 B. 人际技能 C. 技术技能 D. 业务技能 -----------------选择：A 4. 按计划的广度分类，可分为（）。 A. 战略计划和战术计划 B. 长期计划和短期计划 C. 综合计划和专项计划 D. 指导性计划和具体性计划 -----------------选择：A 5. “把工人看成是会说话的机器，在体力和技能上受到最大限度的压榨”的是（）。 A. 泰勒 B. 甘特 C.

D. 吉尔布雷斯夫妇 E. 法约尔 -----------------选择：A 6. 被称为管理过程之父、现代经营管理之父的是（）。 A. 泰勒 B. 法约尔 C. 韦伯 D. 甘特 -----------------选择：B 7. 负责将战略管理者制定的总体目标和计划转化为更具体的目标和任务的管理者是（）。 A. 高层管理者 B. 中层管理者 C. 基层管理者 D. 综合管理者 -----------------选择：B 8. 提出“管理就是设计和保持一种良好环境，使人在群体里高效率地完成既定目标”的是（）。 A. 法约尔 B. 西蒙 C. 泰勒 D. 孔茨 -----------------选择：D 9. 下列哪种理论认为组织和组织成员的行为是复杂的，加上环境的复杂性和不断地变化，使得普遍适用的有效管理方法实质上是不可能存在的？（） A. 经验管理理论 B. 决策理论 C. 系统管理理论 D. 权变管理理论 -----------------选择：D 10. 在波士顿矩阵中，具有低增长、高市场份额特点的是（）。 A. 金牛 B. 明星 C. 问号 D. 瘦狗-----------------选择：A 大工15春《管

大连理工大学优化方法上机大作业

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)

大工15春《电工学》在线作业1满分答案

大工15春《电工学》在线作业1 试卷得分：100 一,单选题 1. 已知某感性负载的阻抗|Z|=7.07Ω，R=5Ω，则其功率因数是（）. A. 0.5 B. 0.6 C. 0.707 D. 0.8 ? 正确答案：C 2. RL电路发生换路时，（）不能跃变. A. uL B. uR C. iL D. IR ? 正确答案：C 3. RL电路的时间常数τ越大，则瞬态（）. A. 越慢 B. 不定 C. 越快 D. 不变 ? 正确答案：A 4. 下列不是变压器绕组参数的是（）. A. 电阻 B. 漏电抗 C. 漏容抗 D. 漏阻抗 ? 正确答案：C 5. 变压器的功率损耗等于（）. A. 磁滞损耗+涡流损耗 B. 铁损耗+铜损耗 C. 铁损耗+涡流损耗 D. 磁滞损耗+铜损耗 ? 正确答案：B

6. R=10Ω，XC=20Ω，XL=10Ω，三者并联后接于220V的交流电源上，电路的总电流为（） A. A. 24.6∠-26.57° B. 24.6 C. 55 D. 56 ? 正确答案：A 7. 某三相电动机，每相绕组的额定电压为220V，对称三相电源的线电压为380V，则三相绕组应作（）. A. 星形连接 B. 星形连接并接中线 C. 三角形连接 D. 星形、三角形连接均可 ? 正确答案：A 8. 一空载变压器，其原绕组电阻为22Ω，当原边加上额定电压220V时，原绕组中电流为（）. A. 10A B. 远大于10A C. 远小于10A D. 0 ? 正确答案：C 9. 有一RLC串联电路，已知R=XL=XC=5Ω，端电压U=10V，则I=（）A. A. 2/3 B. 1/2 C. 2 D. 1 ? 正确答案：C 10. 变压器带负载运行，已知变比N1/N2=10，当一次电流I1=10A时，二次电流I2=（）A. A. 100 B. 1 C. 10 D. 1000 ? 正确答案：A 二,判断题

大工15春在线作业：中国古建筑史100分

一、判断题（共 15 道试题，共 30 分。） V 1. 由于唐朝的黄帝姓李，因此尊崇老子，并在大明宫中建造了道观——三清殿。 A. 错误 B. 正确 满分：2 分 2. 白马寺中保存有宋、金、元、明、清等各个朝代的建筑和文物。 A. 错误 B. 正确 满分：2 分 3. 在古代，只有天子才有资格祭天。 A. 错误 B. 正确 满分：2 分 4. 西域沙漠上的汉长城是用红柳和夯土混合筑成的。 A. 错误 B. 正确 满分：2 分 5. 北京的天、地、日、月四座坛分别位于北京城的东、南、西、北方向上。 A. 错误 B. 正确 满分：2 分 6. 故宫中的御花园规模巨大，仅稍小于大明宫太液池。 A. 错误 B. 正确 满分：2 分 7. 故宫太和殿中屋顶的造型、台阶的层数、面宽和进深的间数以及所有的陈列和摆设都是有象征意义的。 A. 错误 B. 正确 满分：2 分 8. 石结构和木结构各有千秋，并无优劣之分。 A. 错误

B. 正确 满分：2 分 9. 建筑的物质功能取决于需求。 A. 错误 B. 正确 满分：2 分 10. 叠涩的功能与斗栱类似，但斗栱用于木结构，叠涩用于石结构。 A. 错误 B. 正确 满分：2 分 11. 由于木结构不容易保存，古人攻破都城时为了彰显胜利往往也会烧毁敌方的宫殿，因此至今完好保留下来的宫 殿建筑很少。 A. 错误 B. 正确 满分：2 分 12. 包拯在开封去世后就地安葬，合肥包拯墓是后人为纪念他而专门修建的衣冠冢。 A. 错误 B. 正确 满分：2 分 13. 孔庙、孔府和孔林的位置是相邻的。 A. 错误 B. 正确 满分：2 分 14. 因为中国盛产木材，西方盛产石材，所以中国古建筑多为木结构，西方古建筑多为石结构。 A. 错误 B. 正确 满分：2 分 15. 明朝中后期，朝廷决定在北京城外再修一圈新城墙，但由于花费巨大，不久就无力支付，只好改变计划，将南 面已修好的部分向内部弯折，与原来南面的城墙相接，这样形成了“凸”字形。 A. 错误 B. 正确

大连理工大学矩阵与数值分析上机作业

矩阵与数值分析上机作业 学校：大连理工大学 学院： 班级： 姓名： 学号： 授课老师：

注：编程语言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); %第一种原函数

大连理工大学矩阵大作业

2013级工科硕士研究生 《矩阵与数值分析》课程数值实验报告 大连理工大学 Dalian University of Technology

一、设 6 2 2 10 1 N N j S j = = - ∑，分别编制从小到大和从大到小的顺序程序分别计算 100001000000 , S S 并指出两种方法计算结果的有效位数。 程序代码： 从小到大： function f=s(N); %定义函数s f=0; %初始值为0 for j=N:-1:3 %j从3到n循环（从小到大） ft=1000000/(j^2-1); %Sj f=f+ft; %SN end 从大到小： function f=s(N); %定义函数s f=0; %初始值为0 for j=N:-1:3 %j从3到n循环（从小到大） ft=1000000/(j^2-1); %Sj f=f+ft; %SN end 执行结果： 从小到大： s(10000) ans = 4.16566671666167e+05 s(1000000) ans =

4.166656666671731e+05 有效数字：16,16 从大到小： s(10000) ans = 4.165666716661668e+05 s(1000000) ans = 4.166656666671667e+05 有效数字：16,16 分析： 小数和大数相加时，按照从大到小的顺序和按照从小到大的顺序得出的结果不同，前者由 于舍入误差的影响而使结果不准确，所以应避免大数吃小数的现象。 二、解线性方程组 1．分别利用Jacobi 迭代法和Gauss-Seidel 迭代法求解线性方程组Ax b =，其中常向量为()21n -维随机生成的列向量，系数矩阵A 具有如下形式 1111 11 1122n n n n n n n n T I I I A I I T I --------+-?? ?- ?= ? - ? -+? ? ， 其中1 211112n T --?? ? - ?= ?- ? -? ? 为1n -阶矩阵，1n I -为1n -阶单位矩阵，迭代法计算停止的条件为：10 12 10k k x x -+-<，给出10,100,1000n =时的不同迭代步数． 程序代码：

大工15春《工程估价》在线作业1满分答案

大工15春《工程估价》在线作业1 单选题多选题 一、单选题（共 15 道试题，共 60 分。） 1. 根据《建筑安装工程费用项目组成》（建标[2003]206号）文件规定，下列属于规费的是（）。 A. 环境保护费 B. 工程排污费 C. 安全施工费 D. 文明施工费 -----------------选择：B 2. 我国工程造价管理中，一般认为（）影响造价的程度最大。 A. 初步设计阶段 B. 技术设计阶段 C. 施工图设计阶段 D. 施工阶段 -----------------选择：A 3. 以下各项中，属于建筑安装工程直接工程费的是（）。 A. 模板及支架费 B. 大型机械设备进出场及安拆费 C. 二次搬运费 D. 人工费 -----------------选择：D 4. 按我国现行规定，国产设备的原价通常指（）。 A. 设备运至工地仓库的价格 B. 不带有备件的价格 C. 设备制造厂的交货价 D. 运至工地仓库并且安装完毕后的价格 -----------------选择：C 5. 根据“建筑安装工程费用项目组成”(建标[2003]206号)文件的规定，大型机械设备进出场及安拆费应计入（）。 A. 直接工程费 B. 措施费 C. 间接费 D. 施工机械使用费 -----------------选择：B 6. 如果以一所学校为一个建设项目，那么（）可以理解为单项工程。 A. 土建工程 B. 教学楼工程 C. 土石方工程 D. 暖通工程 -----------------选择：B

7. 工程估价的计价过程具有多次性的特征，其中（）应进行的为预算造价。 A. 施工图设计阶段 B. 招投标阶段 C. 合同实施阶段 D. 初步设计阶段 -----------------选择：A 8. 按我国现行规定，下列费用中属于企业管理费的是（）。 A. 劳动保险费 B. 医疗保险费 C. 住房公积金 D. 养老保险费 -----------------选择：A 9. 世界银行工程造价的构成中，（）包括未明确项目的准备金和不可预见准备金两个部分。 A. 项目直接建设成本 B. 项目间接建设成本 C. 应急费 D. 建设成本上升费用 -----------------选择：C 10. 按我国现行规定，以下各项（）不属于直接工程费。 A. 人工费 B. 施工机械使用费 C. 材料费 D. 二次搬运费 -----------------选择：D 11. 根据《建筑安装工程费用项目组成》（建标[2003]206号）文件规定，对构件和建筑安装物进行一般鉴定和检查所发生的费用列入（）。 A. 材料费 B. 规费 C. 措施费 D. 研究试验费 -----------------选择：A 12. 建设项目的工程造价在量上应与（）相等。 A. 建设项目总投资 B. 固定资产投资 C. 建筑安装工程投资 D. 静态投资 -----------------选择：B 13. 按我国现行规定，基本预备费的取费基数为（）。 A. 工程费用 B. 工程建设其他费用 C. 工程费用+工程建设其他费用 D. 安装费用 -----------------选择：C

大连理工c语言第一次上机作业参考答案

第一次上机作业参考答案： 1．大写字母转换成小写字母 从键盘输入一个大写英文字母，输出相应的小写字母。 例：输入G 输出g #include void main() { char c; c=getchar(); if(c>='A' && c<='Z') c+=32; putchar(c); } 2．求平方根 输入1 个实数x，计算并输出其平方根（保留1 位小数）。 例：输入17 输出The square root of 17.0 is 4.1 #include #include void main() { float x,root; scanf("%f",&x); if(x>0) root=sqrt(x); else printf("Input Error!\n"); printf("The square root of %.1f is %.1f\n",x,root); } 3．Temperature Conversion Design a program which converts from degrees Fahrenheit temperature to degrees Celsius temperature. c = 5/9(f-32) a)Input Fahrenheit temperature will be type float. b)Display the temperatures with 2 places of precision #include void main() { float f,c; scanf("%f",&f); c = 5.0/9*(f-32); printf("Fahrenheit %.2f is equal to Celsius %.2f\n",f,c); } 4. 计算旅途时间 输入2 个整数time1 和time2，表示火车的出发时间和到达时间，计算并输出旅途时间。 (有效的时间范围是0000 到2359，不需要考虑出发时间晚于到达时间的情况。) 例：输入712 1411 （出发时间是7：12，到达时间是14：11） 输出The train journey time is 6 hrs 59 mins. #include void main( ) { int time1, time2, hours, mins; scanf("%d%d", &time1, &time2);

大连理工大学机械设计大作业

目录 一、设计任务书及原始数据 (2) 二、根据已知条件计算传动件的作用力 (3) 2.1计算齿轮处转矩T、圆周力F t 、径向力F r及轴向力F a .. 3 2.2计算链轮作用在轴上的压力 (3) 2.3计算支座反力 (4) 三、初选轴的材料，确定材料的机械性能 (4) 四、进行轴的结构设计 (5) 4.1确定最小直径 (5) 4.2设计其余各轴段的直径和长度，且初选轴承型号 (5) 4.3选择连接形式与设计细部结构 (6) 五．轴的疲劳强度校核 (6) 5.1轴的受力图 (6) 5.2绘制弯矩图 (7) 5.3绘制转矩图 (8) 5.4确定危险截面 (9) 5.5计算当量应力，校核轴的疲劳强度 (9) 六、选择轴承型号，计算轴承寿命 (10)

6.1计算轴承所受支反力 (10) 6.2计算轴承寿命 (11) 七、键连接的计算 (11) 八、轴系部件的结构装配图 (12) 一、设计任务书及原始数据 题目二：二级展开式斜齿圆柱齿轮减速器输出轴组合结构设计

表1 设计方案及原始数据 二、根据已知条件计算传动件的作用力 2.1计算齿轮处转矩T、圆周力F t、径向力F r及轴向力F a 已知：轴输入功率P=4.3kW，转速n=130r/(min)。 (1)齿轮上的力 转矩计算公式：T=9.550×106P/n 将数据代入公式中，得：T=315885(N·mm) 圆周力计算公式： F t =2T/，==416(mm) (认为是法面模数) 将转矩T带入其中，得：F t =1519(N) 径向力计算公式：F r =F t ×tanα/cos，= 将圆周力F t 带入其中，得：F r =558(N) 轴向力计算公式：F a = F t ×tan 将圆周力F t 带入其中，得：F a =216(N) 2.2计算链轮作用在轴上的压力 链轮的分度园直径 链速v= 链的圆周力F= 链轮作用在轴上的压力

大连理工大学优化方法上机作业

大连理工大学优化方法上机 作业 -标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

优化方法上机大作业 学院：电子信息与电气工程学部 姓名： 学号： 指导老师：

上机大作业（一） %目标函数 function f=fun(x) f=100*(x(2)-x(1)^2)^2+(1-x(1))^2; end %目标函数梯度 function gf=gfun(x) gf=[-400*x(1)*(x(2)-x(1)^2)-2*(1-x(1));200*(x(2)-x(1)^2)]; End %目标函数Hess矩阵 function He=Hess(x) He=[1200*x(1)^2-400*x(2)+2,-400*x(1); -400*x(1), 200;]; end %线搜索步长 function mk=armijo(xk,dk) beta=0.5; sigma=0.2; m=0; maxm=20; while (m<=maxm) if(fun(xk+beta^m*dk)<=fun(xk)+sigma*beta^m*gfun(xk)'*dk) mk=m; break; end m=m+1; end alpha=beta^mk newxk=xk+alpha*dk fk=fun(xk) newfk=fun(newxk) %最速下降法 function [k,x,val]=grad(fun,gfun,x0,epsilon) %功能：梯度法求解无约束优化问题：minf(x) %输入：fun，gfun分别是目标函数及其梯度，x0是初始点， % epsilon为容许误差 %输出：k是迭代次数，x，val分别是近似最优点和最优值 maxk=5000; %最大迭代次数