北航《先进加工技术及装备》在线作业一
一、单选题（共10 道试题，共30 分。）
1. 在数控机床CNC系统组成中，核心部分为()。
A. 输入输出设备
B. PLC
C. CNC装置
D. 伺服机构
正确答案：
2. 步行机器人的行走机构多为()。
A. 滚轮
B. 履带
C. 连杆机构
D. 齿轮机构
正确答案：
3. 以下不是数控机床伺服系统的形式是()。
A. 开环控制
B. 全闭环控制
C. 半闭环控制
D. 循环控制
正确答案：
4. 用电火花加工冲模时，若火花间隙能保证配合间隙的要求，应选用的工艺方法是()。
A. 直接配合法
B. 修配冲头法
C. 修配电极法
D. 阶梯电极法
正确答案：
5. 在CAD/CAM系统中，()是连接CAD、CAM的纽带。
A. CAE
B. CAG
C. CAPP
D. CAQ
正确答案：
6. 工业机器人的额定负载是指在规定范围内()所能承受的最大负载允许值
A. 手腕机械接口处
B. 手臂
C. 末端执行器
D. 机座
正确答案：
7. CIMS所研究的对象是()。
A. 产品
B. 零件
C. 典型样件
D. 一个车间或一个企业
正确答案：
8. 高功率密度的电子束加工适于()。
A. 钢板上打孔
B. 工件表面合金化
C. 电子束曝光
D. 工件高频热处理
正确答案：
9. FMC是()的简称。
A. 加工中心
B. 计算机控制系统
C. 永磁式伺服系统
D. 柔性制造单元。
正确答案：
10. 粗加工后安排调质处理可提高零件的()。
A. 硬度
B. 耐磨性
C. 强度
D. 综合力学性能
正确答案：
北航《先进加工技术及装备》在线作业一
二、多选题（共10 道试题，共40 分。）
1. CIM技术是传统的制造技术与()系统工程技术等的有机结合。
A. 现代信息技术
B. 管理技术
C. 网络技术
D. 自动化技术、
正确答案：
2. 快速成形技术（RP）最早产生于二十世纪70年代末到80年代初，目前RP技术的主流有()技术。
A. SLA立体光造型
B. LOM薄材叠层制造
C. SLS激光烧结
D. FDM熔融成型
正确答案：
3. 成功企业的关键因素主要有：上市时间、()、灵活性、环境。
A. 产品销售价格
B. 质量
C. 产品成本
D. 服务
正确答案：
4. 数控程序的编制方法有()。
A. 手工编程
B. 自动编程
C. CAD/CAM软件编程
正确答案：
5. CIMS是一种基于CIM哲理构成的计算机化()的制造系统。
A. 信息化
B. 网络化
C. 智能化
D. 集成优化。
正确答案：
6. 在电火花加工中，提高电蚀量和加工效率的电参数途径有()。
A. 提高脉冲频率
B. 增加单个脉冲能量
C. 减少脉冲间隔
D. 增大脉冲间隔
正确答案：
7. 电火花加工的表面质量主要是指被加工零件的()。
A. 表面粗糙度
B. 表面变质层
C. 表面力学性能
D. 表面导电性
正确答案：
8. 典型FMS的子系统有()。
A. 加工系统
B. 运储系统
C. 计算机控制系统
D. 售后服务
正确答案：
9. 工业机器人由()组成。
A. 伺服装置
B. 控制装置
C. 驱动装置
D. 操作机
E. 传感机
正确答案：
10. FMS由()等于系统组成。
A. 物流系统
B. 信息流系统
C. 能源流系统
D. 加工系统
正确答案：
北航《先进加工技术及装备》在线作业一
三、判断题（共10 道试题，共30 分。）
1. 3B代码编程法是最先进的电火花线切割编程方法。
A. 错误
B. 正确
正确答案：
2. 弛张式脉冲电源电能利用率相当高，所以在电火花加工中应用较多。
A. 错误
B. 正确
正确答案：
3. 阶梯形变幅杆振幅放大倍数最高，但受负载阻力时振幅衰减严重，且容易产生应力集中。
A. 错误
B. 正确
正确答案：
4. 法拉第电解定律认为电解加工时电极上溶解或析出物质的量与通过的电量成正比，它也适用于电镀。
A. 错误
B. 正确
正确答案：
5. 离子束加工必须在真空条件下进行。
A. 错误
B. 正确
正确答案：
6. 电子束加工必须在真空条件下进行。
A. 错误
B. 正确
正确答案：
7. 由于离子的质量远大于电子，故离子束加工的生产率比电子束高，但表面粗糙度稍差。
A. 错误
B. 正确
正确答案：
8. 在采取适当的工艺保证后，数控线切割也可以加工盲孔。
A. 错误
B. 正确
正确答案：
9. 当电极丝的进给速度明显超过蚀除速度，则放电间隙会越来越小，以致产生短路。
A. 错误
B. 正确
正确答案：
10. 通常慢走丝线切割加工中广泛使用直径为0.1 mm以上的黄铜丝作为电极丝。
A. 错误
B. 正确
正确答案：
《数值分析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 北航最优化方法大作业参考 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 北航《国际金融》在线作业一 单选题多选题判断题 一、单选题（共 25 道试题，共 75 分。） 1. 汇率作为外汇的（），是重要的经济变量。 A. 价值 B. 价格 C. 变动率 D. 收益率 -----------------选择：B 2. 布雷顿森林协议在对会员国缴纳的资金来源中规定( ) A. 会员国份额必须全部以黄金或可兑换的黄金货币缴纳 B. 会员国份额的50％以黄金或可兑换成黄金的货币缴纳，其余50％的份额以本国货币缴纳 C. 会员国份额的25％以黄金或可兑换成黄金的货币缴纳，其余75％的份额以本国货币缴纳 D. 会员国份额必须全部以本国货币缴纳 -----------------选择：C 3. 进口商与银行订立远期外汇合同，是为了 ( ) A. 防止因外汇汇率上涨而造成的损失 B. 防止因外汇汇率下跌而造成的损失 C. 获得因外汇汇率上涨而带来的收益 D. 获得因外汇汇率下跌而带来的收益 -----------------选择：B 4. 如果一种货币的远期汇率低于即期汇率，称之为( ) A. 升水 B. 贴水 C. 隔水 D. 平水 -----------------选择：B 5. ( ) ，又称为英美式外汇市场，是指无具体交易场所的外汇市场。 A. 有形市场 B. 无形市场 C. 官方外汇市场 D. 平行外汇市场 -----------------选择：B 6. 在浮动汇率制度下，如果外汇市场上的外国货币求过于供时，则( _) A. 外币价格上涨，外汇汇率上升 B. 外币价格下跌，外汇汇率下降 C. 外币价格上涨，外汇汇率下降 D. 外币价格下跌，外汇汇率上升 -----------------选择：A 7. SDR s是( ) A. 欧洲经济货币联盟创设的货币 应用数理统计第一次大作业 学号： 姓名： 班级： 2013年12月 国家财政收入的多元线性回归模型 摘 要 本文以多元线性回归为出发点，选取我国自1990至2008年连续19年的财政收入为因变量，初步选取了7个影响因素,并利用统计软件PASW Statistics 17.0对各影响因素进行了筛选，最终确定了能反映财政收入与各因素之间关系的“最优”回归方程： 46?578.4790.1990.733y x x =++ 从而得出了结论，最后我们用2009年的数据进行了验证，得出的结果在误差范围内，表明这个模型可以正确反映影响财政收入的各因素的情况。 关键词：多元线性回归，逐步回归法，财政收入，SPSS 0符号说明 变 量 符号 财政收入 Y 工 业 X 1 农 业 X 2 受灾面积 X 3 建 筑 业 X 4 人 口 X 5 商品销售额 X 6 进出口总额X7 1 引言 中国作为世界第一大发展中国家，要实现中华民族的伟大复兴，必须把发展放在第一位。近年来，随着国家经济水平的飞速进步，人民生活水平日益提高，综合国力日渐强大。经济上的飞速发展并带动了国家财政收入的飞速增加，国家财政的状况对整个社会的发展影响巨大。政府有了强有力的财政保证才能够对全局进行把握和调控，对于整个国家和社会的健康快速发展有着重要的意义。所以对国家财政的收入状况进行研究是十分必要的。 国家财政收入的增长，宏观上必然与整个国家的经济有着必然的关系，但是具体到各个方面的影响因素又有着十分复杂的相关原因。为了研究影响国家财政收入的因素，我们就很有必要对其财政收入和影响财政收入的因素作必要的认识，如果能对他们之间的关系作一下回归，并利用我们所知道的数据建立起回归模型这对我们很有作用。而影响财政收入的因素有很多，如人口状况、引进的外资总额，第一产业的发展情况，第二产业的发展情况，第三产业的发展情况等等。本文从国家统计信息网上选取了1990-2009年这20年间的年度财政收入及主要影响因素的数据，包括工业，农业，建筑业，批发和零售贸易餐饮业，人口总数等。文中主要应用逐步回归的统计方法，对数据进行分析处理，最终得出能够反映各个因素对财政收入影响的最“优”模型。 2解决问题的方法和计算结果 2.1 样本数据的选取与整理 本文在进行统计时，查阅《中国统计年鉴2010》中收录的1990年至2009年连续20年的全国财政收入为因变量，考虑一些与能源消耗关系密切并且直观上 北航计算机控制系统大作业 ————————————————————————————————作者：————————————————————————————————日期: 计算机控制系统 大作业 姓名：陈启航 学号： 教师:周锐 日期: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 北航《刑事诉讼法》在线作业一单选题多选题 一、单选题（共15道试题，共60分。） 1.在法庭审理过程中，有权申请通知新的证人到庭、调取新的物证，申请重新鉴定或者勘验的人是（） A.被告人 B.被告人及其辩护人 C..当事人 D.当事人和辩护人 ----------------- 选择：D 2.乙市M区人民法院对孙某盗窃罪和抢劫罪作出判决后，人民检察院不抗诉，但孙某不服提岀上诉。市中级人民法院审理后认为，一审判决事实淸楚，证据确实、充分，但疑刑不当。一审对孙某盗窃罪和抢劫罪分别判处2年和9年有期徒刑，决立执行的刑期为10年，而两罪准确量刑分别为5年和7年。二审法院应如何处理？（） A.直接改判两罪刑罚分別为5年和7年有期徒刑，并在7年以上12年以下决立应当执行的刑期 B.直接改判两罪刑罚分别为5年和7年有期徒刑，并决左应当执行的刑期为10年 C.维持一审判决 D.维持一审判处的2年有期徒刑，将一审判处的9年有期徒刑改为7年，并在7年以上9 年以下决定应当执行的刑期 ----------------- 选择：D 3.二审人民法院审理人民检察院提出抗诉的案件，应在开庭多长时间以前通知人民检察院查阅案卷？（） A.7日 B.3日 C.5日 D.10 日 ----------------- 选择：D 4.人民检察院审查公安机关移送起诉的案件，犯罪嫌疑人有刑事诉讼法第十五条规左情形之一的，应当作出（） A..起诉决宦 B..不起诉决泄 C..撤销案件决定 D.?退回公安机关补充侦查的决立 ----------------- 选择：B 5.根据我国刑事诉讼法的规左，下列人员中不属于我国刑事诉讼中的当事人的是：（） A.证人 B.犯罪嫌疑人和被告人 C.被害人 D.附带民事诉讼的原告人 ----------------- 选择：A 6.某市中级人民法院在审理自诉案件原告李某提岀上诉的案件时，该自诉案件的一审被告人张某对李某提出反诉。市中级法院对该反诉应当按照下列哪种方式处理？（） 数理统计第二次大作业材料行业股票的聚类分析与判别分析 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 航空科学与工程学院 《航空工程大型通用软件应用》大作业 机翼结构设计与分析 组号第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 =-， 计算机控制系统 大作业 姓名：陈启航 学号： 教师：周锐 日期：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、一阶向后差分 数值分析第二次大作业 史立峰 SY1505327 一、 方案 （1）利用循环结构将sin(0.50.2)() 1.5cos( 1.2)() {i j i j ij i j i j a +≠+==(i,j=1,2,……,10)进行赋值，得到需要变换的 矩阵A ； （2）然后，对矩阵A 利用Householder 矩阵进行相似变换，把A 化为上三角矩阵A (n-1)。 对A 拟上三角化，得到拟上三角矩阵A (n-1)，具体算法如下： 记A(1)=A ，并记A(r)的第r 列至第n 列的元素为()n r r j n i a r ij ,,1,;,,2,1) ( +==。 对于2,,2,1-=n r 执行 1. 若 ()n r r i a r ir ,,3,2) ( ++=全为零，则令A(r+1) =A(r),转5；否则转2。 2. 计算 () ∑+== n r i r ir r a d 1 2 )( ()( )r r r r r r r r r r d c a d a c ==-=++则取,0sgn ) (,1)(,1若 )(,12r r r r r r a c c h +-= 3. 令 () n T r nr r r r r r r r r R a a c a u ∈-=++) ()(,2)(,1,,,,0,,0 。 4. 计算 r r T r r h u A p /)(= r r r r h u A q /)(= r r T r r h u p t /= r r r r u t q -=ω T r r T r r r r p u u A A --=+ω)()1( 5. 继续。 （3）使用带双步位移的QR 方法计算矩阵A (n-1)的全部特征值，也是A 的全部特征值，具体算法如下： 1. 给定精度水平0>ε和迭代最大次数L 。 2. 记n n ij n a A A ?-==][) 1()1()1(，令n m k ==,1。 目标：使用带双步位移的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) 逐步回归分析 学院：机械工程学院 专业：材料加工工程 日期：2014年12月7日 摘要：本文介绍多元线性回归分析方法以及逐步回归法，然后结合实际，以我国1995-2012年的财政收入为因变量，选取了8个可能的影响因素，选用逐步回归法对各影响因素进行了筛选分析，最终确定了其“最优”回归方程。 关键字：多元线性回归 逐步回归法 财政收入 SPSS 1 引言 自然界中任何事物都是普遍联系的，客观事物之间往往都存在着某种程度的关联关系。为了研究变量之间的相关关系，人们常用回归分析的方法，而回归分析是数理统计中一种常用方法。数理统计作为一种实用有效的工具，广泛应用于国民经济的各个方面，在解决实际问题中发挥了巨大的作用，是一种理论联系实践、指导实践的科学方法。 财政收入，是指政府为履行其职能、实施公共政策和提供公共物品与服务需要而筹集的一切资金的总和。财政收入表现为政府部门在一定时期内（一般为一个财政年度）所取得的货币收入。财政收入是衡量一国政府财力的重要指标，政府在社会经济活动中提供公共物品和服务的范围和数量，在很大程度上决定于财政收入的充裕状况。 本文将以回归分析为方法，运用数理统计工具探求财政收入与各种统计指标之间的关系，总结主要影响因素，并对其作用、前景进行分析和展望。 2 多元线性回归 2.1 多元线性回归简介 在实际问题中，某一因素的变化往往受到许多因素的影响，多元回归分析的任务就是要找出这些因素之间的某种联系。由于许多非线性的情形都可以通过变换转化为线性回归来处理，因此，一般的实际问题都是基于多元线性回归问题进行处理的。对多元线性回归模型简要介绍如下： 如果随机变量y 与m )2(≥m 个普通变量m x x x 21,有关，且满足关系式： εββββ++++=m m x x x y 22110 2,0σεε==D E （2.1） 其中，2210,,,σββββm 是与m x x x 21,无关的未知参数，ε是不可观测的随机变量，),0(~2N I N σε。 学号姓名成绩 《冲击式涡轮和反力式涡轮的设计计算》 总结：对冲击式涡轮和反力式涡轮进行设计计算，得到计算结果，具体见表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 。 二.反力式涡轮参数的选择 北航《人力资源管理》在线作业一 单选题多选题判断题 一、单选题（共 10 道试题，共 40 分。） 1. 人与职业相匹配的职业选择理论是由( )提出的。 A. 美国波士顿大学教授帕森斯 B. 美国约翰霍普金斯大学教授霍兰德 C. 中国人民大学萧鸣政教授 D. 中国清华大学赵平教授 -----------------选择：A 2. 工作分析中方法分析常用的方法是( ) A. 关键事件技术 B. 职能工作分析 C. 问题分析 D. 流程图 -----------------选择：C 3. 认为人力资源是具有劳动能力的全部人口的观点属于( )。 A. 人口观 B. 在岗人口观 C. 人员素质观 D. 人员数量观 -----------------选择：A 4. 目前世界上大多数国家实行的养老保险类型是( )。 A. 强制储蓄型养老保险 B. 国家统筹型社会保险 C. 投保资助型社会保险 D. 子女抚养型养老 -----------------选择：C 5. 劳动争议仲裁委员会对每一起劳动争议实行( )，即行终结的法律制度。 A. 一次裁决 B. 二次裁决 C. 三次裁决 D. 四次裁决 -----------------选择：A 6. （）成为人力资源配置的核心。 A. 知识 B. 市场 C. 人才 D. 竞争 -----------------选择：A 7. 把“员工视为活动主体、公司主人”是哪种人力资源管理模式( )。 A. 自我中心式、非理性化家族管理 B. 以人为中心、非理性化家族管理 C. 以人为中心、理性化团队管理 D. 自我中心式、理性化团队管理 -----------------选择：C 8. （）是建立考核制度和实施考核工作的前提。 A. 调查分析 B. 公平 C. 领导重视 D. 科学的方法 -----------------选择：B 9. 以人性为核心的人本管理主体是何种要素( )。 A. 职工 B. 环境 C. 文化 D. 价值观 -----------------选择：A 10. 现代意义上的人事管理是随着（）的产生而发展起来的。 A. 工厂制度 B. 资本主义制度 C. 工会制度 D. 工业革命 -----------------选择：D 北航《人力资源管理》在线作业一 单选题多选题判断题 二、多选题（共 10 道试题，共 40 分。） 1. 最近几年提出的，越来越多的人把人力看作是人员素质综合发挥的作用力。认为人力资源是劳动生产过程中，可以直接投入的( )总和 A. 体力 B. 脑力 C. 心力 D. 能力 -----------------选择：ABC 2. 人力资源开发与管理的基本环境包括（）。 A. 政治因素 B. 经济因素 C. 社会因素 D. 技术因素 E. 环境因素 -----------------选择：ABCD 3. 人力资源思想演进有四个阶段（）。 A. 传统劳动管理时期 《数值分析A》 一、算法设计方案 整个程序主要分为四个函数，主函数，拟上三角化函数，QR分解函数以及使用双步位移求解矩阵特征值、特征向量的函数。因为在最后一个函数中也存在QR分解，所以我没有采用参考书上把矩阵M进行的QR分解与矩阵Ak的迭代合并的方法，而是在该函数中调用了QR分解函数，这样增强了代码的复用性，减少了程序长度；但由于时间关系，对阵中方法的运算速度没有进行深入研究。 1.为了减少QR分解法应用时的迭代次数，首先对给定矩阵进行拟上三角化处理。 2.对经过拟上三角化处理的矩阵进行QR分解。 3.注意到计算特征值与特征向量的过程首先要应用前面两个函数，于是在拟上三角化矩阵的基础上对QR分解函数进行了调用。计算过程中，没有采用goto语句，而是根据流程图采用其他循环方式完成了设计，通过对迭代过程的合并，简化了程序的循环次数，最后在计算特征向量的时候采用了列主元高斯消去法。 二、源程序代码 #include北航最优化方法大作业参考
15秋北航国际金融在线作业一答案满分
北航数理统计回归分析大作业
北航计算机控制系统大作业
北航数值分析大作业第一题幂法与反幂法
免费15秋北航《刑事诉讼法》在线作业一答案满分
北航数理统计第二次大作业-数据分析模板
北航航空工程大型通用软件应用大作业样本
北航数值分析报告第三次大作业
北航计算机控制系统大作业
北航数值分析大作业第二题
北航数值分析大作业第二题精解
北航数理统计回归分析大作业
北航涡轮泵大作业
15秋北航《人力资源管理》在线作业一答案满分
北航数值分析第二次大作业--QR分解
相关主题
文本预览