北航《操作系统》在线作业一

北航《操作系统》在线作业一
北航《操作系统》在线作业一

北航《操作系统》在线作业一

一、单选题(共15 道试题,共60 分。)

1. 银行家算法中的数据结构包括有可利用资源向量Available、最大需求矩阵Max、分配矩阵Allocation、需求矩阵Need,下列选项正确的是()。

A. Max[i,j]=Allocation[i,j]+Need[i,j]

B. Need[i,j]= Allocation[i,j]+ Max[i,j]

C. Max[i,j]= Available[i,j]+Need[i,j]

D. Need[i,j]= Available[i,j]+ Max[i,j]

正确答案:

2. 用户在一次计算过程中,或者一次事物处理中,要求计算机完成所做的工作的集合,这是指()。

A. 进程

B. 程序

C. 作业

D. 系统调用

正确答案:

3. 在最佳适应算法中是按()顺序形成空闲分区链。

A. 空闲区首址递增

B. 空闲区首址递减

C. 空闲区大小递增

D. 空闲区大小递减

正确答案:

4. ()不是基本的操作系统。

A. 批处理操作系统

B. 分时操作系统

C. 实时操作系统

D. 网络操作系统

正确答案:

5. 位示图方法可用于( )。

A. 盘空间的管理

B. 盘的驱动调度

C. 文件目录的查找

D. 页式虚拟存贮管理中的页面调度

正确答案:

6. 操作系统中的作业管理是一种()。

A. 宏观的高级管理

B. 宏观的低级管理

C. 系统刚开始加电

D. 初始化引导完成

正确答案:

7. 避免死锁的一个著名的算法是()。

A. 先入先出法

B. 银行家算法

C. 优先级算法

D. 资源按序分配法

正确答案:

8. 多道程序设计是指( )。

A. 在实时系统中并发运行多个程序

B. 在分布系统中同一时刻运行多个程序

C. 在一台处理机上同一时刻运行多个程序

D. 在一台处理机上并发运行多个程序

正确答案:

9. 在操作系统中,进程是一个()概念,而程序是一个静态的概念。

A. 运行态

B. 组合态

C. 动态

D. 等待态

正确答案:

10. 特别文件是与()有关的文件。

A. 文本

B. 图象

C. 硬件设备

D. 二进制数据

正确答案:

11. 作业与进程的主要区别为()。

A. 前者是由用户提交,后者是由系统自动生成

B. 两者执行不同的程序段

C. 前者以用户任务为单位,后者是操作系统控制的单位

D. 前者是批处理的,后者是分时处理的

正确答案:

12. 作业调度是从输入井中处于()状态的作业中选取作业调入主存运行。

A. 运行

B. 收容

C. 输入

D. 就绪

正确答案:

13. 在分页存储管理系统中,从页号到物理块号的地址映射是通过()实现的。

A. 段表

B. 页表

C. PCB

D. JCB

正确答案:

14. 通道是一种特殊的()。

A. I/O 设备

B. 设备控制器

C. 处理机

D. I/O控制器

正确答案:

15. 设两个进程共用一个临界资源的互斥信号量mutex,当mutex=1时表示()。

A. 一个进程进入了临界区,另一个进程等待

B. 没有一个进程进入临界区

C. 两个进程都进入了临界区

D. 两个进程都在等待

正确答案:

北航《操作系统》在线作业一

二、判断题(共10 道试题,共40 分。)

1. SPOOLing技术可以解决进程使用设备死锁问题。

A. 错误

B. 正确

正确答案:

2. 有了线程之后,程序只能以线程的身份运行。

A. 错误

B. 正确

正确答案:

3. 有关分页管理的叙述中,正确的是程序和数据是在开始执行前一次性装入的。

A. 错误

B. 正确

正确答案:

4. 虚存容量的扩大是以牺牲CPU工作时间以及内、外存交换时间为代价的。

A. 错误

B. 正确

正确答案:

5. 页式的地址是一维的,段式的地址是二维的。

