各位大侠帮帮忙画庞加莱映射和频谱图！！！急急急

基于MATLAB

dy=[x(3)
    x(4)
    0.1*sin(10*t)-2.5*x(3)+1.5*x(4)-x(1)+0.5*x(2)
    0.15*sin(10*t)-2.5*x(4)+1.5*x(3)-x(2)+0.5*x(1)];
各位我自己编了个程序 但是得到的结果不理想 希望各位能帮帮忙
画poincare映射：
M文件：
function poincare(tstop)
主函数：
global omega;
omega=10;
T=2*pi/omega;%激励的周期
step=T/10;%定义步长为T/10
y0=[0;0;0;0];
tspan=[0:step:100*T];%定义时间范围
[t,y]=ode45('zhendong',tspan,y0);
for i=500:10:1000%稳态过程每个周期取一个点
plot(y(i,1),y(i,3),'b.');
hold on;% 保留上一次的图形
end
xlabel('x1');ylabel('dx1/dt');
上面是其中的一个庞加莱界面程序。
>> tstop=50;
>> y0=[0;0;0;0];
>> [t,y]=ode45('zhendong',tstop,y0,[]);
>> yy=fft(y(end-1000:end,1));
>> N=length(yy);
>> power=abs(yy);
>> freq=(1:N-1)*1/step/N;
>> plot(freq(1:N/2),power(1:N/2));
这是一个频谱的程序。
如果我说的有不清楚的地方请留言，真的很需要大家的帮忙 我是MATLAB新手 刚接触