版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、. -. -可修编- . z.老年人跌倒检测技术研究Research on Fall Detection Technology for the Elderly摘要跌倒重威胁着老年人的安康和生命,提供自由、实时的平安监护对老年人的生活质量和生命保证有重大的应用价值和研究意义,本文正是针对老年人对跌倒自动呼救的需求,围绕面向老人的便携式跌倒监护系统展开应用研究,基于跌倒的运动学原理,设计了跌倒检测算法。该跌倒检测算法通过对人体加速度和姿态的监测,构建了基于不同跌倒阶段的阈值判断的跌倒检测逻辑。经过系统实验验证,该跌倒检测算法对常见的跌倒情况的检测准确率在95%以上,能够满足老年人跌倒的危险情况检测
2、需求。关键词:老年人,穿戴式,跌倒检测,惯性传感器. -. -可修编- . z.1 绪论1.1 研究背景20世纪后,全世界正以无法想象的速度步入全球老龄化时代 1。我国也于20世纪末进入老龄化社会,预计到2050年将进入深度老龄化阶段,成为老龄化程度最严重的国家2,这也将从多方面给我国社会带来巨大的压力和挑战。在这种情况下,老年人的生活照料、出行等需求日益凸显,其中表达在平安保健问题尤为突出,主要表达在两点:第一,老年人身体机能退化,他们极易意外或因病尤其是心脑血管疾病跌倒。据国外统计,约有1/3的65岁老人每年至少跌倒一次3,跌倒在老年人死亡原因中的比例高达25% 4,跌倒非常容易造成老人骨
3、折、脏震荡;如果跌倒后得不到及时救助,会进一步提高致死率和致残率。跌倒不仅给老年人造成生理上的伤害,还会带来心理上的阴影。此外,当前社会老年人跌倒后无人敢搀扶、无人敢救,渐渐成为一种普遍的社会现象,在这种情况下,跌倒对老年人更是致命的。如何在信任危机日益严重的社会维护老年人的生命安康权利,成为日渐重要的议题。因此针对老年人对平安保证的迫切需求,利用电子信息技术解决此问题有着极大的社会价值。1.2 研究现状目前,跌倒检测算法研究的着眼点包括三点:第一是重点分析临界阶段或冲击,第二是重点分析跌倒后阶段,最后是两者同时分析。跌倒的临界阶段的研究Wu等人(University of Vermont,
4、USA)通过视频分析人体在跌倒时特定点的运动发现,这些点在跌倒临界阶段的速度为其他有意识行为的3倍5。Charif等人尝试通过计算机视觉系统追踪头部的运动,并通过粒子滤波和固定水平和速度阈值方法检测跌倒6。Noury等人(University of Grenoblein France)设计了一种佩带在腋下的自主检测装置,包括加速度计、一个倾角计和一个振动传感器。检测算法特征为速度阈值、从直立到水平的体位和跌倒后的静息状态7。Mathie等人使用位于胸部的三轴加速度传感器,Hwang等人使用陀螺仪,Bourke和Lyons等人使用双轴陀螺构建阈值算法,在临界阶段检测跌倒8-10。Tamura等人
5、提出在跌倒早期检测成功后,可以释放平安气囊保护跌倒者,但是该系统只能保护向后跌倒的情况11。跌倒的冲击的研究Williams等人在1998年首次描述了佩戴在腰部的跌倒检测装置,主要原理为通过压电传感器检测冲击,然后激发水银倾角计检测体位12。Doughty et al.继续研究了这种方法并由Tunstall (Whitley lodge, Yorkshire, England)公司推向市场13。Lindemann等人将三轴加速度传感器放入助听器中,通过3个阈值检测跌倒过程中的冲击:*y平面的总加速度为2g;竖直方向的最大速度为0.7m/s;所有方向的总加速度为6g14。跌倒后阶段的研究在跌倒后
6、阶段,人体一般处于平行于地面且静息状态。因此可以通过倾角计测量人体姿态或通过位于鞋底的压力传感器检测跌倒。但是跌倒后阶段的状态和被测者平时卧床休息的状态容易混淆,因此以上的检测方法一般需要和其他阶段的检测方法结合使用。综上所述,跌倒是行走过程中常见的危险事件,目前跌倒检测的研究前沿为跌倒检测准确率为90%左右。1.3 目前该领域存在的问题跌倒检测的准确率和实用性还有待提高。目前大多数跌倒检测研究的实用性不强,检测准确率在90%左右,对于危害严重的跌倒事件来说,还缺乏以完全监护使用者的行走平安。因此,跌倒检测系统的检测准确性和还有一定的提升空间。跌倒检测算法复杂难以产品化由于大局部跌倒检测算法复
7、杂,时效性差,很难工程化,目前市场上也没有相关的老年人监护产品,将算法做到简单实用,便于工程化也是本课题的难点。1.4 研究思路本文首先通过分析跌倒运动的定义及运动学原理,通过对跌倒进展建模,提取特征模式,进而通过样机平台搭建为算法研究提供硬件支持,最后通过系统实验验证算法的有效性和实用性。. -. z.2 跌倒运动学机理2.1 跌倒的定义与原因无意识的倒向地面,或者其他一些未能造成强烈冲击的情况,如失去意识、中风、癫痫发作。这个跌倒的定义仍被许多研究机构沿用至今,扩展后还包含了晕厥和心血管疾病所造成的跌倒 15。跌倒时间段的划分跌倒主要分为四个阶段:跌倒前阶段、临界阶段、跌倒后阶段和恢复阶段
8、图1。跌倒前阶段与日常行为一样;临界阶段以人体倒向地面为起始,以人体承受剧烈冲击为完毕,这个阶段大概持续;跌倒后阶段以人体承受剧烈冲击为起点,该阶段人体处于不活动状态。恢复阶段是跌倒者自行站立或在他人帮助下恢复站立16。图1 跌倒的阶段划分关于跌倒的定义,尽管缺乏统一的标准,但是近来研究人员广泛采用是:人体有意识或是无意识的倒向并躺在地面或是其较低水平面的事件17。跌倒的原因是多种因素相互结合的。主要包括生理因素、疾病因素、药物因素、境因素、性别及社会心理因素等18-21。2.2 跌倒的运动学描述身体重心的变化由跌倒的定义可知,跌倒事件是人体从一个较高平面倒向另一个较低平面的过程,从而必然伴随
9、着人体重心的位置变化一般情况是重心的降低。如图2示,以向前跌倒为例,人体重心图中蓝色圆点随时间呈现出逐渐减低的过程。在理想情况下,跌倒的运动轨迹重心高度的时间变化轨迹应该呈现为抛物线形状,其运动学方程为:2.1其中,h0为人体重心的初始高度,单位为m;g为重力加速度,单位为m/s2;t为时间,单位为s,初始时刻t=0为跌倒开场时刻。th图2 跌倒时重心变化的运动学特征在实际应用中,采用加速度积分的方法描述人体重心变化通常存在较大累积误差,因此,一般结合其他传感量进展修正,譬如检测高度的气压计等。身体姿态的变化在跌倒过程中另一个运动学特点是人体姿态的变化。按照跌倒时的姿态分类,跌倒可以分为向前跌
10、倒、向后跌倒和侧向跌倒。以人体面向*轴正方向为例,向前跌倒可以等效为人体绕Y轴逆时针旋转,向后跌倒可以等效为人体绕Y轴顺时针旋转,同理,侧向跌倒可以等效为人体绕*轴的旋转。因此,假设将人体跌倒简化为刚体的运动,则其姿态的变化可以用欧拉角来表示。如图3所示,设定*yz轴蓝色为地理坐标系的参考轴,*YZ轴红色为捷联在人体上的载体坐标系的坐标轴,则人体的运动可以用欧拉角表示。但是采用欧拉法确定姿态角的方法有一定的局限性,由于微分方程中存在三角函数,实时计算比拟困难,而且当为90时,方程出现奇异点,姿态角无法解算。. -. z.3 跌倒检测监护装置硬件设计跌倒检测监护装置系统分为跌倒检测终端和监护中心
11、两局部。跌倒检测终端由电源模块、处理器模块、传感器模块、无线通信模块、人机交互模块构成,分别用于终端供电、运行跌倒检测算法、物理惯性信息的检测、跌倒检测算法的执行、报警信息的发布和取消以及信息的无线传输。监控中心分为则包括网管和PC,分别用于网络的管理和信息的储存和后处理。如图4所示。图4 跌倒检测监护装置硬件原理1处理器模块该局部主要有AVR单片机ATmega8L最小系统构成,包括电源、晶振、复位和下载电路,其原理图如图5。图5 处理器模块电路图2惯性传感器模块该局部由三轴加速度传感器MMA7260及其外部RC滤波电路构成,其原理图如图6:图6 惯性传感器模块电路图MMA7260是飞思卡尔公
12、司生产的性价比最高的微型电容式加速度传感器,其采用了信号调理、单级低通滤波器和温度补偿技术。可用三轴加速度计利用重力分量换算原理,来测量角度,与其他数字量倾角传感器相比自然要精准许多,因为模拟量的,可将电压值换算对应倾斜角度值。并且可通过G1G2管脚调整加速度的量程,适用于检测人体的加速度。3无线通信模块该局部主要由CC2430 Zigbee SoC构成,采用AC*公司的AT7020芯片天线,以到达节省体积的目的。其原理图如图7:图7 无线模块电路图4电源模块该局部主要有锂电池、低压差调节器LPO和电量计组成,锂电池采用3.7V可充电电池,电池端具有保护电路;低压差调节器采用TPS79633芯
13、片,将3.7V调节到3.3V;电量计采用BQ27510芯片,检测锂电池的电量和预测剩余使用时间。另外,需要注意的是,加速度传感器的电源应该与其他电源隔离,以减少噪声干扰。该局部原理图如图8:图8 电源模块5人机交互模块该局部主要由各种外设包括按键开关、电源接口、LED指示灯和蜂鸣器等构成,其原理图如图9所示:图9 人机交互模块. -. z.4 跌倒检测算法设计跌倒检测的目标是能够将跌倒(Fall)与日常生活的正常动作(Activities of Daily Life,ADL)区分开来,准确地检测跌倒的发生,并智能判断(并执行)是否需要报警求助。从而尽可能地缩短救助时闻,减小跌倒带来的伤害(尤其
14、是长时间晕厥),降低误报率,最终提升被监测者的生活质量。4.1 跌倒运动学模式分析4.1.1 跌倒与日常行为分类无论是跌倒还是日常行为,如果以分解动作来看都是不同静态姿势之间的转化,当然转化过程中的剧烈程度和时间点也是不同跌倒和日常行为的主要区别。人体的常见静态姿势包括:左/右侧躺、仰卧、俯卧、站立、坐、蹲等。在这些静态姿势之间的转化就形成了人体的运动。常见的跌倒主要分为3种:向前跌倒,即从站立到俯卧的剧烈过程;向后跌倒,即从站立到仰卧的剧烈过程;侧向跌倒,包括向左/右跌倒,即从站立到侧卧的剧烈过程。据ONeil等人的研究说明,前向跌倒是最为常见的跌倒,约占整个跌倒事件的60%左右22。当然实
15、际生活中,还可能出现更加复杂的跌倒情况,这些跌倒模式较为复杂且极少发生,将作为下一步研究的对象,不在本文的讨论围。为简化研究,6种容易和跌倒事件混淆的日常行为被作为主要研究对象,其包括:走;跑;跳;上下楼;坐下;躺下。本研究中的跌倒检测即是探讨常见的3种跌倒事件与易混淆的6种日常行为的主要区别和识别方式。由于其主要区别包含两个方面:其一是运动所对应的静态姿势不同,对应的惯性量即是姿态的变化;其二是在不同静态姿势间转换的剧烈程度不同,对应的惯性量即是加速度的变化。下面分别从这两个方面讨论跌倒和日常行为的运动模式。4.1.2 加速度分析数据预处理加速度的时域信号表达了人体运动过程中剧烈程度,以及在
16、不同时刻剧烈程度的分布。然而在实际信号采集过程中,原始数据往往包含了有效信号和各种干扰信号:人体运动加速度;重力加速度;人体抖动;测量噪声;惯性传感器与身体的相对运动其中,人体运动加速度和重力加速度是有效信号,这两者可以表达人体的运动状态,而后三个信号属于干扰信号,应当予以抑制和去除。由于这三种干扰信号都与脉冲噪声相似,因此选用中值滤波的方法对加速度信号进展预处理,其处理效果如图10所示。图10 中值滤波处理加速度信号总加速度强度Acceleration Vector Magnitude总加速度即为三轴加速度的矢量和,它能比拟直观的表达人体运动的剧烈程度,其数学表示为:5.1常见的6种日常行为
17、的总加速度强度AVM与竖直轴加速度A*如图11所示,其中总加速度AVM以黑色曲线表示,竖直轴加速度A*以蓝色曲线表示。分析日常行为的加速度曲线可得以下特征;总加速度AVM的幅值以1g为中值摆动,由于重力加速度的存在。竖直轴加速度A*的幅值以-1g为中值摆动,由于*轴正方向与重力方向相反。坐下起立、走、躺下和上楼的总加速度AVM变化幅值在2g之,说明运动剧烈程度较小。跑和跳的总加速度AVM变化幅值在4g之,说明运动剧烈程度较大。在人体保持直立状态的日常行为中除躺下外的其他5种,总加速度AVM与竖直加速度A*的变化趋势根本一致。躺下时,由于人体从直立状态转化到平躺状态,因此竖直加速度A*曲线存在明
18、显阶跃变化。图11 常见日常行为的总加速度黑色与竖直轴加速度曲线蓝色常见的3种跌倒事件的总加速度强度AVM与竖直轴加速度A*如图12所示,其中总加速度AVM以黑色曲线表示,竖直轴加速度A*以蓝色曲线表示。分析这3种跌倒事件的加速度曲线可得以下特征:不同方向的跌倒事件加速度曲线非常相似,变化趋势根本保持一致。总加速度曲线AVM中存在明显峰值,该峰值是由身体与地面冲击造成的。总加速度曲线AVM中峰值之前时刻,对应的竖直加速度曲线A*存在一个明显的加速度谷值,该谷值是由身体在跌倒过程中的失重造成的。总加速度曲线AVM中峰值之后一段时间,对应的竖直加速度曲线A*存在一个明显阶跃变化,这是由于人体姿态变
19、化造成的。图12 常见跌倒事件的总加速度黑色与竖直轴加速度曲线蓝色比拟日常行为和跌倒事件的和加速曲线AVM特征可知,两者总加速度AVM峰值的大小存在明显差异。经过大量模拟跌倒和日常行为实验,统计汇总后,跌倒和日常行为的总加速度峰值围如表6和表7所示。表1 日常行为总加速度最值ADL行走坐蹲上楼下楼慢跑跳Ma*(g)1.881.541.751.781.984.53.11Min(g)1.551.471.051.121.542.783.61表2 跌倒总加速度最大值围Fall前倒后倒侧倒Ma*(g)5.785.685.78Mean(g)4.984.454.37Min(g)4.764.824.454.1
20、.3 姿态分析人体处于静止状态时,可以通过重力加速度在三个方向的分布估计人体姿态,仅需加速度传感器支持,且计算简单。图13显示了跌倒事件中竖直方向加速度的变化,通过重力加速度的分布变化,同时也反映了人体姿态的变化。在人体直立时,竖直方向加速度为-1g正方向为垂直地面向上,而在人体平躺时,竖直方向加速度为0g左右。从而,该曲线表达了在跌倒过程中人体姿态的典型变化:在失重阶段人体与竖直轴夹角由0左右快速变为90左右,随后在平躺阶段时,人体与竖直轴夹角保持90左右一段时间。虽然跌倒时人体姿态变化与躺下的情况一致,但是足以区分其他日常行为。图15 跌倒时竖直方向加速度曲线4.2 跌倒检测算法设计通过以
21、上对日常行为和跌倒的加速度和人体姿态分析可知,跌倒可以分为三个阶段,第一阶段是失重阶段,此时人体开场无意识的倒向地面;第二阶段是冲击阶段,此时人体其他部位与地面发生冲击;第三阶段是平躺阶段,此时人体平躺在地面上处于静息状态。并且,跌倒过程中突出的特征可以归结为如图16所示:在人体处于失重阶段时,竖直轴加速度Ay先下降再急剧增加。在人体承受冲击阶段时,总加速度的激增。冲击之后,由于人体姿态变化而处于平躺状态,总加速度(AVM)维持在1g左右,Y轴加速度维持在0g附近。图16 跌倒的运动学特征经采样分析日常行为事件几乎不能同时满足以上3个特征,这为该跌倒检测系统具有较高的特异性提供了保障。为有效识
22、别以上3个特征,我们通过总加速度AVM与人体与竖直轴夹角TAy的变化时序来设计跌倒检测算法。总加速度AVM的定义已经在前文加速度分析中给出,而在人体处于静息状态时非运动状态,竖直轴夹角TAy的计算采用重力加速度分布的方法估算,其公式如下:5.2其中TAy表示人体与竖直轴夹角,单位为deg;Ay表示竖直轴加速度,单位为m/s2;G表示重力加速度,单位为m/s2。基于以上分析和大量试验,总结出以下跌倒检测算法,其流程图如图17所示。该算法的核心思想是加速度阈值判断结合人体姿态检测,因此需要定义三个阈值用于评估人体运动剧烈状态和人体姿态:HAT为加速度高阈值,当AVM大于HAT时,说明人体受到强烈冲
23、击。LAT为加速度低阈值,当AVM小于LAT时,说明人体处于不活泼状态。AT为角度阈值,当TAy绝对值小于AT时,说明人体处于平躺状态。此外,还需要定义三个延时用来控制不同阈值检测的时序。所有阈值和延时都经过大量日常行为和模拟跌倒实验调校和验证。图17 跌倒检测算法流程图其主要流程描述如下:三轴加速度传感器的测量采样频率为200Hz。每次采样后有三轴加速度计算总加速度AVM。设定总加速度高阈值HAT,假设检测到AVM HAT,说明存在强烈冲击,继续判断。设定总加速度低阈值LAT,假设3s检测AVM LAT ,说明人体处于静息状态,继续判断。假设以上条件满足,则通过三轴加速度分布进展姿态判断,延
24、时循环两次,假设|TAy| AT,则为接近水平状态,判定跌倒。. -. z.5 实验验证跌倒检测的目的是在于准确区分日常行为和跌倒事件,因此跌倒检测实验也是围绕这两种行为展开的。实验容是:志愿者按照要求执行相应的日常行为和模拟跌倒事件。5.1 跌倒检测实验设计跌倒检测的实验工程也借鉴了之前N.Noury的关于跌倒检测准确性的评估37。综合常见的日常行为和跌倒事件,在该实验中志愿者执行的工程如表3所示。表3 跌倒检测实验工程实验工程工程细分报警种类向后跌倒以坐姿完毕是以平卧完毕是以侧卧完毕是快速恢复否向前跌倒膝盖着地是胳膊着地是以平卧完毕是倒后翻滚以平卧完毕是倒后翻滚以侧卧完毕是快速恢复否左侧跌
25、倒以平卧完毕是快速恢复否右侧跌倒以平卧完毕是快速恢复否晕厥沿墙跌坐是日常行为坐 起立否躺 起立否行走否弯腰拾物否咳嗽或打喷嚏否6.4.2 跌倒检测实验结果经过志愿者总共1000次行为模拟,其中包括400次模拟跌倒和600次模拟日常行为,跌倒检测的检测准确率统计如表13所示。表13 跌倒检测实验结果统计行为种类实验次数报警次数未报警次数准确率日常行为5000500100%向前跌倒10099099%向后跌倒1001000100%向左跌倒10094694%向右跌倒10095595%由于跌倒检测的输出只有跌倒报警和未报警两个值,因此跌倒检测的准确性可以由以下四种情况发生的概率来评价。真阳性事件true
26、 positive, TP,即跌倒确实发生,并且跌倒报警发出。假阳性事件false positive, FP,即跌倒没有发生,而跌倒报警发出。真阴性事件true negative, TN,即跌倒没有发生,并且没有跌倒报警。假阴性事件false negative, FN,即跌倒确实发生,而没有跌倒报警。因此,跌倒检测的准确性可以通过以下两个指标来评价:检测敏感性,即在所有跌倒事件中,成功检测出跌倒的比率,与其相反的为漏警率。其计算公式如下:6.1检测特异性,即在所有日常行为事件中,成功检测出未跌倒的比率,与其相反的为虚警率。其计算公式如下:6.2其中,TP,FP分别真阳性和假阳性事件,TN,FN
27、分别真阴性和假阴性事件。由于在跌倒检测实验中,所有日常行为事件都被准确检测出,因此假阳性事件为零,说明该跌倒检测算法具有很高的特异性在实验中为100%。同时,由于存在小局部跌倒事件没有被成功检测出,即检测存在假阴性事件,但是该情况发生概率很小见表13,因此该跌倒检测算法同样具有很高的敏感性在实验于95%。. -. z.参考文献. 联合国人口司. 联合国老龄化议题-老龄化问题概况EB/OL, ./chinese/esa/ageing/introduction.htm, 2013-03-09国务院新闻办. 2010年第六次全国人口普查主要数据公报EB/OL, .gov./wszb/zh
28、ibo449, 2010-04-28吕筠, 立明. 老年伤害研究回忆J. 疾病控制杂志, 1999, 3(4): 300-303郝晓宁, 胡鞍钢. 中国人口老龄化: 安康不平安及应对政策J. 中国人口. 资源与环境, 2010, 2(3): 73-78Wu G. Distinguishing fall activities from normal activities by velocity characteristicsJ. Journal of Biomechanics, 2000, 33:1497-1500Nait-Charif H., Mckenna S. J. Activity Su
29、mmarisation and Fall Detection in a Supportive Home EnvironmentA. Pattern Recognition, 2004, ICPR 2004, Proceedings of the 17th International Conference on IEEEC, 2004:323Noury, N., Herv, T., Rialle, V., et al. Monitoring behavior in home using a smart fall sensor and position sensorsA. Microtechnol
30、ogies in Medicine and Biology, 1st Annual International Conference on IEEEC, 2000:607Mathie, M., Celler, B., Lovell, N.H., et al. Classification of basic daily movements using a tria*ial accelerometerJ. Medical and Biological Engineering and puting, 2004, vol. 42, no. 5:679-687Hwang, J., Kang, J., J
31、ang, Y., et al. Development of novel algorithm and real-time monitoring ambulatory system using Bluetooth module for fall detection in the elderlyA. Engineering in Medicine and Biology Society, IEMBS04, 26th Annual International Conference of the IEEEC, 2004:2204Bourke, A., Lyons, G. A threshold-bas
32、ed fall-detection algorithm using a bi-a*ial gyroscope sensorJ. Medical engineering & physics, 2008, vol. 30, no. 1:84-90Tamura, T., Yoshimura, T., Sekine, M., et al. A Wearable Airbag to Prevent Fall InjuriesJ. IEEE Transaction on Information Technology in Biomedicine, 2009, vol. 13, no. 6:910-914.
33、Williams, G., Doughty, K., Cameron, K., et al. A smart fall and activity monitor for telecare applicationsA. Engineering in Medicine and Biology Society, Proceedings of the 20th Annual International Conference of the IEEEC, 1998:1151Doughty, K., Lewis, R., McIntosh, A. The design of a practical and reliable fall detector for munity and institutional telecareJ. Journal of telemedicine and telecare, 2000, vol. 6, no. suppl 1:150-154Lindemann, U., Hock, A., Stuber, M., et al. Evaluation of a fall detector based on accelerometers: A pilot studyJ.
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024版沙子买卖合同
- 二零二五年海南二手房买卖及配套设施完善合同3篇
- 西安交通大学《过程分子生物学》2023-2024学年第一学期期末试卷
- 二零二五年度鞋类批发市场购销合同市场地位巩固
- 二零二五年度酒店消防器材维护保养及更换合同3篇
- 二零二五年度水利工程安全评价技术服务合同3篇
- 二零二五年度新能源汽车电池回收利用合伙协议书3篇
- 二零二五年股东股权置换合同参考范本6篇
- 二零二五版生物科技研发技术顾问聘用协议2篇
- 二零二五版物流企业劳动安全及货物保护协议合同3篇
- 2023年保安公司副总经理年终总结 保安公司分公司经理年终总结(5篇)
- 中国华能集团公司风力发电场运行导则(马晋辉20231.1.13)
- 中考语文非连续性文本阅读10篇专项练习及答案
- 2022-2023学年度六年级数学(上册)寒假作业【每日一练】
- 法人不承担责任协议书(3篇)
- 电工工具报价单
- 反歧视程序文件
- 油气藏类型、典型的相图特征和识别实例
- 流体静力学课件
- 顾客忠诚度论文
- 实验室安全检查自查表
评论
0/150
提交评论