北航《金属工艺学》在线作业三
北航《金属工艺学》在线作业三
一、单选题（共 15 道试题，共 30 分。）
1. 纯铁的熔点是（）。
A. 1538℃
B. 1438℃
C. 1338℃
D. 1238℃
正确答案：A
2. 工程上所用的材料，一般要求其屈强比（）。
A. 越大越好
B. 越小越好
C. 大些，但不可过大
D. 小些，但不可过小
正确答案：C
3. 铁碳合金相图中，ACD线是（）。
A. 液相线
B. 固相线
C. 共晶线
D. 共析线
正确答案：A
4. （）是将淬火钢重新加热到A1点以下的某一温度，保温一定时间后，冷却到室温的一种操作。
A. 正火
B. 回火
C. 退火
D. 淬火
正确答案：B
5. 共晶转变的产物是（）。
A. 奥氏体
B. 渗碳体
C. 珠光体
D. 莱氏体
正确答案：D
6. 在钢的整体热处理过程中，常用的作预备热处理的是（）。
A. 正火和淬火
B. 正火和回火
C. 淬火和回火
D. 退火和正火
正确答案：D
7. 造船用的碳素钢属于（）。
A. 普通质量碳素钢
B. 优质碳素钢
C. 特殊质量碳素钢
D. 以上都不是
正确答案：B
8. 塑料的主要成分（）。
A. 稳定剂
B. 填充剂
C. 增强剂
D. 树脂
正确答案：D
9. 金属的韧性通常随加载速度提高、温度降低、应力集中程度加剧而（）。
A. 变好
B. 变差
C. 无影响
D. 难以判断
正确答案：B
10. 钢的含碳量一般在（）。
A. 0.77%以下
B. 2.11%以下
C. 4.3%以下
D. 6.69%以下
正确答案：B
11. 金属疲劳的判断依据是（）。
A. 强度
B. 塑性
C. 抗拉强度
D. 疲劳强度
正确答案：D
12. 用金刚石圆锥体作为压头可以用来测试（）。
A. 布氏硬度
B. 洛氏硬度
C. 维氏硬度
D. 以上都可以
正确答案：B
13. 为了减小淬火内应力和降低脆性，表面淬火后一般要（）。
A. 正火
B. 低温回火
C. 中温回火
D. 高温回火
正确答案：B
14. 玻璃钢是（）。
A. 复合材料
B. 塑料
C. 一种钢
D. 陶瓷
正确答案：A
15. Wc<0.77%铁碳合金冷却至A3线时，将从奥氏体中析出（）。
A. 铁素体
B. 渗碳体
C. 珠光体
D. 莱氏体
正确答案：A
北航《金属工艺学》在线作业三
二、多选题（共 10 道试题，共 20 分。）
1. 金属的晶格类型主要有（）三大类。
A. 体心立方晶格
B. 面心立方晶格
C. 密排八方晶格
D. 密排六方晶格
正确答案：ABD
2. 通常用（）热处理热处理来强化钢的表面层。
A. 表面
B. 化学
C. 机械
D. 物理
正确答案：AB
3. 下列选项属于塑料按树脂性质不同分类的是（）。
A. 通用塑料
B. 工程塑料
C. 热固性塑料
D. 热塑性塑料
正确答案：CD
4. 按所有、用气体的不同，气体保护电弧焊有（）等。
A. 氩弧焊
B. 氮气焊
C. 埋弧焊
D. 二氧化碳气体保护焊
正确答案：AD
5. 下列哪几项属于铸件的特点？（）
A. 适应性广
B. 铸件质量高
C. 生产成本低
D. 形状和尺寸精度高
正确答案：AC
6. 零件的失效形式主要有（）。
A. 过量变形失效
B. 断裂失效
C. 表面损伤失效
D. 裂纹失效
正确答案：ABCD
7. 切削用量的三要素是（）。
A. 切削速度
B. 进给量
C. 背吃刀量
D. 侧吃刀量
正确答案：ABC
8. 常用的硬度表示方法有（）。
A. 布氏硬度
B. 维氏硬度
C. 俚氏硬度
D. 洛氏硬度
正确答案：ABD
9. 钢的热处理工艺曲线包括（）三个阶段。
A. 高温
B. 加热
C. 保温
D. 冷却
正确答案：BCD
10. 合金钢按主要质量等级分为（）。
A. 优质合金钢
B. 特殊质量合金钢
C. 碳素钢
D. 工具钢
正确答案：AB
北航《金属工艺学》在线作业三
三、判断题（共 25 道试题，共 50 分。）
1. 工艺是劳动者利用生产工具对各种原材料、半成品进行增值加工或处理，最终使之成为制成品的方法与过程。（）
A. 错误
B. 正确
正确答案：B
2. 晶体在固态下原子在物质内部做有规则排列，而非晶体在固态下物质内部原子成无序堆积状态。（）
A. 错误
B. 正确
正确答案：B
3. 镦粗是使坯料高度缩小，横截面积增大的锻造工序。（）
A. 错误
B. 正确
正确答案：B
4. 焊条由焊芯和药皮两部分构成。（）
A. 错误
B. 正确
正确答案：B
5. 自由锻造的基本工序主要有镦粗、拔长、冲孔、弯曲、切断等。（）
A. 错误
B. 正确
正确答案：B
6. 金属工艺学是研究工程上常用材料性能和加工工艺的一门综合性的技术基础课。（）
A. 错误
B. 正确
正确答案：B
7. 铸铁的石墨化就是指从铸铁中结晶析出石墨的过程。（）
A. 错误
B. 正确
正确答案：B
8. 零件的使用要求包括零件形状、尺寸、加工精度和表面粗糙度等的外部质量要求，以及具体工作条件下对零件成分、组织、性能的内部质量要求。（）
A. 错误
B. 正确
正确答案：B
9. 细化晶粒的主要措施有增加过冷度和液态金属在结晶前加入变质剂进行变质处理两种。（）
B. 正确
正确答案：B
10. 常用的退火方法有完全退火、球化退火和去应力退火。（）
A. 错误
B. 正确
正确答案：B
11. 金属型铸造比砂型铸造尺寸精度低。（）
A. 错误
B. 正确
正确答案：A
12. 回火的目的是降低淬火钢的脆性，减少或消除内应力，使组织趋于稳定并获得所需要的性能。（）
A. 错误
B. 正确
正确答案：B
13. 粗加工时车刀角度前后角都大。（）
A. 错误
B. 正确
正确答案：A
14. 淬火时，在水中加入盐或碱，可增加在650℃-550℃范围内的冷却速度，避免产生软点。（）
A. 错误
B. 正确
正确答案：B
15. 锻造前坯料加热的目的在于提高塑性和降低变形抗力。（）
A. 错误
B. 正确
正确答案：B
16. 金属在力作用下所显示与弹性和非弹性反应相关或涉及力—应变关系的性能，叫做金属力学性能（机械性能）。（）
A. 错误
B. 正确
正确答案：B
17. 在我国低合金高强度钢主要加入的元素为硅。（）
A. 错误
B. 正确
正确答案：A
18. 测定某种钢的力学性能时，已知试棒的直径是10mm，其标距长度是直径的五倍，Fs=20.68KN，则此钢的σs为263.44Mpa。（）
A. 错误
B. 正确
正确答案：B
19. 铸铁是一系列主要由铁、碳、硅组成的合金的总称。（）
B. 正确
正确答案：B
20. 冒口的主要作用是排气和补缩。（）
A. 错误
B. 正确
正确答案：B
21. 拉伸实验中，试样所受的力为动态力。（）
A. 错误
B. 正确
正确答案：A
22. 合金是指由两种或两种以上化学元素组成的具有金属特性的物质。（）
A. 错误
B. 正确
正确答案：B
23. 锻造前对金属进行加热，目的是提高塑性。（）
A. 错误
B. 正确
正确答案：B
24. 机床的床身一般选锻造。（）
A. 错误
B. 正确
正确答案：A
25. 自由锻造按使用设备不同，又可分为手工锻造和机器锻造。（）
A. 错误
B. 正确
正确答案：B
北航最优化方法大作业参考
北航最优化方法大作业参考
1 流量工程问题 1.1 问题重述 定义一个有向网络G=(N,E)，其中N是节点集，E是弧集。令A是网络G的点弧关联矩阵，即N×E阶矩阵，且第l列与弧里(I,j)对应，仅第i行元素为1，第j行元素为-1，其余元素为0。再令b m=(b m1,…,b mN)T，f m=(f m1,…,f mE)T，则可将等式约束表示成： Af m=b m 本算例为一经典TE算例。算例网络有7个节点和13条弧，每条弧的容量是5个单位。此外有四个需求量均为4个单位的源一目的对，具体的源节点、目的节点信息如图所示。这里为了简单，省区了未用到的弧。此外，弧上的数字表示弧的编号。此时，c=((5,5…,5)1 )T， ×13 根据上述四个约束条件，分别求得四个情况下的最优决策变量x=((x12,x13,…,x75)1× )。 13 图 1 网络拓扑和流量需求
1.2 7节点算例求解 1.2.1 算例1（b1=[4;-4;0;0;0;0;0]T） 转化为线性规划问题： Minimize c T x1 Subject to Ax1=b1 x1>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x1*=[4 0 0 0 0 0 0 0 0 0 0 0 0]T 对应的最优值c T x1=20 1.2.2 算例2（b2=[4;0;-4;0;0;0;0]T） Minimize c T x2 Subject to Ax2=b2 X2>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x2*=[0 4 0 0 0 0 0 0 0 0 0 0 0]T 对应的最优值c T x2=20 1.2.3 算例3（b3=[0;-4;4;0;0;0;0]T） Minimize c T x3 Subject to Ax3=b3 X3>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x3*=[4 0 0 0 4 0 0 0 0 0 0 0 0]T 对应的最优值c T x3=40
北航数值分析大作业一
《数值分析B》大作业一 SY1103120 朱舜杰 一．算法设计方案： 1.矩阵A的存储与检索 将带状线性矩阵A[501][501]转存为一个矩阵MatrixC[5][501] . 由于C语言中数组角标都是从0开始的，所以在数组MatrixC[5][501]中检索A的带内元素a ij的方法是： A的带内元素a ij=C中的元素c i-j+2,j 2.求解λ1，λ501，λs ①首先分别使用幂法和反幂法迭代求出矩阵按摸最大和最小的特征值λmax和λmin。λmin即为λs； 如果λmax>0,则λ501=λmax；如果λmax<0,则λ1=λmax。 ②使用带原点平移的幂法（mifa（）函数），令平移量p=λmax，求 出对应的按摸最大的特征值λ，max， 如果λmax>0,则λ1=λ，max+p；如果λmax<0,则λ501=λ，max+p。 3.求解A的与数μk=λ1+k（λ501-λ1）/40的最接近的特征值λik （k=1,2，…，39）。 使用带原点平移的反幂法，令平移量p=μk，即可求出与μk最接近的特征值λik。 4.求解A的（谱范数）条件数cond（A）2和行列式d etA。 ①cond（A）2=|λ1/λn|，其中λ1和λn分别是矩阵A的模最大和 最小特征值。
②矩阵A的行列式可先对矩阵A进行LU分解后，detA等于U所有对角线上元素的乘积。 二．源程序 #include 机械原理大作业 This model paper was revised by the Standardization Office on December 10, 2020 机械原理大作业三 课程名称：机械原理 设计题目：齿轮传动设计 院系： 班级： 设计者： 学号： 指导教师： 设计时间： 1、设计题目 机构运动简图 机械传动系统原始参数 2、传动比的分配计算 电动机转速min /745r n =，输出转速m in /1201r n =，min /1702r n =， min /2303r n ，带传动的最大传动比5.2max =p i ,滑移齿轮传动的最大传动比4m ax =v i ，定轴齿轮传动的最大传动比4m ax =d i 。 根据传动系统的原始参数可知，传动系统的总传动比为： 传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三部分实现。设带传动的传动比为5.2max =p i ，滑移齿轮的传动比为321v v v i i i 、、，定轴齿轮传动的传动比为f i ，则总传动比 令 4max 1==v v i i 则可得定轴齿轮传动部分的传动比为 滑移齿轮传动的传动比为 设定轴齿轮传动由3对齿轮传动组成，则每对齿轮的传动比为 3、齿轮齿数的确定 根据滑移齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮5、6、7、8、9和10为角度变位齿轮，其齿数： 35,18,39,14,43,111098765======z z z z z z ；它们的齿顶高系数1=* a h ，径向间 隙系数25.0=*c ，分度圆压力角020=α，实际中心距mm a 51'=。 一、齿轮传动的基本概念 渐开线齿轮的啮合特点：(1)渐开线齿廓能够保证定传动比;(2)渐开线齿廓之间的正压力方向不变;(3)渐开线齿廓传动具有可分性。 齿轮机构的特点是：传动平稳、适用范围广、效率高、结构紧凑、工作可靠、寿命长。但制造和安装精度高、制造费用大，且不适合于距离较远的两轴之间的传动。齿轮传动可以用来传递任意轴间的运动和动力。 齿轮传动按照一对齿轮传递的相对运动分为平面齿轮传动和空间齿轮传动，平面齿轮传动又分为直齿圆柱齿轮传动、斜齿圆柱齿轮传动和人字齿轮传动;按照工作条件可以分为开式传动、半开式传动和闭式传动。 齿轮传动的基本要求是：传动准确、平稳;承载能力强。 二、齿轮传动的设计与计算 齿廓曲线与齿廓啮合基本定律：在啮合传动的任一瞬时，两轮齿廓曲线在相应接触点的功法线必须通过按给定传动比确定的该瞬时的节点。 渐开线齿轮啮合的正确条件：啮合轮齿的工作侧齿廓的啮合点必须总是在啮合线上，即两齿轮的模数和压力角应该分别相等。 齿轮传动的无侧隙啮合及标准齿轮的安装：一个齿轮节圆上的齿厚等于另一个齿轮节圆上的齿槽宽是无侧隙啮合的条件;外啮合齿轮的标准中心距为，内啮合是标准中心距为。 齿轮及其变位的相关计算：相关参数为齿数、模数、分度圆压力角、齿顶高系数和顶隙系数及标准直齿轮的几何尺寸计算，包括分度圆直径、齿顶高、齿根高、齿全高、齿顶圆直径、齿根圆直径、基圆直径、齿距、齿厚、齿槽宽、中心距、顶隙以及变位齿轮的变位系数等。 渐开线齿轮的根切现象：用展成法加工齿轮式，若刀具的齿顶线或齿顶圆与啮合线的焦点超过被切齿轮的极限点，则刀具的齿顶会将被切齿轮的齿根的渐开线齿廓切去了一部 分。避免根切的最小齿数，用标准齿条刀具切制标准齿轮时，因为 ，最少齿数为17。 三、机构的组成 构件指独立的运动单元，两个构件直接接触组成仍能产生某些相对运动的连接叫运动副。运动副按照相对运动的范围可以分为平面运动副和空间运动副;按运动副元素分为：低副-面接触、应力低;高副-点接触或线接触，应力高。其中运动副元素是只形成运动副的组建之间直接接触的部分。 四、机构自由度的计算 机构相对于机架所具有的独立运动的数目，叫机构的自由度。设一个平面机构由N个构件组成，其中必定有一个构件为机架，其活动构件数为n=N-1.设机构共有个低副、 个高副，因为在平面机构中每个低副和高副分别限制两个自由度和一个自由度，故平面机构的自由度为。在计算平面机构的自由度时，应该注意三种特殊情况：(1)复合铰链：三个或更多的构件在同一处联接成同轴线的两个或更多个转动副，就构成了复合铰链，计算自由度时应该按照两个或更多个运动副计算。(2)局部自由度：在有些机构中，为了其他一些非运动的原因，设置了附加机构，这种附加机构的运动是完全独立的，对整个 北航数理统计大作业-多元线性回归 应用数理统计多元线性回归分析 （第一次作业） 学院： 姓名： 学号： 2013年12月 交通运输业产值的多元线性回归分析 摘要：本文基于《中国统计年鉴》（2012年版）统计数据，寻找影响交通运输业发展的因素，包括工农业发展水平、能源生产水平、进出口贸易交流以及居民消费水平等，利用统计软件SPSS对各因素进行了筛选分析，采用逐步回归法得到最优多元线性回归模型，并对模型的回归显著性、拟合度以及随机误差的正态性进行了检验，最后可以利用有效的最优回归模型对将来进行预测。 关键字：多元线性回归，逐步回归，交通运输产值，工业产值，进出口总额1，引言 交通运输业指国民经济中专门从事运送货物和旅客的社会生产部门，包括铁路、公路、水运、航空等运输部门。它是国民经济的重要组成部分，是保证人们在政治、经济、文化、军事等方面联系交往的手段，也是衔接生产和消费的一个重要环节。交通运输业在现代社会的各个方面起着十分重要的作用，因此研究交通运输业发展水平与各个影响因素间的关系显得十分重要，建立有效的数学相关模型对于预测交通运输业的发展，制定相关政策方案提供依据。根据经验交通运输业的发展受到工农业发展、能源生产、进出口贸易以及居民消费水平等众因素的影响，故建立一个完整精确的数学模型在理论上基本无法实现，并且在实际运用中也没有必要，一种简单有效的方式就是寻找主要影响因素，分析其与指标变量的相关性，建立多元线性回归模型就是一种有效的方式。 变量与变量之间的关系分为确定性关系和非确定性关系，函数表达确定性关系。研究变量间的非确定性关系，构造变量间经验公式的数理统计方法称为 北航计算机控制系统大作业 ————————————————————————————————作者：————————————————————————————————日期: 计算机控制系统 大作业 姓名：陈启航 学号： 教师:周锐 日期:2016年6月1日 综合习题1 已知: 4 4)(+= s s D , 1） 试用 Z 变换、一阶向后差分、向前差分、零极点匹配、Tus ｔｉn 变换和预修正的Tus ｔin （设关键频率＝4）变换等方法将D （s)离散化，采样周期分别取为０.１s 和 0.4s; ２） 将 Ｄ（z ）的零极点标在Z 平面图上 3） 计算D （j ω)和各个Ｄ(e j ωT )的幅频和相频特性并绘图,ｗ由０~ 20r ad ，计算40 个点,应包括=4 点,每个T 绘一张图（Z 变换方法单画) 4） 计算 D(s)及T=0．1,T=０.4 时D （z ）的单位脉冲响应，运行时间为4 秒 5) 结合所得的结果讨论分析各种离散化方法的特点 6) 写出报告，附上结果。 解： （1） Z 变换法: a.离散化： T =0.1s 时， D (z )= 4z z ?0.6703； T =0.4s 时, D (z )= 4z z ?0.2019 ； ｂ.D (z )的零极点 c. D (jω)和D(e jωT )幅频相频特性曲线 连续系统: -1 -0.8-0.6-0.4-0.200.20.40.60.81 -1-0.8-0.6-0.4-0.200.20.40.60.81 零点 T=0.1s 时极点T=0.4s 时极点 T=0.1s时 T=0.4s时 《数值分析》计算实习题目 第一题： 1. 算法设计方案 （1）1λ，501λ和s λ的值。 1)首先通过幂法求出按模最大的特征值λt1，然后根据λt1进行原点平移求出另一特征值λt2，比较两值大小，数值小的为所求最小特征值λ1，数值大的为是所求最大特征值λ501。 2)使用反幂法求λs ，其中需要解线性方程组。因为A 为带状线性方程组，此处采用LU 分解法解带状方程组。 （2）与140k λλμλ-5011=+k 最接近的特征值λik 。 通过带有原点平移的反幂法求出与数k μ最接近的特征值 λik 。 （3）2cond(A)和det A 。 1）1=n λλ2cond(A)，其中1λ和n λ分别是按模最大和最小特征值。 2）利用步骤（1）中分解矩阵A 得出的LU 矩阵，L 为单位下三角阵,U 为上三角阵，其中U 矩阵的主对角线元素之积即为det A 。 由于A 的元素零元素较多，为节省储存量，将A 的元素存为6×501的数组中，程序中采用get_an_element()函数来从小数组中取出A 中的元素。 2.全部源程序 #include Harbin Institute of Technology 机械原理大作业二 课程名称：机械原理 设计题目：凸轮机构设计 姓名：李清蔚 学号：1140810304 班级：1408103 指导教师：林琳 一.设计题目 设计直动从动件盘形凸轮机构，其原始参数见表 1 表一：凸轮机构原始参数 升程（mm ) 升程 运动 角（o） 升程 运动 规律 升程 许用 压力 角（o） 回程 运动 角（o） 回程 运动 规律 回程 许用 压力 角（o） 远休 止角 （o） 近休 止角 （o） 40 90 等加 等减 速30 50 4-5-6- 7多 项式 60 100 120 二.凸轮推杆运动规律 （1）推程运动规律（等加速等减速运动） 推程F0=90° ①位移方程如下： ②速度方程如下： ③加速度方程如下： （2）回程运动规律（4-5-6-7多项式） 回程，F0=90°，F s=100°，F0’=50°其中回程过程的位移方程，速度方程，加速度方程如下： 三.运动线图及凸轮线图 本题目采用Matlab编程，写出凸轮每一段的运动方程，运用Matlab模拟将凸轮的运动曲线以及凸轮形状表现出来。代码见报告的结尾。 1、程序流程框图 开始 输入凸轮推程回 程的运动方程 输入凸轮基圆偏 距等基本参数 输出ds,dv,da图像 输出压力角、曲率半径图像 输出凸轮的构件形状 结束 2、运动规律ds图像如下： 速度规律dv图像如下： 加速度da规律如下图： 3.凸轮的基圆半径和偏距 以ds/dfψ-s图为基础，可分别作出三条限制线（推程许用压力角的切界限D t d t,回程许用压力角的限制线D t'd t'，起始点压力角许用线B0d'')，以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件。 得图如下：得最小基圆对应的坐标位置O点坐标大约为（13，-50）经计算取偏距e=13mm，r0=51.67mm. 航空科学与工程学院 《航空工程大型通用软件应用》大作业 机翼结构设计与分析 组号第3组 小组成员11051090 赵雅甜 11051093 廉佳 11051100 王守财 11051108 刘哲 11051135 张雄健 11051136 姜南 6月 目录 一 CATIA部分....................................... 错误!未定义书签。( 一) 作业要求..................................... 错误!未定义书签。( 二) 作业报告..................................... 错误!未定义书签。 1、三维模型图................................... 错误!未定义书签。 2、工程图....................................... 错误!未定义书签。 二 FLUENT部分...................................... 错误!未定义书签。( 一) 作业要求..................................... 错误!未定义书签。( 二) 作业报告..................................... 错误!未定义书签。 1、计算方法和流程............................... 错误!未定义书签。 2、网格分布图................................... 错误!未定义书签。 3、气动力系数................................... 错误!未定义书签。 4、翼型表面压力曲线............................. 错误!未定义书签。 5、翼型周围压力云图............................. 错误!未定义书签。 6、翼型周围x方向速度云图....................... 错误!未定义书签。 7、翼型周围y方向速度云图....................... 错误!未定义书签。 8、翼型周围x方向速度矢量图..................... 错误!未定义书签。 9、翼型周围y方向速度矢量图..................... 错误!未定义书签。 10、流线图...................................... 错误!未定义书签。 三 ANSYS部分....................................... 错误!未定义书签。( 一) 作业要求..................................... 错误!未定义书签。( 二) 作业报告..................................... 错误!未定义书签。 1、机翼按第一强度理论计算的应力云图............. 错误!未定义书签。 2、机翼按第二强度理论计算的应力云图............. 错误!未定义书签。 3、机翼按第三强度理论计算的应力云图............. 错误!未定义书签。 4、机翼按第四强度理论计算的应力云图............. 错误!未定义书签。 数值分析第三次大作业 一、算法的设计方案： （一）、总体方案设计： x y当作已知量代入题目给定的非线性方程组，求（1）解非线性方程组。将给定的(,) i i 得与(,)i i x y 相对应的数组t[i][j],u[i][j]。 （2）分片二次代数插值。通过分片二次代数插值运算，得到与数组t[11][21],u[11][21]]对应的数组z[11][21]，得到二元函数z=(,)i i f x y 。 （3）曲面拟合。利用x[i],y[j],z[11][21]建立二维函数表，再根据精度的要求选择适当k 值，并得到曲面拟合的系数矩阵C[r][s]。 （4）观察和(,)i i p x y 的逼近效果。观察逼近效果只需要重复上面（1）和（2）的过程，得到与新的插值节点(,)i i x y 对应的(,)i i f x y ，再与对应的(,)i i p x y 比较即可，这里求解 (,)i i p x y 可以直接使用（3）中的C[r][s]和k 。 （二）具体算法设计： （1）解非线性方程组 牛顿法解方程组()0F x =的解* x ，可采用如下算法： 1）在* x 附近选取(0) x D ∈，给定精度水平0ε>和最大迭代次数M 。 2）对于0,1, k M =执行 ① 计算() ()k F x 和()()k F x '。 ② 求解关于() k x ?的线性方程组 () ()()()()k k k F x x F x '?=- ③ 若() () k k x x ε∞∞ ?≤，则取*()k x x ≈，并停止计算；否则转④。 ④ 计算(1) ()()k k k x x x +=+?。 ⑤ 若k M <，则继续，否则，输出M 次迭代不成功的信息，并停止计算。 （2）分片双二次插值 给定已知数表以及需要插值的节点，进行分片二次插值的算法： 设已知数表中的点为： 00(0,1,,) (0,1,,)i j x x ih i n y y j j m τ=+=???=+=?? ，需要插值的节点为(,)x y 。 1) 根据(,)x y 选择插值节点(,)i j x y ： 若12h x x ≤+ 或12 n h x x ->-，插值节点对应取1i =或1i n =-， 数理统计第二次大作业材料行业股票的聚类分析与判别分析 2015年12月26日 材料行业股票的聚类分析与判别分析摘要 1 引言 2 数据采集及标准化处理 2.1 数据采集 本文选取的数据来自大智慧软件的股票基本资料分析数据，从材料行业的股票中选取了30支股票2015年1月至9月的7项财务指标作为分类的自变量，分别是每股收益（单位：元）、净资产收益率（单位：%）、每股经营现金流（单位：元）、主营业务收入同比增长率（单位：%）、净利润同比增长率（单位：%）、流通股本（单位：万股）、每股净资产（单位：元）。各变量的符号说明见表2.1，整理后的数据如表2.2。 表2.1 各变量的符号说明 自变量符号 每股收益（单位：元）X1 净资产收益率（单位：%）X2 每股经营现金流（单位：元）X3 主营业务收入同比增长率（单位：%）X4 净利润同比增长率（单位：%）X5 流通股本（单位：万股）X6 每股净资产（单位：元）X7 表2.2 30支股票的财务指标 股票代码X1 X2 X3 X4 X5 X6 X7 武钢股份600005-0.0990-2.81-0.0237-35.21-200.231009377.98 3.4444宝钢股份6000190.1400 1.980.9351-14.90-55.011642427.88 6.9197山东钢铁600022-0.11650.060.0938-20.5421.76643629.58 1.8734北方稀土6001110.0830 3.640.652218.33-24.02221920.48 2.2856 杭钢股份600126-0.4900-13.190.4184-36.59-8191.0283893.88 3.4497抚顺特钢6003990.219310.080.1703-14.26714.18112962.28 1.4667盛和资源6003920.0247 1.84-0.2141-5.96-19.3739150.00 1.2796宁夏建材6004490.04000.510.3795-22.15-92.3447818.108.7321宝钛股份600456-0.2090-2.53-0.3313-14.81-6070.2043026.578.1497山东药玻6005290.4404 5.26 1.2013 6.5016.7825738.018.5230国睿科技6005620.410011.53-0.2949 3.3018.9416817.86 3.6765海螺水泥600585 1.15169.05 1.1960-13.06-25.33399970.2612.9100华建集团6006290.224012.75-0.57877.90-6.4034799.98 1.8421福耀玻璃6006600.790014.250.9015 3.6017.27200298.63 6.2419宁波富邦600768-0.2200-35.02-0.5129 3.1217.8813374.720.5188马钢股份600808-0.3344-11.710.3939-21.85-689.22596775.12 2.6854亚泰集团6008810.02000.600.1400-23.63-68.16189473.21 4.5127博闻科技6008830.503516.71-0.1010-10.992612.8023608.80 3.0126新疆众和6008880.0523 1.04-0.910662.64162.0464122.59 5.0385西部黄金6010690.0969 3.940.115115.5125.5712600.00 2.4965中国铝业601600-0.0700-2.920.2066-9.0882.79958052.19 2.3811明泰铝业6016770.2688 4.66-1.09040.8227.8640770.247.4850金隅股份6019920.1989 3.390.3310-10.05-39.01311140.26 6.7772松发股份6032680.35007.00-0.3195-4.43-9.622200.00 6.0244方大集团0000550.0950 5.66-0.480939.2920.6742017.94 1.6961铜陵有色0006300.0200 1.220.6132 3.23-30.74956045.21 1.5443鞍钢股份000898-0.1230-1.870.7067-27.32-196.21614893.17 6.4932中钢国际0009280.572714.45-0.4048-14.33410.2441286.57 4.2449中材科技0020800.684610.27 1.219547.69282.1740000.00 6.8936中南重工0024450.1100 4.300.340518.8445.0950155.00 2.7030 2.2 数据的标准化处理 由于不同的变量之间存在着较大的数量级的差别，因此要对数据变量进行标准化处理。本文采用Z得分值法标准化的方法进行标准化，用x的值减去x的均值再除以样本的方差。也就是把个案转换为样本均值为0、标准差为1的样本。如果不同变量的变量值数值相差太大，会导致计算个案间距离时，由于绝对值较小的数值权数较小，个案距离的大小几乎由大数值决定，标准化过程可以解决此类问题，使不同变量的数值具有同等的重要性。经Z标准化输出结果见表 2.2。 表2.2 经Z标准化后的数据 ZX1ZX2ZX3ZX4ZX5ZX6ZX7 计算机控制系统 大作业 姓名：陈启航 学号： 教师：周锐 日期：2016年6月1日 综合习题1 已知: 4 4 )(+= s s D ， 1） 试用 Z 变换、一阶向后差分、向前差分、零极点匹配、Tustin 变换和 预修正的Tustin （设关键频率＝4）变换等方法将D(s)离散化，采样周期分别取为 和 ； 2） 将 D(z)的零极点标在Z 平面图上 3） 计算D （j ω）和各个D(e j ωT )的幅频和相频特性并绘图,w 由0~ 20ra d ，计算40 个点，应包括＝4 点，每个T 绘一张图（Z 变换方法单画） 4） 计算 D(s)及T=，T= 时D(z)的单位脉冲响应，运行时间为4 秒 5） 结合所得的结果讨论分析各种离散化方法的特点 6） 写出报告，附上结果。 解： （1） Z 变换法： a.离散化： T =0.1s 时， D (z )=4z z ?0.6703 ； T =0.4s 时， D (z )=4z z ?0.2019 ； b.D (z )的零极点 c. D (jω)和D(e jωT )幅频相频特性曲线 连续系统： T =0.1s 时 T =0.4s 时 d. D(s)和D(z)单位脉冲响应 D(s)单位脉冲响应： D(z)单位脉冲响应： T=0.1s时 T=0.4s时 （2）各种离散化方法： a.离散化后的D(z) 1、一阶向后差分： T=0.1s时 D(z)= 0.2857z z?0.7143 T=0.4s时 D(z)= 0.6154z z?0.3846 2、一阶向前差分：T=0.1s时 D(z)= 0.4 z?0.6 T=0.4s时 D(z)= 1.6 z+0.6 3、零极点匹配T=0.1s时 D(z)=0.1648(z+1) z?0.6703 T=0.4s时 D(z)=0.3991(z+1) z?0.2019 4、Tustin变换T=0.1s时 D(z)=0.1667(z+1) z?0.6667 T=0.4s时 D(z)= 0.4444(z+1) 5、预修正的Tustin变换（设关键频率=4） T=0.1s时 D(z)=0.1685(z+1) z?0.6629 T=0.4s时 D(z)=0.5073(z+1) z+0.0146 b.D(z)的零极点 1、一阶向后差分 材料学院研究生会 学术部 2011年12月 2007-2008学年第一学期期末试卷 一、（6分，A 班不做）设x 1，x 2，…，x n 是来自正态总体2(,)N μσ的样本，令 )x x T -= ， 试证明T 服从t -分布t （2） 二、（6分，B 班不做）统计量F-F(n,m)分布，证明 111(,)F F n m αααα-的（0<<1）的分位点x 是。 三、（8分）设总体X 的密度函数为 其中1α>-，是位置参数。x 1，x 2，…，x n 是来自总体X 的简单样本，试求参数α的矩估计和极大似然估计。 四、（12分）设总体X 的密度函数为 1x exp x (;) 0 , p x μμσσσ??-? -≥??? =????? ，其它， 其中,0,μμσσ-∞<<+∞>已知，是未知参数。x 1，x 2，…，x n 是来自总体X 的简单样本。 （1）试求参数σ的一致最小方差无偏估计σ∧ ； （2）σ∧ 是否为σ的有效估计？证明你的结论。 五、（6分，A 班不做）设x 1，x 2，…，x n 是来自正态总体211(,)N μσ的简单样本，y 1，y 2，…，y n 是来自正态总体222(,)N μσ的简单样本，且两样本相互独立，其中221122,,,μσμσ是未知参数，2212σσ≠。为检验假设012112:, :,H H μμμμ=≠可令12, 1,2,..., , ,i i i z x y i n μμμ=-==-则上述假设检验问题等价于0111:0, :0,H H μμ=≠这样双样本检验问题就变为单检验问题。基于变换后样本z 1，z 2，…，z n ，在显著性水平α下，试构造检验上述问题的t-检验统计量及相应的拒绝域。 六、（6分，B 班不做）设x 1，x 2，…，x n 是来自正态总体20(,)N μσ的简单样本，0μ已知，2σ未知，试求假设检验问题 22220010:, :H H σσσσ≥<的水平为α 的UMPT 。 七、（6分）根据大作业情况，试简述你在应用线性回归分析解决实际问题时应该注意哪些方面？ 八、（6分）设方差分析模型为 总离差平方和 试求A E(S )，并根据直观分析给出检验假设012:...0P H ααα====的拒绝域形式。 九、（8分）某个四因素二水平试验，除考察因子A 、B 、C 、D 外，还需考察A B ?，B C ?。今选用表78(2)L ，表头设计及试验数据如表所示。试用极差分析指出因子的主次顺序和较优工艺条件。 目标：使用带双步位移的QR 分解法求矩阵10*10[]ij A a =的全部特征值，并对其中的每一个实特征值求相应的特征向量。已知：sin(0.50.2)() 1.5cos( 1.2)(){i j i j ij i j i j a +≠+== (i,j=1,2, (10) 算法： 以上是程序运作的逻辑，其中具体的函数的算法，大部分都是数值分析课本上的逻辑，在这里特别写出矩阵A 的实特征值对应的一个特征向量的求法： ()[]()() []()[]()111111I 00000 i n n n B A I gause i n Q A I u Bu u λλ-?-?-=-?-?? ?-=????→=??????→= ?? ? 选主元的消元 检查知无重特征值 由于=0i A I λ- ，因此在经过选主元的高斯消元以后，i A I λ- 即B 的最后一行必然为零，左上方变 为n-1阶单位矩阵[]()()11I n n -?-，右上方变为n-1阶向量[]()11n Q ?-，然后令n u 1=-，则 ()1,2,,1j j u Q j n ==???-。 这样即求出所有A所有实特征值对应的一个特征向量。 #include 机械原理课程机构设计 实验报告 题目：建筑垃圾破碎机的设计与分析小组成员与学号： 班级： 第1页 建筑垃圾破碎机的设计与分析 摘要 本文简单介绍了建筑垃圾回收再利用的重要性，与工艺性，并自主设计了将颚式破碎机与反击式破碎机相结合的建筑垃圾破碎机。通过solidworks软件对设计机构进行建模，用adams进行仿真分析，验证所设计的机构均达到设计需要与可行性。 关键词：建筑垃圾破碎机、连杆机构、凸轮廓线设计 第2页 目录 1．机构的引出 (4) 1.1 建筑垃圾及其回收利用价值 (4) 1.2颚式破碎机和反击式破碎机各自的利弊分析 (4) 1.3设计新的建筑垃圾破碎机 (6) 2．机构的结构、功能介绍及建模 (7) 2.1 机构设计简图及各部分功能 (7) 2.2尺寸设计及建模 (8) 2.2.1主动轮和各从动轮的传动比 (8) 2.2.2凸轮廓线设计与挡板行程 ................................... 错误！未定义书签。 3．机构的仿真分析 (12) 3.1颚式破碎机的急回特性 (12) 3.2颚式破碎机的传动角验证 (14) 3.3停歇运动导杆机构所带动的下挡板往复运动的间歇性 (14) 4．总结 (17) 第3页 第4页 1． 机构的引出 1.1 建筑垃圾及其回收利用价值 二十一世纪是一个飞速发展的时代，随着城市人口的增加、新农村建设以及城市地铁的大规模扩建，建筑行业的新陈代谢全面加速，建筑垃圾的排放量也随之增加。然而，传统的方法处理建筑垃圾是将建筑垃圾运往乡村或郊外，露天堆放或掩埋。这样不仅破坏植被，降低土壤的生产能力，而且会让建筑垃圾中的有害物质渗入地下水层，污染环境，给人们的生活带来困扰。因此，如何实现建筑垃圾的高效、环保循环利用成为当今人们所面临的一个难题。 建筑垃圾的主要组成部分是废弃混凝土和砖块，而它们都是由水泥和天然砂石拌合而成的，这些都是砖块等建筑材料的重要组成部分。为了最大程度的利用建筑垃圾，首先应该解决的问题就是对其中的大块物料进行破碎，只有这样，破碎后的小快物料才能很好的还原天然砂石的性能，实现建筑垃圾的循环利用。 1.2颚式破碎机和反击式破碎机各自的利弊分析 目前应用较广的破碎机有颚式破碎机与反击式破碎机两种。 颚式破碎机的主体构造如图 1 图 1 颚式破碎机的主体构造 其工作原理为：轮①通过皮带和电机上的主动轮相连，①的转动带动杆②进而带动构件③的摆动（构件③的上端和机架铰接）。构件③通过摆动将体积较大 北航2010《应用数理统计》考试题及参考解答 09B 一、填空题（每小题3分，共15分） 1，设总体X 服从正态分布(0,4)N ，而12 15(,,)X X X 是来自X 的样本，则22 110 22 11152() X X U X X ++=++服从的分布是_______ . 解：(10,5)F ． 2，?n θ是总体未知参数θ的相合估计量的一个充分条件是_______ . 解：??lim (), lim Var()0n n n n E θθθ→∞ →∞ ==． 3，分布拟合检验方法有_______ 与____ ___. 解：2 χ检验、柯尔莫哥洛夫检验． 4，方差分析的目的是_______ . 解：推断各因素对试验结果影响是否显著． 5，多元线性回归模型=+Y βX ε中，β的最小二乘估计?β 的协方差矩阵?βCov()=_______ . 解：1?σ-'2Cov(β) =()X X ． 二、单项选择题（每小题3分，共15分） 1，设总体~(1,9)X N ，129(,, ,)X X X 是X 的样本，则___B___ . （A ） 1~(0,1)3X N -； （B ）1 ~(0,1)1X N -； （C ） 1 ~(0,1) 9X N -； （D ~(0,1)N ． 2，若总体2(,)X N μσ，其中2σ已知，当样本容量n 保持不变时，如果置信度1α-减小，则μ的 置信区间____B___ . （A ）长度变大； （B ）长度变小； （C ）长度不变； （D ）前述都有可能. 3，在假设检验中，就检验结果而言，以下说法正确的是____B___ . （A ）拒绝和接受原假设的理由都是充分的； （B ）拒绝原假设的理由是充分的，接受原假设的理由是不充分的； （C ）拒绝原假设的理由是不充分的，接受原假设的理由是充分的； （D ）拒绝和接受原假设的理由都是不充分的. 4，对于单因素试验方差分析的数学模型，设T S 为总离差平方和，e S 为误差平方和，A S 为效应平方和，则总有___A___ . 学号姓名成绩 《冲击式涡轮和反力式涡轮的设计计算》 总结：对冲击式涡轮和反力式涡轮进行设计计算，得到计算结果，具体见表1 和表2。 表1 反力式涡轮的计算结果 表2 冲击式涡轮的计算结果 根据计算结果，我们对比可以得到冲击式涡轮和反力式涡轮的相同点 是： 冲击式涡轮和反力式涡轮在计算功率时，均由泵的功率决定，由 T P N N =∑ 计算。 不同点具体见表3. 表3 反力式涡轮和冲击式涡轮的比较 1. 冲击式涡轮出口压力值取决于涡轮排气是直接排入周围环境还是导入辅助喷管，但两种情况下出口压强和反力式相比均很小。而反力式涡轮通常用于补燃式的液体火箭发动机中的涡轮泵中，所 以在不记喷注器压降的条件下，涡轮的出口压力等于燃烧室的压力。 2.在计算反力式涡轮的参数时，由于反力度容易确定，在分析过程 中广泛采用热力反力度。 反力式涡轮的设计计算 一.反力式涡轮参数的选择 在具有冲击式涡轮的供应系统(无补燃发动机系统)中，由燃气发生器产生的富燃燃气驱动涡轮，涡轮不冷却，富燃燃气的温度在1000～1200K 的范围内，比富氧燃气的允许温度(600～800K)高得多。另外，富燃燃气的气体常数比富氧燃气的气体常数大一些，这些都有利于减小需通过涡轮的燃气流量。 涡轮流量m t q 是具有冲击式涡轮的供应系统的主要参数之一。m t q 值越小，发动机的比冲就越高。涡轮流量m t q 可由泵和涡轮的功率平衡： T Pf Po N N N =+ 泵的需用功率降低，可减小通过涡轮的燃气流量，因此应尽量提高泵的效率。选定泵的结构并确定其效率后，可根据功率平衡求出所需的涡轮燃气流量，由此确定涡轮的效率。 涡轮入口压力(燃气发生器压力)取决于氧化剂泵的出口压力。当用燃料冷却推力室时，燃料泵出口压力比氧化剂泵的出口压力高。 涡轮出口压力之值取决于涡轮排气是直接排入周围环境还是导入辅助喷管。 冲击式涡轮计算的原始数据为： (1)涡轮的设计功率：涡轮功率T N 由泵所需的功率决定，由涡轮泵装置设计任务给定： 其中，T N —涡轮的设计功率，又称涡轮的轴功率； Pf N —燃料泵的轴功率； Pf N —氧化剂泵的轴功率。 (2)涡轮的设计角速度：涡轮的设计转速ω由泵不发生汽蚀时允许的最大角速度确定； (3)涡轮工质的物理常数和温度：涡轮进口总压*0P 、进口总温*0T 、和出口静压2P ；涡轮工质的绝热指数k 和气体常数R 。 二.反力式涡轮参数的选择机械原理大作业
专业点题北航机械原理
最新北航数理统计大作业-多元线性回归
北航计算机控制系统大作业
北航数值分析大作业第一题幂法与反幂法
哈工大机械原理大作业凸轮机构第四题
北航航空工程大型通用软件应用大作业样本
北航数值分析报告第三次大作业
北航数理统计第二次大作业-数据分析模板
北航计算机控制系统大作业
北航数理统计期末考试题
北航数值分析大作业第二题精解
北航七系机械学院机械原理大作业
北航应用数理统计考试题及参考解答
北航涡轮泵大作业