作者cancellarame (计时之神)
看板MATLAB
标题Re: [其他] for 回圏问题
时间Fri Apr 13 23:30:31 2012
抱歉 其实我不太懂程式 只是觉得这巢状回圈 J=1:3 应该会跑3次
可市跑出来的ˇ只有J=3
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
point=[c*ones(1,50);T];
ptt=[ptt point];
end
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 120.107.198.122
1F：→ Append:你为什麽认为跑出来的只有J=3? 04/13 23:32
2F：推 vincent59:你看到workspace才这麽说的吧??? 04/13 23:41
3F：→ vincent59:就目前看来 的确已跑了3次 04/13 23:42
4F：→ vincent59:你想验证的话 很简单 你在回圈里 打j 04/13 23:43
5F：→ cancellarame:抱歉 词不达意 应该是说跑了3回 但是取值只看到1回 04/13 23:55
6F：→ cancellarame:应该怎麽写 才能取到3回的值 谢谢 04/13 23:58
7F：推 rex0707:c要用向量写吧 不然每换一个j c就被覆盖一次 04/14 01:10
8F：→ rex0707:阿 抱歉没看清楚 你是哪个要取三回的值? 04/14 01:11
9F：→ rex0707:要取3回的值你有将他设成是矩阵吗? 04/14 01:12
10F：→ rex0707:再依序将值写入 04/14 01:12
11F：→ cancellarame:应该是ptt 要取3回的值 04/15 01:15
12F：→ cancellarame:程式目的是要判别是否为周期解 现在是将C值改变 04/15 01:25
13F：推 vincent59:你能否告诉我，你最後跑出来ptt变数的大小? 04/15 07:40
14F：推 vincent59:ptt的内容??? 04/15 07:44
15F：→ vincent59:还有你贴上来的程式其实不能跑的??? 04/15 07:45
16F：推 vincent59:是哪3回的值???? 04/15 08:21
17F：→ vincent59:不要让我们猜阿?? 04/15 08:22
18F：→ vincent59:如果猜的没错，其实你都抓到那3个值 04/15 08:23
19F：→ vincent59:其中一行改成 point=[c*ones(50,1) T']; 04/15 08:23
20F：→ vincent59:你就知道我想表达的?? 04/15 08:24
21F：→ mp19990920:你如果仔细看看就会发现，i回圈里面没有用到j，所以才 04/15 08:50
22F：→ mp19990920:会有这样的结果。 04/15 08:50
23F：→ cancellarame:49 50栏的值分别是 72.7573 72.8094 04/15 11:59
24F：→ cancellarame:49 50行 vincent59大 改或不改结果都一样耶 04/15 12:02
25F：推 vincent59:看样子 越接近你要问的 04/16 09:15
26F：推 vincent59:请问mp大 i回圈里面没有用到j是指???? 04/16 09:24
27F：推 vincent59:请问原po 您说的一样是指?? 04/16 10:50
28F：→ cancellarame:抱歉 小弟能力还没法独自解决回圈问题 04/16 12:29
29F：→ cancellarame:point=[c*ones(1,50);T]; 换成 point=[c*ones(50,1) 04/16 12:30
30F：→ cancellarame:跑出的 49 50行 的值分别是 72.7573 72.8094 04/16 12:32
31F：→ cancellarame:[c*ones(1,50);T]换成point=[c*ones(50,1) T']; 04/16 12:58
32F：推 vincent59:这叫一样???? 好!咱们来一行一行来研究XD 04/16 16:25
33F：→ diva:v大加油 04/16 19:19
34F：推 vincent59:diva大应该也很强吧 xd 04/16 19:40