xinquan 发表于 2006-12-8 10:48

非线性微分方程的求解

方程如下：
x''+x'+kx+kz=0;
z'=x'-abs(x')*abs(z)^(n-1)*z-x'*abs(z)^n;
其中x‘’，z‘，x’都是指对时间的导数
这样的方程应该如何求解呢？
论文急用，请高手指点，谢谢！

tammy 发表于 2006-12-8 16:50

将方程写成标准形式

然后用ode45之类的求解，具体看帮助

yuanmm427 发表于 2006-12-9 16:15

function fff()
options=odeset('Reltol',1e-8,'absTol',);
=ode45(@rigid,,,options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
function dy=rigid(t,y)
k2=0.5;
k3=1.2;
dy=zeros(3,1);
dy(1)=-(y(1)+k2*y(2)+k3*y(3));
dy(2)=y(1)-abs(y(1))*abs(y(2)).^2*y(2)-y(1)*abs(y(2)).^3;
dy(3)=y(1);
%y(1)=x';y(2)=z;y(3)=x;
你照着试一下，我运行过了，可行，中间的系数你自己改把

xinquan 发表于 2006-12-11 13:39

谢谢楼上的，我去试一下
请问高手，如果振动方程不是齐次的，如果是在地震波激励下，那么地震波记录该如何输入呢？紧急求助啊！

[ 本帖最后由 ChaChing 于 2009-7-12 12:20 编辑 ]

handb 发表于 2007-1-12 09:08

这个程序是不是先把两个Function函数存为m文件,然后从k2=0.5开始运行程序啊,我这样做提示出错,这是怎么回事?谢谢
页: [1]
查看完整版本: 非线性微分方程的求解