位移传递率matlab编程,各种谱计算,频响函数,传递率
A.信号与谱的分类
由于时域信号有不同的分类, 变换后对应的频域也有不同的谱
信号可分为模拟(连续)信号和数字(离散)信号,
连续信号变换后称为谱密度, 离散信号变换
后称为谱.
连续信号又可分为绝对可积,平方可积(能量有限),均方可积(功率有限)
绝对可积信号有傅里叶谱(线性谱)和傅里叶谱密度(线性谱密度),如时域信号单位为电压V,
则前者单位为V,后者单位为V/Hz.
均方可积信号有功率谱PS(单位为V2)和功率谱密度PSD(单位为V2/
Hz.).
平方可积信号有能量谱密度ESD(单位为V2 s /
Hz.).
注1
平方量称为功率,平方量乘秒称为能量,谱分量除以频率称为谱密度.
注2
功率谱密度另一定义(离散信号的功率谱密度)见下述,
连续信号的功率谱密度.
为连续(光滑)曲线,
离散信号的功率谱密度为不连续的阶梯形..
注3
随机信号求功率谱密度时为减少方差,可采用平均,重叠和加窗处理(Welch法).
数字信号又可分为绝对可和,平方可和,均方可和.
B.各种谱计算
1.
线性谱Linear Spectrum:
对时域离散信号作DFT(离散傅里叶变换)得到,
采用方法为FFT(快速傅里叶变换)法.X(f)=FFT(x(t))
2.
自功率谱APS=Auto Power Spectrum:
离散信号的线性谱乘其共轭线性谱
APS(f)=X(f)*conj(X(f)), conj=conjugate共轭(实部不变,虚部变符号).
3.
互功率谱 CPS=Cross Power
Spectrum::x(t)的线性谱乘y(t)的共轭线性谱
互功率谱是复数,可表示为幅值和相位或实部和虚部等.
CPS(f)=X(f)
*conj(Y(f)) Y(f)=FFT(y(t))
4.
(自)功率谱密度PSD(=Power Spectrum
Density):
PSD(f)=APS(f)/Δf Δf—频率分辨率(Hz),
自功率谱密度与自相关函数成傅立叶对应关系
故功率谱密度也称为规一化的功率谱.
5.互功率谱密度CSD=CPS(f)/Δf
A.频响函数FRF, 传递率
A1.频响函数.FRF为响应的傅里叶变换与力的傅里叶变换之比或力和响应的互谱与力
的自谱之比 后者可通过平均减少噪声,故较常用.
H(f)=X(f ) / F(f)=X(f)*conj(F(f)) / F(f)*conj(F(f))=CPS /
APS.
A2. 频响函数有三种表达形式
频响函数表达成分子多项式与分母多项式(特征多项式)之比,也称有理分式.
(两多项式求根后)
频响函数表达成极点,零点和增益ZPK形式.
频响函数表达成部分分式,也称极点留数形式,(
部分分式的分子项称为留数.),
例如:最常见的单自由度(位移)频响函数H(ω)=X(ω)/F(ω)
H = 1 / (k+(jω)
2*m+jωC) 有理分式(多项式之比)
= (1 /m )*
1/(jω-p1)(jω-p2) 极点,零点和增益ZPK形式
= R1/(jω-p1)
+ R2/(jω-p2). 部分分式(极点和留数形式)
本例特殊,
分子非多项式,无根(无零点),留数为共轭虚数(一般为共轭复数)
a.共轭极点( 分母多项式的根 )
p: p1=σ+jωd, p2=σ-jωd, J=√-1
ωd--有阻尼固有频率, ωd=ωn
*√1-ζ2
b.共轭留数R: R1=1/2j*ωd R2=
-1/2j*ωd
c.增益K: K = 1/m
计算留数可用待定系数法或(复变函数中的)留数定理.
多自由度系统中留数包含振型信息.
A3.
频响矩阵: 当N点测力,N点测响应时, 频响函数为N x
N矩阵,但独立元素只有N个,
测试时既可只测一行(如H11,H12,H13,…H1N,
即激多点,测一点);也可只测一列
(如H11,H21,H31,…HN1,即激一点,测多点)
B.
传递率(Transmissibility)
传递率为同量纲物理量傅里叶变换之比,如电压传递率,力传递率,位移传递率等,
以位移传递率为例:
Tij=Xi(f)/Xj(f)= Xj(f)*conj(Xj(f)) / Xj(f)*conj(Xj(f))=CPSij /
APSjj
式中:
Xi(f)--
位移xi的傅里叶变换, Xj(f)--
位移xj的傅里叶变换,(不测力法无频响函数,只能用传递率求振型,此时xj位置保持不变,称为参考(基准)位移.
位移传递率matlab编程,各种谱计算,频响函数,传递率相关推荐
- 位移传递率matlab编程,机械振动设计计算与VB编程实例
1 编写机械振动计算程序的一般步骤1 1.1 明确编写程序的目的和了解清楚相关的计算公式.物理量纲1 1.2 在简单的界面上用BASIC编写计算程序并完成调试2 1.3 设计一个较为完备的用户界面4 ...
- matlab多自由度弹簧响应,用matlab编程实现法计算多自由度体系的动力响应..doc
用matlab编程实现法计算多自由度体系的动力响应. 用matlab编程实现法计算多自由度体系的动力响应 姓名:王 辉 学号 班级:硕1205班 专业:桥梁与隧道工程 2012年11月27日 用mat ...
- matlab编程刀尖频响,用半理论法预测主轴系统刀尖点频响函数
1问题的提出主轴系统刀尖点频响函数是衡量和评价数控机床切削稳定性的重要指标,是优化机床加工工艺.评价机床动态特性以及综合性能的依据.其中,机床的自激振动(又称再生颤振)对切削稳定的影响极大.目前,主要 ...
- matlab麦克斯韦速率,MATLAB编程首先建立计算麦克斯韦分布律的子程序mxwm-Read.PPT...
MATLAB编程首先建立计算麦克斯韦分布律的子程序mxwm-Read * * 第5章 综合实例 用 MATLAB 求解问题时,一般要经历建模和编程两个过程,只有在建模正确的前提下,方能得出正确的结果. ...
- 【MATLAB编程实战】计算图片的信息熵
欢迎关注,本专栏主要更新MATLAB仿真.界面.基础编程.画图.算法.矩阵处理等操作,拥有丰富的实例练习代码,欢迎订阅该专栏!(等该专栏建设成熟后将开始收费,快快上车吧~~) [MATLAB编程实战] ...
- ANSYS APDL谐响应分析——悬臂梁的频响函数计算以及幅值、角度(相位)、分贝计算
问题描述 研究一根悬臂梁,材质为钢材.长度 L = 2 L=2 L=2 米:截面为矩形,矩形的长度为 H = 5 c m H = 5cm H=5cm,宽度为 B = 2 c m B = 2cm B=2 ...
- 用步进正弦激励法测量计算机声卡输出输入幅值频响函数
一. 实验目的 (1) 通过对计算机声卡系统的实测加深对频响函数的理解: (2) 对声音振动及测量有直观了解: (3) 对计算机声卡系统的频响特性有直观的认识. 二. 实验内容 (1) 用声卡发出的步 ...
- 汽车振动响应分析-频响函数法(附程序)
文章目录 汽车的四自由度振动模型 频响函数法 MATLAB程序 汽车的四自由度振动模型 m:车身的等效质量: m_1:前轴的等效质量(包含电机的质量在内): m_2:后轴的等效质量: m_3:人体和座 ...
- 深入理解傅里叶变换的性质:实函数、卷积、相关、功率谱、频响函数
深入理解傅里叶变换的性质:实函数.卷积.相关.功率谱.频响函数 1实函数傅里叶变换的性质 1.1实函数傅里叶变换的性质 1.2实偶函数傅里叶变换的性质 1.3实奇函数傅里叶变换的性质 2傅里叶变换的基 ...
最新文章
- vue-route动态路由
- UVa11646 - Athletics Track(水题)
- Python列表的用法和基本操作
- Vim技能修炼教程(16) - 浮点数计算函数
- wxWidgets:wxSpinCtrlDouble类用法
- windows7 安装IIS没有default web site 解决方法
- Spring Cloud Consul 基础使用介绍
- 【高并发】java中的CAS,你需要知道的东西
- html真实雾效果图,HTML5 Canvas图片上云雾弥漫驱散特效
- 从标准输入读取几行输入。每行输入都要打印到标准输出上,前面加上行号。...
- window系统怎么解压tar.bz2文件
- `node2vec` `TSNE` 待解决问题
- HC3i论坛电子病历系统热门资料30篇(二)(
- 方志朋-2018年文章汇总
- 小学计算机网络教室管理制度,计算机网络教室使用管理制度[大全五篇]
- 智能计算机与应用期刊,智能计算机与应用期刊是省级期刊吗
- 【程序员如何买基金 序】基金学习地图
- c语言float m1 m2什么意思,M0、M1、M2的涵义及其作用
- PyQt6 使用 QAxWidget 打开 IE/Word
- contiki编程(一)
热门文章
- 天池训练营——基于人脸的常见表情识别(3)——模型搭建、训练与测试
- 任正非在荣耀送别会上的讲话(泪奔)
- 纯前端实现—抽奖小游戏
- 目前最新版AndroidStudio详细安装教程,纯萌新安装
- 知道python程序设计基础期末答案_Python程序设计基础_知到答案章节作业期末答案...
- 经典日内策略:ORB突破策略(期货)
- 想要搭建自己的腾讯云服务器怎么做
- net-java-php-python-校园二手图书销售网站计算机毕业设计程序
- Mac系统中QQ记录的存储路径
- UnityVR--小程序4--第一人称控制器