大连理工优化方法大作业MATLAB编程

fun ctio n [x,dk,k]=fjqx(x,s) flag=0;

a=0;

b=0;

k=0;

d=1;

while (flag==0)

[p,q]=getpq(x,d,s);

if (P<0)

b=d;

d=(d+a)/2;

end

if(p>=0) &&( q>=0)

dk=d;

x=x+d*s;

flag=1;

end

k=k+1;

MATLAB编程作业

《Matlab 编程训练》作业专业学生姓名班级学号指导教师完成日期

实训一 MATLAB 语言介绍和数值计算 1．先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存变量。 12 2sin851z e =+ . 2. 已知 1234413134787,2033657327A B --???? ????==???? ????-???? ，求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） A+6*B:

A-B+I: (2)A*B和A.*B A*B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] c=A*B 结果： A.*B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=A.*B 结果：

(3)A^3和A.^3 A^3程序： A=[12 34 -4;34 7 87;3 65 7] E=A^3 结果： A．^3程序： A=[12 34 -4;34 7 87;3 65 7] C=A.^3 (4)A/B及B\A A/B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] C=A/B 结果：

B\A程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=B\A 结果： (5)将矩阵C=B\A的右下角2*2子矩阵赋给D, 并（3）保存变量（mat文件）程序： A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; C=B*inv(A); D=C(2:3,2:3) 结果：

MATLAB程序设计作业

Matlab程序设计班级姓名学号

《MATLAB程序设计》作业 1、考虑如下x-y 一组实验数据： x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot的原始数据、一次拟合曲线和三次拟合曲线，给出MATLAB代码和运行结果。代码如下： x=[1,2,3,4,5,6,7,8,9,10]; y=[1.2,3,4,4,5,4.7,5,5.2,6,7.2]; plot(x,y); title('原始数据'); p=polyfit(x,y,1); q=polyval(p,x); figure,plot(x,q); title('一次拟合'); p=polyfit(x,y,2); q=polyval(p,x); figure,plot(x,q); title('二次拟合'); 运行结果如下：

1 2 3 4 5 6 7 8 9 10 12 3 4 5 6 7 8 原始数据 123 456789 102 2.5 3 3.54 4.5 55.56 6.57一次拟合 123456789 101 2 3 4 5 6 7 二次拟合 2、在[0，3π]区间，绘制y=sin(x)曲线（要求消去负半波，即(π，2π)区间内的函数值置零），求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。代码如下： clear clc x=(0:0.01:3*pi); y=sin(x); plot(x,y); y1=(y>=0).*y; figure,plot(x,y1);

实例matlab-非线性规划-作业

现代设计方法-工程优化理论、方法与设计姓名学号班级研问题：某厂向用户提供发动机，合同规定，第一、二、三季度末分别交货40台、60台、80台。每季度的生产费用为（元），其中x 是该季生产的台数。若交货后有剩余，可用于下季度交货，但需支付存储费，每台每季度c 元。已知工厂每季度最大生产能力为100台，第一季度开始时无存货，设a=50、b=0.2、c=4，问工厂应如何安排生产计划，才能既满足合同又使总费用最低。讨论a 、b 、c 变化对计划的影响，并作出合理的解释。问题的分析和假设：问题分析：本题是一个有约束条件的二次规划问题。决策变量是工厂每季度生产的台数，目标函数是总费用（包括生产费用和存储费）。约束条件是生产合同，生产能力的限制。在这些条件下需要如何安排生产计划，才能既满足合同又使总费用最低。问题假设： 1、工厂最大生产能力不会发生变化； 2、合同不会发生变更； 3、第一季度开始时工厂无存货； 4、生产总量达到180台时，不在进行生产； 5、工厂生产处的发动机质量有保证，不考虑退货等因素； 6、不考虑产品运输费用是否有厂家承担等和生产无关的因素。符号规定： x1——第一季度生产的台数； x2——第二季度生产的台数； 180-x1-x2——第三季度生产的台数； y1——第一季度总费用； y2——第二季度总费用； y3——第三季度总费用； y ——总费用（包括生产费用和存储费）。 ()2bx ax x f +=

