主要内容

modalfrf

频率特性函数进行模态分析

描述

例子

降维= modalfrf (x,y,fs,窗口)估计频率响应函数的矩阵,降维,从激励信号,x,和响应信号,y,所有的采样率fs。输出,降维,是一个H1估计使用韦尔奇的计算方法窗口窗口的信号。xy必须有相同数量的行。如果xy是一个矩阵,每一列表示一个信号。

系统响应,y,被认为包含加速度测量。计算一个频率特性函数从位移或速度测量,使用“传感器”论点。modalfrf总输出频率特性函数动态灵活性(敏感性)格式无关的传感器类型。

降维= modalfrf (x,y,fs,窗口,noverlap)指定noverlap毗邻的重叠部分的样本。

例子

降维= modalfrf (___,名称,值)指定选项使用名称-值参数,从先前的语法使用输入的任意组合。选项包括估计量,测量配置和传感器测量系统响应的类型。

例子

(降维,f,)= modalfrf (___)还返回频率向量对应于每个频率特性函数,以及多个相干矩阵。

(降维,f)= modalfrf (sys)计算模型确定的频率特性的函数sys。使用估计命令等党卫军(系统辨识工具箱),n4sid(系统辨识工具箱),或特遣部队(系统辨识工具箱)创建sys从时域输入和输出信号。这个语法允许只使用的“传感器”名称-值参数。你必须有一个系统辨识工具箱™许可使用这种语法。

降维= modalfrf (sys,f)指定的频率来计算降维。这个语法允许只使用的“传感器”名称-值参数。你必须有一个系统辨识工具箱许可使用这种语法。

例子

modalfrf (___)没有输出参数块在当前图的频率响应函数。故事情节是有限的前四个激励,和四个反应。

例子

全部折叠

对于可视化的频率特性函数输入/锤激发。

加载数据文件,其中包含:

  • Xhammer - - - - - - 输入激励信号组成的五个锤打击定期交付。

  • Yhammer - - - - - - 响应系统的输入。Yhammer测量位移。

在4千赫采样的信号。情节激励和输出信号。

