专利名称：基于阶跃响应的传感器动态误差频域修正技术的制作方法
技术领域：
本发明涉及传感器动态误差的频域修正技术，特别是一种适用于根据传感器阶跃 响应实验数据求取传感器动态误差频域修正函数并以此对传感器动态响应进行连续频域 修正的技术，解决了现有频域修正技术在对传感器阶跃响应实验数据进行动态误差修正时 存在的问题。
背景技术：
动态误差是动态测试中首先必须解决的一个问题。测量系统无失真动态测试的条 件是，测量系统的传递函数等于一个放大环节的传递函数K，也就是等于静态传递函数。换 句话说，测量系统各个环节均无惯性与阻尼等。在这种情况下，测量系统的频率特性曲线是 一条与横坐标平行的直线。但是，在实际中传感器以及整个测量系统的频率响应总是有限 的，即不可能是一条与频率轴平行的直线，特别是传感器的频率响应范围往往成为限制整 个测量系统频率响应的主要环节。由于传感器幅频特性的不平坦，因而被测信号中的各次 谐波，有的被放大，有的被衰减。由于传感器相频特性不是理想的直线，因而各次谐波之间 的相位差也改变了。这些使得传感器输出信号与输入信号之间存在着畸变。这就是所谓的 “动态误差”。它与静态测试中的系统误差不同。传感器动态误差无法用一个系数去修正。 这是因为整个响应曲线都发生了畸变，并且与被测信号的频谱有关。黄俊钦简要介绍了三种传感器动态误差的修正方法频率域修正法、数值微分法 和叠加积分法(黄俊钦，测试系统动力学，北京国防工业出版社，1996年，pp. 51-55，第2 章第6节动态误差及无失真动态测试)。徐科军介绍了传感器动态误差修正方法，包括 频域修正法、数字微分法、叠加积分法、反滤波动态误差修正法、实时数字微分修正法，并对 频域修正方法进行了改进，还做到了频域修正方法的在线实现(徐科军，传感器动态特性 的实用研究方法，合肥中国科学技术大学出版社，1999年，pp. 101-112，第7章动态误差 修正方法)。传感器动态误差频域修正的基本步骤是(1)对传感器进行动态标定实验， 得到传感器输入和输出的实验数据；( 对传感器的输入和输出数据做快速傅里叶变换， 得到输入和输出信号的频谱，再据此求取传感器的频率响应函数；C3)根据传感器的频率 响应函数，构造相应的修正传递函数；(4)在实际测量中，得到传感器的动态响应的输出数 据后，对输出数据做傅里叶变换；( 将变换结果与修正传递函数相乘，再进行傅里叶反变 换，就得到修正后的传感器输出信号。这个经过修正后的传感器输出信号，与原来的输出信 号相比，将更接近于输入信号(被测信号)，即较为准确地反映被测信号的情况。将现有的传感器动态误差频域修正的技术应用于传感器脉冲响应的动态误差修 正，可以取得较好的效果。但是，脉冲响应实验的幅值和落点控制精度不高、重复性较差。 而阶跃响应实验具有较高的幅值控制精度和较好的重复性，且其能标定的频带更宽，还适 用于多维力传感器的力矩实验，因而在传感器的动态标定实验中应用较多。然而，将现有的 传感器动态误差频域修正技术应用于阶跃响应的动态误差修正，在求解动态误差频域修正 函数时，出现了问题。具体表现为对于同一个传感器，根据一组传感器阶跃响应实验数据
4求取的传感器动态误差频域修正函数，用于另外一组阶跃响应实验数据的误差修正时，无 法得到满意的结果，有时甚至是错误的结果。而这两组数据可能仅仅是阶跃响应的起跳点 (起跳时刻)不同。这样，求取的传感器动态误差频域修正函数就没有普适性，也就无法应 用于实际传感器动态响应的误差修正。因为，在实际应用中，传感器实际响应数据与其动态 阶跃响应实验数据一般都是不一样的。针对上述问题深入分析，发现原因如下一、对传感器阶跃响应实验中的输出数据 进行基于FFT的频谱分析时，不可能也没有必要采集无限长数据，一般是采集一段数据进 行频谱分析。由于FFT机理的原因，在做频谱分析时，相当于对实验数据进行了周期延拓， 即将原来由低到高(或者由高到低)的一个阶跃波形视为一个矩形波(由低到高，再由高 到低；或者由高到低，再由低到高)。而矩形波的频谱中会出现零值，即它的幅频特性曲线 中有的频率点处的幅值为零。这样，动态误差频域修正函数就会出现无效值(这里有两种 情况，一是零比零，即分子和分母均为零；二是零值)。在要修正的频带范围内由于无效值 的存在使得动态误差频域修正函数有效修正带宽缩小为无效值对应的最小频率以内的范 围，所以，该动态误差频域修正函数无法对其它阶跃响应数据或实际应用中的动态响应进 行正确的动态误差修正。二、原本的阶跃激励信号在做基于FFT的频谱分析时，由于周期延 拓被视为矩形波信号，而响应输出信号的周期延拓又与矩形波激励的传感器实际响应输出 不符。所以，根据传感器阶跃响应实验数据直接求取动态误差频域修正函数会引入较大的 误差，无法用于传感器实际动态响应的误差修正。另外，在获得传感器动态误差频域修正函数之后，现有频域修正技术对传感器的 实际响应输出数据进行误差修正时，由于截取数据的窗函数泄漏效应导致每次修正后数据 段的两端的修正误差较大，从而造成对传感器动态响应进行连续的误差修正时误差较大。

发明内容
本发明要解决现有传感器动态误差频域修正技术无法直接应用于动态阶跃响应 实验数据的误差修正问题，提供一种有效的动态误差频域修正函数的求解方法以及对传感 器动态响应输出数据进行连续的动态误差频域修正方法。本发明所采用的技术方案是先对传感器进行动态阶跃响应实验，采集传感器的 激励(输入)数据和响应(输出)数据；再在所采集的输入数据和输出数据中分别截取待 研究的阶跃激励和阶跃响应过程的数据段，并保证一定的数据长度和阶跃高低占空比，为 消除频域修正函数求解结果中的无效值点提供条件；接着，对截取的阶跃响应数据段的前 端数据进行拼接处理，形成单周期矩形波激励信号对应的传感器响应信号(即输入信号为 一个单周期矩形波，输出信号为一个既有正阶跃响应又有负阶跃响应的信号)，以大幅度地 消除阶跃激励及其响应在加窗后进行频谱分析时由于周期延拓引起的估计误差；然后，分 别对截取的阶跃激励数据和拼接处理后的阶跃响应数据加特定的窗函数(矩形窗或海明 窗)，为消除频域修正函数求解结果中的无效值点提供条件；最后，对加特定窗后的阶跃激 励和阶跃响应数据，采用功率谱估计法求取频域修正函数，插补其要修正频带范围内的无 效值点，并将其要修正频带范围外的函数值置零，得到最终的动态误差频域修正函数。采用 连续、滑动的误差修正方式，即采用对修正结果的数据片段连续拼接的方式，对传感器实际 响应输出进行动态误差修正，以取得良好的修正效果。
本发明的技术流程为动态阶跃响应实验1 —数据截取2 —数据拼接处理3 —信 号加窗4 —动态误差频域修正函数求解5 —传感器动态响应的误差修正6，如图1所示。所述图1中的动态阶跃响应实验1设备主要由阶跃激励源7、传感器8、数据采集 仪9、上位机10组成，如图3所示。阶跃激励源7产生阶跃激励信号x(t)，激励信号x(t) 和传感器8的响应信号y (t)由数据采集仪9以采样率fs同步采集，并传送给上位机10进 行存储。所述数据截取2即为采用时间长度为T(T = N/fs，N为截取的数据点数)的矩形 窗对传感器阶跃激励数据X(t)和响应数据y(t)从同一时刻开始进行同步截取，分别得到 数据h*yN。截取得到的数据&和71^的占空比σ决定了数据中的零谱值分量的频率点。所述数据拼接处理3是将矩形窗所截取的数据yN中反映阶跃响应动态过程的数 据段&进行倒相并替代yN数据前端与&时间长度相同的S1段数据，进而得到拼接处理后 的数据y/。所述信号加窗4是将截取的激励数据&和拼接处理后的响应数据yN'分别乘以相 同的窗函数W( ·)，进而得到加窗的数据和λν’。窗函数选择矩形窗或海明窗(Hamming)。所述动态误差频域修正函数求解5通过计算％和^讷自功率谱Gxx (f)、Gyy (f)及 其互功率谱Gxy (f)、Gyx (f)，然后，据此计算出有偏估计频域修正函数钱CO和应2(/)，并在要 修正的频带
内分别插补钱(/)和式(/)中的无效值点，得到HJf)和吐⑴，对！^⑴ 和H2 (f)相加和取平均，得到频域修正函数々(/)。当成(/)和应2(/)在要修正的频带
之外的频率点函数值置零，得到最终的传感器动态误差 频域修正函数H (f)。所述传感器动态响应的误差修正6采取对传感器的动态响应进行连续、滑动的动 态误差频域修正方式，每次修正数据段i~N(i)为采用矩形窗对传感器的动态响应r(t)进行 滑动截取所得，矩形窗每次滑动时间τ小于ι·Ν(υ的时间长度；再对ι·Ν(υ进行傅里叶变换 得R(f)；然后，将R(f)与所求的动态误差频域修正函数H(f)进行点乘得修正后的频域数 据c(f)；接着，对C(f)进行逆傅里叶变换，得到修正后的数据cN(i)；取每次滑动截取的数 据段的修正结果％(1)的中间段时间长度为τ的数据gcji)进行拼接，得到动态误差修 正后的反映实际被测量的真实值。本发明方法的优点是能够直接针对幅值控制精度高、标定频带宽的传感器动态 阶跃响应实验数据求取有效修正频带宽、动态误差修正精度高的传感器动态误差频域修正 函数，实现对传感器实际动态响应进行连续的动态误差修正。


图1是本发明技术方案流程框图，即基于阶跃响应的传感器动态误差频域修正技 术方案流程框图；图2是本发明具体实施例一的技术方案流程框图；图3是本发明具体实施例中动态标定实验设备组成框图；图4是本发明具体实施例中杆式风洞应变天平阶跃响应实验示意图5是本发明具体实施例中数据截取示意图；图6是本发明具体实施例中数据拼接处理示意图；图7是本发明具体实施例一中动态误差频域修正函数求解流程框图；图8是本发明具体实施例中传感器频域修正前后的频率特性曲线示意图；图9是本发明具体实施例中传感器动态误差的频域修正流程框图；图10是本发明具体实施例中连续、滑动的修正传感器动态响应误差的示意图；图11是本发明具体实施例二的技术方案流程框图；图12是本发明具体实施例三的技术方案流程框图；图13是本发明具体实施例二 /三中动态误差频域修正函数求解流程框图。
具体实施例方式下面结合附图对本发明做进一步说明本发明的设计思想是基于传感器动态阶跃响应实验数据，分别采取频域修正函 数零点插补、利用特定窗函数的泄漏特性、设定截取的阶跃激励数据的占空比三种方案来 克服动态误差频域修正函数求解中存在无效值的问题，采取数据拼接处理的方式来克服由 于阶跃激励和传感器响应截取造成的周期延拓效应引起的误差，再通过对阶跃激励和拼接 处理后的响应信号进行功率谱估计求取两个有偏估计频域修正函数，然后将这两个有偏估 计的频域修正函数进行相加和求平均，并将要修正频带范围之外的函数值置零，最终得到 动态误差频域修正函数；在实际测量中，采用该动态误差频域修正函数对传感器实际响应 数据进行连续、滑动的动态误差修正。本发明专利以杆式风洞应变天平为例。杆式风洞应变天平是一种6维力传感器， 是风洞实验中测量飞行器模型所受空气动力/力矩的重要设备，为飞行器的飞行角度、流 线型外形设计、气流对飞行控制的影响等诸多方面的研究提供重要的数据资料。但是，由于 它阻尼比过小、固有频率较低，从而动态响应速度慢、阶跃响应超调量大、调节时间长，即动 态误差大，无法满足动态测试的要求。本发明技术方案流程框图如图1所示。动态阶跃响应实验1所采集的实验数据 依次通过数据截取2、数据拼接处理3、信号加窗4和动态误差频域修正函数求解5求得传 感器的动态误差频域修正函数，再采用该频域修正函数，通过传感器动态响应的误差修正6 实现对传感器实际动态响应进行连续的动态误差修正。依据传感器动态误差频域修正函数求解过程中修正频带内无效值问题的克服方 案的不同，本发明方法提出三种具体实施例。实施例一图2所示为本发明方法的第一种具体实施例(简称为具体实施例一)的技术方案 流程框图。进行动态阶跃响应实验1的动态标定实验系统组成框图如图3所示，包括阶跃激 励源7、传感器8、数据采集仪9及上位机10。阶跃激励源7为传感器8提供阶跃激励信号 X (t)，传感器8在阶跃激励信号X (t)的激励下，产生响应信号y (t)，数据采集仪9以采样频 率仁同步采集X(t)和y(t)，并上传至上位机10进行存储。具体实施例中杆式风洞应变天 平阶跃响应实验示意图如图4所示。杆式风洞应变天平安装在加载台上，钢丝挂上砝码绕过定位滑轮连接至杆式风洞应变天平，对其某一方向施加一稳定的载荷，然后通过快速剪 断钢丝的方法进行卸载，从而实现对杆式风洞应变天平施加负阶跃载荷。实验中，钢丝剪断 点设置在靠近杆式风洞应变天平处。数据截取2即为采用时间长度为T(T = N/fs, N为截取的数据点数)的矩形窗对 实验数据X(t)和y(t)进行同步截取。数据截取示意图如图5所示。截取的数据段需分别 包含阶跃激励信号x(t)和响应信号y(t)的动态过程。数据x(t)和y(t)的截取保证同 步，即两组数据截取的起始时刻和结束时刻相同。截取的数据段占空比σ均为 0. 5(即50% )，以使&的周期延拓为周期为T、频率Af为1/T的方波。从而，&的频谱中 频率为Af的偶次倍频的频率分量为零。Af为对做傅里叶变换时的频率分辨率。数据拼接处理3的拼接处理示意图如图6所示，将所述数据截取2所截取的数据 yN中反映阶跃响应动态过程的数据段&进行倒相并替代yN数据前端与&时间长度相同的 S1S数据，进而得到拼接处理后的数据y/。图6 (a)所示为数据拼接处理前的数据7,，图 6(b)所示为数据拼接处理后的数据y/。yN'中的数据段S1计算如下S1 = -S2+A式中，A为阶跃响应yN的静态阶跃幅值。信号加窗4采用N(N = T · fs)点矩形窗函数点乘所述数据截取2截取的数据％ 和所述拼接处理3处理得到的数据y/，得加窗后的数据％和5^ ’。动态误差频域修正函数求解5的求解流程框图如图7所示。求解步骤为步骤一对信号加窗4加窗处理后的数据乓和扎’做傅里叶变换，得X(f)和Y(f)。步骤二 依据X (f)和Y (f)求取数据％和5^ ’的自功率谱Gxx (f)、Gyy (f)和互功率 tb (;xy(f)、Gyx(f)。步骤三求有偏估计频域修正函数钱(/) ,H2(f)。步骤四对钱(/)和应2(/)中在要修正频带
内由于％的频谱中频率为Δ f 的偶次倍频的频率分量为零导致的无效值频率点分别进行插补得HJf)和吐江)。H^f)和 H2 (f)计算式如下
H^n-Af) ^ H,(n-Af), η为奇数
!//,(π·Δ/) = [^1 ((η-1).Δ/) + H1 {{η +1)·Δ/)]/2，η%偶数
//2( ·Δ/) = ^2( ·Δ/),"为奇数
{Η2(η·Δ/) = [H2((η-1)·Δ/) + H2{{η +1)·Δ/)]/2，η为偶数步骤五对所求的H1 (f)和H2 (f)进行相加和取平均，得到频域修正函数应(/)。步骤六将频域修正函数应(/)中位于要修正频带W，fmax]之外的频率点的函数值 全部置零，得到所要求的传感器动态误差频域修正函数H(f)。图8所示为传感器频域修正前后的频率特性曲线示意图。图8(a)为幅频增益特 性曲线，IA代表原传感器的频响函数的幅频增益特性，2A代表所述动态误差频域修正函数 求解5所求得的传感器动态误差频域修正函数H(f)的幅频增益特性，3A代表传感器幅频增 益特性IA经H(f)修正后的幅频增益特性。图8(b)为相移频率特性曲线，IB代表原传感器的频响函数的相移频率特性，2B代表所述动态误差频域修正函数求解5所求得的传感器动 态误差频域修正函数H(f)的相移频率特性，:3B代表传感器相移频率特性IB经H(f)修正后 的相移频率特性。修正后，传感器的频带被拓宽到fmax。传感器动态响应的误差修正6采用所述动态误差频域修正函数求解5求取的传感 器动态误差频域修正函数H(f)对传感器实际动态响应r(t)进行连续、滑动的动态误差修 正。传感器动态误差的频域修正流程框图如图9所示。先采用与H(f)点数相同的矩 形窗函数截取传感器实际动态响应r (t)，得到第i个N点数据段i~N(i)，再对ι·Ν (i)做傅里 叶变换得R(f)，然后将R(f)与H(f)进行点乘得C(f)，最后再对C(f)做逆傅里叶变换，得 到第i段数据的动态误差修正结果cN (i)。为避免数据i~N(i)加窗的泄漏效应导致的修正结果cN(i)数据两端较大的误差对 数据r(t)连续修正结果的影响，传感器动态响应的误差修正采取连续、滑动的频域修正方 式。连续、滑动的修正传感器动态响应误差的示意图如图10所示。图10(a)所示为误差修 正前的实际数据，每次待修正的数据段i~N(i)采用矩形窗进行连续、滑动截取，滑动时间为 τ，τ小于i^i)的时间长度。图10(b)所示为误差修正后的数据段的连续拼接示意图，取 每次修正结果cN(i)的中间段时间长度为τ的数据gcji)进行拼接，得到连续、滑动频 域修正后的反映实际被测量的真实值。具体实施例一的优点在于(1)数据截取时保证占空比为50%，从而使得动态误差频域修正函数中无效值点 规律性地出现在频域修正函数频率分辨率Δι的偶次倍频的频率点上，便于插补；(2)数据拼接处理降低了数据在频谱分析时由于周期延拓带来的较大误差；(3)将动态误差频域修正函数H(f)中位于要修正频带W，fmJ之外的频率点的修 正函数值置零，从而避免传感器动态响应在进行频域修正时高频噪声被放大；(4)采取连续、滑动的方式对传感器实际动态响应进行动态误差频域修正，降低了 数据在误差修正时由于加窗泄漏效应引起的连续修正结果的误差。实施例二图11所示为本发明方法的第二种具体实施例(简称为具体实施例二 )的技术方 案流程框图。具体实施例二中，动态阶跃响应实验1、数据截取2、数据拼接处理3、传感器动态 响应的误差修正6与具体实施例一相同。信号加窗4采用N(N = T · fs)点海明窗(Hamming)函数点乘所述数据截取2截 取的数据％和所述拼接处理3处理得到的数据yN'，得到加窗后的数据知和^/。％加海明 窗后得到的数据、在做傅里叶变换时由于频谱泄漏的特性，使得其变换结果中不含零值分 量，从而保证了频域修正函数咸(/)和氛(/)中不含无效值点。由于y/加同样的海明窗后 得到的数据Λ '与、具有同样的频谱泄漏特性，从而对频域修正函数A1C/)和应2(/)计算结 果的精度影响甚微。动态误差频域修正函数求解5的求解流程框图如图13所示。图13与图7不同 的是，由于攻(/)和应2(/)的计算结果中不含无效值点，所以，不需要图7中所示的步骤四对成(/)和应2(/)进行无效值点插补。图13所示具体求解步骤为步骤一对信号加窗4加窗处理后的数据％和扎’做傅里叶变换，得X(f)和Y(f)。步骤二 依据X (f)和Y(f)求取数据％和^/的自功率谱Gxx(f)、Gyy(f)和互功率 tb (;xy(f)、Gyx(f)。步骤三求有偏估计频域修正函数攻(/) ,H2(f)。步骤四对钱(/)和应2(/)进行相加和取平均，得到频域修正函数应(/)。步骤五将频域修正函数应(/)中位于要修正频带
之外的频率点的函数值 全部置零，得到最终所要求的传感器动态误差频域修正函数H(f)。在实际测量中，采用所述动态误差频域修正函数求解5求取的传感器动态误差频 域修正函数H(f)对传感器动态响应进行连续、滑动的动态误差修正。具体实施例二的优点在于除拥有具体实施例一优点外，利用了海明窗的频谱泄 漏特性，克服了动态误差频域修正函数计算中存在的无效值问题，所以，在传感器动态误差 频域修正函数求解5中无需进行插补，精度较高。实施例三图12所示为本发明方法的第三种具体实施例(简称为具体实施例三)的技术方 案流程框图。具体实施例三中，动态阶跃响应实验1、数据拼接处理3、信号加窗4、传感器动态 响应的误差修正6与第一种具体实施例相同，动态误差频域修正函数求解5与第二种具体 实施例相同。数据截取2则为采用时间长度为T(T = N/fs，N为截取的数据点数)的矩形窗同步 截取阶跃激励X(t)和传感器响应y(t)的动态过程数据段％和7,，％和7,的占空比为σ。 设Z表示整数域，Z+表示正整数域。％和7,傅里叶变换时的频率分辨率Af = 1/Τ，频域分 析的频率点f (η) = η· Δ ·(η彡0且nez)。由于傅里叶变换时视％为周期为Τ、频率为 Af的单周期矩形波，所以，当no e Ζ(η彡0且η e Ζ)时，信号％基频的η倍频谐波分量 的幅值IX (η · Af) I为零。因此，为避免％做傅里叶变换时零幅值分量的频率点位于要修 正的频带
(fmax ( fs/2)范围内，数据截取2采取下述步骤步骤一令零分量频率点η · Δ f > fmax，即得η > fmJ Δ f，η e Z ；步骤二 选取σ满足条件① N0 = N · σ e Z+，0 < σ < 1 ；②当且仅当η e (fmax/ 厶乜 + ~)且11￡2时，11.0 GZ。步骤三从阶跃信号的起跳零幅值点位置向前取N-N。，从其后面一点向后取N。 点，进而确定矩形窗截取X(t)数据的位置。矩形窗截取y(t)的位置与其截取x(t)的位置 在时间上同步。经过上述步骤后，数据截取2截取的实验数据的动态过程数据段％和yN的傅里叶 变换的零幅值分量频率点均位于要修正频带
之外，从而保证了动态误差频域修正 函数求解5所求的动态误差频域修正函数H(f)在要修正频带
范围内不存在无效 值，提高了精度。具体实施例三的优点在于除拥有具体实施例一的⑵、⑶、⑷优点外，在数据截取2截取数据时，就将最终动态误差频域修正函数的无效值排除在要修正频带范围之 外，免除了动态误差频域修正函数求解5中的插补运算，且该实施例动态误差频域修正函 数求解5求得的传感器动态误差频域修正函数H(f)的理论精度要比前两种具体实施例高。
权利要求
1.基于阶跃响应的传感器动态误差频域修正技术，用于根据传感器动态阶跃响应实验 数据求解传感器的动态误差频域修正函数，并将所求的频域修正函数用于传感器实际动态 响应的动态误差修正，技术流程包括动态阶跃响应实验、数据截取、数据拼接处理、信号加 窗、动态误差频域修正函数求解、传感器动态响应的误差修正，其特征在于由动态阶跃响应实验获取传感器阶跃激励和阶跃响应实验数据；然后，分别截取一段 实验数据；对截取的阶跃响应数据进行拼接处理；再对截取的阶跃激励数据和拼接处理后 的阶跃响应数据加上特定的窗函数；再根据加窗后的数据，求解有效修正频带宽、动态误差 修正精度高的传感器动态误差频域修正函数；在实际测量中，据此频域修正函数对传感器 实际动态响应进行连续、滑动的动态误差修正。
2.如权利要求1所述的基于阶跃响应的传感器动态误差频域修正技术，其特征在于 数据截取为采用时间长度为T(T = N/fs，N为截取的数据点数，fs为信号采样频率)的矩形 窗对传感器阶跃激励数据x(t)和响应数据y(t)从同一时刻开始进行同步截取，分别得到 数据％和yN，xN的周期延拓为周期为T、频率Δ f为1/T的方波，Δ f亦为对&和yN做傅里 叶变换的频率分辨率；针对具体实施例一和具体实施例二，截取的数据占空比σ均为0.5，以使 xN的频谱中频率为Δι的偶次倍频的频率分量为零；针对具体实施例三，设ζ表示整数域，Z+表示正整数域，数据截取采取下述步骤步骤一令零分量频率点η · Δ f > fmax，即得η > fmJ Δ f，n e Z ；步骤二 选取σ满足条件(DN0 = N · σ e Z+, 0 < σ < 1 ；②当且仅当η e (fmax/八乜+ ^)且1!￡2时，11.0 e Ζ；步骤三从阶跃信号的起跳零幅值点位置向前取Ν-Ν。，从其后面一点向后取N。点，进 而确定矩形窗截取x(t)数据的位置；矩形窗截取y(t)的位置与其截取x(t)的位置在时间 上同步。
3.如权利要求1所述的基于阶跃响应的传感器动态误差频域修正技术，其特征在于 数据拼接处理是将矩形窗所截取的数据yN中反映阶跃响应动态过程的数据段&进行倒相 并替代yN数据前端与&时间长度相同的S1段数据，进而得到拼接处理后的数据y/ ；yN' 中S1段数据计算如下S1 = -S2+A式中，A为阶跃响应yN的静态阶跃幅值。
4.如权利要求1所述的基于阶跃响应的传感器动态误差频域修正技术，其特征在于 信号加窗是将截取的激励数据&和拼接处理后的响应数据y/分别乘以相同的窗函数 K ·)，进而得到加窗的数据知和夕；/ ；针对具体实施例一和具体实施例三，窗函数W( ·)采 用N(N = T · fs)点矩形窗函数；针对具体实施例二，窗函数W( ·)采用N(N = T · fs)点海 明窗函数(Hamming)。
5.如权利要求1所述的基于阶跃响应的传感器动态误差频域修正技术，其特征在 于动态误差频域修正函数求解通过计算的自功率谱(ixx(f)、Gyy(f·)及其互功率谱 Gxy (f)、Gyx (f)，然后据此计算出有偏估计频域修正函数钱(/)和应2(/)；针对具体实施例一， 在要修正的频带
/2, η为偶数H2(n-Af) = H2(n-Af),"为奇数H1 (η. Af) = [H2 {{η -1). Af) + H2 ((η +1). Δ/)] / 2，η 为偶数对氏(《相加和取平均即得频域修正函数泠(/)，将点(/)中位于要修正的频带
之外的频率点的函数值置零，得到最终的传感器动态误差频域修正函数H(f)；针对实施例二和实施例三，直接对^(乃和众“/)相加和取平均，得到频域修正函数 应(/)，将应(/)中位于要修正频带[0，ffflaJ之外的频率点的函数值置零，得到最终的传感器 动态误差频域修正函数H(f)。
6.如权利要求1所述的基于阶跃响应的传感器动态误差频域修正技术，其特征在于 对传感器实际动态响应数据采取连续、滑动的动态误差频域修正方式，每次修正数据段 rN(i)为采用矩形窗对传感器动态响应r(t)进行滑动截取所得，矩形窗每次滑动时间τ小 于i~N(i)的时间长度，再对i~N(i)进行傅里叶变换得R(f)，然后，将R(f)与所求的动态误差 频域修正函数H(f)进行点乘得修正后的频域数据C(f)，接着，对C(f)进行逆傅里叶变换 即得修正后的数据cN(i)，取每次滑动截取的数据段的修正结果cN(i)的中间段时间长度为 τ的数据gcji)进行拼接，得到动态误差修正后的反映被测量真实情况的数值。
全文摘要
本发明为一种基于阶跃响应的传感器动态误差频域修正技术。针对传感器动态阶跃响应实验数据，截取实验数据中的一段并保证一定的占空比，再进行拼接处理；然后，将截取的阶跃激励数据和拼接处理后的阶跃响应数据加上矩形窗或海明窗；针对加窗后的数据，采用功率谱估计、在修正频带内进行无效值点插补和将修正频带外函数值置零的方法，得到传感器动态误差频域修正函数。采用此频域修正函数对传感器实际动态响应数据进行连续、滑动的动态误差修正，实现高精度动态测量。
文档编号G01D3/02GK102095430SQ201010548738
公开日2011年6月15日 申请日期2010年11月18日 优先权日2010年11月18日
发明者周全, 徐科军, 杨双龙 申请人:合肥工业大学