建模： 1、第一、二、三季度末分别交货40台、60台、80台； 2、每季度的生产费用为（元）； 3、每季度生产数量满足40 ≤x1≤100，0≤x2≤100，100≤x1+x2 ≤180； 4、要求总费用最低，这是一个目标规划模型。目标函数: y1 2111x b x a Z ?+?= y2()4012222-?+?+?=x c x b x a Z y3()()()10018018021221213 -+?+--?+--?=x x c x x b x x a Z y x x x x x x Z Z Z Z 68644.04.04.0149201 212221321--+++=++= 40≤x1≤100 0≤x2≤100 100≤x1+x2≤180 ()2 bx ax x f +=

大连理工大学网络教育学院《管理学》课程大作业满分

网络教育学院《管理学》课程大作业学习中心：层次：专业：年级：学号：姓名：完成日期：

大工20春《管理学》大作业及要求第一部分：注意：请从以下题目中任选其一作答！题目一：谈谈如何正确理解管理既是一门科学又是一门艺术。在实践工作中如何运用这一基本原理？题目二：谈谈现代管理理论中具有代表性的管理理论学派的主要思想。题目三：不同层次的管理者在应具备的技能上有何侧重?请举例说明。题目四：试述影响集权与分权的因素。题目五：结合实际论述领导者应具备的用人艺术。题目三：不同层次的管理者在应具备的技能上有何侧重?请举例说明。答：管理者分为高层管理者、中层管理者、基础管理者。不同层次的管理者都应该具备技术技能、人际技能和概念技能。只是有不同的侧重点。技术技能：对于基层管理者最重要，对于中层管理者较重要，对于高层管理者不重要。人际技能：对于任何层次的管理者都重要。概念技能：对于高层管理者最重要，对于

中层管理者较重要，对于基层管理者不重要。比如一个房地产企业，高层管理者为总裁、副总裁、股东等等。他们制定和实施公司总体战略，完成董事会下达的年度经营目标，按照发展战略开展具体的经营工作，负责建设高效的组织团队等；中层管理者为项目经理，区域经理等，他们按照高层管理者战略要求，负责或协助基础管理者工作，发挥着承上启下作用。房产开发项目经理，房产销售经理等都要保证各个项目顺利进行，努力完成高层的要求。基层管理者为房地产开发包工头，销售主管主要负责管理他们的团队，让作业人员能顺利开展工作。本身要求自己要熟悉这块业务，才能给底下员工更多的指导与帮助。包工头对于建设商品房的每个环节都要很熟悉，能控制成本，及时完工。销售主管管理好销售团队，做好每天日常考勤、仪表、销售报表等。第二部分：学习心得通过管理学这个课程，我深刻地意识到一个企业的成功离不开每个管理者，而每个管理者必须具备相应的管理技能。认识了管理在企业中的重要性。一个好的管理者能让企业迅速发展，管理层制定的管理决策影响整个企业的未来。比如华为集团，他们凭什么在手机行业瑶瑶领先？不管是高层的决策，中层的实施，基层的管理都很到位。领先专业技术管理、狠抓业务，带好团队。一个不称职的管理者会让企业走向末路，比如10年前的“三鹿奶粉事件”，他们为了利益，不顾产品质量管理。作为管理者必须要加大对企业内管质量人员的教育力度，使他们认识到质量就是企业的生命，质量问题是企业最大的灭亡隐患。杜绝不合格的奶制品在商业腐败中流向市场。管理学同样与我们息息相关，管理是一切组织的根本，管理工作适用于各种大小规模的组织；盈利与非盈利的企事业单位、制造业以及服务性行业；因此，学好管理学对于我们现在的工作岗位都有其非常重要的意义。目前我们公司绩效管理和有效的激励机制很符合管理者的要求，我一定要学好管理学这个课程。

matlab程序设计第三章课后习题答案

1. p138 第6题在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点。 >> t=0:0.01:pi; >> x1=t; >> y1=2*x1-0.5; >> x2=sin(3*t).*cos(t); >> y2=sin(3*t).*sin(t); >> plot(x1,y1,'r-',x2,y2,'g-') >> axis([-1,2,-1.5,1]) >> hold on >> s=solve('y=2*x-0.5','x=sin(3*t)*cos(t)','y=sin(3*t)*sin(t)'); >> plot(double(s.x),double(s.y),'*'); 截图：

