matlab编程中的问题

来源:百度知道 编辑:UC知道 时间:2024/05/16 06:44:27
function hit
clear all;
clc;
y0 = [1;0;0;1;
1;0;1;-1;]

tn=[0];
for i=1:2

opts = odeset('events',@events1);
[t,y,te]=ode45(@f,[tn(1) Inf],y0,opts);
comet(y(:,3),y(:,1));
plot(y(:,3),y(:,1));
hold on;
comet(y(:,7),y(:,5));
plot(y(:,7),y(:,5),'m');
hold on;

if i<3
tn=[te(i)]
end
if i>=3
tn=[te(2)]
end

y01=[]
for n=1:4:8
if y(end,n)<0.0001
y01=[y01;y(end,n);-y(end,n+1)]
else
y01=[y01;y(end,n);y(end,n+1)]
end
if abs(y(end,n+2)-2)<0.0001|abs(y(end,n+2)+2)<0.0001
y01=[y01;y(end,n+2);-y(end,n+3)]
else
y01=[y01;y(end,n+2);y(end,n+3)]
end
[color=Red] for m=(n+4):4:8

问问老师吧

纯计算类问题,算法中有问题,希望经过你仔细检查可以找出错误,说实话,我看了2个小时才看懂,太丢人了,呵呵

虽然我学过,但看了半个多小时还是根本看不懂,碰上这个问题你可有的想了!
还是问问老师吧