A. 错误

B. 正确

正确答案:

6. 为了实现对临界区的共享,在每个进程中的临界区前面应设置V操作,在临界区之后应设置P操作。

A. 错误

B. 正确

正确答案:

7. 对临界资源应采取互斥访问方式来实现共享。

A. 错误

B. 正确

正确答案:

8. 文件系统中分配存储空间的基本单位是记录。

A. 错误

B. 正确

正确答案:

9. 在多道程序系统中,内存中作业的总数越多,系统利用率越高。

A. 错误

B. 正确

正确答案:

10. 单级目录结构能够解决文件重名问题。

A. 错误

B. 正确

正确答案:

惯性导航作业

惯性导航作业

一、数据说明: 1:惯导系统为指北方位的捷连系统。初始经度为116.344695283度、纬度为39.975172度,高度h为30米。初速度 v0=[-9.993908270;0.000000000;0.348994967]。 2:jlfw中为600秒的数据,陀螺仪和加速度计采样周期分别为为1/100秒和1/100秒。 3:初始姿态角为[2 1 90](俯仰,横滚,航向,单位为度),jlfw.mat中保存的为比力信息f_INSc(单位m/s^2)、陀螺仪角速率信息wib_INSc(单位rad/s),排列顺序为一~三行分别为X、Y、Z向信息. 4: 航向角以逆时针为正。 5:地球椭球长半径re=6378245;地球自转角速度wie=7.292115147e-5;重力加速度g=g0*(1+gk1*c33^2)*(1-2*h/re)/sqrt(1-gk2*c33^2); g0=9.7803267714;gk1=0.00193185138639;gk2=0.00669437999013;c33=sin(lat纬度); 二、作业要求: 1:可使用MATLAB语言编程,用MATLAB编程时可使用如下形式的语句读取数据:load D:\...文件路径...\jlfw,便可得到比力信息和陀螺仪角速率信息。用角增量法。 2:(1) 以系统经度为横轴,纬度为纵轴(单位均要转换为:度)做出系统位置曲线图; (2) 做出系统东向速度和北向速度随时间变化曲线图(速度单位:m/s,时间单位:s); (3) 分别做出系统姿态角随时间变化曲线图(俯仰,横滚,航向,单位转换为:度,时间单位:s); 以上结果均要附在作业报告中。 3:在作业报告中要写出“程序流程图、现阶段学习小结”,写明联系方式。

北航数值分析大作业一

《数值分析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 #include #include #include #include #include #include #define E 1.0e-12 /*定义全局变量相对误差限*/ int max2(int a,int b) /*求两个整型数最大值的子程序*/ { if(a>b) return a; else return b; } int min2(int a,int b) /*求两个整型数最小值的子程序*/ { if(a>b) return b; else return a; } int max3(int a,int b,int c) /*求三整型数最大值的子程序*/ { int t; if(a>b) t=a; else t=b; if(t

机械原理大作业

机械原理大作业 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'=。

北航惯性导航综合实验五实验报告

惯性导航技术综合实验 实验五惯性基组合导航及应用技术实验

惯性/卫星组合导航系统车载实验 一、实验目的 ①掌握捷联惯导/GPS组合导航系统的构成和基本工作原理; ②掌握采用卡尔曼滤波方法进行捷联惯导/GPS组合的基本原理; ③掌握捷联惯导 /GPS组合导航系统静态性能; ④掌握动态情况下捷联惯导 /GPS组合导航系统的性能。 二、实验内容 ①复习卡尔曼滤波的基本原理(参考《卡尔曼滤波与组合导航原理》第二、五章); ②复习捷联惯导/GPS组合导航系统的基本工作原理(参考以光衢编著的《惯性导航原理》第七章); 三、实验系统组成 ①捷联惯导/GPS组合导航实验系统一套; ②监控计算机一台。 ③差分 GPS接收机一套; ④实验车一辆; ⑤车载大理石平台; ⑥车载电源系统。 四、实验内容 1)实验准备 ①将IMU紧固在车载大理石减振平台上,确认IMU的安装基准面紧靠实验平台; ②将IMU与导航计算机、导航计算机与车载电源、导航计算机与监控计算