p366 第4题绘制极坐标曲线，并分析对曲线形状的影响。 function [ output_args ] = Untitled2( input_args ) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here theta=0:0.01:2*pi; a=input('请输入a的值：'); b=input('请输入b的值：'); n=input('请输入n的值：'); rho=a*sin(b+n*theta); polar(theta,rho,'k'); end 下面以a=1，b=1，n=1的极坐标图形为基础来分析a、b、n的影响。

对a的值进行改变：对比发现a只影响半径值的整倍变化对b的值进行改变：对比发现b的值使这个圆转换了一定的角度

对n的值进行改变：对比发现当n>=2时有如下规律 1、当n为整数时，图形变为2n个花瓣状的图形 2、当n为奇数时，图形变为n个花瓣状的图形分别让n为2、3、4、5

大连理工大学结构优化复习总结

结构优化设计－基于结构分析技术，在给定的设计空间实现满足使用要求且具有最佳性能或最低成本的工程结构设计的技术优化设计的三要素:设计变量;约束条件;目标函数凸域:基于n维空间的区域s里，如果取任意两点x1和x2，连接这两点的线段也属于s，该区域称凸域（=αx1+(1-α)x2 ）凸函数:如果函数f(x)定义在n维空间的凸域s上，而且对s中的任意两点x1和x2和任意常数α,0.0<=α<=1.0,有f[αx1+(1- α)x2]<=αf(x1)+(1- α)f(x2),则f(x)称为s上的凸函数严格凸函数:上式小于严格成立凸规划:如果可行域是凸域，目标函数是凸函数，这样构成的数学规划问题为凸规划问题。准则设计法:依靠工程经验;效率高;缺乏严格数学基础最优准则法基于库塔克（K-T）条件:需构造迭代求解算法;通用性不强数学规划方法:有严格的数学基础，有较好的通用性，计算效率要考虑。结构优化问题的求解布骤 I. 建立优化模型。给定初始设计方案。 II. 结构分析(有限元） III.优化（收敛性）检验。满足则结束程序，否则继续IV IV. 灵敏度分析 V. 求解优化问题，修改结构模型，返回II。优化求解的两大类方法：准则法；数学规划法准则设计方法：用优化准则代替原来的优化问题同步失效准则设计的评价： {优点：简单、方便，特别是独立约束个数n=m时；工程实用；适合于构件设计。缺点：只能处理简单构件设计；缩小了设计空间，不能保证最优解；若n < m ，可能无解；当n > m时，确定哪些破坏模式应同时发生比较困难。改进：为了弥补等式约束代替不等式约束的缺陷，引入松弛因子ψi σi (X ) =ψiσip , 0 ≤ψi ≤1, i =1,2,......n 启发：用准则代替原来的优化问题，准则法的基本思想；如果将桁架的每根杆看作一种可能的破坏模式，桁架看作一个元件。可以得到满应力准则满应力方法的缺点：完全无视重量会漏掉最轻设计；中间点一般是不可行设计，对工程实际不利。希望得到可行的中间设计点。齿形法：采用射线步进行可行性调整，适用于桁架一类刚度与设计变量成正比的结构。将所有设计变量同时乘以一个常数ξ：A n i=ξA i o} 线性函数都是凸函数，线性规划是凸规划。

BP神经网络matlab实例(简单而经典)

p=p1';t=t1'; [pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化 net=newff(minmax(pn),[5,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络 net.trainParam.show=2000; % 训练网络 net.trainParam.lr=0.01; net.trainParam.epochs=100000; net.trainParam.goal=1e-5; [net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络 pnew=pnew1'; pnewn=tramnmx(pnew,minp,maxp); anewn=sim(net,pnewn); %对BP网络进行仿真 anew=postmnmx(anewn,mint,maxt); %还原数据 y=anew'; 1、BP网络构建（1）生成BP网络 = net newff PR S S SNl TF TF TFNl BTF BLF PF (,[1 2...],{ 1 2...},,,) R?维矩阵。 PR：由R维的输入样本最小最大值构成的2 S S SNl：各层的神经元个数。 [ 1 2...] { 1 2...} TF TF TFNl：各层的神经元传递函数。 BTF：训练用函数的名称。（2）网络训练 [,,,,,] (,,,,,,) = net tr Y E Pf Af train net P T Pi Ai VV TV （3）网络仿真 = [,,,,] (,,,,) Y Pf Af E perf sim net P Pi Ai T {'tansig','purelin'},'trainrp'

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

大连理工大学《桥涵水文》大作业及要求

网络教学学院《桥涵水文》离线作业学习中间：层次：专业：年级：年春/秋季学号：学生：教导老师：杨颖完结日期：年月日大工20春《桥涵水文》大作业及要求留意：从以下五个题目中任选两个进行回答（留意：从题目一、二中挑选一道计算题，并从题目三、四、五中挑选一道问答题，别离进行回答，不可以一起挑选两道计算题或两道问答题）；回答前，需将所选的题目进行仿制（使教师清晰你所选的题目）。题目一：计算题某水文站有22年的年最大流量观测材料，并已知其统计参数，，计算、、和的抽样差错，以及榜首项的经历频率的抽样差错。（，B=2.9）题目二：计算题已知某海湾22年的年最高潮水位的观测材料如表1所示，要求计算重现期为100年和50年的最高潮位。（计算过程中如有需求可直接在表1中增加计算内容）表1 某海湾22年的年最高潮水位观测数据 m 年最高潮位（cm） m 年最高潮位（cm） 1 330 108900 13 275 75625 2 315 99225 14 274 75076 3 310 96100 15 270 72900

4 303 91809 16 268 71824 5 297 88209 17 264 69696 6 286 81706 18 263 69169 7 285 81225 19 261 68121 8 284 80656 20 256 65536 9 284 80656 21 254 64516 10 282 79524 22 254 64515 11 280 78400 6173 1740763 12 278 77284 题目三：桥梁工程师有必要知道桥位河段的水文和河槽演化特性,试述桥位河段如何分类。不一样河段应如何布设桥孔及墩台？（可罗列工程实例或简略算例，主张必要处附图像阐明。）答复内容应与题目要求相共同，字数不少于1000，字体宋体小四字，1.5倍行距离。题目四：查阅有关材料，论述桥下河槽的冲刷表象和过程如何？它们构成的缘由是啥？桥梁冲刷计算时选用了哪些处理方法？（主张必要处附图像阐明。）答复内容应与题目要求相共同，字数不少于1000，字体宋体小四字，1.5倍行距离。题目五：请谈谈当前公路和桥梁工程面对哪些新方式？呈现哪些新理念？（可罗列工程实例，主张必要处附图像阐明。）答复内容应与题目要求相共同，字数不少于1000，字体宋体小四字，1.5倍行距离。作业详细要求： 1．以附件方式上交离线作业（附件的巨细约束在10M以内），挑选已完结的作业（留意命名），点提交即可。如下图所示。 2．封面格局封面称号：大连理工大学桥涵水文大作业，字体为宋体加黑，字号为小一；名字、奥鹏卡号、学习中间等字体为宋体，字号为小三号。 3．正文格局作业正文内容一致选用宋体，字号为小四。留意：作业大概独立完结，禁绝抄袭其他网站或许请人代做，如有相同作业，分数以零分计。

实验二--MATLAB程序的设计(含实验报告)

实验二 MATLAB 程序设计一、实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。二、实验的设备及条件计算机一台（带有MATLAB7.0以上的软件环境）。 M 文件的编写：启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正三、实验容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。（提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。要求：（1）用switch 语句实现。（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。（提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2→1 3→10→5→16→8→4→2→1 6→3→10→5→16→8→4→2→1

大连理工大学大作业

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

题目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,3…n) ∑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 >+-

Matlab编程与应用习题和一些参考答案

Matlab 上机实验一、二 3.求下列联立方程的解???????=+-+-=-+=++-=--+4 1025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; >> b=[4;4;9;4]; >> c=a\b 4．设???? ??????------=81272956313841A ，??????????-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。 >> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8]; >> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; >> C1=A*B' >> C2=A'*B >> C3=A.*B >> inv(C1) >> inv(C2) >> inv(C3) 5．设 ?? ????++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。 >> x=linspace(0,2*pi,101); >> y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); >> plot(x,y,'r') 6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。（mean var ） a=randn(8,6) mean(a) var(a) k=mean(a) k1=mean(k) i=ones(8,6) i1=i*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g)

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

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

(完整版)大连理工大学《高层建筑结构》大作业.doc

大连理工大学《高层建筑结构》大作业学习中心：姓名：学号：题目二：底部剪力法计算题钢筋混凝土 4 层框架经质量集中后计算简图如下图所示，各层高均为4m，集中于各楼层的重力荷载代表值分别为： G1 435kN ，G2 440kN ，G3 430kN ，G4380kN 。结构阻尼比0.05 ，自振周期为 T10.383s ，Ⅰ1类场地类别，设计地震分组为第一组，抗震设防烈度为8 度（设计基本地震加速度为0.30g）。按底部剪力法计算结构在多遇地震时的水平地震作用及地震剪力。 G4380kN G3430kN G2440kN G1435kN （ a）计算简图解：（1）计算结构等效总重力荷载代表值 G eq0.85 G 0.85 G1G2G3G4 0.85 380 430 440435

1432.25kN （2）计算结构总水平地震作用 F EK1 G eq0.139 5997.6833.7kN （3）计算顶部附加水平地震作用 1.4T g 1.4 0.40 5.6s T10.467s 故，不需考虑顶部附加水平地震作用，即n0 。（4）计算各楼层水平地震作用（如下图所示） G i H i 1 n F EK 分别计算各楼层水平地震作用，如根据公式 F in G j H j j 1 下： F1 270 9.8 3.5 833.7 166.7kN 270 9.8 3.5 270 9.8 7 180 9.8 10.5 F2 270 9.8 7.0 833.7 333.5kN 270 9.8 3.5 270 9.8 7 180 9.8 10.5 F3 270 9.8 10.5 833.7 333.5kN 270 9.8 3.5 270 9.8 7 180 9.8 10.5 （5）计算各楼层的层间地震剪力（如下图所示） V1 F1 F2 F3 833.7kN V2F2F3667.0kN V3F3333.5kN

matlab程序设计作业

Matlab程序设计作业姓名: 学号: 专业:

? MATLAB 程序设计》作业 1、考虑如下x-y 一组实验数据: x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y 二[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot 的原始数据、一次拟合曲线和三次拟合曲线，给出原始曲线 MATLAB 代码和运行结果。 7 6 5 4 3 2 2 3 4 5 6 7 8 9 10

7 6.5 6 5.5 5 4.5 4 3.5 3 2.5 10 一次拟合三次拟合

x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]; y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2]; figure; plot(x,y) p1=polyfit(x,y,1); y1=polyval(p1,x); figure; plot(x,y1) p2=polyfit(x,y,3); y2=polyval(p2,x); figure; plot(x,y2) 2、在［0, 3n区间，绘制y二Sin(x)曲线(要求消去负半波，即(n 2n)区间内的函数值置零)，求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 x=0:pi/1000:3*pi; y=Sin(x); y1=(y>=0).*y; %消去负半波figure(1); plot(x,y1, 'b' ); a=mean(y1) %求出y1 的平均值 b=max(y1) %求出y1 的最大值b, 以及最大值在矩阵中的位置; d=x(find(y1==b)) >> ex1 a = 0.4243 b = 1 d = 1.5708 7.8540 >>

