24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

田山东

捐助贵宾 (著名写手)

应助: 5 (幼儿园)
金币: 77
散金: 739
红花: 8
帖子: 1577
在线: 138.1小时
虫号: 800307
注册: 2009-06-28
性别: GG
专业: 基础物理学

[求助] matlab矢量写法，这样可以吗？急求啊！！

大家好，在论坛上和网上也学到了一些将for语句写成矢量的写法，但是我的这个里面有if句法，如何来实现呢？我试了下，没有成功
就这样几句话：
ny0=64;
fi=2*pi/16;
t=2.8;
I=1i;
H0=zeros(ny0,ny0);
for j=1:ny0;
if(mod(jl,2)==1)
m=(j(l)+1)/2;
else
m=j(l)/2;
end
Ax=-(m-1)*0.5*fi;
   if (mod(j,4)==1)
   H0(j,j+1)=t*exp(I*Ax);
   end
end
我试着这样去做：

ny0=64;
fi=2*pi/16;
t=2.8;
I=1i;
H0=zeros(ny0,ny0);
j=1:ny0;
l=1:1:length(j);
if(mod(j(l),2)==1)
m=(j(l)+1)/2;
else
m=j(l)/2;
end
Ax=-(m-1)*0.5*fi;
   if (mod(j(l),4)==1)
   H0(j(l),j(l)+1)=t*exp(I*Ax);
   end
但是很奇怪啊，里面没数据啊。
有没有人可以帮我下啊？

回复此楼

» 猜你喜欢

实验室接单子已经有4人回复
全日制（定向）博士已经有4人回复
假如你的研究生提出不合理要求已经有6人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复
26申博（荧光探针方向，有机合成）已经有4人回复

everythinghasitsseason.enjoyyourlife.

1楼 2012-12-05 19:58:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

【答案】应助回帖

引用回帖:

7楼: Originally posted by 田山东 at 2012-12-06 09:03:01
谢谢你了。你说的是上面那几行语句，你看下下面那几行，l=1:1：length（j）表示j的长度...

j=1:ny0;
循环标志l = 1:length(j),所以j(l)==l
你的代码就是

CODE:

 ny0=64;

fi=2*pi/16;

t=2.8;

I=1i;

H0=zeros(ny0,ny0);

j=1:ny0;

for l=1:length(j)

    if(mod(l,2)==1)

        m=(l+1)/2;

    else

        m=l/2;

    end

    Ax=-(m-1)*0.5*fi;

    if (mod(l,4)==1)

        H0(l,l+1)=t*exp(I*Ax);

    end

end

然后分析一下:
j = 1:64;
开始循环j的元素l
如果是奇数,自加1除以2,如果是偶数,除以2,赋值给m
根据m得到Ax
如果当前元素除4余1,给H0的(l,l+1)元素赋值.
j中所有除以4余1的元素是

CODE:

 1     5     9    13    17    21    25    29    33    37    41    45    49    53    57    61

一共16个,等于其实你只需要给H0的这16个元素赋值

CODE:

 (1,2)     (5,6)     (9,10)    (13,14)    (17,18)    (21,22)    (25,26)    (29,30)    (33,34)    (37,38)    (41,42)    (45,46)    (49,50)    (53,54)    (57,58)    (61,62)

既然给H0赋值的元素个数只有16个,为什么要循环64次呢?等于有48次循环计算的m和Ax根本没有用处.
那么修改一下:revise 1

CODE:

for l=1:length(j)

    if (mod(l,4)==1)

        m=(l+1)/2;

        Ax=-(m-1)*0.5*fi;

        H0(l,l+1)=t*exp(I*Ax);

    end

end

既然mod(l,4)==1,就有l=4*n+1=2*2n+1,说明l除以2也余1,直接计算m然后给H0元素赋值.

revise2,matlab的for可以指定步长,连if判断都不需要了

CODE:

for l=1:4:length(j)

    m=(l+1)/2;

    Ax=-(m-1)*0.5*fi;

    H0(l,l+1)=t*exp(I*Ax);

end

向量化的办法:
revise3,既然是给全0矩阵H0的16个元素赋值,那你直接运用向量化计算出来结果,然后赋值就行了.不过可能有些绕,需要用好几个函数来实现
可以通过sparse来实现,不过sparse之后用full只能还原到61*62维,因为其他元素都是0
所以需要用局部矩阵赋值的办法

CODE:

 ny0=64;

fi=2*pi/16;

t=2.8;

I=1i;

H0=zeros(ny0,ny0);

j=1:ny0;

l = j(mod(j,4)==1); % 行下标,16个

n = l+1; % 列下标,16个

m = (l+1)/2; % m,16个

Ax = -(m-1)*0.5*fi; % Ax,16个

expAx = t*exp(I*Ax); % H0的被赋值元素,16个

h = sparse(l,n,expAx); % 根据下标和元素生成稀疏矩阵

H0(1:max(l),1:max(n)) = full(h); % 给H0对应位置赋值

revise 4,终极向量化办法,得到行列号后,直接得到下标序列,将H0看做一维向量赋值

CODE:

 ny0=64;

fi=2*pi/16;

t=2.8;

I=1i;

H0=zeros(ny0,ny0);

j=1:ny0;

l = j(mod(j,4)==1); % 行下标,16个

n = l+1; % 列下标,16个

m = (l+1)/2; % m,16个

Ax = -(m-1)*0.5*fi; % Ax,16个

expAx = t*exp(I*Ax); % H0的被赋值元素,16个

H0(sub2ind([ny0,ny0],l,n)) = expAx; % 直接给H0对应位置赋值