用matlab 求解非线性方程周期解

lnzj0224 · 发表于 2012-10-16 14:49

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

本帖最后由 lnzj0224 于 2012-10-16 15:24 编辑

对于如下方程

d^2y/dt^2=-sin(t);
已知y(0)=y(2*pi)，y'(0)=y'(2*pi)

显然方程的解析解可以是y=sin(t);
但我用打靶法加newton迭代计算总是发散，请问是什么原因。
哪位高手可以指点一下具体点的算法。谢谢

zqyadam · 发表于 2012-10-17 19:10

syms y(t)
Dy = diff(y);
dsolve(diff(y,2) == -sin(t),y(0)== y(2*pi),Dy(0) == Dy(2*pi))

复制代码

用matlab直接求解，版本2012b

账号		自动登录	找回密码
密码			我要加入

[非线性振动] 用matlab 求解非线性方程周期解

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。