matlab函数计算的一些简单例子1

MATLAB作业一1、试求出如下极限。（1） 23 25 (2)(3) lim (5) x x x x x x x ++ + →∞ ++ + ，（2） 23 3 1 2 lim () x y x y xy x y →- → + + ，（3） 22 22 22 1cos() lim ()x y x y x y x y e+ → → -+ + 解：（1）syms x； f=((x+2)^(x+2))*((x+3)^(x+3))/((x+5)^(2*x+5)) limit（f，x，inf） =exp（-5）（2）syms x y； f=(x^2*y+x*y^3)/(x+y)^3； limit(limit(f,x,-1),y,2) =-6；（3）syms x y； f=(1-cos(x^2+y^2))/(x^2+y^2)*exp(x^2+y^2)； limit(limit(f,x,0),y,0) =0 2、试求出下面函数的导数。（1 ）() y x=, (2)22 atan ln() y x y x =+ 解; （1）syms x; f=sqrt(x*sin(x)*sqrt(1-exp(x))); g= diff(f,x); g== (sin(x)*(1 - exp(x))^(1/2) + x*cos(x)*(1 - exp(x))^(1/2) - (x*exp(x)*sin(x))/(2*(1 - exp(x))^(1/2)))/(2*(x*sin(x)*(1 - exp(x))^(1/2))^(1/2)) pretty(g)= (2)syms x y; f=atan(y/x)-log(x^2+y^2) pretty(-simple(diff(f,x)/diff(f,y)))= 2 x + y =------- x - 2 y (3) 假设1 cos u- =，试验证 22 u u x y y x ?? = ???? 。解：syms x y; u=1/cos(sqrt(x/y)); diff(diff(u,x),y)-diff(diff(u,y),x)=0; 所以： 22 u u x y y x ?? = ????

