如果非线性微分方程中含有自变量的分段函数,如何用ode函数求解

eatnche594

這次我已經降階好的運動方程式

function  yp=ivp13(t,y)

yp=zeros(1,16);  % note: column vector
yp(1)=y(2);
yp(2)=((6.15*10.^6*(y(3)-y(1)))-(6.15*10.^6*y(1)))/(13.2)-y(16);
yp(3)=y(4);
yp(4)=((8.1*10.^7*(y(5)-y(3)))-(6.15*10.^6*(y(3)-y(1))))/(1.3)-y(16);
yp(5)=y(6);
yp(6)=(6.03*10.^6*(y(7)-y(5))-8.1*10.^7*(y(5)-y(3)))/(2.28)-y(16);
yp(7)=y(8);
yp(8)=(8.1*10.^7*(y(9)-y(7))-6.03*10.^6*(y(7)-y(5)))/(1.12)-y(16);
yp(9)=y(10);
yp(10)=((V+0.5*A*y(9)+0.5*J*y(9).^2)*8.1*10.^7.*(D+y(11)-(D+V*y(9)+0.5*A*y(9).^2)))-2.01*10.^6*(y(9)-y(7))/(3.66)-y(16);
yp(11)=y(12);
yp(12)=(2.01*10.^6*(y(13)-y(11))-2.83*10.^8.*(D+y(11)-(D+V*y(9)+0.5*A*y(9).^2)))/(0.35)-y(16)*V-A*y(15).^2;
yp(13)=y(14);
yp(14)=(1000-2.01*10.^6*(y(13)-y(11)))/(5)-y(16)*V-A*y(15).^2;
yp(15)=y(16);
y(16)=(107.257-D*(0.25+0.104+5)*(V*y(16).^2+1000))/(0.096+13.2+1.3+2.28+1.12+3.66+((0.25+0.104+5)*D.^2));
yp=yp';

其中D,V,A,J是有自变量的分段函数,如附件 Appendix_A.doc (56.5 KB, 下载次数: 23)

2007-8-21 16:16 上传

点击文件名下载附件

想請問大俠,俠女我該如何用ode函数求解

給定條件為

tspan=[0 0.4];
x0=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
求解yp(16),yp(14)

麻煩了 感激不盡

无水1324

里面的X是自变量，相当于t吗？

eatnche594

x0是初始值全部都為0

无水1324

哦，不好意思我#2写错了，
是你附件中的X

eatnche594

θa 是凸輪作動角,是指凸輪旋轉角 θi 在此轉位角內,將驅動輸出軸做分割旋轉,在轉位角外,輸出軸並不跟輸入軸連動而處於暫停狀態  0~50度dwell  60~310度是active 320~360度dwell

无水1324

%           x-b, x>b
%      f=   0  , abs(x)<b
%          x+b,x<-b

x''+c*x‘+f=0

funcion dy=fun(t,x)
% b为常数
if x>b
f=x(1)-b;
elseif x<-b
f=x(1)+b;
else
f=0
end
dy=[x(2);-c*x(2)-f]

这是一个简单的分段的例子，我对你的物理模型及参数还是没有看懂，所以只能给你写一个典型的单自由度的例子

eatnche594

掰謝 讓您看不懂了 會照這例題 試著解看看

咕噜噜

建议楼主用简体字，繁体字看不懂