目录
? 1 本质模态函数(IMF)
? 2 经验模态分解(EMD)
? 3 结论
? 4 相关条目
? 5 参考文献
? 6 外部链接
[编辑]本质模态函数(IMF)
任何一个资料，满足下列两个条件即可称作本质模态函数。
⒈局部极大值(local maxima)以及局部极小值(local minima)的数目之和必须与零交越点(zero crossing)的数目相等或是最多只能差1，也就是说一个极值后面必需马上接一个零交越点。
⒉在任何时间点，局部最大值所定义的上包络线(upper envelope)与局部极小值所定义的下包络线，取平均要接近为零。
因此，一个函数若属于IMF，代表其波形局部对称于零平均值。此类函数类似于弦波（sinusoid-like），但是这些类似于弦波的部分其周期与振幅可以不是固定。因为，可以直接使用希尔伯特转换，求得有意义的瞬时频率。
[编辑]经验模态分解(EMD)
EMD算法流程图
建立IMF是为了满足希尔伯特转换对于瞬时频率的限制条件之前置处理，也是一种转换的过程。我们将IMF来做希尔伯特转换可以得到较良好的特性，不幸的是大部分的资料并不是IMF，而是由许多弦波所合成的一个组合。如此一来，希尔伯特转换并不能得到正确的瞬时频率，我们便无法准确的分析资料。为了解决非线性（non-linear）与非稳态（non-stationary）资料在分解成IMF时所遇到的困难，便发展出EMD。
经验模态分解是将讯号分解成IMF的组合。经验模态分解是借着不断重复的筛选程序来逐步找出IMF。
以讯号为例，筛选程序的流程概述如下:
步骤 1 : 找出中的所有局部极大值以及局部极小值，接着利用三次样条
(cubic spline)，分别将局部极大值串连成上包络线与局部极小值串连成下包络线。
步骤 2 : 求出上下包络线之平均，得到均值包络线。
步骤 3 : 原始信号与均值包络线相减，得到第一个分量。
步骤 4 : 检查是否符合IMF的条件。如果不符合，则回到步骤1并且将
当作原始讯号，进行第二次的筛选。亦即
重复筛选次
直到符合IMF的条件，即得到第一个IMF分量，亦即
步骤 5 : 原始讯号减去可得到剩余量，表示如下式
步骤 6 : 将当作新的资料，重新执行步骤1至步骤5，得到新的剩余量。如此重复次
.
.
.
当第个剩余量已成为单调函数(monotonic function)，将无法再分解IMF 时，整个EMD的分解过程完成。原始讯号可以表示成个IMF分量与一个平均趋势(mean trend)分量的组合，亦即
如此一来，原始资料便分解成n个IMF和一个趋势函数，我们便可将IMF做希尔伯特转换来进行瞬时频率的分析。
[编辑]结论
傅立叶变换是将一个讯号分解成无限多个弦波来分析资料，但是希尔伯特-黄转换则是将一个讯号分解成数个近似于弦波的讯号（周期、振幅不固定）和一个趋势函数来做分析。
两者各有其优缺点，整理如下
优点：
1.避免复杂的数学运算
2.可分析频率会随时间变化的讯号
3.较适于分析气候、经济等具有趋势的资料
4.可以找出一个函数的趋势
缺点：
1.缺乏严谨的物理意义
2.需要复杂的递回，运算时间反而比短时距傅立叶变换要长
3.希尔伯特转换未必能正确计算出本质模态函数之瞬时频率
4.无法使用快速傅立叶变换
5.只有在特例（组合较简单的资料）时使用希尔伯特-黄转换较快
[编辑]相关条目
?Norden E. Huang, et al. "The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series
analysis." Proc. R. Soc. Lond. A (1998) 454, 903–995 (Link)
?陈韦佑, "以希尔伯特-黄转换法为GPS接收机抑制调频干扰", 国立台湾大学电机工程研究所硕士论文, 2007.
?Jian-Jiun Ding, Time frequency analysis and wavelet transform class ppt,the Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2009.
?Jian-Jiun Ding,"时频分析与小波转换",available in
https://www.doczj.com/doc/718710298.html,.tw/TFW.htm
2.Matlab代码
https://www.doczj.com/doc/718710298.html,/matlabcentral/fileexchange/19681
2.1 EMD分解代码（emd.m）
function imf = emd(x)
% Empiricial Mode Decomposition (Hilbert-Huang Transform) % imf = emd(x)
% Func : findpeaks
x = transpose(x(:));
imf = [];
while ~ismonotonic(x)
x1 = x;
sd = Inf;
while (sd > 0.1) | ~isimf(x1)
s1 = getspline(x1);
s2 = -getspline(-x1);
x2 = x1-(s1+s2)/2;
sd = sum((x1-x2).^2)/sum(x1.^2);
x1 = x2;
end
imf{end+1} = x1;
x = x-x1;
end
imf{end+1} = x;
% FUNCTIONS
function u = ismonotonic(x)
u1 = length(findpeaks(x))*length(findpeaks(-x));
if u1 > 0, u = 0;
else, u = 1; end
function u = isimf(x)
N = length(x);
u1 = sum(x(1:N-1).*x(2:N) < 0);
u2 = length(findpeaks(x))+length(findpeaks(-x));
if abs(u1-u2) > 1, u = 0;
else, u = 1; end
function s = getspline(x)
N = length(x);
p = findpeaks(x);
s = spline([0 p N+1],[0 x(p) 0],1:N);
2.2 找极值代码（findpeaks.m）
function n = findpeaks(x)
% Find peaks.
% n = findpeaks(x)
n = find(diff(diff(x) > 0) < 0);
u = find(x(n+1) > x(n));
n(u) = n(u)+1;
2.3 绘制时-频曲线以及尺度分解代码（plot_hht.m）
function plot_hht(x,Ts)
% Plot the HHT.
% plot_hht(x,Ts)
% Ts ：time interval （sec）
% :: Syntax
% The array x is the input signal and Ts is the sampling period.
% Func : emd
% Get HHT.
imf = emd(x);
for k = 1:length(imf)
b(k) = sum(imf{k}.*imf{k});
th = angle(hilbert(imf{k}));
d{k} = diff(th)/Ts/(2*pi);
end
[u,v] = sort(-b);
b = 1-b/max(b);
% Set time-frequency plots.
N = length(x);
c = linspace(0,(N-2)*Ts,N-1);
for k = v(1:2)
figure, plot(c,d{k},'k.','Color',b([k k k]),'MarkerSize',3); set(gca,'FontSize',8,'XLim',[0 c(end)],'YLim',[0 1/2/Ts]); xlabel('Time'), ylabel('Frequency');
end
% Set IMF plots.
M = length(imf);
N = length(x);
c = linspace(0,(N-1)*Ts,N);
for k1 = 0:4:M-1
figure
for k2 = 1:min(4,M-k1), subplot(4,1,k2), plot(c,imf{k1+k2}); set(gca,'FontSize',8,'XLim',[0 c(end)]); end
xlabel('Time');
end
金翔公司体系管理程序文件
为使新产品及变更的产品之工程文件能得到及时供给生产现场使用，为客户提供合格样板，确保变更后的产品、材料得到有效管理，特编制本文件。 2.范围： 适用于本公司所有产品的工程段文件编制，ECN管控、客户样板制作及相关工程技术问题处理方法。 3.职责： 3.1.销售部负责：提供客户信息及传达市场信息予研究所，并对样板制作的期限完成要求 做出跟进； 3.2.研究所负责： 3.2.1.样板制作资料： 321.1.将样板制作资料要求下达到采购部、品质部、生产部； 3.2.1.2.样板的制作要求及进度跟踪、最终确认样品能否达到要求； 3.2.2.新产品及ECN后工程资料提供： 3.2.2.1.负责在新产品量产前将图纸资料/机箱配置一览表提供予生产现场及作业场
3222 负责所有ECN产品投入生产前，将更改的工程资料提供生产现场及作业场 所； 3.3.采购部负责样板外购物料的采购及跟催； 34 生产部负责对要求的样板进行生产、组装，并保证与要求的相符负责，在新产品量产前将 各工序作业指导书提供予生产现场及作业场所； 3.5.品质部： 3.5.1.负责对制作好的样板进行检查合格后放行； 3.5.2.负责编制用于来料、制程、出货检查的各类规范性文件； 4、定义 “工程变更”定义：本厂的产品，由客户提出的更改或由本厂提出并经得客 户同意的更改，标准、结构或工艺。 4.程序： 4.1.客户样板的制作：（指已通过设计开发评测后的产品，设计开发段的产品样板制作 按《设计开发控制流程》执行） 4.1.1. 销售部应做到： 4.1.1.1.接到客户通知或根据市场调查需制作产品样板予客户确认时，由销售部开出 《样板制作需求单》分发研究所。 4.1.1.2.《样板制作需求单》应注明产品结构、颜色、外观、交期等需求；销售部负
工程变更管理制度 1目的和范围 1.1建设过程中的工程变更，对工程质量、工期和造价都可能产生影响。为了规范工程变更的申报、审查、批准等工作程序，特制订本制度。 1.2 本制度所指的工程变更，是指在施工过程中由某种原因引起的，需要对原工程设计做出局部调整或修改的事项。 1.3 本制度适用于由建设单位或其它参建单位（包括设计单位、施工单位、调试单位、监理单位、设备制造厂家等）提出的工程变更要求。 1.4在施工图审查过程中向设计单位提出的设计调整或修改意见，设计院直接修改图纸,可不统计为工程变更。 2 职责 2.1 所有工程变更都必须经过建设单位批准。在审批工程变更时,建设单位应确认因变更引起的工期和工程造价等方面的影响. 2.2 所有经建设单位批准的工程变更，都必须由设计单位发出设计变更通知单才可付诸实施。 2.3项目监理部负责审查工程变更的申请；监督和检查经建设单位批准的工程变更的实施。 2.4施工单位负责工程变更的实施。 3．管理要求 3.1 工程变更的提出 3.1.1 提出工程变更要求的原因主要有以下几方面： 1 设计图纸有差错或设计内容深度不够； 2 设计与实际情况不符合，或者设计条件(地质、设备等)有变化； 3 由于现场条件所限，设计采用的材料规格、品种、质量不能完全符合要求； 4 上级（项目审批）单位提出变更要求； 5 因施工问题需要做出变更； 6 为节约投资、改善运行条件需要做出变更；
7 根据技术改进和合理化建议需要做出变更。 3.1.2 工程变更的要求可以由设计单位、建设单位、施工单位、调试单位、设备制造单位等提出。 1 设计单位的工程变更要求，可以使用施工图升版的方式或者直接向建设单位发出设计变更通知单的方式提出。设计变更通知单应说明变更原因、内容、工程量增减及其相关预算费用等内容。 2 建设单位的工程变更要求，可使用工程联系单直接向设计单位提出，比较重大的变更要求应向设计单位发出正式公函。上级（项目审批）单位和运行单位的工程变更要求，都应当通过建设单位提出。 3项目监理部的工程变更要求,可使用工程联系单直接向设计单位提出, 建设单位 审批 4 施工单位或设备供应商的工程变更要求，应填写工程变更申请单，报项目监理部审查、建设单位审批。 3.2 工程变更的分类 工程变更可分为四大类： 3.2.1小型变更：不改变设计原则、不影响质量和安全经济运行、不影响外观形象，而且不增减预算费用的变更事项。例如图纸尺寸的差错更正、材料等强换算代用、图纸细部增补详图、图纸间矛盾问题的处理等。这类变更不引起工程费用变更或变更甚少。 3.2.2 一般变更：工程内容与工程量有少量变化，但不涉及可研或初步设计已审定的原则，或对局部施工计划与施工进度有一定影响，但不影响工程总进度。变更引起的工程费用的增減,按施工合同约定无需调整。 3.2.3 较大变更：施工图的设计范围、工艺流程、设备布置有一定变化，但未违反初步设计审定的设计原则，不影响工程质量和建设总工期，或者变更引起的工程费用增減,按施工合同约定需要调整。 3.2.4 重大变更：涉及可研或初步设计审定的设计原则、方案或规模、主要设备换型，工程费用增加超限从而将导至原审定的概算调整。 3.3 工程变更的审批
变更控制管理程序 目的 规定对本厂已批准的各类文件、设备、设施、物料供应商及包装材料发 生变更时控制的程序, 以保证产品的生产过程始终处于受控制的状态。 适用范围 范围适用于全厂所有产品的上述变更。 职责 各个部门或个人可根据工作职责提出变更申请。 受变更影响的各部门对变更申请进行评估、审核、列出相关的实施计划。 并对经批准的变更申请和行动计划负责实施，负责将实施情况书面报告给质量管理部门。 4.2.3质量管理部部负责变更的管理，指定专人负责变更控制工作，界定变更分类，组织变更评估和审核，制订变更实施计划，跟踪变更的实施， 对变更效果进行评价，及时反馈变更信息。?? ?? 企业负责人对所有变更申请和实施计划进行批准以及对变更进行批准。 4 程序 变更类型 4.1.1所有已批准的SOP文件及各种记录的变更 4.1.2技术文件的变更技术文件包括：工艺规程、质量标准、分析方法。 4.1.3关键设备仪器、设施的变更： 4.131关键设备仪器：直接用于生产和QC测试，对产品质量直接造成影 响的各类设备仪器。 4.1.3.2QC 关键设施：对洁净室环境直接造成影响，对产品质量直接造成影响的各 类公用设施。
4.1.4 物料供应商变更：包括原料、辅料、包装材料供应商的变更。 4.1.5 标签、说明书、包装材料的变更：包括标签、说明书、单盒、中盒和 外箱上印刷的文字、颜色、图案和尺寸、材质等的变化。 生产用物料的贮存条件和（或）有效期的变更，包括原辅料、包装材 料、中间体和成品。 4.1.10 在日常生产过程中发生的非计划性变更 变更程序 421在需要变更时，由需要发生变更的岗位责任人员向 QA 口头提出申请, QA 下发已编号的变更申请表，并在变更登记表中登记。 4.2.2 提出变更的申请人按以下要求填写 “变更申请表”， 交本部门负责人 评估、签署意见： （ 修改后版本号 ） ，文件应写完整编号。 4.2.2.4 申请人、申请日期及所属部门。 4.2.3 申请变更部门负责人按照以下要求对变更进行评估： 4.2.3.1 变更是否违背了政策法规、法定标准。 4.2.3.2 变更是否属于需验证的范围。 判断标准参见：验证主计划、 子系统 验证 主计划及验证计划；设施、设备和仪器的确认；工艺验证；清洁验证； 4.1.6 委托生产商要求的变更。 4.1.7 政府部门要求的变更。 4.1.8 人员组织机构图及其他变更。 4.1.9 4.2.2.1 变更名称： 明确变更的主题及变更类型。 4.2.2.2 变更理由： 描写该变更提出的原因。 4.2.2.3 变更内容： 说明 “原来内容” （ 原来版本号 ）、 “ 修改后内容 ”
工程变更管理实施细则 第一章总则 第一条为了加强对政府投资非经营性建设项目的全过程监督管理，规范工程建设中的工程变更行为，合理有效控制工程投资规模，根据有关法律法规和《XXX政府投资非经营性建设项目代建管理办法（试行）》、《XXX重大建设工程项目监督检查办法》、《XXX重点项目工程招投标及工程质量管理实施意见》等有关规定，制订本实施细则。 第二条本实施细则所称工程变更指设计变更和现场签证。 第三条根据《中华人民共和国招标投标法》规定，必须进行施工招标的工程，原则上应实行无工程变更的管理制度，可变可不变的项目不办理变更，但发生特殊情况时，由于非施工单位原因造成的工程内容及工程量的增减，可以办理工程变更。 第四条项目建设管理过程中涉及的工程变更必须根据本实施细则规定的原则和程序办理相关手续，并在招标文件及合同条款中明确写明执行本实施细则。 第五条工程变更必须遵循合法性、真实性、科学性、全面性和时效性的原则。设计文件一经审核批准，原则上不允许随意变更。因特殊情况确需变更的项目，应加强规范管理，严格控制投资规模，不得以任何方式肢解或重复工程变更。 第六条在工程建设过程中，需要变更设计图纸和增减工程项目，应在确保工程技术标准、质量标准保持不变和安全保障的情况下变更，并需按程序审批和签证。 第二章设计变更 第七条设计变更是指在工程实施过程中，因工程项目自身的性质和特点，或因设计文件深度不够，或因政策法规调整，或因不可预见因素与环境情况变化，需要变更原有设计文件时，由设计单位充分论证后，对工程项目标准、功能、材料、工艺、质量、构造、尺寸、指标及数量等做出修改或补充设计文件的行为。 第八条设计变更的程序 一、设计变更的提出
质量管理体系文件Quality Management System File 产品工程变更控制程序 文件编号：**-QP-051 版本号：A/0版 **公司20**-**-**实施 ** 发布
文件发行/修改履历表 1目的
规范产品工程变更，确保所有的变更均经过验证和授权人员的批准，避免变更的影响没有被充分识别而影响产品的品质。 2范围 本文件适用于**公司产业群整车产品工程变更的管理。 3引用标准 3.1ISO/TS16949: 2009 质量管理体系——汽车生产件及相关服务件组织应用 ISO9001：2008的特殊要求 3.2TSG Z0004-2007 特种设备制造、安装、改造、维修质量保证体系基本要求 （适用于整车生产**公司车用气瓶的安装，即CNG车型）。 3.3强制性产品认证《**公司质量保证能力要求》，**公司产业群所涉及的强制性产品 认证和所引用的CCC实施细则包括： CNCA-02C-023:2008 机动车辆类强制性认证实施规则 **产品 CNCA-02C-056:2005机动车辆产品强制性认证实施规则机动车回复反射器产品CNCA-02C-059:2005机动车辆产品强制性认证实施规则汽车后视镜产品 CNCA-02C-061:2005机动车辆产品强制性认证实施规则汽车门锁及车门保持件产品 CNCA-02C-060:2005机动车辆产品强制性认证实施规则汽车内饰件产品 CNCA-02C-062:2005机动车辆产品强制性认证实施规则汽车燃油箱产品 CNCA-02C-058:2005机动车辆产品强制性认证实施规则汽车外部照明及光信号装置产品 CNCA-02C-063:2005 机动车辆产品强制性认证实施规则汽车座椅及座椅头枕 产品 CNCA-02C-026:2005机动车辆产品强制性认证实施规则汽车安全带产品 CNCA-02C-055:2005机动车辆产品强制性认证实施规则机动车用喇叭产品 CNCA-02C-066:2005机动车辆产品强制性认证实施规则汽车行驶记录仪产品 4术语 4.1工程变更： 来自内部或外部对产品设计、制造工艺、物料、供应商等的变更。包括产品设计变更和工艺变更。 4.2产品设计变更： 来自内部或外部原因对产品设计的更改，包括对产品及其零部件尺寸、结构、功能、材料等以及由此而引起的数模、图纸、BOM、技术条件、供应商等方面的变 更。 4.3工艺变更：
几种时频分析方法综述2——希尔伯特—黄变换 夏巨伟 (浙江大学空间结构研究中心) 摘要：希尔伯特—黄变换由经验模态分解（empirical mode decomposition ，简称EMD ）和Hilbert 谱分析两部分组成。经验模态分解方法是一种自适应的、高效的数据分解方法。由于这种分解是以局部时间尺度为基础，因此，它适应于非线性、非平稳过程。通过经验模型分解，任何复杂的数据集都可以被分解为个数有限的、而且常常是为数不多的几个固有模函数(intrinsic mode functions ，简称IMF)的线性叠加。通过分解得到IMF 后，就可以对每一个分量做希尔伯特变换，得到其瞬时频率和幅度。本文详细对Hilbert-Huang Transform 的过程进行了阐述，并用算例分析指出了其优势所在。 关键词：希尔伯特—黄变换；时频分析技术； 1 希尔伯特—黄变换（Hilbert-Huang Transform ） 1.1 希尔伯特变换与瞬时频率（Hilbert Transform and instantaneous frequency ） 对于任意一个时间序列X(t)，它的希尔伯特变换具有如下形式： -1 ()(t)=,-X Y P d t ττπτ ∞∞? 其中，P ——积分的柯西主值； 希尔伯特变换对于任何属于L p 空间中的函数都存立，即上式中X(t)∈L p （— ∞，+∞）。 通过上述定义，X(t)和Y(t)成为一组复共轭对，同时能够构造一个实部和虚部分为X(t)和Y(t)的解析信号(Analytic Signal)Z(t)，Z(t)表示为： ()()(t)=(t)(t)=a ,i t Z X iY t e θ+ 其中， ()()1/222 (t)a =(t)+(t),arctan .X(t)Y t X Y t θ????= ????? 理论上讲有无数种方式去定义虚部，但是希尔伯特变换是唯一能够得到解析 信号结果的方法。 X(t)的Hilbert 变换实质上是将X(t)与函数1/t 在时域上做卷积，这就决定了通过X(t)的Hilbert 变换能够考察其局部特性。得到X(t)的瞬时相位函数后，其瞬时频率为： ()() (t).d w t dt θ= 1.2 经验模态分解与固有模态函数（Empiricalmode decomposition/EMD and Intrinsic mode function/IMF ） 固有模态函数需要满足两个条件：（1）极值与零点的数量必须相等或最多相差一个；（2）由局部极大值包络和局部极小值包络定义的平均包络曲线上任何一点的值为0；
山西中部引黄工程 施工监理Ⅶ标 工程变更管理监理实施细则 编号： ZBYH/Ⅶ/010 山西北龙工程监理有限公司 中部引黄工程施工监理Ⅶ标监理部 二○一二年七月
山西中部引黄工程 施工监理Ⅶ标 工程变更管理监理实施细则 编号： ZBYH/Ⅶ/010 批准：王仁龙 审核：陈继斌 编制：王仁龙丰巨武 关明杰杨光宇 山西北龙工程监理有限公司 中部引黄工程施工监理Ⅶ标监理部 二○一二年七月
目录 1目的与依据 (1) 1.1目的 (1) 1.2编制依据 (1) 2范围 . (1) 3监理职责 (1) 4工程变更条件与内容 (1) 5工程变更处理程序与要求 (2) 6工程变更的申报内容与要求 (5) 7工程变更的合同支付 (5)
山西省中部引黄工程施工监理Ⅶ标 工程变更监理实施细则 1目的与依据 1.1 目的 为规范山西省中部引黄工程施工监理Ⅶ标工程变更监理工作，使工程变更监 理工作程序化、规范化，编制本实施细则。 1.2 编制依据 （1）山西省中部引黄工程施工监理Ⅶ标建设监理合同(合同编 号:SXSZBYHGC-QT-JL-011 （ 2012）)； （2）山西省中部引黄工程施工监理Ⅶ标工程建设施工合同及其合同文件； （3）山西省中部引黄工程施工监理Ⅶ标工程设计文件；（4） 《水利工程建设项目施工监理规范》（SL288-2003）；（5） 山西省中部引黄工程施工监理Ⅶ标工程建设监理规划。 2适用范围 本实施细则适用于监理范围内山西省中部引黄工程施工监理Ⅶ标工程的工程 变更监理工作。 3监理职责 3.1 总监理工程师负责组织实施工程变更审查，签署合同文件授权范围内的 一切工程变更指令。 3.2 监理工程师在总监理工程师领导下完成工程变更监理具体工作。 4工程变更条件与内容 4.1 工程变更包括设计变更和施工变更，是指因设计条件、施工现场条件、 设计方案、施工方案发生变化，或项目法人与监理机构认为必要时，为实现合同 目标而对设计文件或施工状态所作出的变更与修改。 4.2 监理机构仅接受下列工程变更指示、通知或建议： (1)执行项目法人发出的工程变更指示。 (2)设计单位根据施工现场条件变化所提出的工程变更建议或设计修改通知。 (3)施工单位根据施工现场条件、施工方案或施工状态发生变化，依照合同文
云南睿城建设项目管理有限公司工程 变更管理办法及流程 第一条、目的 1、为了加强变更管理，规范工作流程，有效地控制成本，确保工程质量和工程进度，特制定本变更管理办法及流程。 2、通过对变更申报资料进行审查、审批，确保变更的及时性、合理性和经济性，消除变更对工程成本和进度带来的消极影响。 第二条、变更是对原设计内容进行完善、修改及优化，变更共分为三类： 1、一般变更：不改变设计原则，不影响使用功能，不影响工程的质量和安全，不影响美观；变更发生费用在2万元（含）以下的； 2、较大变更：不改变设计原则，不影响使用功能，不影响工程的质量和安全，不影响美观；变更发生费用在2万元至10万元（含）以下的； 3、重大变更：对原方案、原系统、主要结构布置、主要尺寸、坐标、主要标高、主要设备及主要使用功能改变及变更发生费用在10万元以上的。
第三条、变更的体现形式分为四类： 1、由建设单位(业主单位)提出的变更； 2、由监理单位提出的项目变更； 3、由设计单位提出的项目变更； 4、由施工单位提出的项目变更。 第四条对上述提出的工程变更，提出部门备齐相关原始资料，按本变更管理办法中图一及图二进行逐级上报审批。 第五条变更应将工程变更内容描述清楚。如：工程名称、变更原因、变更时间、变更部位、图纸比例、图示尺寸、规格型号、材料材质等，应达到根据变更单可准确计算工程量。 第六条变更单由项目部分专业依发生先后顺序进行编号。 第七条变更的控制 1、变更控制原则： 1.1 符合国家规范：变更应是对原设计中不满足国家规范、法规的部分进行变更，使之满足国家相关规范、法规； 1.2 保证使用功能：变更应是对原设计中不合理的部分进行变更，变更后应比原设计更合理、更满足使用功能；
（一）简单的EMD程序 function imf = emd(x) % Empiricial Mode Decomposition (Hilbert-Huang Transform) % imf = emd(x) % Func : findpeaks x = transpose(x(:));%转置 imf = []; while ~ismonotonic(x) %当x不是单调函数，分解终止条件 x1 = x; sd = Inf;%均值 %直到x1满足IMF条件，得c1 while (sd > 0.1) | ~isimf(x1) %当标准偏差系数sd大于0.1或x1不是固有模态函数时，分量终止条件 s1 = getspline(x1);%上包络线 s2 = -getspline(-x1);%下包络线 x2 = x1-(s1+s2)/2;%此处的x2为文章中的h sd = sum((x1-x2).^2)/sum(x1.^2); x1 = x2; end imf{end+1} = x1; x = x-x1; end imf{end+1} = x; % FUNCTIONS function u = ismonotonic(x) %u=0表示x不是单调函数，u=1表示x为单调的 u1 = length(findpeaks(x))*length(findpeaks(-x)); if u1 > 0, u = 0; else, u = 1; end function u = isimf(x) %u=0表示x不是固有模式函数，u=1表示x是固有模式函数 N = length(x); u1 = sum(x(1:N-1).*x(2:N) < 0); u2 = length(findpeaks(x))+length(findpeaks(-x)); if abs(u1-u2) > 1, u = 0; else, u = 1; end function s = getspline(x) %三次样条函数拟合成元数据包络线 N = length(x); p = findpeaks(x); s = spline([0 p N+1],[0 x(p) 0],1:N);
第34卷第2期福州大学学报(自然科学版)Vol.34No.2 2006年4月Journal of Fuzhou University(Natural Science)Apr.2006 文章编号:1000-2243(2006)02-0260-05希尔伯特-黄变换谱及其在地震信号分析中的应用 陈子雄,吴琛,周瑞忠 (福州大学土木建筑工程学院,福建福州350002) 摘要:介绍了希尔伯特-黄变换(HHT)这一非线性、非平稳信号处理方法,并利用HHT处理了地震工程中 常用的El Centro地震波,得到了该信号的Hilbert谱、边际谱和能量谱,提取了该信号的主要动力特性,并与 该信号的Fourier分析结果进行了对比,显示出HHT这一方法的优越性. 关键词:希尔伯特-黄变换;经验模态分解;固有模态函数;地震信号 中图分类号:TU311.3文献标识码:A Hilbert-Huang transform spectru m and its application in seismic signal analysis CHEN Zi-xiong,W U Chen,ZHOU Rui-zhong (College of Civil Engineering and Architecture,Fuzhou University,Fuzhou,Fujian350002,China) Abstract:HHT is a ne w method to deal with non-linear and non-stationary data.El Centro earth- quake wave is analyzed by HHT.Through the way,Hilbert spectrum,marginal spectrum and energy spec trum are got and dynamic property is extrac ted.The comparison between HHT spectrum and Fourier spec trum is made and the superiority of HHT is demonstrated. Keyw ords:Hilbert-Huang transform;empirical mode decomposition;intrinsic mode function;seismic signal 地震信号具有短时、突变等特点,是一种典型的非平稳随机信号,必须对其进行分析与处理,才可以提取信号的主要特征.传统的Fourie r变换能够表述信号的频率特性,但不提供任何时域信息[1],而小波分析虽然在时域和频域都具有很好的局部化性质,但本质上仍是一种窗口可调的Fourier变换,在小波窗内的信号必须是平稳的,因而没有根本摆脱Fourier分析的局限[2].小波基的选择也是信号分析中的一个重要问题,另外,小波基的有限长会造成信号能量的泄漏,使信号的能量-频率-时间分布很难定量表述. Hilbert-Huang变换(HH T)的信号处理方法被认为是近年来对以Fourier变换为基础对线性和稳态谱分析的一个重大突破[2].它由经验模态分解(E mpirical Mode Decomposition,E MD)方法和Hilbert变换(H T)两部分组成,其核心是E MD分解.该方法采用了固有模态函数(Intrinsic Mode Function,I MF)概念以及将任意信号分解为I MF组成的思想,即E MD法,使得瞬时频率具有实际的物理意义[3].它不受Fourier分析的局限,可依据数据本身的时间尺度特征进行模态分解,分解过程中保留了数据本身的特性,再对各I MF分量进行Hilbert变换,得到信号能量在时间尺度上的分布规律,实现地震动力特性的提取. 1Hilbert-Huang变换 1.1经验模态分解和固有模态函数 经验模态分解(EMD)的目的是通过对非线性非平稳信号的分解获得一系列表征信号特征时间尺度的固有模态函数(I MF),使得各个I MF是窄带信号,可以进行Hilbert分析.首先设定两个条件:1整个时间序列的极大极小值数目与过零点数目相等或最多相差一个;o时间序列的任意点上,由极大值确 收稿日期:2005-07-27 作者简介:陈子雄(1981-),男,硕士研究生;通讯联系人:周瑞忠,教授. 基金项目:教育部博士点专项科研基金资助项目(20040386004)
XXXXXXXXX公司 工程变更管理办法 第一章总则 第一条为规范XXXXXXXXX公司的工程变更管理工作，依据政府有关法规，结合工程建设管理实际情况，制定本办法。 第二条本办法所称工程变更，系指工程签订合同至工程正式验收前，变更设计图纸或变更工程现场工程施工的活动。经确定的设计施工图，任何单位和个人不得擅自变更，确需变更的，依据本办法进行。 第三条本办法适用于公司及下属企业单个合同金额在50万元以上的工程项目变更管理。 公司是指XXXXXXXXX公司。下属企业是指公司出资的全资、控股企业以及拥有实际控制权的各级企业。项目公司是指具体负责建设工程项目的下属企业。 第四条本办法中工程变更限额范围为累计工程变更价款在单个合同总价的0%－10%（含）之间，超出10%以上的需重新申报追加合同金额并上报公司招标领导小组审批。 第五条工程变更应从有利于工程出发，或加快竣工，或降低费用，或提高工程效率和价值、或提高安全性以及满足要求。 第六条审核工程变更应以合同为准绳、以事实为依据，从技术、经济、法律等方面进行全面审核。 第七条变更指令是实施变更的唯一途径，有指令才能实施变更，
无指令不能实施任何变更。 第八条监理单位和跟踪审计控制工程变更的外部管理。监理单位接受工程变更申报，全面审核工程变更，确定变更的原因和责任，发布变更指令，监督变更的实施，做好变更的管理台账；跟踪审计负责审核工程变更的真实性、合理性、必要性，并审核与确定承包商提交的工程变更价款。 第九条项目公司控制工程变更的内部管理，负责资料完善和传递，组织相关单位、部门洽商、初审、审核审查论证，负责跟踪、检查、督办、监督变更实施，建立变更管理内部台帐等。 第十条如工程项目未配备监理单位，由项目公司发布变更指令，并履行监理单位的变更管理职责。 第十一条工程变更必须严格坚持“先批准，后变更；先设计，后施工”的原则。承包商未经公司批准而擅自变更和施工的，自行承担全部责任，并不予验收和计量；项目公司未按本办法规定的程序报经批准而擅自发出变更指令的，应按公司相关制度追究责任。 第十二条对现场应急、抢险、突发的、急需的、紧急或特殊情况下的重大变更，项目公司可事先口头或书面向工程技术部报告，经同意后先行紧急变更，事后及时完备图纸等工程变更资料，并补办审批手续。 第十三条工程变更的各项管理工作应确定具体的时限要求，各部门都应遵照执行并承担超限责任。 第二章工程变更管理
用希尔伯特黄变换（HHT）求时频谱和边际谱 1.什么是HHT？ HHT就是先将信号进行经验模态分解（EMD分解），然后将分解后的每个IMF分量进行Hilbert变换，得到信号的时频属性的一种时频分析方法。 2.EMD分解的步骤。
EMD分解的流程图如下：
3.实例演示。 给定频率分别为10Hz和35Hz的两个正弦信号相叠加的复合信号，采样频率fs=2048Hz的信号，表达式如下：y=5sin(2*pi*10t)+5*sin(2*pi*35t) (1)为了对比，先用fft对求上述信号的幅频和相频曲线。 代码: function fftfenxi clear;clc; N=2048; %fft默认计算的信号是从0开始的
t=linspace(1,2,N);deta=t(2)-t(1);1/deta x=5*sin(2*pi*10*t)+5*sin(2*pi*35*t); % N1=256;N2=512;w1=0.2*2*pi;w2=0.3*2*pi;w3=0.4*2*pi; % x=(t>=-200&t<=-200+N1*deta).*sin(w1*t)+(t>-200+N1*deta&t<=-200+N2*det a).*sin(w2*t)+(t>-200+N2*deta&t<=200).*sin(w3*t); y = x; m=0:N-1; f=1./(N*deta)*m;%可以查看课本就是这样定义横坐标频率范围的 %下面计算的Y就是x(t)的傅里叶变换数值 %Y=exp(i*4*pi*f).*fft(y)%将计算出来的频谱乘以exp(i*4*pi*f)得到频移后[-2,2]之间的频谱值 Y=fft(y); z=sqrt(Y.*conj(Y)); plot(f(1:100),z(1:100)); title('幅频曲线') xiangwei=angle(Y); figure(2) plot(f,xiangwei) title('相频曲线') figure(3) plot(t,y,'r') %axis([-2,2,0,1.2]) title('原始信号')
工程变更控制程序 （IATF16949-2016/ISO9001-2015） 1.0目的 明确工程变更过程、各相关部门的职责和工作流程，使本公司各项工程变更在受控状态下进行，以确保变更后产品的质量能够满足客户需求和降低工程变更造成的质量风险。 2.0适用范围 本程序文件适用于制造过程中所有工程变更。 3.0职责 3.1生管部负责涉及产品结构、物料规格及工艺标准等方面工程变更的提出和编制，评估并验证其可行性。 3.2生产部负责与生产作业方式或方法相关的工程变更的提出和编制，评估并验证其可行性。 3.3 行政部保证工程变更单的受控下发、登记、文件保存。 3.4 生产部负责实施工程变更，必要时需对工程变更实施后的产品进行标识。 3.5 品质部负责对工程变更实施监督并对变更实施的效果进行现场验证及相关反馈。
3.6 生管部负责与客户联络，将客户的要求传达公司内部。 4.0定义 4.1工程变更：指因生产过程、产品品质及客户需求等原因而对公司目前采 用的生产工艺、生产设备、生产线变更或转移、产品结构、物料规格、标准 规范等方面实施的变更。不仅包括公司内的相应变更，也包括供应商及其重 要分供方的相应变更。 4.2 ECN：Engineering Change Notice的缩写，即工程变更通知。 4.3 变更等级： 4.3.1 一级变更：变更可能会影响产品安全性能。 4.3.2 二级变更：变更可能会影响产品基本性能。 4.3.3 三级变更：不会影响产品基本性能。 4.4 风险工程变更：没有时间完成必要的验证，但为了产品出货需要，必须 马上实施的工程变更。 4.5 临时工程变更：产品生产过程中发生的制造过程、过程控制、材料控制 等各种行为的临时性的工程变更。 4.6 紧急过程更改通知单 AVO：Avoid Verbal Orders. 紧急情况（夜间或节假日，并且工程变更签批负责人不在公司时）下过程发生
希尔伯特-黄变换（Hilbert-Huang Transform，HHT） 0 前言 传统的数据分析方法都是基于线性和平稳信号的假设，然而对实际系统，无论是自然的还是人为建立的，数据最有可能是非线性、非平稳的。 希尔伯特-黄变换（Hilbert-Huang Transform，HHT）是一种经验数据分析方法，其扩展是自适应性的，所以它可以描述非线性、非平稳过程数据的物理意义。 1 HHT简介[贺礼平.希尔伯特-黄变换在电力谐波分析中的应用研究[D].湖南:中南大学,2009]HHT的发展。 1995年，Norden E.Huang为研究水表面波构思出一种所谓“EMD--HSA”的时间序列分析法，通过这种方法他发现水波的演化不是连续的，而是突变、离散、局部的。 1998年，Norden E.Huang等人提出了经验模态分解方法，并引入了Hilbert谱的概念和Hilbert谱分析的方法，美国国家航空和宇航局（NASA）将这一方法命名为Hilbert-Huang Transform，简称HHT，即希尔伯特-黄变换。 HHT是一种新的分析非线性非平稳信号的时频分析方法，由两部分组成： 第一部分为经验模态分解（Empirical Mode Decomposition，EMD）（the sifting process，筛选过程），它是由Huang提出的,基于一个假设：任何复杂信号都可以分解为有限数目且具有一定物理定义的固有模态函数（Intrinsic Mode Function，IMF；也称作本征模态函数）；EMD方法能根据信号的特点，自适应地将信号分解成从高到低不同频率的一系列IMF；该方法直接从信号本身获取基函数，因此具有自适应性，同时也存在计算量大和模态混叠的缺点。 第二部分为Hilbert谱分析（Hilbert Spectrum Analysis，HSA），利用Hilbert变换求解每一阶IMF 的瞬时频率，从而得到信号的时频表示，即Hilbert谱。 简单说来，HHT处理非平稳信号的基本过程是：首先，利用EMD方法将给定的信号分解为若干IMF，这些IMF是满足一定条件的分量；然后，对每一个IMF进行Hilbert变换，得到相应的Hilbert谱，即将每个IMF表示在联合的时频域中；最后，汇总所有IMF的Hilbert谱就会得到原始信号的时间-频率-能量分布，即Hilbert谱。 在HHT中，为了能把复杂的信号分解为简单的单分量信号的组合，在进行EMD方法时，所获得的IMF 必须满足下列两个条件： 1）在整个信号长度上，一个IMF的极值点和过零点数目必须相等或至多只相差一点。 2）在任意时刻，由极大值点定义的上包络线和由极小值点定义的下包络线的平均值为零，也就是说IMF的上下包络线对称于时间轴。
1.0目的： 明确工程更改的管理程序及步骤,用以管理客户要求的产品更改和内部的产品及过程更改,确保产品的工程更改水平得到有效控制。 2.0范围: 适用于本公司的产品和过程更改。 3.0职责权限 3.1项目部是本程序归口管理部门，负责所有工程更改的实施和协调。 3.2业务部负责顾客工程更改的内部信息传递和工程更改后顾客的批准事宜。 3.3相关部门负责提出更改的申请和配合工程更改的实施。 3.4品保部负责工程更改资料的登记分发及保存。 4.0定义
6.1变更的提出: 6.1.1业务部及相关部门在接获客户方提出的产品方面变更和量产时接获顾客品质改善的 要求，由业务部提出顾客的格式变更联络/通知或开立「顾客设变通知单」并附上设变资料转交开发部。 6.1.2对公司内部产品方面和过程方面的变更需求，由提出部门填写「变更申请单」和相 关资料通知开发部。 6.1.3供应商提出的变更由供应商或开发部填写「变更申请单」。
6.1.4产品设变时涉及的数模版本由开发进行控制，具体见「外来文件及数据管理程序」6.2可行性评审： 6.2.1开发部在接到「顾客设变通知单」和「变更申请单」时需根据变更内容发联络单通 知技术部、质量部、制造厂、模具厂必要时包括外协厂商进行可行性评审，结果记 录在「变更评审表」上。 6.2.2开发部在接到「顾客设变通知单」和「变更申请单」时，需确认变更对现有产品的 影响，包括库存品和原料、模具等，并组织技术部对需要变更的产品进行比对，确 立变更方案并估计成本增加数额，以及变更对现有文件的影响和确定费用承担的责 任方。 6.2.3对变更需要顾客承担的费用，由业务部按「产品报价作业办法」重新向客户提出报 价，由客户确认后实施。 6.3变更的准备： 6.3.1变更涉及到模具工装的变更由技术部针对产品变更或模具品质改进对于模具结构加 工方式进行评估，完成后确定模具变更方案按照「设计图设计变更管理办法」执行 修订。 6.3.2产品图纸的变更由技术部按照「工程文件管理程序」执行工程文件的更改。 6.4变更的实施： 6.4.1模具产品的设计变更由技术部发出「设计变更通知单（模具）」、零部件产品的设计变 更由工程部发出「设计变更通知单（部品）」通知相关部门设变的实施。 6.4.2设变执行时开发需要发出「4M变更记录表」以控制设变方案得到有效执行。各相关部 门根据「设计变更通知单（模具/部品）」的方案实施，各部门对库存品处理，模具修 改，文件修改等必须记录在「4M变更记录表」上。 6.4.3过程4M变更中列“A”项者为影响品质较重要者，在变更前需经各部门评审并经开发 部主管核准后才可执行变更作业。 6.4.4对内部提出的变更,会影响客户专用权的设计、外形、配合和功能（包括性能和/或耐 久度）时，必须由业务向客户提出申请，核准后才可执行变更作业。 6.5变更的验证和确认： 6.5.1开发部组织生产部门对变更品进行试制。品保部在接获「4M变更记录表」时依厂内 试作规定排定的计划时程，对变更的产品品质确认。 6.5.2需重新提交PPAP文件的情况执行「生产件批准程序」的相关要求。 6.5.3产品属B项者需经厂内质量部进行品质确认且合格后，即可准备量产；但须保存各 项记录文件，待客户相关人员监查时，主动提出接受查验。 6.5.4由工程变更后之各相关文件、制程作业标准书、质量管理计划表、检查基准书------- 等，因工程变更而需修订者，应于量产前修订完成，并根据「工程文件管理程序」 的相关要求执行，并将变更的相关文件登记在「产品履历表」上。 6.6变更前产品的处理： 6.6.1相关部门在接到「4M变更记录表」时，制造及相关部门对变更前的产品视同不合格 品进行标识并按照「不合格品管理程序」的相关要求执行，结果记录在「4M变更记 录表」上。 6.6.2变更的产品依「产品标识与追溯管理程序」的标示使用办法进行标识识别，产品属 〝B〞项者可免标示。 6.6.3分包商工程变更，不论是工程变更定义之A类或B类，均需提出申请，经品保部确 认及其部门主管核准后，方可量产交货。 6.6.3客户有特殊要求时，依客户的要求执行。 6.8变更的记录按照「记录管理程序」的相关要求进行归档储存。 7.0支持性文件: 7.1工程文件管理程序（文件编号：） 7.2记录管理程序（文件编号：） 7.3生产件批准程序（文件编号：）
工程设计变更/技术核定和现场签证 管理流程 一、目的 为了加强在建项目工程设计变更、技术核定及现场签证的管理，规范其工作管理流程，确保工程质量和工程进度，有效的控制工程成本，特制定本管理办法。 二、适用范围 本管理流程适用于公司在建各项目的工程设计变更、技术核定及现场签证的管理。 三、工程设计变更、技术核定及现场签证的实施原则 1．时间限制原则：公司对工程设计变更、技术核定、现场签证事项以及其结算实行严格的时间限制，保持与实施进度同步，严禁事后补办。 2．一单一算原则：工程设计变更单、技术核定单或现场签证单应遵循一单一算原则，且对应于相应的工程合同。 3．量价分离原则：工程设计变更单、技术核定单及现场签证单所涉及项目工程由项目工程部主管责任工程师确认完成状况和发生的工程量，由专业造价工程师负责复核工程量、确认其单价和造价变量额。 4．先估算后施工原则：现场签证在施工前应先由主管造价工程师明确估价，按相应审批权限审批执行。 5．一月一清原则：每月26日前，造价部就截止至本月25日已完成且手续完备的工程设计变更单、技术核定单及现场签证单进行汇总，建立台账。 6. 完工确认原则：当工程设计变更单、技术核定单及现场签证单涉及工程完工后，主管责任工程师和造价工程师、监理单位、施工单位必须在完工后5日内共同签
字确认，如属隐蔽工程，必须在其隐蔽前签字确认。 7．原件结算原则：工程设计变更单、技术核定单及现场签证的结算必须要有齐备的、有效的资料原件作为结算依据。 8．标准表格原则：所有的工程设计变更单、技术核定单及现场签证都必须使用规定的标准表格，统一编号。 9．权利限制原则：对工程设计变更单、技术核定单及现场签证管理实行严格的权限规定，严禁越权签字。 10．签证避让原则：项目工程部或产品研发中心、设计院等可以用工程变更或技术核定单表示清楚和量化的工程项目内容，不得再现场签证。 11．坚持优化原则：工程变更、技术核定单及现场签证必须执行严格审批、规范程序的原则，坚持有利于工程项目在功能、质量、投资方面实施优化的原则。 四、工程设计变更及技术核定单的管理审批流程 一> 工程变更产生的原因及种类 1.修改工艺技术，包括设备的改变； 2.增减工程内容； 3.改变使用功能； 4.设计错误、遗漏； 5.合理化建议，在保证产品品质的条件下为节约工程造价发生的变更； 6.施工中产生错误； 7.使用的材料品种的改变； 8.工程地质勘察资料不准确而引起的修改，如基础加深等； 9.为满足客户要求的变更；
HHT-希尔伯特·黄变换 1998年，Norden E. Huang等人提出了经验模态分解方法，并引入了Hilbert谱的概念和Hilbert谱分析的方法，美国国家航空和宇航局（NASA）将这一方法命名为Hilbert-Huang Transform，简称HHT，即希尔伯特-黄变换。 HHT主要内容包含两部分，第一部分为经验模态分解（Empirical Mode Decomposition，简称EMD），它是由Huang提出的；第二部分为Hilbert谱分析（Hilbert Spectrum Analysis，简称HAS）。简单说来，HHT处理非平稳信号的基本过程是：首先利用EMD方法将给定的信号分解为若干固有模态函数（以Intrinsic Mode Function或IMF 表示，也称作本征模态函数），这些IMF是满足一定条件的分量；然后，对每一个 IMF进行Hilbert变换，得到相应的Hilbert谱，即将每个IMF表示在联合的时频域中；最后，汇总所有IMF的Hilbert谱就会得到原始信号的Hilbert谱。 与传统的信号或数据处理方法相比，HHT具有如下特点： (1)HHT能分析非线性非平稳信号。 传统的数据处理方法，如傅立叶变换只能处理线性非平稳的信号，小波变换虽然在理 论上能处理非线性非平稳信号，但在实际算法实现中却只能处理线性非平稳信号。历 史上还出现过不少信号处理方法，然而它们不是受线性束缚，就是受平稳性束缚，并 不能完全意义上处理非线性非平稳信号。HHT则不同于这些传统方法，它彻底摆脱了 线性和平稳性束缚，其适用于分析非线性非平稳信号。 (2)HHT具有完全自适应性。 HHT能够自适应产生“基”，即由“筛选”过程产生的IMF。这点不同于傅立叶变换 和小波变换。傅立叶变换的基是三角函数，小波变换的基是满足“可容性条件”的小 波基，小波基也是预先选定的。在实际工程中，如何选择小波基不是一件容易的事， 选择不同的小波基可能产生不同的处理结果。我们也没有理由认为所选的小波基能够 反映被分析数据或信号的特性。 (3)HHT不受Heisenberg测不准原理制约——适合突变信号。 傅立叶变换、短时傅立叶变换、小波变换都受Heisenberg测不准原理制约，即时间窗口与频率窗口的乘积为一个常数。这就意味着如果要提高时间精度就得牺牲频率精度，反之亦然，故不能在时间和频率同时达到很高的精度，这就给信号分析处理带来一定 的不便。而HHT不受Heisenberg测不准原理制约，它可以在时间和频率同时达到很 高的精度，这使它非常适用于分析突变信号。 (4)HHT的瞬时频率是采用求导得到的。