这个方程怎么编程序（关于悬臂梁受到一冲击）

tiange

在得到悬臂梁的运动方程，利用振型函数展开后，
如果给悬臂梁施加一个冲击力，怎么用MATLAB 编写程序求其广义坐标Y（t）?

dsolve('D2y+2*m*n*Dy+n*n*y=f(t)’）
但是这个f(t)是一个冲击函数 只是在t=0时刻的一个冲击力
在t=0时刻为一值，而在其他时刻均为零
怎么体现在函数里面的

笑石头

在simulink下达建模型很简单的就解出来了。D2y+2*m*n*Dy+n*n*y=f(t)经拉普拉斯变换（假设初始条件为零）得
s^2*y+2*m*n*s*y+n^2*y=f(t),上式可以变换为y=f(t)/(s^2+2*m*n*s+n^2)。在simulink下信号源为f(t)，即一个阶越信号（step time舍得很小，如0.01，initial value为t=0时刻的冲击力的值，final value为0），中间为传递函数1/(s^2+2*m*n*s+n^2),右边接一个示波器，示波器的输出就是所求的y了。

tiange

谢谢