负载modaldata次要情节(2,1,1)情节(tham, Xhammer (:)) ylabel (“力(N)”次要情节(2,1,2)情节(tham, Yhammer (:)) ylabel (“位移(m)”)包含(“时间(s)”)

图包含2轴对象。坐标轴对象1包含一个类型的对象。坐标轴对象2包含一个类型的对象。

计算和显示频率特性函数。使用矩形窗窗的信号。指定的窗口覆盖锤打击之间的时期。

clf winlen =大小(Xhammer, 1);modalfrf (Xhammer (:), Yhammer (:), fs, winlen,“传感器”,“说”)

图包含2轴对象。坐标轴对象1润扬悬索桥与标题11包含一个类型的对象。坐标轴对象2包含一个类型的对象。

计算两个输入和两个输出系统的频率特性函数随机噪声而兴奋不已。

加载一个文件,其中包含的数据Xrand,输入激励信号Yrand,系统响应。计算频率特性函数使用一个5000 -样损害窗口和50%重叠相邻的数据段。指定输出测量位移。

负载modaldatawinlen = 5000;频= modalfrf (Xrand Yrand fs,损害(winlen), 0.5 * winlen,“传感器”,“说”);

使用的绘图功能modalfrf可视化的反应。

modalfrf (Xrand Yrand fs,损害(winlen), 0.5 * winlen,“传感器”,“说”)

图包含8轴对象。坐标轴对象1润扬悬索桥与标题11包含一个类型的对象。坐标轴对象2包含一个类型的对象。轴3润扬悬索桥与标题12包含一个对象类型的线。坐标轴对象4包含一个类型的对象。润扬悬索桥坐标轴对象与标题5 21包含一个类型的对象。6轴对象包含一个类型的对象。坐标轴对象7与标题频22包含一个类型的对象。8轴对象包含一个类型的对象。

估计对于一个简单的单输入/系统频率特性函数和比较它的定义。

一维离散振动系统由一个单位质量, 附在墙上,弹簧的弹性常数 k = 1 。一个传感器样品质量的位移 F 年代 = 1 赫兹。质量阻尼阻碍的运动通过施加一个力与速度成正比,与阻尼常数 b = 0 01

产生3000次样本。定义采样间隔 Δ t = 1 / F 年代

Fs = 1;dt = 1 / f;N = 3000;t = dt * (0: n - 1);b = 0.01;

系统状态空间模型可以描述的

x ( k + 1 ) = 一个 x ( k ) + B u ( k ) , y ( k ) = C x ( k ) + D u ( k ) ,

在哪里 x = ( r v ] T 状态向量, r v 分别的位移和速度,质量, u 是驱动力, y = r 是测量的输出。了状态矩阵

一个 = 经验值 ( 一个 c Δ t ) , B = 一个 c - - - - - - 1 ( 一个 - - - - - - ) B c , C = ( 1 0 ] , D = 0 ,

2 × 2 身份,和连续时间状态空间矩阵

一个 c = ( 0 1 - - - - - - 1 - - - - - - b ] , B c = ( 0 1 ]

Ac = [0 1; 1 - b];一个= expm (Ac * dt);公元前= [0,1];Ac \ B =(。(2)公元前*;C = 0 [1];D = 0;

质量是由随机输入第一2000秒然后离开回到休息。用状态方程模型计算的时间演化系统从一个零初始状态。情节的位移作为时间的函数。

rng默认的u = randn (1, N) / 2;u(2001:结束)= 0;y = 0;x = (0, 0);k = 1: N y (k) = x + D C * * u (k);x = x + B * * u (k);结束情节(t, y)

图包含一个坐标轴对象。坐标轴对象包含一个类型的对象。

估计系统的模态频率特性函数。使用损害窗口一半只要测量信号。指定输出质量的位移。

风=损害(N / 2);润扬悬索桥[f] = modalfrf (u ' y ', Fs,风,“传感器”,“说”);

一个离散时间系统的频率特性函数可以表示为时域的z变换系统的传递函数,评价单位圆。比较modalfrf估计的定义。

[b] = ss2tf (a, b, C, D);nfs = 2048;fz = 0:1 / nfs: 1/2-1 / nfs;z = exp (2 j *π* fz);ztf = polyval (b, z)。/ polyval (a, z);情节(f, 20 * log10 (abs(降维)))情节(fz * Fs, 20 * log10 (abs (ztf)))网格ylim (40 [-60])

图包含一个坐标轴对象。坐标轴对象包含2线类型的对象。

估计固有频率和振型的阻尼比。

博士(fn) = modalfit(润扬悬索桥,f, f, 1“FitMethod”,“页”)
fn = 0.1593
博士= 0.0043

比较的固有频率 1 / 2 π ,这是无阻尼系统的理论价值。

西奥= 1 /(2 *π)
西奥= 0.1592

估计频率特性函数和模态参数的一个简单的多输入/多输出系统。

一个理想的一维振动系统由两种质量, 1 2 ,在两堵墙之间。单位是这样 1 = 1 2 = μ 。每个质量是附加到最近的墙由一个弹簧弹性常数 k 。一个相同的弹簧连接两个群众。三种阻尼器阻碍运动的群众对他们施加力量与速度成正比,与阻尼常数 b 。传感器样品 r 1 r 2 ,群众的位移 F 年代 = 50 赫兹。

产生30000次样本,相当于600秒。定义采样间隔 Δ t = 1 / F 年代

Fs = 50;dt = 1 / f;N = 30000;t = dt * (0: n - 1);

系统状态空间模型可以描述的

x ( k + 1 ) = 一个 x ( k ) + B u ( k ) , y ( k ) = C x ( k ) + D u ( k ) ,

在哪里 x = ( r 1 v 1 r 2 v 2 ] T 状态向量, r v 分别的位置和速度 th质量, u = ( u 1 u 2 ] T 是向量的输入驱动力, y = ( r 1 r 2 ] T 是输出向量。了状态矩阵

一个 = 经验值 ( 一个 c Δ t ) , B = 一个 c - - - - - - 1 ( 一个 - - - - - - ) B c , C = ( 1 0 0 0 0 0 1 0 ] , D = ( 0 0 0 0 ] ,

4 × 4 身份,和连续时间状态空间矩阵

一个 c = ( 0 1 0 0 - - - - - - 2 k - - - - - - 2 b k b 0 0 0 1 k / μ b / μ - - - - - - 2 k / μ - - - - - - 2 b / μ ] , B c = ( 0 0 1 0 0 0 0 1 / μ ]

k = 400年 , b = 0 1 , μ = 1 / 10

k = 400;b = 0.1;m = 1/10;Ac = [0 1 0 0; 2 * 2 k * b k b; 0 0 0 1; k k / m / m b / m 2 * 2 * b / m];一个= expm (Ac * dt);公元前= [0 0;1 0;0 0,0 1 / m];Ac \ B =(。(4)公元前*;C = [1 0 0 0, 0 0 0 1);D = 0 (2);

群众是由随机输入整个测量。用状态方程模型计算的时间演化系统从一个零初始状态。

rng默认的u = randn (2 N);y = (0, 0);x = (0, 0, 0, 0);kk = 1: N y (:, kk) = C * x + D * u (:, kk);x = x + B * * u (:, kk);结束

使用输入和输出数据来估计系统的传递函数作为频率的函数。使用15000 -样本损害窗口9000个样本的重叠相邻段。指定输出位移测量。

风=损害(15000);小说= 9000;润扬悬索桥[f] = modalfrf (u ' y ', Fs,风,小说,“传感器”,“说”);

计算理论传递函数作为时域的z变换传递函数,评价单位圆。

nfs = 2048;fz = 0:1 / nfs: 1/2-1 / nfs;z = exp (2 j *π* fz);(b1, a1) = ss2tf (A, B, C, D, 1);(b2 a2) = ss2tf (A, B, C, D, 2);频(1:1)= polyval (b1 (1:), z)。/ polyval (a1, z);频(1:2)= polyval (b1 (2:), z)。/ polyval (a1, z);频(2:1)= polyval (b2 (1:), z)。/ polyval (a2, z);频(2:2)= polyval (b2 (2:), z)。/ polyval (a2, z);

估计和叠加理论预测的阴谋。

jk = 1:2kj = 1:2次要情节(2,2,2 * (jk-1) + kj)图(20 * log10 (abs(润扬悬索桥(kj jk:,))))情节(fz * Fs, 20 * log10 (abs(润扬悬索桥(jk: kj))))轴([0 f / 2 -100 0])标题(sprintf (“输入% d、输出% d '、jk kj))结束结束

图包含4轴对象。坐标轴对象1标题输入1,输出1包含2线类型的对象。坐标轴对象2标题输入1,输出2包含2线类型的对象。坐标轴对象3标题输入2,输出1包含2线类型的对象。坐标轴对象4标题输入2,输出2包含2线类型的对象。

情节估计使用的语法modalfrf没有输出参数。

图modalfrf (u ' y ', Fs,风,小说,“传感器”,“说”)

图包含8轴对象。坐标轴对象1润扬悬索桥与标题11包含一个类型的对象。坐标轴对象2包含一个类型的对象。轴3润扬悬索桥与标题12包含一个对象类型的线。坐标轴对象4包含一个类型的对象。润扬悬索桥坐标轴对象与标题5 21包含一个类型的对象。6轴对象包含一个类型的对象。坐标轴对象7与标题频22包含一个类型的对象。8轴对象包含一个类型的对象。

估计固有频率、阻尼比和模式系统的形状。使用peak-picking方法计算。

(fn,博士、女士)= modalfit(润扬悬索桥,f, f 2“FitMethod”,“页”);fn
fn = fn (:: 1) = 3.8466 3.8466 3.8495 3.8495 fn (:,: 2) = 3.8492 3.8490 3.8552 14.4684

比较无阻尼的固有频率的理论预测系统。

无阻尼=√eig [2 * k - k, k - k / m 2 * / m))) / 2π
无阻尼=2×13.8470 - 14.4259

计算两个输入的频率特性函数/ six-output数据集对应一个钢架。

加载包含输入激励和输出的结构加速度计测量。系统抽样的1024 Hz约3.9秒。

负载modaldataSteelFrameX = SteelFrame.Input;Y = SteelFrame.Output;fs = SteelFrame.Fs;

使用子空间方法来计算频率特性的功能。输入和输出信号划分为不重叠的,1000 -样本段。窗口每一段用一个矩形窗口。指定一个模型的36。

润扬悬索桥[f] = modalfrf (X, Y, fs, 1000,“估计”,“子”,“秩序”、36);

可视化系统的稳定图。确定15物理模式。

modalsd(润扬悬索桥,f, f,“MaxModes”15)

图包含一个坐标轴对象。标题为稳定图的坐标轴对象包含4线类型的对象。这些对象代表频率稳定,稳定频率和阻尼,在频率不稳定,平均响应函数。

输入参数

全部折叠

激励信号,指定为一个向量或矩阵。

数据类型:|

响应信号,指定为一个向量或矩阵。

数据类型:|

采样率,指定为一个积极的标量用赫兹。

数据类型:|

窗口中,指定为一个整数或一个行或列向量。使用窗口把信号分成部分:

  • 如果窗口是一个整数,然后呢modalfrfxy成段的长度窗口和窗户每一段用一个矩形窗口的长度。

  • 如果窗口是一个矢量,然后呢modalfrfxy成段长度相同的向量和windows每一部分使用窗口

  • 如果“估计”被指定为“子”,然后modalfrf忽略的形状窗口和使用它的长度来确定数量的频率点返回的频率特性函数。

如果的长度xy不能准确划分为一个整数的片段noverlap重叠的样本,那么相应的信号截断。

可用窗口的列表,请参阅窗户

例子:损害(N + 1)(1-cos(2 *π* (0:N) / N)) / 2都指定一个损害窗口长度N+ 1。

数据类型:|

重叠的样本数量,指定为一个正整数。

  • 如果窗口是一个标量,然后呢noverlap必须小于窗口

  • 如果窗口是一个矢量,然后呢noverlap的长度必须小于窗口

数据类型:|

确定系统指定为一个模型,确定参数。使用估计命令等党卫军(系统辨识工具箱),n4sid(系统辨识工具箱),或特遣部队(系统辨识工具箱)创建sys从时域输入和输出信号。看到确定模型的模态分析了一个例子。语法使用sys通常需要更少的数据比使用非参数方法的语法。你必须有一个系统辨识工具箱许可使用这个输入参数。

例子:ids ([0.5418 0.8373; -0.8373 - 0.5334], [0.4852, 0.8373], [1 0], 0 (0, 0), (0, 0), 1)生成一个对应于一个单位质量状态空间模型确定附加到一堵墙单元弹性常数的春天和阻尼常数为0.01。质量是采样的位移在1赫兹。

例子:idtf ([0 0.4582 - 0.4566), -1.0752 - 0.99 [1], 1)生成一个确定的传输函数模型对应于一个单位质量附加到一堵墙单元弹性常数的春天和阻尼常数为0.01。质量是采样的位移在1赫兹。

赫兹频率,指定为一个向量。

数据类型:|

名称-值参数

指定可选的双参数作为Name1 = Value1,…,以=家,在那里的名字参数名称和吗价值相应的价值。名称-值参数必须出现在其他参数,但对的顺序无关紧要。

R2021a之前,用逗号来分隔每一个名称和值,并附上的名字在报价。

例子:“传感器”、“或者”,“是”,“标题”指定的输入信号由速度测量和选择是H1的估计量。

估计量,指定为“标题”,“氢气”,“高压”,或“子”。看到传递函数更多的信息H1H2估计。

  • 使用“标题”当噪声不相关的激励信号。

  • 使用“氢气”当响应信号的噪声是不相关的。在这种情况下,激励信号的数量必须等于响应信号的数量。

  • 使用“高压”最小化建模和估计响应数据之间的差异,通过最小化错误的跟踪矩阵。Hv的几何平均数吗H1H2:Hv= (H1H2)1/2

    对于测量必须输入/输出)。

  • 使用“子”状态空间模型来计算频率特性的功能使用。在这种情况下,noverlap参数被忽略。这种方法比非参数方法通常需要更少的数据。看到n4sid(系统辨识工具箱)为更多的信息。

状态空间模型中存在的引线,指定为一个逻辑值。这个论点是可用的前提“估计”被指定为“子”

数据类型:逻辑

测量配置相同数量的激励和响应渠道,指定为“固定”,“rovinginput”,或“rovingoutput”

  • 使用“固定”当激励源和传感器系统的固定地点。激励有助于每一反应。

  • 使用“rovinginput”当从粗纱激发(或测量结果粗纱锤)测试。一个传感器系统的保持在一个固定的位置。一个激励源放置在多个位置和产生一个传感器响应/位置。函数的输出频(:,:我)= modalfrf (x (:, i), y(我):,)

  • 使用“rovingoutput”当从一个测量结果粗纱传感器测试。一个激励源系统的保持在一个固定的位置。一个传感器被放置在多个位置和响应每个位置一个激励。函数的输出频(:,i) = modalfrf (x (:, i), y(我):,)

订单状态空间模型,指定为一个整数或行向量的整数。如果你指定一个整数向量,那么该函数选择一个最优值指定范围。这个论点是可用的前提“估计”被指定为“子”

数据类型:|

传感器类型,指定为acc的,”或者“,或“说”

  • acc的——指定系统的响应信号与加速度成正比。

  • ”或者“——指定系统的响应信号与速度成比例。

  • “说”——指定系统的响应信号与位移成正比。

modalfrf总输出频率特性函数动态灵活性(敏感性)格式无关的传感器类型。

例如:无阻尼谐振子

一个简单的无阻尼谐振子的运动单位质量和弹性常数的采样率 f 年代 = 1 / Δ t 由传递函数描述

H ( z ) = N 年代 e n 年代 o r ( z ) 1 - - - - - - 2 z - - - - - - 1 因为 Δ t + z - - - - - - 2 ,

分子的大小取决于测量:

  • 位移: N d 年代 ( z ) = ( z - - - - - - 1 + z - - - - - - 2 ) ( 1 - - - - - - 因为 Δ t )

  • 速度: N v e l ( z ) = ( z - - - - - - 1 - - - - - - z - - - - - - 2 ) Δ t

  • 加速度: N 一个 c c ( z ) = ( 1 - - - - - - z - - - - - - 1 ) - - - - - - ( z - - - - - - 1 - - - - - - z - - - - - - 2 ) 因为 Δ t

计算三种可能的系统响应传感器的频率特性函数类型。使用一个示例2赫兹和30000个白噪声样本作为输入。

fs = 2;dt = 1 / f;N = 30000;u = randn (N, 1);ydis =过滤器((1-cos (dt)) * (0 1 1), (2 * cos (dt) 1), u);[frfd, fd] = modalfrf (u, ydis fs,损害(N / 2),传感器=“说”);yvel =过滤器(罪(dt) * (0 1 1), (2 * cos (dt) 1), u);(frfv阵线)= modalfrf (u, yvel fs,损害(N / 2),传感器=”或者“);yacc =过滤器([1 - (1 + cos (dt))因为(dt)], [1 2 * cos (dt) 1], u);[frfa, fa] = modalfrf (u, yacc fs,损害(N / 2),传感器=“acc”);重对数(fd、abs (frfd)阵线、abs (frfv),英足总,abs (frfa))网格传奇([“说””或者““acc”),位置=“最佳”)

图包含一个坐标轴对象。坐标轴对象包含3线类型的对象。这些对象代表说,或者,acc。

在所有情况下,生成的频率特性函数的格式对应位移。速度和加速度测量第一和第二次衍生品,分别的位移测量。范围的频率特性函数是等价系统的固有频率。远离固有频率,频率特性的功能有所不同。

输出参数

全部折叠

频率特性函数,返回为一个向量,矩阵,或三维数组。降维有大小p——- - - - - -——- - - - - -n,在那里p是频率垃圾箱的数量,反应的数量,n是激励信号的数量。

modalfrf总输出频率特性函数动态灵活性(敏感性)格式无关的传感器类型。

频率,作为一个向量返回。

多个相干矩阵,作为一个矩阵返回。有一个列每个响应信号。

引用

[1]“动态刚度、合规、移动性和更多…”年代我e米en年代,l一个年代t米od我f我ed 2019,https://community.sw.siemens.com/s/article/dynamic-stiffness-compliance-mobility-and-more

[2]布兰德,安德斯。噪音和振动分析:信号分析和实验过程。英国奇切斯特:约翰威利& Sons, 2011。

[3]欧文,汤姆。“介绍频率响应函数,Vibrationdata, 2000年,https://vibrationdata.com/tutorials2/frf.pdf

[4]Vold,哈佛,约翰·克罗利和g·托马斯Rocklin。“估计频率响应函数的新方法。”声音和振动。18卷,1984年11月,品种马非常页。

版本历史

介绍了R2017a

另请参阅

||(系统辨识工具箱)|

主题