matlab求单自由度振动方程为,单自由度振动方程与Matlab/Simulink求解

1.问题

引用1：质量-弹簧-阻尼系统

引用2：模型推导

2.运动方程

Step1: 将微分方程最高阶变量移到等式左边

式1

Step2: 为每一阶微分式选择状态变量，最高阶除外

2.1

$x_{1} = x$

$x_{2} = x$ '

$x_{3} = x$ ''

...

式2: 通项

2.2: 同时求导

$\dot{x_{1} } =\dot{x}$

$\dot{x_{2} } =\ddot{x}$

2.3：整理为

$\frac{dy}{dt}$ 形式，各阶微分式替换为状态变量

$\frac{dx(1)}{dt} = \dot{x} =x_{2}$

$\frac{dx(2)}{dt} = \ddot{x} =\frac{1}{m} * (\mu -c*x(2)-k*x(1))$

3.1 Matlab 求解代码与结果

clear;clc;

tic

options = odeset('reltol',1e-13);

tspan = [0,80];

[tspan,x]=ode45(@vibration,tspan,[1 1],options);

toc

figure

plot(tspan,x(:,1),tspan,x(:,2));

legend('Position','Acceleration')

%Time Domain

fig = figure('Name','Position');

plot(x(:,1))

%Purpose Function

function dxdt = vibration(~ ,x)

%for t= 0:0.01:40

%f = cos(10*t);

f=0;

c = 0.3;

k = 5;

m = 1;

dxdt = [0;0];

dxdt(1) = x(2);

dxdt(2) =(1/m)*( f - c*x(2) - k*x(1));

end

图1: 加速度与位移

3.2 Simulink程序图与结果

图2：初值x(0) = 0, x(1) = 1, m = 1, k = 4, c = 0.1

Scope显示为：

matlab求单自由度振动方程为,单自由度振动方程与Matlab/Simulink求解相关推荐

matlab求二元函数极值算法_最优化计算与matlab实现（3）——进退法
参考资料 <精通MATLAB最优化计算(第二版)> 数值实现 Matlab 2019a 目录石中居士:最优化计算与Matlab实现--目录zhuanlan.zhihu.com 进退法 ...
亮度均匀性 matlab,求：亮度保持的夜景图像直方图均衡算法 matlab程序
该楼层疑似违规已被系统折叠隐藏此楼查看此楼求:亮度保持的夜景图像直方图均衡算法 matlab程序我是大四学生,最近在做一个论文,头疼死了,不知道这个论文的程序怎么写.这里是matlab论坛我想 ...
matlab求线性方程组的符号解,符号代数线性方程组的MATLAB求解.PPT
符号代数线性方程组的MATLAB求解建模常见问题MATLAB求解 Content 符号代数线性方程(组)的MATLAB求解最优化方法的MATLAB实现一维搜索问题线性规划无约束非线性最优化问 ...
matlab求连个书叫大叔,科学计算与MATLAB语言尔雅答案
[单选题]避雷器与被保护设备(). A. 串联于电源侧 B. 并联于电源侧 C. 串联或并联于电源侧 D. 之间无连接 [单选题]排气式避雷器是()的另一种叫法. A. 保护间隙 B. 管型避雷器 C ...
matlab求曲线每点的斜率_如何用matlab求出图中各条直线的斜率
展开全部 1.瞬时斜率62616964757a686964616fe4b893e5b19e31333431373236:求导. 方法:d(k)=(u(k)-u(k-1))/T(k);d(0)=0; d ...
matlab求洛伦兹方程的解,[转载]用Matlab求解洛伦兹方程
1. 洛伦兹方程求解本文说明用Matlab工具箱求解洛伦兹方程的过程,并给出吸引子的三维动态图象.洛伦兹方程如下: (1)这是一个自洽的方程组,求解过程如下: (1) 建立自定义函数 functio ...
matlab中对于xf(x)的积分,[matlab 积分]MATLAB求积分？
MATLAB求积分? 问题补充:我想用MATLAB求一下这个式子的积分,谁能帮我一下?g(f)=(e^4kxf)/sinh(kctf)未知量是f,最后求出来的式子里保留其他几个字母. ●matcom ...
matlab计算绝对值的导数,如何用matlab求带绝对值函数的导数？你值得一看的技巧...
用matlab求函数的导数:如果一个函数表达式中有绝对值,那对该函数求导就需要注意啦!该函数可能在他的零点处导数不存在,所以我们不能用传统的方法求导,应该在零点处分别求左导数和右导数,如果左导数和右导 ...
SLAM基础- 题目：基础矩阵F、本质矩阵E和单应矩阵H的自由度和秩
主要分析:基础矩阵F.本质矩阵E和单应矩阵H的自由度和秩,首先两个问题,尺度归一化和秩尺度归一化: 以本质矩阵为例,表达两帧的相机归一化坐标之间的对应关系将矩阵写成向量,转化为下式: 由于等式右侧 ...
matlab 求单/多元函数极值
matlab 求单/多元函数极值单元函数极值: 平时如果手算的话,就会先求导数,再求驻点,最终代值算出极值,如果用matlab代码求的话,就可以减少很多不必要的计算. fun=inline('0.5 ...

matlab求单自由度振动方程为,单自由度振动方程与Matlab/Simulink求解

matlab求单自由度振动方程为,单自由度振动方程与Matlab/Simulink求解相关推荐

最新文章

热门文章