matlab
文章平均质量分 78
lijil168
这个作者很懒,什么都没留下…
展开
-
键合图-多体及多物理系统建模利器
键合图-多体及多物理系统建模利器1、 键合图的广义状态变量1.1势是流的原因1.2流是势的原因1.3动量是势对时间的累积效应1.4位移是流对时间的累积效应1.5功率是势和流的相干效应1.6动能(磁能)是流对冲量(磁通量)的累积效应1.7势能(电能)是势对位移(电荷)的累积效应2、键合图到方框图的转换3、键合图列写系统状态方程4、Matlab Simulink方框图建模5、AMESIM键合图建模6、Adams微分方程建模1、 键合图的广义状态变量1.1势是流的原因1.2流是势的原因1.3动量是势对时原创 2023-12-28 10:04:17 · 1121 阅读 · 2 评论 -
通过Matlab编程分析微分方程、SS模型、TF模型、ZPK模型的关系
对系统求微分,传递函数乘s,对系统求积分 ,传递函数乘1/s。原创 2023-08-21 15:15:59 · 1268 阅读 · 0 评论 -
挠性航天器姿态机动动力学模型及PD鲁棒控制
【代码】挠性航天器姿态机动动力学模型及PD鲁棒控制。原创 2023-05-08 21:03:10 · 1041 阅读 · 2 评论 -
PID整定二:基于Ziegler-Nichols的频域响应
*当用鼠标选取与虚轴相交的点时,找到系统振荡开始的增益Km=14,对应的极点有三个(红色的点),振荡频率wm=14rad/s,根据整定公式可得Kp=8.8371;基于稳定性分析的频域响应PID整定方法:对于给定的被控对象传递函数,可以得到其根轨迹,对应穿越jw轴的点,增益Km为系统开始振荡时的增益K值,wn为振荡频率。整定前系统根轨迹图见Figure1,系统有3个极点,根轨迹图有3根,从极点开始,终止与无限远(因为系统无0点,否则有终止于0点的线),**整定后的系统根轨迹图见Figure3,原创 2023-05-03 18:15:02 · 1802 阅读 · 0 评论 -
用matlab编写了一个DSP数据处理小软件
因为当系统出现故障时,故障通常会对系统信号产生调制作用,使系统的响应信号振幅产生周期性的变化(调幅),或者频率产生周期性的变化(调频),这种周期性变化的频率对应于故障发生的频率,变化的大小对应于故障的强弱,经过本软件希尔伯特变换得到瞬时频率和瞬时振幅,即可找到故障发生的频率和故障强弱,希尔伯特变换在故障诊断中具有广泛的运用。5、将调制后的信号进行FFT变换:频分析->频谱,可看出以30HZ的载波频率为中心,间隔5HZ及10HZ的边频(调制信号的频率),通过边频分析定位故障源也是NVH故障分析的一种方法。原创 2022-09-21 19:52:57 · 3354 阅读 · 1 评论 -
matlab simulink响应谱计算
用Abaqus进行仿真计算的加速度冲击响应的结果与SPECTR部分不一致,尝试找一下原因。结果与Abaqus仿真结果基本一致,加速度响应从14HZ到18Hz中途有下降,而。原创 2022-08-09 16:32:14 · 1257 阅读 · 2 评论 -
基于强跟踪卡尔曼滤波的隔振系统故障诊断——matab simulink仿真
一、参考资料1、基于强跟踪滤波器的在线故障诊断方法2、线性连续系统的离散化二、理论推导画重点:当考虑系统受到载荷激励时,微分方程右边为载荷激励;当考虑系统受到位移和速度的激励时,右边为位移速度激励以刚度k和阻尼c为系数的线性合作。这是两个不同的动力学系统。如果要卡尔曼滤波监测k和c,则需要扩展k和c两个状态变量,这时只能用第二种系统模型,否则卡尔曼滤波在状态预测时,不准确的k和c的估计值...原创 2020-01-08 21:20:34 · 2570 阅读 · 23 评论 -
强跟踪卡尔曼滤波STF估算车辆质量——matab simulink仿真
强跟踪卡尔曼滤波估算车辆质量——matab simulink仿真目录参考资料强跟踪滤波器算法解决的问题Simulink模型车辆运动状态仿真代码强跟踪卡尔曼滤波代码目录参考资料强跟踪滤波器理论参考:https://www.docin.com/p-1465349030.html基于强跟踪滤波器的故障诊断方法研究matab simulink仿真参考:卡尔曼滤波估算车辆质量——matab si...原创 2020-01-04 14:27:18 · 7029 阅读 · 81 评论 -
信号调制产生边频的原理及希尔伯特解调
文章目录1、定义1.1 调频1.2 调幅2、希尔伯特解调示例参考资料:https://wenku.baidu.com/view/20b6b3ebb9f67c1cfad6195f312b3169a451ea30.htmlhttps://blog.csdn.net/lijil168/article/details/1028768901、定义调幅、调频与调相1.1 调频1.2 调...原创 2019-11-05 10:21:39 · 6592 阅读 · 0 评论 -
从欧拉公式看希尔伯特变换
从欧拉公式看希尔伯特变换1、定义2、欧拉公式3、希尔伯特变换的意义1、定义希尔伯特变换可看做是一个线性滤波器:希尔伯特变换将x(t)的正频谱翻反转90度,负频谱正转90度,0谱线置0,得到x^(t),2次谱线转到180度位置,四次转回自身,可见,逆变换加个符号即可。jx^(t)的正频率谱线与x(t)相同,负频率谱线相反,x~(t)=x(t)+jx^(t)正好把x(t)的双边谱变成...原创 2019-11-02 23:39:18 · 2566 阅读 · 1 评论 -
从欧拉公式看三角波的单边谱与双边谱
1、在有限区间上,凡满足狄利赫里条件的周期函数都可以展开为傅里叶级数。其三角函数形式的频谱为单边谱,复指数函数形式的频谱为双边谱参考:https://wenku.baidu.com/view/af2ea95b7dd184254b35eefdc8d376eeaeaa17de.html实数系列的FFT后的X(k)是共轭对称的,即实部偶对称,虚部奇对称,matlab fft()得到的是双边谱,前部...原创 2019-10-22 23:17:24 · 7327 阅读 · 0 评论 -
时域补零对于DFT谱的影响
1、提高采样频率、时域补0无法提高物理分辨率提高采样频率不能提高频率分辨率的原因其实很简单,因为提高了采样频率,虽然在相同的观察时长那的点数增多了,但与此同时采样频率也变大了,点数增加几倍采样频率增加几倍,所以不改变观察时长而仅仅提高采样频率并不能提高DFT谱的频率分辨率。但是时域补零呢?采样频率没有变化,而点数增加无疑会减小DFT谱的相邻谱线间隔,相邻谱线间隔的缩小为什么不能提高频率的分辨率...原创 2019-10-17 11:20:35 · 5568 阅读 · 0 评论 -
(转)小波的分解和重构
转:天津大学 小波分析 宗婧 1015202078原理可参考:https://wenku.baidu.com/view/73439a6d5901020207409cd5.html1、单层小波分解%1. 单层小波分解 %读入信号 load leleccum; s=leleccum(1:4000); %通过db4小波基进行离散小波变换 [cA1,cD1]=dwt(s,'d...转载 2019-05-14 11:08:19 · 29574 阅读 · 12 评论 -
卡尔曼滤波simulink例子,位移和速度2变量估计
1、数学模型如下:2、simulink建模3、滤波结果如下图4、系统仿真程序与滤波程序%Plantfunction [u]=kalman_sim(u1,u2,u3,u4)%u1:输入加速度dv=u1;%u2:系统噪声WW=u2;%u3:测量噪声VV=u3;%u4:时间tt=u4;persistent v s Tsif t==0 Ts...原创 2019-07-10 17:15:24 · 6425 阅读 · 3 评论 -
单自由度振动全解:matlab理论计算&& virtual.lab motion仿真&&Excite PU仿真
1、单自由度振动理论推导见:https://blog.csdn.net/lijil168/article/details/88932542其全解为:2、利用matlab编程如下(计算发动机飞轮位移激励下离合器从动盘的位移响应):k1=7.5/(pi/180);m1=0.22;wn=sqrt(k1/m1);s_=0.5;wd=sqrt(1-s_^2)*wn;c1=...原创 2019-06-30 09:03:17 · 7508 阅读 · 1 评论 -
matlab编程进行卡尔曼滤波的简单例子
1、matlab编程进行卡尔曼滤波的简单例子clearN=50;x(1)=0; %理论速度初值v(1)=0ut=-270; %加速度值 F=1; %状态转移矩阵B=0.01; %控制矩阵 步长值H=1; %观测矩阵v=randn(1,N)*5; %观测噪声 均值=0,方差=5R=[5];%观测噪声协方差矩阵w=randn(1,N)*1;%预测噪声,均值为0,方差=1Q...原创 2019-07-04 16:12:29 · 2578 阅读 · 2 评论 -
利用matlab小波GUI工具箱对空挡噪声进行滤波
利用matlab小波GUI工具箱对空挡噪声进行滤波:1、命令wavemenu调出matlab小波GUI工具箱,2、File->load加载噪声音频文件,wavelet 1-D,选择小波函数,设置阀值对细节部分(高频)进行过滤,3、 File->Save->De-noised Signal保存滤掉白噪声后的信号到.mat文件,matlab加载.mat文件到x...原创 2019-06-02 11:35:55 · 2142 阅读 · 2 评论 -
matlab计算曲线形心,并将y轴移动到形心上,使图形居中布置
m=[-4.1 -4 -3.7 -3.4 -2.7 -2 -1.6 0 1.6 2 2.9 3.5 3.9 4 4.2];n=[10 9.1 6 4 2 0.8 0.5 0 0.5 0.8 2 4 6 8 10];plot(m...原创 2019-05-24 09:49:52 · 2560 阅读 · 0 评论 -
matlab小波分析常用函数
参考:https://wenku.baidu.com/view/1973af1a580102020740be1e650e52ea5518ceb8.html1、连续小波变换scale=1:200;cw2=cwt(y,scale,'morl');[X,Y]=meshgrid(t,centfrq('morl')*Fs./scale);%小波的中心频率*采样频率/尺度=》频率%scal2frq...原创 2019-05-14 16:25:16 · 6811 阅读 · 1 评论 -
单缸发动机扭矩动力学计算:理论计算&&virtual.lab motion仿真
一、单缸发动机曲轴连杆活塞系统动力学分析参考:汽车动力总成系统/(伊朗)贝赫鲁兹.马沙迪,(英)戴维.克罗拉著;白先旭,刘勇强译。——北京:机械工业出版社,2018.6(汽车先进技术译丛.汽车创新与开发系列)书名原文:Vehicle Powertrain SystemsISBN 978-7-111-60005-3燃烧室压力如下:X_Colum...原创 2019-07-05 14:05:07 · 4174 阅读 · 0 评论 -
matlab设计声音听听
用不同幅值、频率和相位的信号构成声音信号,看声音的音质是否与波形有关。1、单频率声音幅值不同%%F=[494 440 392 440 494 494];fs=8000;tend=6;t=linspace(0,tend,fs*tend+1);t(end)=[];y=exp(t-6).*sin(2*pi*F(floor(t)+1).*t);plot(t,y);sound(y,...原创 2019-07-15 10:19:20 · 358 阅读 · 0 评论 -
二自由度振动仿真:matlab直接解微分方程&&virtual.lab motion仿真
一、二自由度振动解析法与多种数值算法见:https://blog.csdn.net/lijil168/article/details/67647924二、用matlab符号建立微分方程,并用龙格库塔法 ode45()直接进行数据计算%SH760小轿车空载主要参数clear;m=1340;a=1.54;b=1.29;l=a+b;Ic=2395;%绕质...原创 2019-06-30 23:02:58 · 7559 阅读 · 0 评论 -
matlab对文本文件的读写
https://blog.csdn.net/u012366767/article/details/81565427介绍Matlab读写文本文件的方法,包括load、dlmread、csvread、importdata、textscan、dlmwrite、csvwrite、save等函数的用法。读文件1. 纯数字load test.txt % 导入文本数据,并保存在 test 变量;...转载 2019-08-13 15:21:46 · 6628 阅读 · 0 评论 -
振动系统零频率和重频率问题探讨
振动系统零频率和重频率问题特征方程有零根和重根问题,对应的振动系统,其模态是什么样的呢?1、零频率(0特征根)设质量为m的两个小球用刚度为k的弹簧相连,求该2自由度系统的模态。syms m kM=[m,0;0,m]K=[k,-k;-k,k][eig_vec,eig_val]=eig(inv(M)*K)eig_vec = [ 1, -1][ 1, 1]eig_val =...原创 2019-09-29 13:36:45 · 1287 阅读 · 1 评论 -
matlab虚拟现实仿真
clear;%SH760小轿车空载主要参数m=1340;a=1.54;b=1.29;Ic=2395; %绕质心的转动惯量k1=40000;k2=44000;M=[m,0;0,Ic];K=[k1+k2,-(k1*a-k2*b);-(k1*a-k2*b),k1*a^2+k2*b^2]; [eig_vec,eig_val] = eig(inv(M)*K);原创 2017-03-28 20:05:11 · 1968 阅读 · 0 评论 -
用simulink分析傅里叶变换、短时傅里叶变换和小波变换
参考:https://zhuanlan.zhihu.com/p/22450818?refer=dong51、FFT,分析基波的参与时间对傅里叶变换的影响基波为4个余弦波,用matlab程序控制其参与输入信号的增益。当增益为1,4基波全时参与时,输入信号为:注意:FFT后在Vector Scope显示每帧的频率图时,应设置采样频率,否则,频率计算可能有误!比如,Buffer设...原创 2019-04-15 11:14:19 · 19284 阅读 · 3 评论 -
matlab绘3d图
tn=0.1:0.1:0.5;%x轴:横向时间fn=1:8;%y轴:纵向频率A=ones(8,5);%列向-横向-时间x; %行向-纵向-频率yA(4,2)=0; %Z(y,x)[T,F]=meshgrid(tn,fn);mesh(T,F,A); %三维绘图axis tight;title('三维时频图'); ylabel('Frequency(Hz)'...原创 2019-04-16 22:56:23 · 810 阅读 · 0 评论 -
单自由度阻尼强迫振动通解求导及simulink验证(修正了网上常见的错误)
参考资料见:https://wenku.baidu.com/view/009915e19b89680203d82503.html注意:文中对响应滞后的相位角计算有误! 注意: 公式9推导到10时有误。由公式9,响应滞后的相位角phai应在0到pi之间,所以,当频率比>1时,由10计算的相位角应该+pi ...原创 2019-03-31 18:33:10 · 6748 阅读 · 9 评论 -
SH760模态分析-多种解析与数字计算方法
详细分析参考:https://wenku.baidu.com/view/6915c300f08583d049649b6648d7c1c708a10b62 主振型与固有频率;偏频(假设分配系数=1) 模态分析:可以清楚地看到模态和主模态之间的关系!!!!!! 当取初值x10=1时,仿真的结果为:x1=A11sin(ω1t+φ1)+A12sin(ω2原创 2017-03-28 21:34:03 · 2135 阅读 · 0 评论 -
振型矩阵与正则振型矩阵
%由振型矩阵mode_vec(特征列向量组成)、质量矩阵M求正则振型矩阵mode_vec_nfunction mode_vec_n= mode_n(mode_vec,M)M_p=mode_vec'*M*mode_vec; %求主质量阵n=length(mode_vec);mode_vec_n=zeros(n);for i=1:n mode_vec_n(:,i)=m原创 2017-03-28 20:22:36 · 6928 阅读 · 0 评论 -
非正交坐标系-模态正交的物理及几何意义
模态正交不是指模态向量正交,而是指模态向量关于质量矩阵与刚度矩阵正交,即 Ф(i)TMФ(j) =0Ф(i)TKФ(j) =0物理上表示一模态方向上的力无法影响另一模态方向上的位移及加速度等,故各模态能量独立,不存在能量流动。关于模态的几何意义见下图。 系统在主坐标系中绕主坐标轴做主运动,为轨迹在坐标轴上的投影 (非垂直投影)。原创 2017-03-28 20:20:07 · 4629 阅读 · 0 评论 -
SH760二自由度系统自由振动解析法求解
解析方法:先求系统主模态,再将模态向量正则化,将系统转换到正则坐标系进行解耦,求解完成再转换回来。 %解析法求通解:零输入响应。%SH760小轿车空载主要参数clear;m=1340;a=1.54;b=1.29;l=a+b;Ic=2395; %绕质心的转动惯量rou=sqrt(Ic/m);k1=40*1000;k2=44*1000;M=[m*(b原创 2017-03-28 20:16:15 · 1118 阅读 · 0 评论 -
SH760二自由度系统解析法求全解-过渡及稳态过程
%解析法求特解及全解:简谐力激励的稳态受迫响应与过渡过程。%SH760小轿车空载主要参数clear;m=1340;a=1.54;b=1.29;l=a+b;Ic=2395; %绕质心的转动惯量rou=sqrt(Ic/m);k1=40*1000;k2=44*1000;M=[m*(b^2+rou^2)/l^2,m*(a*b-rou^2)/l^2;m*(a*b-原创 2017-03-28 20:13:59 · 651 阅读 · 0 评论 -
SH760二自由度系统非周期激励解析法求解
%解析法求非周期激励的响应。采样、卷积积分法。%SH760小轿车空载主要参数clear;m=1340;a=1.54;b=1.29;l=a+b;Ic=2395; %绕质心的转动惯量rou=sqrt(Ic/m);k1=40*1000;k2=44*1000;M=[m*(b^2+rou^2)/l^2,m*(a*b-rou^2)/l^2;m*(a*b-rou^2)原创 2017-03-28 20:10:04 · 412 阅读 · 0 评论 -
matlab几何计算程序集
一、求矢量长度(模) function x=leng(a)x=sqrt(sum(a.^2)); 二、求矢量的夹角 function ang=ang(a,b)ang=acos((a'*b)/leng(a)*leng(b)); 三、求点到线上的垂足 %求点到线上的垂足;P0线外点,P1线上点,n线的方向;function P=point_po2l(P原创 2017-03-28 20:00:54 · 1106 阅读 · 0 评论 -
ss状态空间方程matlab计算
k1=1000;k2=2000;k3=3000;m1=1;m2=2;alpha=1;beta=0.005;M=[m1,0;0,m2];K=[k1+k2,-k2;-k2,k2+k3];n=size(M,1);C=alpha*M+beta*K;ssA=[zeros(n),eye(n);-M\K,-M\C];ssB=[C,M;M,zeros(n)]\[eye(n)原创 2017-03-28 19:57:27 · 11265 阅读 · 0 评论 -
ss模型复模态的物理意义及adams复模态振型求解
adams可以线性化输出SS状态空间ABCD矩阵,根据A可以进行系统复模态计算。A可看作是角频率张量矩阵,由模态矩阵进行坐标变换解耦后可以得到所有的模态参数,特征向量(模态系数)为复数不仅影响幅值,还影响相位角,而特征值为复数,其虚部为阻尼频率,而实部会使振动随时间衰减。 GX˙+HX=0ФTGФ Ф-1X˙+ФTHФ Ф-1X=0将X˙转换到主模态基下,G转换原创 2017-03-28 19:48:37 · 1933 阅读 · 0 评论 -
机器学习:利用核函数进行非线性分类
%Preceptron4.m%感知机算法(对偶形式)利用核函数完成分类算法%f(X)=W'*X+b%f([X;1])=[W;b]'*[X;1] 增加维数,简化计算%设W初始为0时,权重W是训练点的线性组合:W=Σaj*yj*Xj%类别yi决定了Xi的系数,ai>=0,正比于被误分后权重更新的次数,可看作Xi信息量的指示%误分判据 yi*((W'*Xi)+b)%%%%%%原创 2017-03-27 22:55:48 · 2417 阅读 · 0 评论 -
机器学习:空间变换进行非线性分类
%Preceptron3.m%感知机算法(对偶形式)%f(X)=W'*X+b%f([X;1])=[W;b]'*[X;1] 增加维数,简化计算%设W初始为0时,权重W是训练点的线性组合:W=Σaj*yj*Xj%类别yi决定了Xi的系数,ai>=0,正比于被误分后权重更新的次数,可看作Xi信息量的指示%误分判据 yi*((W'*Xi)+b)%%%%%%%%%%%%%%%%%原创 2017-03-27 22:09:54 · 1881 阅读 · 0 评论 -
Matlab数据插值-内插、外插
插值就是已知一组离散的数据点集,在集合内部某两个点之间预测函数值的方法。一、一维插值插值运算是根据数据的分布规律,找到一个函数表达式可以连接已知的各点,并用此函数表达式预测两点之间任意位置上的函数值。插值运算在信号处理和图像处理领域应用十分广泛。1.一维插值函数的使用若已知的数据集是平面上的一组离散点集(x,y),则其相应的插值就是一维插值。MATLAB中一维插转载 2013-10-04 18:50:55 · 66027 阅读 · 0 评论