📄 new1.m
字号:
function OUT=new(t,X)
a=0.04;b=2;c=0.16;
x(1)=X(1);x(2)=X(2);x(3)=X(3);x(4)=X(4);
Q=[ X(5), X(9), X(13), X(17);
X(6), X(10), X(14), X(18);
X(7), X(11), X(15), X(19);
X(8), X(12), X(16), X(20);];
K=(abs(x(2)+x(3)+1)-abs(x(2)+x(3)-1))-(abs(x(2)+x(3)+3)-abs(x(2)+x(3)-3))-0.25*(abs(x(2)+x(3)+5)-abs(x(2)+x(3)-5));
dx(1)=-a*x(1)+b*x(2)+K*b;
dx(2)=-b*x(1)+c*x(2)+c*K;
dx(3)=-c*x(3)-x(4)-c*K;
dx(4)=x(3)-c*x(4)+K;
DX1=[dx(1);dx(2);dx(3);dx(4)];
K=(abs(x(2)+x(3)+1)-abs(x(2)+x(3)-1))-(abs(x(2)+x(3)+3)-abs(x(2)+x(3)-3))-0.25*(abs(x(2)+x(3)+5)-abs(x(2)+x(3)-5));
J=[ -a b+b*K b*K 0;
-b c+c*K c*K 0;
0 -c*K -c-c*K -1;
0 K 1 -c];
F=J*Q;
OUT=[DX1; F(:)];
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -