作者jeffppp (烦..)
看板MATLAB
标题Re: [问题] 回圈
时间Fri Nov 16 17:46:39 2012
※ 引述《james1022jk (雁)》之铭言：
: ※ 引述《icrtstyle (嘿嘿)》之铭言：
: : 我想请问一下,要怎麽把以下3个回圈合并在一起????
: : ----------
: : for i=1:n;
: : t=1;
: : x(i,t+1)=x(i,t)+g(i,t);
: : for j=1:n;
: : x(i,t+1)=x(i,t+1)+e*aa(i,j)*x(j,t+1);
: : end
: : end
: : for i=1:n;
: : t=2;
: : x(i,t+1)=x(i,t)+g(i,t);
: : for j=1:n;
: : x(i,t+1)=x(i,t+1)+e*aa(i,j)*x(j,t+1)+e*(exp(-1))*aa(i,j)*(x(j,t-1)+g(j,t-1));
: : end
: : end
: : for i=1:n;
: : t=3;
: : x(i,t+1)=x(i,t)+g(i,t);
: : for j=1:n;
: : x(i,t+1)=x(i,t+1)+e*aa(i,j)*x(j,t+1)+e*(exp(-1))*aa(i,j)*(x(j,t-1)+g(j,t-1))+
: : e*(exp(-2))*aa(i,j)*(x(j,t-2)+g(j,t-2));
: : end
: : end
: : -----------------------------
: for t = 1:3
: for i = 1:n
: x(i,t+1) = x(i,t)+g(i,t);
: for j = 1:n
: switch t
: case 1
: x(i,t+1)=x(i,t+1)+e*aa(i,j)*x(j,t+1);
: case 2
: x(i,t+1)=x(i,t+1)+e*aa(i,j)*x(j,t+1)+...
: e*(exp(-1))*aa(i,j)*(x(j,t-1)+g(j,t-1));
: case 3
: x(i,t+1)=x(i,t+1)+e*aa(i,j)*x(j,t+1)+...
: e*(exp(-1))*aa(i,j)*(x(j,t-1)+g(j,t-1))+...
: e*(exp(-2))*aa(i,j)*(x(j,t-2)+g(j,t-2));
: end
: end
: end
: end
: 我果然是太无聊..........................................................
for i=1:n;
for t=1:100
x(i,t+1)=x(i,t)+g(i,t);
for j=1:n;
temp=x(i,t+1);
for k=1:t
if k==1
temp=temp+x(i,t+1)+e*aa(i,j)*x(j,t+1);
else
temp=temp+e*(exp(1-k))*aa(i,j)*(x(j,t+(1-k))+g(j,t+(1-k)))
end
end
x(i,t+1)=temp;
end
end
end
不确定对不对 XDD
太多层了
可以把100 带小一点的数字试试看
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 163.25.95.251
1F：推 icrtstyle:太感谢!!!几乎都是对的,只有一个小小的瑕疵而已~~~ 11/17 02:05
2F：→ icrtstyle:就是..temp=temp+x(i,t+1)+e*aa(i,j)*x(j,t+1);这行 11/17 02:06
3F：→ icrtstyle:要改成..temp=temp+e*aa(i,j)*x(j,t+1);就这一点点问题 11/17 02:07
4F：→ icrtstyle:感谢你让我学到这种修改回圈方式~~~^^ 11/17 02:08