大连理工优化方法大作业MATLAB编程

function [x,dk,k]=fjqx(x,s) flag=0; a=0; b=0; k=0; d=1; while(flag==0) [p,q]=getpq(x,d,s); if (p<0) b=d; d=(d+a)/2; end if(p>=0)&&(q>=0) dk=d; x=x+d*s; flag=1; end k=k+1;

if(p>=0)&&(q<0) a=d; d=min{2*d,(d+b)/2}; end end %定义求函数值的函数fun，当输入为x0=（x1，x2）时，输出为f function f=fun(x) f=(x(2)-x(1)^2)^2+(1-x(1))^2; function gf=gfun(x) gf=[-4*x(1)*(x(2)-x(1)^2)+2*(x(1)-1),2*(x(2)-x(1)^2)]; function [p,q]=getpq(x,d,s) p=fun(x)-fun(x+d*s)+0.20*d*gfun(x)*s'; q=gfun(x+d*s)*s'-0.60*gfun(x)*s'; 结果： x=[0,1]; s=[-1,1]; [x,dk,k]=fjqx(x,s) x =-0.0000 1.0000 dk =1.1102e-016 k =54

function f= fun( X ) %所求问题目标函数 f=X(1)^2-2*X(1)*X(2)+2*X(2)^2+X(3)^2+ X(4)^2- X(2)*X(3)+2*X(1)+3*X(2)-X(3); end function g= gfun( X ) %所求问题目标函数梯度 g=[2*X(1)-2*X(2)+2,-2*X(1)+4*X(2)-X(3)+3,2*X(3)-X(2)-1,2*X(4)]; end function [ x,val,k ] = frcg( fun,gfun,x0 ) %功能：用FR共轭梯度法求无约束问题最小值 %输入：x0是初始点，fun和gfun分别是目标函数和梯度 %输出：x、val分别是最优点和最优值，k是迭代次数 maxk=5000;%最大迭代次数 rho=0.5;sigma=0.4;