机、GPS 接收机与导航计算机、GPS 天线与GPS 接收机、GPS 接收机与GPS 电池之间的连接线正确连接; ③ 打开GPS 接收机电源,确认可以接收到4颗以上卫星; ④ 打开电源,启动实验系统。 2) 捷联惯导/GPS 组合导航实验 ① 进入捷联惯导初始对准状态,记录IMU 的原始输出,注意5分钟内严禁移动实验车和IMU ; ② 实验系统经过5分钟初始对准之后,进入导航状态; ③ 移动实验车,按设计实验路线行驶; ④ 利用监控计算机中的导航软件进行导航解算,并显示导航结果。 五、 实验结果及分析 (一) 理论推导捷联惯导短时段(1分钟)位置误差,并用1分钟惯导实验数据验证。 1、一分钟惯导位置误差理论推导: 短时段内(t<5min ),忽略地球自转0ie ω=,运动轨迹近似为平面1/0R =,此时的位置误差分析可简化为: (1) 加速度计零偏?引起的位置误差:2 10.88022t x δ?==m (2) 失准角0φ引起的误差:2 02 0.92182g t x φδ==m (3) 陀螺漂移ε引起的误差:3 30.01376 g t x εδ==m 可得1min 后的位置误差值123 1.8157m x x x x δδδδ=++= 2、一分钟惯导实验数据验证结果: (1)纯惯导解算1min 的位置及位置误差图:

专业点题北航机械原理

一、齿轮传动的基本概念 渐开线齿轮的啮合特点:(1)渐开线齿廓能够保证定传动比;(2)渐开线齿廓之间的正压力方向不变;(3)渐开线齿廓传动具有可分性。 齿轮机构的特点是:传动平稳、适用范围广、效率高、结构紧凑、工作可靠、寿命长。但制造和安装精度高、制造费用大,且不适合于距离较远的两轴之间的传动。齿轮传动可以用来传递任意轴间的运动和动力。 齿轮传动按照一对齿轮传递的相对运动分为平面齿轮传动和空间齿轮传动,平面齿轮传动又分为直齿圆柱齿轮传动、斜齿圆柱齿轮传动和人字齿轮传动;按照工作条件可以分为开式传动、半开式传动和闭式传动。 齿轮传动的基本要求是:传动准确、平稳;承载能力强。 二、齿轮传动的设计与计算 齿廓曲线与齿廓啮合基本定律:在啮合传动的任一瞬时,两轮齿廓曲线在相应接触点的功法线必须通过按给定传动比确定的该瞬时的节点。 渐开线齿轮啮合的正确条件:啮合轮齿的工作侧齿廓的啮合点必须总是在啮合线上,即两齿轮的模数和压力角应该分别相等。 齿轮传动的无侧隙啮合及标准齿轮的安装:一个齿轮节圆上的齿厚等于另一个齿轮节圆上的齿槽宽是无侧隙啮合的条件;外啮合齿轮的标准中心距为,内啮合是标准中心距为。

