作者Schulman (Schulman)
看板MATLAB
标题Re: [其他] for 回圏问题
时间Wed Apr 18 23:39:47 2012
※ 引述《cancellarame (计时之神)》之铭言：
: 一直无法解决回圈只跑出1组值的问题 50
: 应该会有50*3 个值
: 请教大家
: 感谢大家解惑!!^^
function ptt=checkpoincarevarycccfinal
^^^
output似乎应该是ptt?
: s=0.005;
: e=0.25;
: time=[];
: ptt=[];
: for j=1:3
: c=0.18+0.02*(j-1);
: x=[1e-6,1e-6,1e-6,1-1e-6];
: [t, y]=ode45(@diflv4dhet,[0 8888],x,[],c,s,e);
: for i=2:length(t)
: if y(i-1,3) >0.1 && y(i,3) <0.1
: w =t(i)+((0.1-y(i,3))*(t(i+1)-t(i))/(y(i+1,3)-y(i,3))) ;
: %time temp
: time=[time w];
: end
: end
T=time(2:end-1)-time(1:end-2) ;
T=T(end-49 : end);
%end <=========把原来的end mark起来
point=[c*ones(1,50);T];
ptt=[ptt point];
end %<==============把end移到这里来
end
这样的话输出的ptt会是2x150的矩阵。应该就是你要的3次的总共结果。
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.109.73.193
1F：推 vincent59:这样真的是原po要的结果吗?? 如果是就太好了。 04/19 08:43