Re: [其他] for 迴圏問題
※ 引述《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
推
04/19 08:43, , 1F
04/19 08:43, 1F
討論串 (同標題文章)
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