齿轮及其变位的相关计算:相关参数为齿数、模数、分度圆压力角、齿顶高系数和顶隙系数及标准直齿轮的几何尺寸计算,包括分度圆直径、齿顶高、齿根高、齿全高、齿顶圆直径、齿根圆直径、基圆直径、齿距、齿厚、齿槽宽、中心距、顶隙以及变位齿轮的变位系数等。 渐开线齿轮的根切现象:用展成法加工齿轮式,若刀具的齿顶线或齿顶圆与啮合线的焦点超过被切齿轮的极限点,则刀具的齿顶会将被切齿轮的齿根的渐开线齿廓切去了一部 分。避免根切的最小齿数,用标准齿条刀具切制标准齿轮时,因为 ,最少齿数为17。 三、机构的组成 构件指独立的运动单元,两个构件直接接触组成仍能产生某些相对运动的连接叫运动副。运动副按照相对运动的范围可以分为平面运动副和空间运动副;按运动副元素分为:低副-面接触、应力低;高副-点接触或线接触,应力高。其中运动副元素是只形成运动副的组建之间直接接触的部分。 四、机构自由度的计算 机构相对于机架所具有的独立运动的数目,叫机构的自由度。设一个平面机构由N个构件组成,其中必定有一个构件为机架,其活动构件数为n=N-1.设机构共有个低副、 个高副,因为在平面机构中每个低副和高副分别限制两个自由度和一个自由度,故平面机构的自由度为。在计算平面机构的自由度时,应该注意三种特殊情况:(1)复合铰链:三个或更多的构件在同一处联接成同轴线的两个或更多个转动副,就构成了复合铰链,计算自由度时应该按照两个或更多个运动副计算。(2)局部自由度:在有些机构中,为了其他一些非运动的原因,设置了附加机构,这种附加机构的运动是完全独立的,对整个

北航数值分析大作业第一题幂法与反幂法

《数值分析》计算实习题目 第一题: 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 #include void init_a();//初始化A double get_an_element(int,int);//取A 中的元素函数 double powermethod(double);//原点平移的幂法 double inversepowermethod(double);//原点平移的反幂法 int presolve(double);//三角LU 分解 int solve(double [],double []);//解方程组 int max(int,int); int min(int,int); double (*u)[502]=new double[502][502];//上三角U 数组 double (*l)[502]=new double[502][502];//单位下三角L 数组 double a[6][502];//矩阵A int main() { int i,k; double lambdat1,lambdat2,lambda1,lambda501,lambdas,mu[40],det;

北航惯性导航综合实验四实验报告

基于运动规划的惯性导航系统动态实验 二零一三年六月十日

实验4.1 惯性导航系统运动轨迹规划与设计实验 一、实验目的 为进行动态下简化惯性导航算法的实验研究,进行路径和运动状态规划,以验证不同运动状态下惯导系统的性能。通过实验掌握步进电机控制方法,并产生不同运动路径和运动状态。 二、实验内容 学习利用6045B 控制板对步进电机进行控制的方法,并控制电机使运动滑轨产生定长运动和不同加速度下的定长运动。 三、实验系统组成 USB_PCL6045B 控制板(评估板)、运动滑轨和控制计算机组成。 四、实验原理 IMU安装误差系数的计算方法 USB_PCL6045B 控制板采用了USB 串行总线接口通信方式,不必拆卸计算机箱就可以在台式机或笔记本电脑上进行运动控制芯片PCL6045B 的学习和评估。 USB_PCL6045B 评估板采用USB 串行总线方式实现评估板同计算机的数据交换,由评估板的FIFO 控制回路完成步进电机以及伺服电机的高速脉冲控制,任意2 轴的圆弧插补,2-4 轴的直线插补等运动控制功能。USB_PCL6045B 评估板上配置了全部PCL6045B 芯片的外部信号接口和增量编码器信号输入接口。由 USB_PCL6045B 评估测试软件可以进行PCL6045B 芯片的主要功能的评估测试。

图4-1-1USB_PCL6045B 评估板原理框图 如图4-1-1 所示,CN11 接口主要用于外部电源连接,可以选择DC5V 单一电源或DC5V/24V 电源。CN12 接口是USB 信号接口,用于USB_PCL6045B 评估板同计算机的数据交换。 USB_PCL6045B 评估板已经完成对PCL6045B 芯片的底层程序开发和硬件资源与端口的驱动,并封装成156 个API 接口函数。用户可直接在VC 环境下利用API 接口函数进行编程。 五、实验内容 1、操作步骤 1)检查电机驱动电源(24V) 2)检查USB_PCL6045B 控制板与上位机及电机驱动器间的连接电缆 3)启动USB_PCL6045B 控制板评估测试系统检查系统是否正常工作。 4)运行编写的定长运动程序,并比较实际位移与设定位移。