matlab仿真实例

matlab 仿真实例实验五MATLAB 及仿真实验一、控制系统的时域分析 (一)稳定性 1、系统传递函数为G(s)，试判断其稳定性。程序： >> nu m=[3,2,5,4,6]; >> den=[1,3,4,2,7,2]; >> sys=tf( nu m,de n); >> figure(1); >> pzmap(sys); >> title(' 零极点图') 由图可知：在S 右半平面有极点，因此可知系统是不稳定的。 2、用MATLA 求出 G(s)=(s A 2+2*s+2)/(s A 4+7*s A 3+5*s+2) 的极点。程序及结果： >> sys=tf([1,2,2],[1,7,3,5,2]); >> p=pole(sys) 矿'. 赳 _ ■ —

-6.6553 0.0327 + 0.8555i 0.0327 - 0.8555i -0.4100 (二)阶跃响应 1、二阶系统G(s)=10/s A2+2*s+10 1)键入程序，观察并记录单位阶跃响应曲线：程序： >> sys=tf(10,[1,2,10]); >> step(sys); >> title('G(s)=10/sA2+2*s+10 单位阶跃响应曲线') 2）计算系统闭环跟、阻尼比、无阻尼振荡频率，并记录程序及结果： >> sys=tf(10,[1,2,10]); >> p=pole(sys)

p = -1.0000 + 3.0000i -1.0000 - 3.0000i >> [wn,z]=damp(sys) wn = 3.1623 3.1623 z = 0.3162 0.3162 3)记录实际测取的峰值大小，峰值时间和过渡过程时间，并填表实际值理论值峰值Cmax 1.35s 峰值时间tp 1.05s 过渡时间+5% 3.54s ts +2% 3.18s 程序: >> sys=tf(10,[1,2,10]); >> step(sys); >> title('G(s)=10/sA2+2*s+10 单位阶跃响应曲线')

