求解微分方程组【x1，y1】
m*D2x1+c*Dx1+(k1-F*k2*(cos(2*pi*n*t1))^2)*x1-F*k2*sin(2*pi*n*t1)*cos(2*pi*n*t1)*y1=m*e*(2*pi*n)^2*cos(2*pi*n*t1)+m*g
m*D2y1 +c*Dy1 +(-F*k2*sin(2*pi*n*t1)*cos(2*pi*n*t1)*x1+(k1-F*k2*(sin(2*pi*n*t1))^2)*y1)=m*e*(2*pi*n)^2*sin(2*pi*n*t1)
其中F是与t1有关的分段函数，如下：
F=1.*(-pi/2+a2<=a1&pi/2-a2>a1)+0.5*(1+cos((a1-pi/2+a2)*pi/(2*a2))).*(pi/2-a2<=a1&pi/2+a2>a1)+0.*(pi/2+a2<=a1& 3*pi/2-a2>a1)+0.5*(1+cos((a1-3*pi/2-a2)*pi/(2*a2))).*(3*pi/2-a2<=a1&3*pi/2+a2>a1)
其中：
a1=2*pi*n*t1;
m=1.2; %单位：g
n=800;%转速
fs=10000;%采频
c=300;
g=9.5;
a2=1.3694;
k1=3.7688e+7;
k2=5.8406+7;
e=0.006;
方程组中，除t1，x1，y1外均已知或与t1有关，
初始条件：'x1(0)=1','Dx1(0)=0','y1(0)=1','Dy1(0)=0',
请大佬看看能不能用matlab写出来，我是想得到一组关于时间t1的（x1，y1）信号，
要是能求出关于t1的解就更好了
