[数学][笔记] 傅里叶变换公式的几种形式 (Several Forms of Fourier Transform Formulas)

实干、实践、积累、思考、创新。 # 傅里叶变换公式的几种形式 研究振动控制、随机振动等,离不开傅里叶变换(Fourier transforms)。对于傅里叶变换,不同文献、书籍有时候会采用不同形式的公式,刚开始看的时候有点凌乱,后面才理清楚,其实不同公式形式本质上都是等价。为了便于后续学习,以下总结几种常见的傅里叶变换公式形式。PS.这里傅里叶变换,指的是非周期函数(周期可以理解为无限长)的傅里叶变换,不是傅里叶级数(Fourier series),傅里叶级数是针对周期函数的说法。 ## 第一种形式 一种通用的形式,\(\Phi \left( p \right)\)和\(F\left( x \right)\)互为傅里叶变换对: $$\Phi \left( p \right) \Leftrightarrow F\left( x \right)$$ $$\Phi \left( w \right) = {1 \over {2\pi …

[数学][软件] FOUR_TRAN Example 4: Stock Periodicity and Volatility Analysis [FOUR_TRAN傅里叶分析工具 案例4: 股票周期性和波动性分析]

实干、实践、积累、思考、创新。 2020 年写的题目,现在已经是2022年了…… 这两天补上笔记                       案例 ( Examples )   [01] [数学][软件] FOUR_TRAN Example 1: Filtering [FOUR_TRAN傅里叶分析工具使用案例1: 滤波]   [02] …

[软件][数学][地震动] FPSA: Fourier and Power Spectra Analysis [地震波频谱分析工具]

实干、实践、积累、思考、创新。 程序图标( Program Icon ) 程序介绍 ( Program Introduction) 有小伙伴找我们,建议我们做一个这样的工具。断断续续花了许多天研究,做了出来。以下建安介绍这个小软件。这是一个采用快速傅里叶变换(FFT)进行傅里叶谱及功率谱密度分析的工具。功率谱密度(PSD)计算采用周期图法(Periodogram),并同时提供了多个不同的窗函数用于加窗处理。具体窗函数包括 Rectwin, Hann, Haming,  Bartlett。 软件可输出数据的傅里叶幅值谱 Fourier Amplitude、相位谱 Fourier Phase及功率谱 Power Specturm (PSD, Power Spectral Amplitude)。 软件支持批量导入信号,批量分析,并可将分析的傅里叶谱及功率谱导入程序,进行平均值计算。软件可以用于地震波功率谱的批量计算、分析及平均处理。 此外,根据网友建议,软件也增加了多种地震波的格式导入,具体包括PEER,YJK,PKPM,Midas等的地震波数据格式快速导入。 程序界面 ( Program Interface ) 20221123 …

[笔记] 用Matlab进行功率谱分析为什么会出现负的功率谱密度?

实干、实践、积累、思考,创新。 最近闲暇的时候翻翻 随机振动的书,又搞到傅里叶分析,功率谱等。 在测试网上的MATLAB功率谱分析程序的时候,发现功率谱居然是负数,(又暴露了知识面的短缺, 😥 ) 后仔思考发现,主要是因为功率谱使用了 分贝 为计量单位,比如通过公式 10*log10(Pxx) 将原本的功率谱Pxx进行转换,可见当Pxx处于0~1之间的时候,通过log10(Pxx)转换出来的公式就是负数了。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[数学][软件] FOUR_TRAN Example 3: Earthquake Ground Acceleration Frequency Spectrum Analysis [FOUR_TRAN傅里叶分析工具 案例3: 地震波频谱分析]

实干、实践、积累、思考、创新。 2020 年写的题目,现在已经是2022年了…… 这两天补上笔记                       关于软件( About the Program)   [01] [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具] 案例 ( Examples …

[振动台][试验][软件] 振动台(Shaking Table Test)白噪声试验数据处理——求结构自振特性

实干、实践、积累、思考、创新。 小伙伴问,振动台试验如何获得结构的自振特性。提取了振动台试验的数据,如何进行处理。 结构进行地震振动台试验前,均会进行结构动力特性试验。 自振特性的测试有很多种方法,如自由振动法、正弦波扫频法,白噪声扫频法。 其中白噪声扫频法的大概意思是,将模型安装在振动台后,进行地震波加载前,在振动台上输入小振幅的白噪声,进行激振试验,测量台面和结构的加速度反应。通过传递函数、功率谱等频谱分析方法,获得结构模型的自振频率、阻尼比、振型等参数。 于是小伙伴随手扔来一个白噪声扫频后测点的响应结果,按上面的思路,试试处理一下。 将测点响应导入本站的 FOUR_TRAN ( [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具] ) 软件,并进行傅里叶分析 (Fourier Analysis),如下图: 可以发现在频率3~3.5Hz位置,幅值谱很大。该位置很可能就是结构的基频。 将数据导出,并进一步导入本站的DataSmoothing ( [工具][试验][编程] DataSmoothing: A Program for Data Smoothing [试验数据曲线平滑+降噪工具] )软件,进行平滑处理。如下图所示: 可以较为清晰的看到结构的基频在3.2HZ左右,其他凸起是结构的其他阶频率。 …

[数学][软件] FOUR_TRAN Example 2: Square Wave Signal Decomposition [FOUR_TRAN傅里叶分析工具使用案例2: 方波信号分解]

实干、实践、积累、思考、创新。 20200205放假期间做的,现在整理出来。整个暑假,因为疫情都待在家中。 简单的傅里叶分析例子,方波信号分解。用前面自己写的 FOUR_TRAN 傅里叶分析工具进行分析。 准备如下的方波信号数据。 打开FOUR_TRAN,导入方波数据,进行傅里叶分析。 将前2阶频率的数据累加,获得下图粉色部分曲线。曲线类似一个sin函数,与方波差异较大。 将前5阶频率的数据累加,获得下图灰色部分曲线,曲线相对有之前多,与方波接近了。 将前20阶频率的数据累加,获得下图蓝色部分曲线,发现曲线以方波为中心进行震荡,并且与方波十分接近了。 接着将前50阶频率的数据累加。 以下是,将前300阶频率的数据累加的结果,可见曲线在离散点上与方波吻合了。方波可以由这么多三角函数累加而成。 绘制成EXCEL图形如下。 关于软件( About the Program)   [01] [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具] 案例 ( Examples ) [01] …

[数学][软件] FOUR_TRAN Example 1: Filtering [FOUR_TRAN傅里叶分析工具使用案例1: 滤波]

实干、实践、积累、思考、创新。 滤波,傅里叶分析的一个应用。用前面自己写的 FOUR_TRAN 傅里叶分析工具进行分析。 函数 y1: y1=10*SIN(2*PI()*0.5*x),振幅为10,频率为0.5; 函数2: y2=0.5*SIN(2*PI()*10*x),振幅为0.5,频率为10; 函数yw 为函数y1与函数y2的叠加: y=y1+y2 从图可见,由于y2振幅小,频率大,y1振幅大,频率小。叠加后y整体趋势与y1类似,相当于在y1的基础上加上了一些波动。 我么记下来通过傅里叶变化,想办法把从y中把这些波动的y2过滤掉。 打开FOUR_TRAN,导入离散后的函数y,并进行傅里叶分析。 从上图可见,傅里叶变换,可准确解析出两个主要频率成分,一个是0.5,一个时10,其中频率为0.5的傅里叶振幅谱大,而频率为10的傅里叶振幅谱小,能量成分低。与我们前面的函数一致。 接着我们把频率为0.5附近的成分累加起来(逆变换),如下图红色部分,红色曲线基本就是最初的y1,y2被成功过滤掉了。 如下,红儿曲线即是过滤掉频率10后的曲线,曲线与y1一致。 接着我们把频率为10附近的成分累加起来(逆变换),得到下图的红色部分,此时红色部分基本是函数y2,y1成功过滤掉了。 这就是简单的滤波的例子。 关于软件( About the Program)   [01] [数学][地震动][软件] FOUR_TRAN: Fourier Analysis …

[结构设计][地震作用][规范]振型分解反应谱法的一些概念总结 (Basic Concepts of Response Spectra Method)

实干、实践、积累、思考、创新。 温故而知新,理论指导实践,实践检验理论。 (1)振型型分解法,首先是进行模态分析,有多少个动力自由度,理论上就有多少个模态,相应的有多少个周期(频率),及振型。 (2)振型向量关于质量矩阵及刚度矩阵正交。因此,无阻尼运动方程可以实现解耦,将耦合的运动方程,解耦为多个广义单自由度运动方程。 (3)如果阻尼矩阵也满足于振型的正交性条件(如,瑞丽阻尼),则有阻尼结构的运动方程也可以解耦,解耦为多个有阻尼的广义单自由度运动方程。 (4)解耦后的单自由度方程的频率就是振型的频率。即,看是错综复杂的多自由度的震动过程其实是多个规则的不同频率的三角函数组成的。(PS. 自然界就是这么神奇,就像傅里叶变换一样,看是动态的,实则背后是静态的,是死的,太可怕了,无规律的东西,从频率来看,背后却是规律的… 这里不扯这个。 (5)振型无绝对大小,只是表示结构按某个具体频率振动时,各个动力自由度的振幅的相对大小。 (6)利用振型将多自由度方程解耦后分,若对解耦的单自由度方程进行时程分析,该方法常称为模态时成分析方法。若对解耦的单自由度方程进行反应谱分析,则称为振型分解反应谱法,这是目前结构设计规范的主流设计方法。 (7)振型分解法依靠振型对运动方程进行解耦,而振型是与弹性刚度及质量相关的,因此,机遇固定的振型对运动方程解耦,也意味着结构必须是弹性,该方法仅适用于弹性分析。 (8)振型分解反应谱法,由于引入了反应谱,使得结构工程师主要关注最大值,查看结果简便了,但是简便也带来了问题,因为反应谱丢掉了时程结果的许多信息。 (9)由于反应谱法只能获得最大值,因此振型分解反应谱涉及多个层次的组合问题。首先,各振型的极大值怎么叠加组合为最后的响应,该部分组合是所谓的“振型组合”,如常见的ABS组合方式,SRSS组合方式,及CQC组合方式等。另外,还有一个组合问题是多个方向的地震响应的组合问题,由于不同方向的地震动严格来说是不同的,所谓的不同,是说具体的时程肯定是不同的,响应的反应谱也是不同的,不同就会导致不同步,不同步那不同方向的结果也需要组合。直接时程分析法考虑多个方向的地震同时作用,直接就把多个方向的地震波加上同时进行分析即可,无非是动力方程的右边项将不同方向的地震波叠加即可,而振型分解反应谱法不行,不同方向的地震响应结果,也需要组合,先进行单个方向的效应分析,然后再把这些单个方向的极大值效应进行组合,该组合即所谓的“方向组合”。 (10)由于振型分解反应谱法的概念是,先计算单个振型的某个效应(如剪力,弯矩等)的最大值(正值),然后将单个振型的结果按一定的方法叠加起来,因此,振型分解再用反应谱分析再叠加的过程,丢掉了方向性。或者说,这些响应量,如剪力,只有一个统一的方向。结果都只有一个方向,那使用起来不方便,不直观,因此,在应用的时候,为了给出方向,又有研究者给出一些建议方法,判定响应方向,比如按主振型的方向,来确定响应的方向。但该方法也仅是对于一些简单结构,给出一个响应的参考方向,对于复杂结构,依然存在问题一些问题。比如,多塔连体结构,由于振型分解反应谱法,给出的不同塔楼的力都是同一个方向的,那振型分解反应谱可能就丢失了塔楼的反向运动,有可能存在隐患。因此,振型分解反应谱法虽然简便好用,但是也有不足,这个时候就需要补充弹性时程分析。这就是为何规范要求对复杂结构进行补充的弹性时程分析的一个重要原因。这个振型分解反应谱法的方向问题,还会引起其他相关的问题,具体工程的时候具体思考和分析。 (11)振型分解反应谱法的振型组合是非线性的,因此会出现诸如振型分解反应谱法的楼层剪力与楼层地震力(外力)不平衡的问题。因为,楼层剪力是多个振型的楼层剪力组合而得到的,单个振型下的楼层剪力是由于楼层地震力根据平衡求解的,满足平衡关系,但是经过振型组合后(如,SRSS,CQC),又不不平衡关系了,因为这些振型组合的方法都不是线性的。因为不能是线性的,为何?简单说,以为各个振型的极大值不是同时出现的。 (12)振型分解反应谱法,实际上是一个等效静力分析,为何这么说,因为运动方程经过解耦,再套上反应谱法,对于每一个振型,相当于在各个动力自由度上加上了一个等效惯性力,然后用这个惯性力进行静力分析,得到该振型下相关的响应量,如构件剪力,弯矩,轴力等,然后再进行振型组合。因此,在有限元求解上,其实是一个静力的求解分析过程。 (13)说到振型分解反应谱法,《高规》及《抗规》,又要扯到“扭转耦联”这个四个字,规范也给出了,两个方法,其中第一个是 a.不考虑扭转耦联的振型分解反应谱法,及b.考虑扭转耦联的振型分解反应谱法。其中,不考虑扭转耦联的振型分解反应谱法采用的是 SRSS组合,仅考虑一个水平方向的振型,即仅进行一个方向的振型分析,不考虑另一个方向质量或扭转惯量的耦合作用。考虑扭转耦联的振型分解反应谱法采用的是 CQC组合(CQC,组合过程中各个振型也是耦联的,需要通过两两振型的周期比及阻尼比参数来计算),分析过程中每一个楼层考虑水平方向及扭转方向3个自由度,振型也包含三个方向的分量。 (14)关于“扭转耦联”,必须说的是,由于一般结果,质量中心及刚度中心很难完全重合,因此,结构的扭转振动总是存在的,因此,进行“考虑扭转耦联的振型分解反应谱法”是相对更精确的。 (15)另外,进行“考虑扭转耦联的振型分解反应谱法”分析与是否考虑多向地震作用或者考虑哪个方向地震作用无关。不要将扭转耦联等同于双向或者三向,不考虑耦联则等同于单向地震。考虑扭转耦联,本质上说的是模态分析的时候,需要考虑平动与扭转自由度的耦联,模态需要能反应扭转的成分。是否考虑多向地震作用,只是振型分解后,方向组合的问题。是否考虑不同角度的地震只是涉及到振型参与系数的计算方法。振型分析,是否考虑扭转耦联或者不考虑扭转耦联,仅仅是结构固有特性的反映。 PS. 最后几点对SRSS和CQC及“扭转耦联”的表述还不是太清楚,借筑信达 李楚舒李总 的话补充一下:完全对称(没有扭转)的SRSS和CQC的结果也有较大区别,SRSS会在地震方向低估作用,而在另一个方向高估(见Wilson一书)。所以用CQC与结构是否扭转没关系,而是振型间存在耦合这一客观存在,所以必须用CQC。所以抗规的“扭转耦联”不对,应该是“振型耦联”——这误导了很多工程师. 相关博文( Related Posts ) [01] …