大连理工大学关于2012-2016级本科生培养方案

大连理工大学关于制订2012-2016级本科生培养方案的指导性意见培养方案是人才培养的基本“蓝图”，是学校实现人才培养目标、保证人才培养质量的基础性文件，也是学生进行学习规划、学校进行教学组织与管理、学部和院系进行课程建设规划的主要依据。 2008年，学校召开了第十四次本科教育教学研讨会，确立了实施精英教育，培养精英人才的目标定位。为了进一步加强本科专业内涵建设，探索特色鲜明的人才培养途径，培养高素质拔尖创新人才，将我校近年来在教育教学改革，特别是在大类招生、大类培养、创新教育、研究型教学等方面取得的成效固化到培养方案中，理顺通识教育与专业教育、人文教育与科学教育、理论教学与实践教学、共性要求与个性培养等关系，优化本科专业的知识结构和课程体系，学校决定制订2012-2016级本科大类和专业培养方案，特提出以下指导性意见。一、指导思想新一轮培养方案的制订要坚持党的教育方针，遵循高等教育的发展规律，充分借鉴国内外一流大学的人才培养经验，以精英人才培养为目标，以促进学生全面发展和适应社会发展需求为基本定位，注重通识教育，确立适合理工科学生和人文与经管类学生选修的两类通识教育核心课程；注重大类培养，构建公共基础和学科基础平台课程；注重个性发展，灵活设置特色专业方向和各类选修课程。以学生实践能力与创新能力培养为导向，注重实践与创新教育，加强创新性实践环节；注重科研与教学紧密结合，将最新科研成果融入教学之中，设置学科前沿课程、学科体验实验等。以优化知识结构为重点，注重本研贯通，设置本科生和研究生层次递进的课程和通选互认课程，科学构建与研究型大学精英人才培养相适应的课程体系，最终形成以能力培养为导向，体现知识、能力、素质协调发展的人才培养方案。二、培养方案构成本科专业培养方案和指导性教学计划的主要内容包括：（1）类别或专业名称 5 1

matlab函数计算的一些简单例子2

MATLAB 作业二 1、请将下面给出的矩阵A 和B 输入到MATLAB 环境中，并将它们转换成符号矩阵。若某一矩阵为数值矩阵，另以矩阵为符号矩阵，两矩阵相乘是符号矩阵还是数值矩阵。 57651653 5501232310014325462564206441211346,3 9636623 51521210760077410120172440773 473 78 124867217110 7 681 5A B ???? ?????????? ????? ?==?????? ????? ?---????????--??? ? 解：A 转换为符号矩阵；a=sym（A） a=[5,7,6,5,1,6,5] [2,3,1,0,0,1,4][6,4,2,0,6,4,4][3,9,6,3,6,6,2][10,7,6,0,0,7,7][7,2,4,4,0,7,0][4,8,6,7,2,1,7]B 转换为符号矩阵；b=sym（B）b = [3,5,5,0,1,2,3][3,2,5,4,6,2,5][1,2,1,1,3,4,6][3,5,1,5,2,1,2][4,1,0,1,2,0,1][-3,-4,-7,3,7,8,12][1,-10,7,-6,8,1,5] 若某一矩阵为数值矩阵，另以矩阵为符号矩阵，两矩阵相乘是符号矩阵例；a*B= [48,3,64,48,159,106,194][17,-26,47,-8,62,26,59][48,-8,52,12,108,64,124][59,22,41,69,151,101,184][43,-22,91,13,175,121,220][22,39,4,53,88,94,147][75,11,115,36,151,70,151] 2、利用MATLAB 语言提供的现成函数对习题1中给出的两个矩阵进行分析，判定它们是否为奇异矩阵，得出矩阵的秩、行列式、迹和逆矩阵，检验得出的逆矩阵是否正确。解：由于a=det（A）=3.7396e+04；故A 是非奇异矩阵。B=det（B）=0，故B 是奇异矩阵；由于a=rank(A)=7，故A 的秩为7；由于b=rank(B)=5，故B 的秩为5；由于a=trace(A)=27，b=trace(B)=26，故A，B 的迹为27,26；由a=inv(A)得A 的逆矩阵如下；