| 查看: 407 | 回复: 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; |
回复此楼» 猜你喜欢
医学类期刊求推荐
已经有5人回复
-
生活琐事由它去
已经有4人回复
-
提交了我也来说说感想
已经有12人回复
-
青B发送上会通知了吗
已经有9人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有6人回复
-
论文撤稿了
已经有8人回复
-
化学专业申博
已经有4人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
-
26/27申博自荐
已经有9人回复
-
博士申请
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 高维数组插值 矢量化 避免for循环
已经有5人回复
- 【活动or资源】vasp晶体结构优化
已经有142人回复
2楼2012-12-05 15:15:35