| 查看: 599 | 回复: 2 | ||
身心天地铁杆木虫 (正式写手)
|
[求助]
索引超出矩阵维度，for循环i值的设置
|
|
t=(1:1:200) counter=zeros(1,200); x(1)=0; y(1)=0; z(1)=0; x(2)=0; y(2)=0; z(2)=0.01; n=0; m=0; while(n<100000) s=0; h=0; for i=3:1:80 l=(-333)*log(rand(1)); b=acos(2*rand(1)-1); e=2*pi*rand(1); x(i)=x(i-1)+l*sin(b)*sin(e); y(i)=y(i-1)+l*sin(b)*cos(e); z(i)=z(i-1)+l*cos(b); a=z(i)-z(i-1); s=s+a; h=h+l; m=m+1 n; if(s)>900 x(i); y(i) m; break end end t=round(abs(h)/60); if abs(x(i))<25000 if abs(y(i))<25000 counter(t)=counter(t)+1; p=x(i); q=y(i) end end n=n+1; end plot(t, counter(t)); 代码如上，问题：就是这个for循环里面的i设置for i=3:1:80，我设置i等于不超过五十左右时，这个程序跑的通，当超过这个数值时就跳出索引超出矩阵维度。求解决，新手摸索两个月了，感激不尽，做不完，都没法回家过年了，唉，跪谢！ |
回复此楼» 猜你喜欢
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有188人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有11人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有6人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
465090264
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4904.3
- 红花: 2
- 帖子: 754
- 在线: 52.5小时
- 虫号: 2070656
- 注册: 2012-10-18
- 性别: GG
- 专业: 应用高分子化学与物理
|
把下面几行连起来看看，你的t超出了counter的范围（这代码没格式，读起来好费劲）counter=zeros(1,200); l=(-333)*log(rand(1));h=h+l;t=round(abs(h)/60);counter(t)=counter(t)+1; 发自小木虫Android客户端 |
2楼2019-01-27 00:44:51
身心天地
铁杆木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 7643.3
- 散金: 56
- 红花: 1
- 帖子: 386
- 在线: 342小时
- 虫号: 2655143
- 注册: 2013-09-14
- 性别: GG
- 专业: 机械摩擦学与表面技术
3楼2019-01-28 21:02:29
