振动联盟浪子 发表于 2012-4-5 17:52

求助用数值法解非线性振动方程的源程序！！！？？？

本帖最后由 振动联盟浪子 于 2012-4-5 18:17 编辑

求助用数值法求解非线性微分方程的源程序！，主要目的是分析系统（连续系统）的非线性动力学特性（分叉或混沌等现象）。我的原始微分方程为这个样子的

最后离散化后的方程是这个样子的

最后转化为求解关于q的非线性方程组
希望各位能给予帮助或给予意见！！！先谢谢了.....

无水1324 发表于 2012-4-6 18:22

在论坛上有很多类似的，虽然不一定是4维的，你还是可以参考着写一个

振动联盟浪子 发表于 2012-4-6 18:50

回复 2 # 无水1324 的帖子

怎么找啊 一个一个看吗？我怎么就没看到呢 希望你们谁有看到的给我个衔接，谢谢啦

无水1324 发表于 2012-4-8 09:46

Duffing.m
function dy = Duffing( t,y )
dy=zeros(2,1);
dy(1)=y(2);
%f=0.35*cos(2*pi*0.1*t);
f=0;
dy(2)=f-(0.2*y(2)+y(1)+0.5*y(1)^3);
end

Duffing_simulation.m
clear;
clc;
=ode45('Duffing',,);
H1=figure;
%相图
plot(Y(:,1),Y(:,2));

这个是二维的，你类似的写成4维的就好了。
http://forum.vibunion.com/thread-110583-1-1.html

振动联盟浪子 发表于 2012-4-8 15:21

回复 4 # 无水1324 的帖子

谢谢....
页: [1]
查看完整版本: 求助用数值法解非线性振动方程的源程序！！！？？？