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

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

本帖最后由 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
页: [1]
查看完整版本: 用matlab 求解非线性方程周期解