哈工大机械原理大作业凸轮机构第四题

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.

北航惯性导航大作业

惯性导航基础课程大作业报告(一)光纤陀螺误差建模与分析 班级:111514 姓名: 学号 2014年5月26日

一.系统误差原理图 二.系统误差的分析 (一)漂移引起的系统误差 1. εx ,εy ,εz 对东向速度误差δVx 的影响 clc;clear all; t=1:0.01:25; g=9.8; L=pi/180*39; Ws=2*pi/84.4*60; Wie=2*pi/24; R=g/(Ws)^2; e=0.1*180/pi; mcVx1=e*g*sin(L)/(Ws^2-Wie^2)*(sin(Wie*t)-Wie*sin(Ws*t)/Ws); mcVx2=e*((Ws^2-(Wie^2)*((cos(L))^2))/(Ws^2-Wie^2)*cos(Ws*t)-(Ws^2)*((sin(L))^2)*cos(Wi e*t)/(Ws^2-Wie^2)-(cos(L))^2); mcVx3=(sin(L))*(cos(L))*R*e*((Ws^2)*cos(Wie*t)/(Ws^2-Wie^2)-(Wie^2)*cos(Ws*t)/(Ws^2-Wi e^2)-1); plot(t,[mcVx1',mcVx2',mcVx3']); title('Ex,Ey,Ez 对Vx 的影响'); xlabel('时间t'); ylabel('Vx(t)'); 0,δλδL ,v v δδ

legend('Ex-mcVx1','Ey-mcVx2','Ez-mcVx3'); grid; axis square; 分析:εx,εy,εz对东向速度误差δVx均有地球自转周期的影响,εx,εy还会有舒勒周期分量的影响,其中,εy对δVx的影响较大。 2.εx,εy,εz对东向速度误差δVy的影响 clc;clear all; t=1:0.01:25; g=9.8; L=pi/180*39; Ws=2*pi/84.4*60; Wie=2*pi/24; R=g/(Ws)^2; e=0.1*180/pi; mcVy1=e*g*(cos(Wie*t)-cos(Ws*t))/(Ws^2-Wie^2); mcVy2=g*sin(L)*e/(Ws^2-Wie^2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t)); mcVy3=g*cos(L)*e/(Ws^2-Wie^2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t)); plot(t,[mcVy1',mcVy2',mcVy3']); title('Ex,Ey,Ez对Vy的影响'); xlabel('时间t'); ylabel('Vy(t)'); legend('Ex-mcVy1','Ey-mcVy2','Ez-mcVy3'); grid; axis square;

北航数值分析报告第三次大作业

数值分析第三次大作业 一、算法的设计方案: (一)、总体方案设计: 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 =-,

北航惯性导航综合实验四实验报告

基于运动规划的惯性导航系统动态实验 GAGGAGAGGAFFFFAFAF

二零一三年六月十日 实验4.1 惯性导航系统运动轨迹规划与设计实验一、实验目的 为进行动态下简化惯性导航算法的实验研究,进行路径和运动状态规划,以验证不同运动状态下惯导系统的性能。通过实验掌握步进电机控制方法,并产生不同运动路径和运动状态。 二、实验内容 学习利用6045B 控制板对步进电机进行控制的方法,并控制电机使运动滑轨产生定长运动和不同加速度下的定长运动。 三、实验系统组成 USB_PCL6045B 控制板(评估板)、运动滑轨和控制计算机组成。 四、实验原理 IMU安装误差系数的计算方法 GAGGAGAGGAFFFFAFAF

USB_PCL6045B 控制板采用了USB 串行总线接口通信方式,不必拆卸计算机箱就可以在台式机或笔记本电脑上进行运动控制芯片PCL6045B 的学习和评估。 USB_PCL6045B 评估板采用USB 串行总线方式实现评估板同计算机的数据交换,由评估板的FIFO 控制回路完成步进电机以及伺服电机的高速脉冲控制,任意 2 轴的圆弧插补,2-4 轴的直线插补等运动控制功能。USB_PCL6045B 评估板上配置了全部PCL6045B 芯片的外部信号接口和增量编码器信号输入接口。由 USB_PCL6045B 评估测试软件可以进行PCL6045B 芯片的主要功能的评估测试。 GAGGAGAGGAFFFFAFAF

图4-1-1USB_PCL6045B 评估板原理框图如图4-1-1 所示,CN11 接口主要用于外部电源连接,可以选择DC5V 单一电源或DC5V/24V 电源。CN12 接口是USB 信号接口,用于USB_PCL6045B 评估板同计算机的数据交换。 USB_PCL6045B 评估板已经完成对PCL6045B 芯片的底层程序开发和硬件资源与端口的驱动,并封装成156 个API 接口函数。用户可直接在VC 环境下利用API 接口函数进行编程。 五、实验内容 GAGGAGAGGAFFFFAFAF

北航数值分析大作业第二题精解

目标:使用带双步位移的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 #include #include #define N 10 #define E 1.0e-12 #define MAX 10000 //以下是符号函数 double sgn(double a) { double z; if(a>E) z=1; else z=-1; return z; } //以下是矩阵的拟三角分解 void nishangsanjiaodiv(double A[N][N]) { int i,j,k; int m=0; double d,c,h,t; double u[N],p[N],q[N],w[N]; for(i=0;i

北航惯性导航综合实验一实验报告

实 验一 陀螺仪关键参数测试与分析实验 加速度计关键参数测试与分析实验 二零一三年五月十二日 实验一陀螺仪关键参数测试与分析实验 一、实验目得 通过在速率转台上得测试实验,增强动手能力与对惯性测试设备得感性认识;通过对陀螺仪测试数据得分析,对陀螺漂移等参数得物理意义有清晰得认识,同时为在实际工程中应用陀螺仪与对陀螺仪进行误差建模与补偿奠定基础。 二、实验内容 利用单轴速率转台,进行陀螺仪标度因数测试、零偏测试、零偏重复性测试、零漂测试实验与陀螺仪标度因数与零偏建模、误差补偿实验。 三、实验系统组成 单轴速率转台、MEMS 陀螺仪(或光纤陀螺仪)、稳压电源、数据采集系统与分析系统。

四、实验原理 1.陀螺仪原理 陀螺仪就是角速率传感器,用来测量载体相对惯性空间得角速度,通常输出与角速率对应得电压信号。也有得陀螺输出频率信号(如激光陀螺)与数字信号(把模拟电压数字化)。以电压表示得陀螺输出信号可表示为: (1-1)式中就是与比力有关得陀螺输出误差项,反映了陀螺输出受比力得影响,本实验不考虑此项误差。因此,式(1-1)简化为 (1-2)由(1-2)式得陀螺输出值所对应得角速度测量值: (1-3) 对于数字输出得陀螺仪,传感器内部已经利用标度因数对陀螺仪模拟输出进行了量化,直接输出角速度值,即: (1-4)就是就是陀螺仪得零偏,物理意义就是输入角速度为零时,陀螺仪输出值所对应得角速度。且 (1-5) 精度受陀螺仪标度因数、随机漂移、陀螺输出信号得检测精度与得影响。通常与表现为有规律性,可通过建模与补偿方法消除,表现为随机特性,可通过信号滤波方法抵制。因此,准确标定与就是实现角速度准确测量得基础。 五、陀螺仪测试实验步骤 1)标度因数与零偏测试实验 a、接通电源,预热一定时间; b、陀螺工作稳定后,测量静止情况下陀螺输出并保存数据;

北航七系机械学院机械原理大作业

机械原理课程机构设计 实验报告 题目:建筑垃圾破碎机的设计与分析小组成员与学号: 班级: 第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 颚式破碎机的主体构造 其工作原理为:轮①通过皮带和电机上的主动轮相连,①的转动带动杆②进而带动构件③的摆动(构件③的上端和机架铰接)。构件③通过摆动将体积较大

(完整)北航惯性导航作业二.

(完整)北航惯性导航作业二. 编辑整理: 尊敬的读者朋友们: 这里是精品文档编辑中心,本文档内容是由我和我的同事精心编辑整理后发布的,发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方,但是任然希望((完整)北航惯性导航作业二.)的内容能够给您的工作和学习带来便利。同时也真诚的希望收到您的建议和反馈,这将是我们进步的源泉,前进的动力。 本文可编辑可修改,如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步,以下为(完整)北航惯性导航作业二.的全部内容。

惯性导航作业

一、数据说明: 1:惯导系统为指北方位的捷连系统.初始经度为116。344695283度、纬度为 39.975172度,高度h为30米。初速度v0=[—9。993908270;0.000000000; 0.348994967]。 2:jlfw中为600秒的数据,陀螺仪和加速度计采样周期分别为为1/100秒和1/100秒。 3:初始姿态角为[2 1 90](俯仰,横滚,航向,单位为度),jlfw。mat中保存的 为比力信息f_INSc(单位m/s^2)、陀螺仪角速率信息wib_INSc(单位rad/s),排列 顺序为一~三行分别为X、Y、Z向信息. 4:航向角以逆时针为正. 5:地球椭球长半径re=6378245;地球自转角速度wie=7。292115147e-5;重力加速度 g=g0*(1+gk1*c33^2)*(1-2*h/re)/sqrt(1—gk2*c33^2);g0=9.7803267714; gk1=0。00193185138639;gk2=0。00669437999013;c33=sin(lat纬度); 二、作业要求: 1:可使用 MATLAB语言编程,用MATLAB编程时可使用如下形式的语句读取数据: load D:\..。文件路径。。.\jlfw,便可得到比力信息和陀螺仪角速率信息。用角增 量法。 2:(1) 以系统经度为横轴,纬度为纵轴(单位均要转换为:度)做出系统位置曲线图; (2)做出系统东向速度和北向速度随时间变化曲线图(速度单位:m/s,时间单位:s); (3) 分别做出系统姿态角随时间变化曲线图(俯仰,横滚,航向,单位转换为:度,时间单位:s); 以上结果均要附在作业报告中. 3:在作业报告中要写出“程序流程图、现阶段学习小结”,写明联系方式。 (注意程序流程图不是课本上的惯导解算流程,而是你程序分为哪几个模块、是 怎样一步步执行的,什么位置循环等,让别人根据该流程图能够编出相应程序)

北航数值分析大作业第二题

数值分析第二次大作业 史立峰 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。

哈工大机械原理大作业凸轮

机械原理大作业二 课程名称: _______ 设计题目: 凸轮机构设计 院 系: ------------------------- 班 级: _________________________ 设计者: ________________________ 学 号: _________________________ 指导教师: ______________________ 哈尔滨工业大学 Harbin I nstituteof Techndogy

设计题目 如右图所示直动从动件盘形凸轮机构,选择一组凸轮机构的原始参数, 据此设计该凸轮机构。 凸轮机构原始参数 二.凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图 凸轮推杆升程运动方程:冷3唱—亦(中] 156 12 .. v 」1 - cos()] 兀1 5 374.4 2 12 ? a 1si n( ) 兀 1 5 % t 表示转角, s 表示位移 t=0:0.01:5*pi/6; %升程阶段 s= [(6*t)/(5*pi)- 1/(2*pi)*si n(12*t/5)]*130; hold on plot(t,s);

