Maple求解常系数非齐次线性微分方程

对于单自由度简谐激励下的有阻尼强迫振动,其运动微分方程是一个二阶常系数非齐次微分方程。m\ddot{y}+c\dot{y}+ky=F_0cos\theta t

打开maple

1.文件---新建----工作表模式

2.第一行,输入微分方程,常数与变量之间用空格隔开。回车enter

3.选择第一行,在右边找到求解微分方程,点击y(t)

就会自动得到第二行和第三行,第三行就是方程的通解

假设初始条件为:y(0)=0; y'(0)=0;

4.在微积分里面,找到一阶微分,然后将f黏贴为y(t)的表达式,把x换成t。然后回车,就得到了y(t)的一阶导数

5.然后利用dsolve函数,来求解c1和c2

这里要注意括号的使用!

对于单自由度简谐激励下的有阻尼强迫振动系统来说,这个通解还需要进行化简,从而才方便分析系统有关的振动特性,以下给出振动力学教材的通解。

\begin{aligned}y=e^{-\xi\omega t}\left(\frac{\dot{y}_0+\xi\omega y_0}{\omega_r}\sin \omega_rt+y_0\cos \omega_rt\right)+Ae^{-\xi\omega t}\left(\frac{\xi\omega\sin \alpha-\theta\mathrm{cos}\alpha}{\omega_r}\mathrm{sin} \omega_rt+\mathrm{sin} \alpha\mathrm{cos}\omega_rt\right)\\&+A\sin{(\theta t-\alpha)}\end{aligned}

其中:

\begin{gathered} A=\beta y_{st} \\ \beta=\frac{1}{\sqrt{(1-\gamma^{2})^{2}+(2\xi\gamma)^{2}}} \\ \alpha=\arctan \frac{2\xi\gamma}{1-\gamma^{2}} \end{gathered}

这里放上公式识别的地址,以供读者copy公式

公式识别 (simpletex.net)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值