结构动力学大作业 对于如下结构，是研究质量块的质量变化和在简支梁上位置的变化对整个系统模态的影响。 1 以上为一个简支梁结构。集中质量块放于梁上，质量块距简支梁的左端点距离为L. 将该简支梁简化为欧拉伯努利梁，并离散为N 个单元。每个单元有两个节点，四个自由度。 单元的节点位移可表示为： ]1122,,,e v v δθθ?=? 则单元内一点的挠度可计作： 带入边界条件： 1 3 32210)(x a x a x a a x v +++=0 1)0(a v x v ===3 322102)(L a L a L a a v L x v +++===1 10 d d a x v x ===θ2 321232d d L a L a a x v L x ++===θ1 0v a =
[]12 3 4N N N N N = 建立了单元位移模式后，其动能势能均可用节点位移表示。单元的动能为： 00111()222 l l T T T ke e e e e y E dx q N Ndxq q mq t ρρ?===??? 其中m 为单元质量阵，并有： l T m N Ndx ρ=? 带入公式后积分可得： 222215622541322413354 1315622420133224l l l l l l l m l l l l l l ρ-?? ??-??= ?? -?? ---? ? 单元势能可表示为 22 200 11()()22 2 T l l T T e pe e e e q y E EI dx EI N N dxq q Kq x ?''''== =??? 其中K 为单元刚度矩阵，并有 ()l T K EI N N dx ''''=? 2 23 2212 612664621261266264l l l l l l EI k l l l l l l l -????-??=??---??-?? 以上为单元类型矩阵，通过定义全局位移矩阵，可以得到系统刚度矩阵和系统质量矩 1 1θ=a )2(1)(3211222θθ+--=L v v L a )(1)(22122133θθ++-= L v v L a 1232133222231)(θ???? ??+-+???? ??+-=L x L x x v L x L x x v 2 2232332223θ??? ? ??-+???? ??-+L x L x v L x L x 2 4231211)()()()()(θθx N v x N x N v x N x v +++=
结构动力学作业 姓名： 学号：
目录 1.力插值法 (1) 1.1分段常数插值法 (1) 1.2分段线性插值法 (4) 2.加速度插值法 (7) 2.1常加速度法 (7) 2.2线加速度法 (9) 附录 (12) 分段常数插值法源程序 (12) 分段线性插值法源程序 (12) 常加速度法源程序 (13) 线加速度法源程序 (13)
1.力插值法 力插值法对结构的外荷载进行插值，分为分段常数插值法和分段线性插值法，这两种方法均适用于线性结构的动力反应计算。 1.1分段常数插值法 图1-1为一个单自由度无阻尼系统，结构的刚度为k ，质量为m ，位移为y (t )，施加的外力为P (t )。图1-2为矩形脉冲荷载的示意图，图中t d 表示作用的时间，P 0表示脉冲荷载的大小。 图1-1 单自由度无阻尼系统示意图 图1-2 矩形脉冲荷载示意图 对于一个满足静止初始条件的无阻尼单自由度体系来说，当施加一个t d 时间的矩形脉冲荷载，此时结构在t d 时间内的位移反应可以用杜哈梅积分得到： 0()sin ()2 (1cos )(1cos ) (0) t st st d P y t t d m t y t y t t T ωττω πω=-=-=-≤≤? (1-1) 如果结构本身有初始的位移和速度，那么叠加上结构自由振动的部分，结构的位移反应为： 02()cos sin (1cos ) (0 )st d y t y t y t t y t t T πωωω =+ +-≤≤ (1-2)
图1-3 分段常数插值法微段示意图 对于施加于结构任意大小的力，将其划分为Δt 的微段，每一段的荷载都为一个常数(每段相当于一个矩形的脉冲荷载)，如图1-3所示，则将每一段的位移和速度写成增量的形式为： 1cos t sin t (1cos t)i i i i y P y y k ωωωω +=?+ ?+-? (1-3) i+1/sin t cos t sin t i i i y P y y k ωωωωω =-?+ ?+ ? (1-4) 程序流程图如下
结构风工程课后思考题参考答案 二、大气边界层风特性 1 对地表粗糙度的两种描述方式：指数律和对数律（将公式写上）。 2 非标准地貌下的风速换算原则（P）和方法（P公式）。1514 3 脉动风的生成： 近地风在流动过程中由于受到地表因素的干扰，产生大小不同的涡旋，这些涡旋的迭加作用在宏观上表现为速度的随机脉动。在接近地面时，由于受到地表阻力的影响，导致风速减慢并逐步发展为混乱无规则的湍流。 脉动风的能量及耗散机制：而湍流运动可以看做是能量由低频脉动向高频脉动过渡，并最终被流体粘性所耗散的过程。在低频区漩涡尺度较大，向中频区（惯性子区）、高频区（耗散区）漩涡尺度逐渐减小，小尺度涡吸收由惯性子区传递过来的能量，能量最终被流体粘性所耗散。 4 Davenport谱的特点：先写出公式 通过不同水平脉动风速谱的比较： （1）D谱不随高度变化，而其他谱（如Kaimal谱、Solari谱、Karman谱）则考虑了近地湍流随高度变化的特点；（D谱不随高度变化，在高频区符合-5/3律，没有考虑近地湍流随高度变化的特点；） （2）D谱的谱值比其它谱值偏大，会高估结构的动力反应，计算结果偏于保守。（3）S(0)=0,意味着L=0，与实际不符。uu5 湍流度随高度及地面粗糙程度的变化规律：随地面粗糙度的增大而增大，随高度的增加而减小。 积分尺度随高度及地面粗糙程度的变化规律：大量观测结果表明，大气边界层中的湍流积分尺度是地面粗糙度的减函数，而且随着高度的增加而增加。 功率谱随高度及地面粗糙程度的变化规律：随着高度增大和粗糙度的减小，能量在频率上的分布趋于集中，谱形显得高瘦；随着高度减小和粗糙度的增大，能量在频率上的分布趋于分散，谱形显得扁平。 相干函数随高度及地面粗糙程度的变化规律：随地面粗糙度的增大而减小，随高度的增加而增大。 6 阵风因子与峰值因子的区别：阵风因子G=U'/U,是最大风速与平均风速的比/ σ是最大脉动风速与脉动风速均方根的比值。g=u 值；峰值因子umax联系：二者可以相互换算：G=(U'+gσ)/U'=1+gσ/U'=1+gI。Uuu 三、钝体空气动力学理论 1 钝体绕流的主要特征有： ）粘性效应：气体粘性随温度升高而增大，液体粘性随温度升高而减小。1（（2）边界层的形成：由于粘性效应，使靠近物体表面的空气流动速度减慢，形 成气流速度从表面等于零逐渐增大到与外层气流速度相等，形成近壁面流动现象。 （3）边界层分离：如果边界层内的流体微粒速度因惯性力减小到使靠近表面的气流倒流，便出现了边界层分离。 （4）再附：在一定条件下，自建筑物前缘分离的边界层会偶然再附到建筑物表面，这时附面层下会形成不通气的空腔，即分离泡。每隔一段时间分离泡破裂产生较大的风吸值，产生一个风压脉冲。 （5）钝体尾流：对于细长钝体，漩涡脱落是在其两侧交替形成的。漩涡脱落时导致建筑物出现横向振动的主要原因。
哈尔滨理工大学毕业设计（论文）任务书毕业设计(论文)是学生毕业前最后一个重要学习环节，是学习深化与升华的重要过程。它既是学生学习、研究与实践成果的全面总结，又是对学生素质与能力的一次全面检验，而且还是对学生的毕业资格及学位资格认证的重要依据。为了保证我校本科生毕业设计(论文)质量，特制定“哈尔滨理工大学本科生毕业设计(论文)撰写规范”。 一、毕业设计(论文)资料的组成A.毕业设计(论文)任务书;B.毕业设计(论文)成绩评定书;C.毕业论文或毕业设计说明书(包括:封面、中外文摘要或设计总说明(包括关键词)、目录、正文、谢辞、参考文献、附录);D.译文及原文复印件;E.图纸、软盘等。 二、毕业设计(论文)资料的填写及有关资料的装订毕业设计(论文)统一使用学校印制的毕业设计(论文)资料袋、毕业设计(论文)任务书、毕业设计(论文)成绩评定书、毕业设计(论文)封面、稿纸(在教务处网上下载用，学校统一纸面格式，使用A4打印纸)。毕业设计(论文)资料按要求认真填写，字体要工整，卷面要整洁，手写一律用黑或蓝黑墨水;任务书由指导教师填写并签字，经院长(系主任)签字后发出。毕业论文或设计说明书要按顺序装订:封面、中外文摘要或设计总说明(包括关键词)、目录、正文、谢辞、参考文献、附录装订在一起，然后与毕业设计(论文)任务书、毕业设计(论文)成绩评定书、译文及原文复印件(订在一起)、工程图纸(按国家标准折叠装订)、软盘等一起放入填写好的资料袋内交指导教师查收，经审阅评定后归档。 三、毕业设计说明书(论文)撰写的内容与要求一份完整的毕业设计(论文)应包括以下几个方面: 1.标题标题应该简短、明确、有概括性。标题字数要适当，
哈工大2009年春季学期数据结构与算法 试 卷一、填空题（每空2分，共20分）1. 在 情况下，等长编码是最优前缀码。2．设有两个算法在同一机器上运行，其执行时间分别为100n 2和2n ，要使前者快于后者，n 至少为 。3．采用堆排序、快速排序、冒泡排序，对初态有序的表，最省时间的是_ 。4． 设二叉树结点的先根序列为ABDECFGH ，中根序列为DEBAFCHG,则二叉树中叶结点是_________.5. 用下标从0开始的N 个元素的数组实现循环队列时,为实现下标变量m 加1后在数组有效下标范围内循环,可采用的表达式是m= 。6. 由带权为3，9，4，2，5的5个叶子结点构成一棵哈夫曼树，则带权路径长度为 。7. 对n 个记录的表进行选择排序，在最坏情况下所需要进行的关键字的比较次数为 。8. 任意一个有n 个结点的二叉树，已知它有m 个叶结点，则度数为2的结点有 。9. n 个顶点的连通图用邻接矩阵表示时，该矩阵至少有 个非零元素10. 举出两种磁带文件的分类方法： 。二、选择题（每题1分，共10 分） 注意 行为 规范 遵守 考场 纪律
1．设一组初始记录关键字序列为(45，80，55，40，42，85)，则以第一个记录关键字45为基准而得到一趟快速排序的结果是（ ）。 (A) 40，42，45，55，80，83(B) 42，40，45，80，85，88 (C) 42，40，55，80，45，85(D) 42，40，45，85，55，80 2．数据的最小单位是（ ）。 (A) 数据项(B) 数据类型(C) 数据元素 (D) 数据变量 3．关键路径是AOE网中( ) 。 A.从始点到终点的最短路径 B.从始点到终点的最长路径 C.从始点到终点的边数最多的路径 D.从始点到终点的边数最少的路径 4．下列说法正确的是（）。 A．最小生成树也是哈夫曼树 B．最小生成树是唯一的 Ｃ．对于n 个顶点的连通无向图，Prim算法的时间复杂性为O(n2) D．Kruskal 算法比Prim算法更适合边稠密的图 5．设栈S和队列Q的初始状态为空，元素E1、E2、E3、E4、E5和E6依次通过栈S，一个元素出栈后即进入队列Q，若6个元素出列的 顺序为E2、E4、E3、E6、E5和E1，则栈S的容量至少应该是（ ）。 (A) 6(B) 4(C) 3(D) 2 6. 将10阶对称矩阵压缩存储到一维数组A中，则数组A的长度最 少为（ ）。 (A) 100 (B) 40(C) 55 (D) 80 7.若数据元素序列11，12，13，7，8，9，23，4，5是采用下列排序 方法之一得到的第二趟排序结果，则该排序算法只能是( )。 A. 插入排序 B.冒泡排序 C. 选择排序 D. 二路归并排序 8.设哈希表长m＝14，哈希函数H（key）＝key％11。表中已有4个结 点：addr(15)=4，addr(38)=5 ， addr(61)=6 ， addr(84)=7 其余地址 为空。如果用二次探测再散列处理冲突，关键字为49的结点的地址 是（） A．8 B .3 C. 5 D. 9 9. 有组记录的输入顺序为（46，79，56，38，40，84），则利用堆排 序方法建立的初始堆为( ) A.79，46，56，38，40，80 B .38，40，56，79，46，84 C. 84，79，56，46，40，38 D. 84，56，79，40，46，38 10. 下列叙述中，不符合m阶B树定义要求的是（）
哈尔滨工业大学计算机科学与技术学院 实验报告 课程名称：数据结构 课程类型：必修 实验项目名称：线性结构及其应用 实验题目：线性结构及其应用 一、实验目的
二、实验要求及实验环境 三、设计思想（本程序中的用到的所有数据类型的定义，主程序的流程图及各程序模块之间的调用关系） 1．逻辑设计 2．物理设计 四、测试结果 五、系统不足与经验体会 六、附录：源代码（带注释） 一、实验目的 输入中缀表达式保存并显示，之后转换为后缀表达式，并且求出表达式的结果。 二、实验要求及实验环境 实验要求 (1)从键盘输入任意一个语法正确的（中缀）表达式，显示并保存该表达式。 (2)利用栈结构，把上述（中缀）表达式转换成后缀表达式，并显示栈的状态变化过程 和所得到的后缀表达式。 (3)利用栈结构，对上述后缀表达式进行求值，并显示栈的状态变化过程和最终结果。 实验环境 Dev-C++软件中运行 Win7系统 三、设计思想 本实验中定义了int 型，char型，struct 型，char *型，struct型
逻辑设计：应用栈后进先出的规律，在转换为后缀表达式时，
将操作运算符压入栈中，碰见更高级运算符时栈中元素出栈，继续比较；否则压栈。这样可以完成表达式的转换。在利用得到的后缀表达式计算结果时，将操作数压栈，遇见符号直接计算，这是后缀表达式的特点。 物理设计：建立一个结构体数组的栈，数组中存放运算符。数组的添加和减少都在数组末尾元素进行。可以视为一个栈。 四、测试结果 样例1. 输入1+2*（3-4/2） 输出为1+2*（3-4/2） ->此为保存并输出的中缀表达式 12342/-*+ ->此为输出后缀表达式
第七章影响线 一判断题 1. 图示梁AB与A0B0，其截面C与C0弯矩影响线和剪力影响线完全相同。（X） 题1图题2图 2. 图示结构Q E影响线的AC段纵标不为零。（X） 3. 图示梁K截面的M K影响线、Q K影响线形状如图a、b所示。 4. 图示梁的M C影响线、Q C影响线形状如图a、b所示。 5. 图示梁的M C影响线、M B影响线形状如图a、b所示。 6. 图示结构M B影响线的AB段纵标为零。 7. 图示梁跨中C截面弯矩影响线的物理意义是荷载P=1作用在截面C的弯矩图形。（X） 8. 用静力法作静定结构某量值的影响线与用机动法作该结构同一量值的影响线是不等价 的。（X） 9. 求某量值影响线方程的方法，与恒载作用下计算该量值的方法在原理上是相同的。（√） 10. 影响线是用于解决活载作用下结构的计算问题，它不能用于恒载作用下的计算。（X） 11. 移动荷载是指大小，指向不变，作用位置不断变化的荷载，所以不是静力荷载。（X） 12. 用静力法作影响线，影响线方程中的变量x代表截面位置的横坐标。（X） 13. 表示单位移动荷载作用下某指定截面的内力变化规律的图形称为内力影响线。（√） 14. 简支梁跨中截面弯矩的影响线与跨中有集中力P时的M图相同。（X） 15. 简支梁跨中C截面剪力影响线在C截面处有突变。（√） 16. 绝对最大弯矩是移动荷载下梁的各截面上最大的弯矩。（√） 17. 静定结构及超静定结构的内力影响线都是由直线组成。（X） 18. 图示结构Q C影响线的CD段为斜直线。 19. 图示结构K断面的剪力影响线如图b所示。（√） 题19图 20. 用机动法作得图a所示Q B左结构影响线如图b。 题20图题21图 21. 图示结构a杆的内力影响线如图b所示 22. 荷载处于某一最不利位置时，按梁内各截面得弯矩值竖标画出得图形，称为简支梁的弯
第一章 单自由度系统 1.1 总结求单自由度系统固有频率的方法和步骤。 单自由度系统固有频率求法有：牛顿第二定律法、动量距定理法、拉格朗日方程法和能量守恒定理法。 1、 牛顿第二定律法 适用围：所有的单自由度系统的振动。 解题步骤：（1） 对系统进行受力分析,得到系统所受的合力； （2） 利用牛顿第二定律∑=F x m ,得到系统的运动微分方程； （3） 求解该方程所对应的特征方程的特征根，得到该系统的固有频率。 2、 动量距定理法 适用围：绕定轴转动的单自由度系统的振动。 解题步骤：（1） 对系统进行受力分析和动量距分析； （2） 利用动量距定理J ∑=M θ ,得到系统的运动微分方程； （3） 求解该方程所对应的特征方程的特征根，得到该系统的固有频率。 3、 拉格朗日方程法： 适用围：所有的单自由度系统的振动。 解题步骤：（1）设系统的广义坐标为θ，写出系统对于坐标θ的动能T 和势能U 的表达式；进一步写求出拉格朗日函数的表达式：L=T-U ； （2）由格朗日方程 θ θ??-???L L dt )( =0，得到系统的运动微分方程； （3） 求解该方程所对应的特征方程的特征根，得到该系统的固有频率。 4、 能量守恒定理法 适用围：所有无阻尼的单自由度保守系统的振动。
解题步骤：（1）对系统进行运动分析、选广义坐标、写出在该坐标下系统的动能T 和势能U 的表达式；进一步写出机械能守恒定理的表达式 T+U=Const （2）将能量守恒定理T+U=Const 对时间求导得零,即0) (=+dt U T d ，进一步得到系统的运动微分方程； （3） 求解该方程所对应的特征方程的特征根，得到该系统的固有频率。 1.2 叙述用衰减法求单自由度系统阻尼比的方法和步骤。 用衰减法求单自由度系统阻尼比的方法有两个：衰减曲线法和共振法。 方法一：衰减曲线法。 求解步骤：（1）利用试验测得单自由度系统的衰减振动曲线，并测得周期和相邻波峰和波谷的幅值i A 、1+i A 。 （2）由对数衰减率定义 )ln( 1 +=i i A A δ， 进一步推导有 2 12ζ πζδ-= ， 因为ζ较小， 所以有 π δζ2= 。 方法二：共振法求单自由度系统的阻尼比。 （1）通过实验，绘出系统的幅频曲线， 如下图：
Decoupling Scheme of Tracking Loop of Seeker Based on Disturbance Compensation XXXXXX1, XXXXXXX 2, XXXXXXXX 2 （1.XXXXXXXXXXXXXXXXXXXXXXX; 2XXXXXXXXXXXXXXXXXXXXXX） Abstract：To get line-of-sight rate of high-accuracy for guidance, finishing aerial attack and intercept when the missile is controlled by lateral jets which will cause high frequency disturbance to missile attitude, a new decoupling scheme is proposed based on disturbance compensation in this paper. Decoupling-ability, tracking performance and noise rejection ability are analyzed and compared with the two existing method. Simulation and analysis show the effectiveness of the scheme proposed. Keywords：Guidance; Lateral jets; Line-of-sight rate; Decoupling; high-accuracy; compensation Introduction Contemporary military strategy has developed from underlining quantity advantage to technical quality advantage. Precision guided weapon is the main way of physical killing, and perform an important role in informatization local wars [1]-[3]. Precision guidance technology is critical for the precision guided weapon, so it is always studied by experts and researchers from many countries. Seeker is a kind of measuring component in controlling loop of homing guidance [4]. It is required to output information of line-of-sight (LOS for short) rate for guidance. So the precision of its output, LOS rate, affects controlling precision of terminal guidance directly. While a seeker is in the state of tracking, the antenna on it should track the direction of the target. A seeker works on a moving missile. So the disturbance caused by the missile body movement must be separated from the antenna to keep it stable in the inertial space [4]-[5]. The traditional design idea is to use a stabilization loop to attenuate the LOS output caused by body disturbance. We call it traditional multi-loop (TML for short) scheme [4]-[6]. There are also some schemes based on TML scheme, for example, decoupling-loop scheme which changes the position of the loop’s output to improve the decoupling-ability and at the same time have some advantages on design of guidance [7]-[8]. However, it is difficult to decouple disturbance by using the schemes mentioned above. Line-of-sight reconstruction(LOSR for short) scheme, by F.William and Paul Zarchan, can decouple completely in ideal conditions in which the transfer function of the receiver and the gyro are both equal to 1 [9]. In fact, these two conditions can’t be satisfied in the presence of high-frequency disturbance. In modern war, missiles need larger maneuverability to intercept aerial target with high speed and large maneuverability. New types of physical actuators bring high-frequency characteristics to missile movement. The high-frequency characteristics require more effective decoupling method for the seeker. Considering weight and cost, the open-loop gain can’t be too high while the time constant of the servo motor and the gyro fixed on the antenna can’t be too small. These restrictions limit further enhance of decoupling-ability. The existing schemes can only reject disturbance in low frequency segment and very high frequency segment effectively. In the frequency segment which missile moving disturbance can actually reach, decoupling-ability of seekers becomes weaker when missile movement frequency turns high. In fact, missile disturbance is able to be measured. Gyro fixed on missile for guidance and attitude
结构动力学试题 2016年4月 重庆交通大学结构工程硕士研究生考试 1.试述结构动力问题和静力问题的主要区别（10分） 答：结构静力学相比，动力学的复杂性表现在： （1）动力问题具有随时间而变化的性质； （2）数学解答不是单一的数值，而是时间的函数； （3）惯性力是结构内部弹性力所平衡的全部荷载的一个重要部分； （4）引入惯性力后涉及到二阶微分方程的求解； （5）需考虑结构本身的动力特性：刚度分布、质量分布、阻尼特性分布的影响。 2.什么是结构动力系统的阻尼？一般结构系统的阻尼有何特性？在结构分析中 阻尼问题的处理方法有哪些？（20分） 答：（1）结构在震动过程中的能量耗散作用称为阻尼； （2）阻尼的特性：a、阻尼耗能与质量（反映附属部分大小）和刚度（反映位移大小）有关。b、难以采用精确的理论分析方法； （3）对于多自由度体系：在结构动力分析中，通常从系统响应这个角度来考虑阻尼，而且能量的损耗是由外界激励来平衡的。一个振动系统可能存在多种不同类型的阻尼，一般来说，要用数学的方法来精确描述阻尼目前是比较困难的。因此，人们根据经验提出了一些简化模型，常用的阻尼模型有黏性阻尼和结构阻尼。黏性阻尼系统：黏性阻尼的特点是阻尼力和运动速度成真封闭。 在用振型叠加法进行分析时，能否将联立的运动方程化为解耦的一系列单自由度运动方程，将取决于阻尼矩阵的性质，即结构的振型是否关于阻尼阵满足正交条件。如果满足阻尼阵的正交条件，则采用振型叠加法分析时，就可以把多自由度体系的动力反应问题化为一系列单自由度问题求解；如果不满足阻尼阵的正交条件，则对位移向量用振型展开后，关于振型坐标的运动方程成为耦联的，必须联立求解，与解耦方程相比，增加了难度和计算量。 3.试述多自由度体系振型矩阵关于质量矩阵和刚度矩阵的正交性的意义，并写出广义正交性的表达式且加以证明。（20分） 答：（1）由振型关于质量、刚度正交性公式可知，i振型上的惯性力在j振型上作的虚功为0。由此可知，既然每一主振型相应的惯性力在其他主振型上不做功，那么它的振动能量就不会转移到别的主振型上去。换句话说，当一个体系只按某一主振型振动时，不会激起其他主振型的振动。这说明各个主振型都能单独出现，彼此线性无关。这就是振型正交的物理意义。一是可用于校核振型的正确性；二是在已知振型的条件下，可以通过折算质量与折算刚度计算对应的频率。而更主要的是任一同阶向量均可用振型的线性组合来表示，在受迫振动分析中，利用振型的正交性，在阻尼矩阵正交的假设下可使运动方程解藕. （2）振型正交性的证明在Clough书中应用的是Betti互易定理，就像D’Alember 原理一样考虑了惯性力，是运动学中功的互等定理。实际振型正交性的证明可
春季学期MATLAB期末作业 学院：机电工程学院 专业：机械制造设计及其自动化 学号： 班号： 姓名：
2013年春季学期 MATLAB 课程考查题 姓名: 学号： 学院：机电学院 专业：机械制造 一、必答题： 1.matlab常见的数据类型有哪些？各有什么特点? 常量：具体不变的数字 变量：会根据已知条件变化的数字 字符串：由单引号括起来的简单文本 复数：含有复数的数据 2.MATLAB中有几种帮助的途径？ （1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器； （2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息； （3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数 （4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息。 3.Matlab常见的哪三种程序控制结构及包括的相应的语句？ 1.顺序结构：数据输入A=input(提示信息,选项) 数据输出disp(X) 数据输出fprintf(fid,format,variables) 暂停pause 或 pause(n) 2.选择结构： If语句： if expression （条件）
statements1（语句组1） else statements2（语句组2） End Switch 语句： switch expression （表达式） case value1 （表达式1） statement1（语句组1） case value2 （表达式2） statement2（语句组2） ... ... case valuem （表达式m） statementm（语句组m） otherwise statement （语句组） end 3.循环结构： For循环： for variable=expression(循环变量) statement（循环体） end While循环： while expression （条件<循环判断语句>） statement（循环体） end 4.命令文件与函数文件的主要区别是什么？
结构动力学大作业（威尔逊- 法） ： 学号： 班级： 专业：
威尔逊-θ法原理及应用 【摘要】在求解单自由度体系振动方程时我们用了常加速度法及线加速度法等数值分析方法。在多自由度体系中，也有类似求解方法，即中心差分法及威尔逊-θ法。实际上后两种方法也能求解单自由度体系振动方程。对于数值方法，有三个重要要求：收敛性、稳定性及精度。本文推导了威尔逊-θ法的公式，并利用MATLAB 编程来研究单自由度体系的动力特性。 【关键词】威尔逊-θ法 冲击荷载 阻尼比 【正文】威尔逊-θ法可以很方便的求解任意荷载作用下单自由度体系振动问题。实际上，当 1.37θ>时，威尔逊-θ法是无条件收敛的。 一、威尔逊-θ法的原理 威尔逊-θ法是线性加速度法的一种拓展（当1θ=时，两者相同），其基本思路和实现方法是求出在时间段[],t t t θ+?时刻的运动，其中1θ≥，然后通过插得到i t t +?时刻的运动（见图 1.1）。 图 1.1 1、公式推导 推导由t 时刻的状态求t t θ+?时刻的状态的递推公式： 对τ积分
{}{}{}{}{}{})(623 2 t t t t t t t y y t y y y y &&&&&&&-?+++=?++θτ θτττ {}{}{}{}{})2(6)(2t t t t t t t y y t y t y y &&&&&+?+?+=?+?+θθθθ {}{}{}{}{}t t t t t t t y y t y y t y &&&&&26 )()(62-?--?=?+?+θθθθ []{}{} {}[]{}{}{}[]{}{}{})223()26)(6( )(2t t t t t t t t t t y t y y t c y y t y t m P P P R &&&&&&?++?++?+?+-+=?+θθθθθ 2、MA TLAB 源程序： clc;clear; K=input('请输入结构刚度k(N/m)'); M=input('请输入质量(kg)'); C=input('请输入阻尼(N*s/m)'); t=sym('t');%产生符号对象t Pt=input('请输入荷载); Tp=input('请输入荷载加载时长(s)'); Tu=input('请输入需要计算的时间长度(s) '); dt=input('请输入积分步长(s)'); Sita=input('请输入θ'); uds=0:dt:Tu;%确定各积分步时刻 pds=0:dt:Tp; Lu=length(uds); Lp=length(pds); if isa(Pt,'sym')%荷载为函数 P=subs(Pt,t,uds); %将荷载在各时间步离散 if Lu>Lp P(Lp+1:Lu)=0; end elseif isnumeric(Pt)%荷载为散点 if Lu<=Lp
XX大学 硕士学位论文开题报告 题目： 院（系）____________________________ 学科____________________________ 导师____________________________ 研究生 __________________________________ 学号____________________________ 开题报告日期 ______________________________ 研究生院制 二?一二年三月
一、开题报告应包括下列主要内容： 目录 1．课题来源及研究的背景和意义 1.1．课题的来源 1.2．课题研究的背景和意义（不少于 500 字） 2．国内外在该方向的研究现状及分析 2.1．国外研究现状 2.2．国内研究现状 （注意对所引用国内外文献的准确标注） 2.3．国内外文献综述的简析（不少于 500 字）（综合评述：国内外研究取得的成果，存在的不足或有待深入研究的问题）3．主要研究内容（不少于 1000 字） （撰写宜使用将来时态，切忌将论文目录直接作为研究内容，要突出本人研究内容）4．研究方案及进度安排，预期达到的目标和取得的研究成果 4.1．研究方案（不少于 500 字） 4.2．预期达到的目标和取得的研究成果 4.3．进度安排（建议从进入研究课题时间开始） 5．为完成课题已具备和所需的条件和经费 6．预计研究过程中可能遇到的困难和问题，以及解决的措施 7．主要参考文献 二、对开题报告的要求 1．开题报告的字数应在 5000 字以上； 2．参考文献应在___篇以上，其中外文资料应不少于三分之一。硕士研究生应在导师的指导下着重查阅近年内发表的中、外文期刊文章，参考的近五年内（从开题时间算起）文献一般不少于三分之一。本学科的基础和专业课教材一般不应列为参考文献。 三、开题报告时间应最迟不得超过第三学期的第三周末。 四、如硕士生首次开题报告未通过，需在一个月内再进行一次。若仍不通过，则停止硕士论文工作。 五、开题报告进行后，此表同硕士学位论文开题报告评议结果存各系（院）研究生秘书书处，以备 研究生院及所属学院进行检查。 六、字体、字号及其他规定 论文中所用中文字体（除各级标题外）为宋体，各级标题用黑体；论文中所用数字、英文为新罗马字体。 节标题小__号字，建议段前0.5行，段后0.5行； 条标题—号字，建议段前0.5行，段后0.5行；
软件工程复习题 一：问答题： 1.结构化系统开发方法的基本思想是什么？ 结构化分析建模是根据结构化分析准则，在需求分析过程中建立的3种模型： （1）数据模型：用实体-联系图描绘数据对象及数据对象之间的关系，以建立系统的数据模型。 （2）功能模型：用数据流图描绘当数据在软件系统中移动时被变换的逻辑过程，以指名系统具有的变换数据的功能，因此，数据流图即建立系统的功能模型。 （3）行为模型：用状态转换图指名作为外部事件结果的系统行为，即状态转换图描绘了系统的各种行为模式和在不同状态间转换的方式，因此，状态转换图即建立系统的行为模型。2.名词解释：软件工程，软件危机，软件开发模型，软件测试，软件项目的质量管理，软件； 成本/效益分析，软件生命周期； 软件工程:1：将系统的，规范的，可定量的方法用到软件的开发、运行、维护之中；2对1中的方法的研究。三要素：工具、过程、方法 软件危机：在软件的开发和维护过程中遇到的一系列严重问题 软件开发模型：瀑布，快速原型，增量，螺旋，形式化方法，基于组件的开发模型 软件测试：广义：对软件生存周期内的所有检查；狭义：对软件产品的质量的检查和评价。就是找出软件的错误 软件项目的质量管理: 制定软件质量保证计划，按照质量评价体系控制软件质量要素，对阶段性的软件产品进行评审，对最终软件产品进行确认，确保软件质量。。 软件：按照特定顺序组织的计算机指令和数据的集合 成本/效益分析：通过成本效益分析来评估项目价值的方法 软件生命周期：从软件的问题定义一直到该软件被淘汰 3.什么是软件危机？软件危机的表现是什么？其产生的原因是什么？ 软件危机：计算机软件在它的开发和维护过程中所遇到的一系列严重问题
数据结构试卷（一） 一、单选题（每题2 分，共20分） 1.栈和队列的共同特点是( )。 A.只允许在端点处插入和删除元素 B.都是先进后出 C.都是先进先出 D.没有共同点 2.用链接方式存储的队列，在进行插入运算时( ). A. 仅修改头指针 B. 头、尾指针都要修改 C. 仅修改尾指针 D.头、尾指针可能都要修改 3.以下数据结构中哪一个是非线性结构？( ) A. 队列 B. 栈 C. 线性表 D. 二叉树 4.设有一个二维数组A[m][n]，假设A[0][0]存放位置在644(10)，A[2][2]存放位置在 676(10)，每个元素占一个空间，问A[3][3](10)存放在什么位置？脚注(10)表示用10进制表示。 A．688 B．678 C．692 D．696 5.树最适合用来表示( )。 A.有序数据元素 B.无序数据元素 C.元素之间具有分支层次关系的数据 D.元素之间无联系的数据 6.二叉树的第k层的结点数最多为( ). A．2k-1 B.2K+1 C.2K-1 D. 2k-1 7.若有18个元素的有序表存放在一维数组A[19]中，第一个元素放A[1]中，现进行二 分查找，则查找A［3］的比较序列的下标依次为( ) A. 1，2，3 B. 9，5，2，3 C. 9，5，3 D. 9，4，2，3 8.对n个记录的文件进行快速排序，所需要的辅助存储空间大致为 A. O（1） B. O（n） C. O（1og2n） D. O（n2） 9.对于线性表（7，34，55，25，64，46，20，10）进行散列存储时，若选用H（K） =K %9作为散列函数，则散列地址为1的元素有（）个， A．1 B．2 C．3 D．4 10.设有6个结点的无向图，该图至少应有( )条边才能确保是一个连通图。 A.5 B.6 C.7 D.8 二、填空题（每空1分，共26分） 1.通常从四个方面评价算法的质量：_________、_________、_________和_________。 2.一个算法的时间复杂度为(n3+n2log2n+14n)/n2，其数量级表示为________。 3.假定一棵树的广义表表示为A（C，D（E，F，G），H（I，J）），则树中所含的结点数 为__________个，树的深度为___________，树的度为_________。 4.后缀算式9 2 3 +- 10 2 / -的值为__________。中缀算式（3+4X）-2Y/3对应的后缀算式 为_______________________________。 5.若用链表存储一棵二叉树时，每个结点除数据域外，还有指向左孩子和右孩子的两个指 针。在这种存储结构中，n个结点的二叉树共有________个指针域，其中有________个指针域是存放了地址，有________________个指针是空指针。 6.对于一个具有n个顶点和e条边的有向图和无向图，在其对应的邻接表中，所含边结点 分别有_______个和________个。 7.AOV网是一种___________________的图。 8.在一个具有n个顶点的无向完全图中，包含有________条边，在一个具有n个顶点的有 向完全图中，包含有________条边。
第四章 树与二元树 填空题 1.假定一棵树的广义表表示为A(B(E)，C(F(H，I，J)，G)，D)，则该树的度为 ① ，树 高度为 ② ，终端结点的个数为 ③ ，单分支节点的个数为 ④ ，双分支结点的个数为 ⑤ ，三分支结点的个数为 ⑥ ，C结点的双亲结点为 ⑦ ，其孩子结点 ⑧ 和 ⑨ 结。该树先根、中根和后根遍历序列分别为 ⑽ 、⑾ 和⑿。该树对应的 二元树为 ⒀ ，此二元树的先根、中根和后根遍历顺序序列分别为⒁、⒂和⒃。 2.由带权为3，9，6，2，5的5个叶子结点构成一棵哈夫曼树，则带权路径长度为 ① ， 该最优二元树共有 ② 个结点，度数为0、1、2的结点的个数分别为③ ，④ 和 ⑤ 个。 3.已知字符集{A、B、C、D、E} 的字符出现的概率分别为{ 3/25 ，9/25，6/25，2/25， 5/25}。画出该字符集的Huffman编码树② ， 字符A、B、C、D、E的编码分别为 ③， ④ ，⑤ ，⑥ ，⑦ ，该字符集的Huffman编码的平均编码长度为⑧ 。若采用二进制 等长编码方案，该字符集的编码长度为 ⑨ 。读该字符集而言，Huffman编码比等长编码平均压缩了 ⑽ %。 4.对于一棵具有n个结点的二元树，当进行链接存储时，其左右链存储结构中的指针域的 总数为 ①个，其中，② 个用于链接孩子结点， ③个空闲着。 5.在一棵二叉树中，度为0的结点个数为n0，度为1的结点个数为n1，度为2的结点个 数为n2，则有n0= ① 。 6.由a，b，c 三个结点构成的二叉树，共有 ① 种不同结构。 7.一棵高度为K的完全二叉树的结点总数最少为 ① 个，最多为 ② 个；第K层最多有 ③ 个结点，最少有 ④ 个结点。 选择题 8.假定在一棵二元树中，双分支结点数为15，单分支结点数为30，则叶子结点数为( ) 个。 A．15 B．16 C．17 D．47 9.在一棵二叉树上第5层的结点数最多为( ) 。 A．8 B．16 C．15 D．32 10.用顺序存储的方式将完全二叉树中的所有结点逐层存放在数组R[ 1…n]中，结点R[i] 若有子树，则左子树是结点( )。