| 查看: 426 | 回复: 1 | ||
[求助]
matlab优化去掉for，用矢量化运行，这样为什么不行呢？
|
|
clc; clear; ny0=64; fi=2*pi/16; t=2.8; I=1i; H0=zeros(ny0,ny0); H1=zeros(ny0,ny0); for j=1:ny0 if(mod(j,2)==1) m=(j+1)/2; else m=j/2; end Ax=-(m-1)*fi*0.5; if (mod(j,4)==1) H0(j,j+1)=t*exp(I*Ax); if(j==1) H0(j,ny0)=t; else H0(j,j-1)=t; end end end 这里for换成矢量运算，有个哥们给个算法，但是没看懂，请指教： tmpcond=mod((1:ny0),4)==1; H(tmpcond,[false;tmpcond(1:end-1)])=t*exp(I*Ax); 我试着运行一个简单的： tmpcond=mod((1:ny0),4)==1; H1(tmpcond,[false;tmpcond(1:end-1)])=t; 结果出错了： Error using vertcat CAT arguments dimensions are not consistent. Error in youhuaxuexi (line 30) H1(tmpcond,[false;tmpcond(1:end-1)])=t; |
回复此楼» 猜你喜欢
fileCode有新解读？
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 高维数组插值 矢量化 避免for循环
已经有5人回复
- 【活动or资源】vasp晶体结构优化
已经有142人回复
2楼2012-12-05 15:15:35
100