t= 5*pi/6:0.01:pi; %远休止阶段 s=130; hold on plot(t,s); t=pi:0.01:14*pi/9; %回程阶段 s=65*[1+cos(9*(t-pi)/5)]; hold on plot(t,s); t=14*pi/9:0.01:2*pi; %近休止阶段 s=0; hold on plot(t,s); grid on % t表示转角,令3 1=1 t=0:0.01:5*pi/6; %升程阶段v=156*1*[1-cos(12*t/5)]/pi hold on plot(t,v); t= 5*pi/6:0.01:pi; %远休止阶段

北航数值分析第二次大作业--QR分解

《数值分析A》

一、算法设计方案 整个程序主要分为四个函数,主函数,拟上三角化函数,QR分解函数以及使用双步位移求解矩阵特征值、特征向量的函数。因为在最后一个函数中也存在QR分解,所以我没有采用参考书上把矩阵M进行的QR分解与矩阵Ak的迭代合并的方法,而是在该函数中调用了QR分解函数,这样增强了代码的复用性,减少了程序长度;但由于时间关系,对阵中方法的运算速度没有进行深入研究。 1.为了减少QR分解法应用时的迭代次数,首先对给定矩阵进行拟上三角化处理。 2.对经过拟上三角化处理的矩阵进行QR分解。 3.注意到计算特征值与特征向量的过程首先要应用前面两个函数,于是在拟上三角化矩阵的基础上对QR分解函数进行了调用。计算过程中,没有采用goto语句,而是根据流程图采用其他循环方式完成了设计,通过对迭代过程的合并,简化了程序的循环次数,最后在计算特征向量的时候采用了列主元高斯消去法。

二、源程序代码 #include #include #include int i,j,k,l,m; //定义外部变量double d,h,b,c,t,s; double A[10][10],AA[10][10],R[10][10],Q[10][10],RQ[10][10]; double X[10][10],Y[10][10],Qt[10][10],M[10][10]; double U[10],P[10],T[10],W[10],Re[10]={0},Im[10]={0}; double epsilon=1e-12; void main() { void Quasiuppertriangular(double A[][10]); void QRdecomposition(double A[][10]); void DoublestepsQR(double A[][10]); int i,j; for(i=0;i<10;i++) { for(j=0;j<10;j++) { A[i][j]=sin(0.5*(i+1)+0.2*(j+1)); Q[i][j]=0; AA[i][j]=A[i][j]; } A[i][i]=1.5*cos(2.2*(i+1)); AA[i][i]=A[i][i];

北航数值分析报告大作业第八题

北京航空航天大学 数值分析大作业八 学院名称自动化 专业方向控制工程 学号 学生姓名许阳 教师孙玉泉 日期2014 年11月26 日

一.题目 关于x , y , t , u , v , w 的方程组(A.3) ???? ?? ?=-+++=-+++=-+++=-+++79 .0sin 5.074.3cos 5.007.1cos sin 5.067.2cos 5.0y w v u t x w v u t y w v u t x w v u t (A.3) 以及关于z , t , u 的二维数表(见表A-1)确定了一个二元函数z =f (x , y )。 表A-1 二维数表 t z u 0 0.4 0.8 1.2 1.6 2 0 -0.5 -0.34 0.14 0.94 2.06 3.5 0.2 -0.42 -0.5 -0.26 0.3 1.18 2.38 0.4 -0.18 -0.5 -0.5 -0.18 0.46 1.42 0.6 0.22 -0.34 -0.58 -0.5 -0.1 0.62 0.8 0.78 -0.02 -0.5 -0.66 -0.5 -0.02 1.0 1.5 0.46 -0.26 -0.66 -0.74 -0.5 1. 试用数值方法求出f (x , y ) 在区域}5.15.0,8.00|), {≤≤≤≤=y x y x D (上的近似表达式 ∑∑===k i k j s r rs y x c y x p 00 ),( 要求p (x , y )以最小的k 值达到以下的精度 ∑∑==-≤-=10020 7210)],(),([i j i i i i y x p y x f σ 其中j y i x i i 05.05.0,08.0+==。 2. 计算),(),,(* ***j i j i y x p y x f (i =1,2,…,8 ; j =1,2,…,5) 的值,以观察p (x , y ) 逼 近f (x , y )的效果,其中j y i x j i 2.05.0,1.0**+==。

相关文档
最新文档