版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

田山東

捐助貴賓 (著名寫(xiě)手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時(shí)
蟲(chóng)號(hào): 800307
注冊(cè): 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

[求助] matlab矢量寫(xiě)法，這樣可以嗎？急求啊�。�

大家好，在論壇上和網(wǎng)上也學(xué)到了一些將for語(yǔ)句寫(xiě)成矢量的寫(xiě)法，但是我的這個(gè)里面有if句法，如何來(lái)實(shí)現(xiàn)呢？我試了下，沒(méi)有成功
就這樣幾句話：
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
但是很奇怪啊，里面沒(méi)數(shù)據(jù)啊。
有沒(méi)有人可以幫我下啊？

回復(fù)此樓

» 猜你喜歡

一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑已經(jīng)有6人回復(fù)
材料學(xué)求調(diào)劑已經(jīng)有5人回復(fù)
考研調(diào)劑已經(jīng)有4人回復(fù)
281求調(diào)劑已經(jīng)有4人回復(fù)
0805 316求調(diào)劑已經(jīng)有6人回復(fù)
085601求調(diào)劑總分293英一數(shù)二已經(jīng)有3人回復(fù)
08工學(xué)調(diào)劑已經(jīng)有17人回復(fù)
340求調(diào)劑已經(jīng)有4人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)
食品專碩一志愿雙一流 328 已經(jīng)有4人回復(fù)

everythinghasitsseason.enjoyyourlife.

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

田山東

捐助貴賓 (著名寫(xiě)手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時(shí)
蟲(chóng)號(hào): 800307
注冊(cè): 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

給定j，得到m,得到Ax，然后對(duì)H0賦值

贊一下

回復(fù)此樓

everythinghasitsseason.enjoyyourlife.

2樓2012-12-05 20:18:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

田山東

捐助貴賓 (著名寫(xiě)手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時(shí)
蟲(chóng)號(hào): 800307
注冊(cè): 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

看了網(wǎng)上的討論，新版本中向量化跟for循環(huán)差距已經(jīng)不是很大，看來(lái)我是多慮了。
同意的頂下，也是對(duì)我摸索這么久的一個(gè)獎(jiǎng)勵(lì)

看帖不頂是小狗

回帖是一種美德

不一定要答案。

要人氣。。

贊一下

回復(fù)此樓

everythinghasitsseason.enjoyyourlife.

3樓2012-12-05 20:28:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

田山東

捐助貴賓 (著名寫(xiě)手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時(shí)
蟲(chóng)號(hào): 800307
注冊(cè): 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

可以花錢(qián)買(mǎi)BB 嗎？花點(diǎn)錢(qián)買(mǎi)點(diǎn)，給大家送點(diǎn)，攢點(diǎn)人氣啊。。

贊一下

回復(fù)此樓

everythinghasitsseason.enjoyyourlife.

4樓2012-12-05 20:29:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

田山東

捐助貴賓 (著名寫(xiě)手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時(shí)
蟲(chóng)號(hào): 800307
注冊(cè): 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

誰(shuí)能幫忙解決了，還是用下矢量化方法吧。還是希望比較一下速度。

贊一下

回復(fù)此樓

everythinghasitsseason.enjoyyourlife.

5樓2012-12-05 21:09:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

libralibra

至尊木蟲(chóng) (著名寫(xiě)手)

驃騎將軍

程序強(qiáng)帖: 40
應(yīng)助: 817 (博后)
金幣: 12914.1
紅花: 64
帖子: 2238
在線: 287.3小時(shí)
蟲(chóng)號(hào): 696514
注冊(cè): 2009-02-05
專業(yè): 計(jì)算機(jī)軟件

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1
xzhdty: 謝謝驃騎將軍 2012-12-05 21:59:11

你能仔細(xì)檢查一下的你的代碼嗎?
第一個(gè)if(mod(j(l),2)==1)少了j后面的括號(hào),但是l又是什么?前面沒(méi)有定義啊

贊一下

回復(fù)此樓

matlab/VB/python/c++/Java寫(xiě)程序請(qǐng)發(fā)QQ郵件:790404545@qq.com

6樓2012-12-05 21:33:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

田山東

捐助貴賓 (著名寫(xiě)手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時(shí)
蟲(chóng)號(hào): 800307
注冊(cè): 2009-06-28
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

引用回帖:

6樓: Originally posted by libralibra at 2012-12-05 21:33:22
你能仔細(xì)檢查一下的你的代碼嗎?
第一個(gè)if(mod(j(l),2)==1)少了j后面的括號(hào),但是l又是什么?前面沒(méi)有定義啊

謝謝你了。你說(shuō)的是上面那幾行語(yǔ)句，你看下下面那幾行，l=1:1：length（j）表示j的長(zhǎng)度

贊一下

回復(fù)此樓

everythinghasitsseason.enjoyyourlife.

7樓2012-12-06 09:03:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

baobiao007

木蟲(chóng) (職業(yè)作家)

中國(guó)特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時(shí)
蟲(chóng)號(hào): 505962
注冊(cè): 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

來(lái)捧個(gè)場(chǎng)咯

回復(fù)此樓

我同意叔本華的觀點(diǎn)，人們投身藝術(shù)和科學(xué)領(lǐng)域的強(qiáng)烈愿望之一就是逃離痛苦、殘酷和枯燥無(wú)味的現(xiàn)實(shí)生活，逃離自己飄忽不定的七情六欲的桎梏。--愛(ài)因斯坦

8樓2012-12-06 17:02:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

libralibra

至尊木蟲(chóng) (著名寫(xiě)手)

驃騎將軍

程序強(qiáng)帖: 40
應(yīng)助: 817 (博后)
金幣: 12914.1
紅花: 64
帖子: 2238
在線: 287.3小時(shí)
蟲(chóng)號(hào): 696514
注冊(cè): 2009-02-05
專業(yè): 計(jì)算機(jī)軟件

【答案】應(yīng)助回帖

引用回帖:

7樓: Originally posted by 田山東 at 2012-12-06 09:03:01
謝謝你了。你說(shuō)的是上面那幾行語(yǔ)句，你看下下面那幾行，l=1:1：length（j）表示j的長(zhǎng)度...

j=1:ny0;
循環(huán)標(biāo)志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;
開(kāi)始循環(huán)j的元素l
如果是奇數(shù),自加1除以2,如果是偶數(shù),除以2,賦值給m
根據(jù)m得到Ax
如果當(dāng)前元素除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個(gè),等于其實(shí)你只需要給H0的這16個(gè)元素賦值

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賦值的元素個(gè)數(shù)只有16個(gè),為什么要循環(huán)64次呢?等于有48次循環(huán)計(jì)算的m和Ax根本沒(méi)有用處.
那么修改一下: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(wèi)=4*n+1=2*2n+1,說(shuō)明l除以2也余1,直接計(jì)算m然后給H0元素賦值.

revise2,matlab的for可以指定步長(zhǎng),連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個(gè)元素賦值,那你直接運(yùn)用向量化計(jì)算出來(lái)結(jié)果,然后賦值就行了.不過(guò)可能有些繞,需要用好幾個(gè)函數(shù)來(lái)實(shí)現(xiàn)
可以通過(guò)sparse來(lái)實(shí)現(xiàn),不過(guò)sparse之后用full只能還原到61*62維,因?yàn)槠渌囟际?
所以需要用局部矩陣賦值的辦法

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); % 行下標(biāo),16個(gè)

n = l+1; % 列下標(biāo),16個(gè)

m = (l+1)/2; % m,16個(gè)

Ax = -(m-1)*0.5*fi; % Ax,16個(gè)

expAx = t*exp(I*Ax); % H0的被賦值元素,16個(gè)

h = sparse(l,n,expAx); % 根據(jù)下標(biāo)和元素生成稀疏矩陣

H0(1:max(l),1:max(n)) = full(h); % 給H0對(duì)應(yīng)位置賦值

revise 4,終極向量化辦法,得到行列號(hào)后,直接得到下標(biāo)序列,將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); % 行下標(biāo),16個(gè)

n = l+1; % 列下標(biāo),16個(gè)

m = (l+1)/2; % m,16個(gè)

Ax = -(m-1)*0.5*fi; % Ax,16個(gè)

expAx = t*exp(I*Ax); % H0的被賦值元素,16個(gè)

H0(sub2ind([ny0,ny0],l,n)) = expAx; % 直接給H0對(duì)應(yīng)位置賦值

贊一下

回復(fù)此樓

matlab/VB/python/c++/Java寫(xiě)程序請(qǐng)發(fā)QQ郵件:790404545@qq.com

9樓2012-12-06 22:24:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主田山東的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +5	@taotao 2026-03-19	5/250	2026-03-24 21:07 by greychen00
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +7	天天好運(yùn)來(lái)上岸?/a> 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +6	小羊睡著了? 2026-03-23	9/450	2026-03-24 20:21 by 小羊睡著了?
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑 +15	吃吃吃才有意義 2026-03-19	16/800	2026-03-24 19:51 by 了了了了。。
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 070300化學(xué)求調(diào)劑 +9	苑豆豆 2026-03-20	9/450	2026-03-24 17:15 by licg0208
[考研] 一志愿南航材料專317分求調(diào)劑 +5	炸呀炸呀炸薯?xiàng)l 2026-03-23	5/250	2026-03-24 16:52 by 星空星月
[考研] 材料與化工328分調(diào)劑 +4	。，。，。，。i 2026-03-23	4/200	2026-03-24 11:03 by 544594351
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[考研] 291求調(diào)劑 +8	hhhhxn.. 2026-03-23	8/400	2026-03-23 23:15 by peike
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[考研] 316求調(diào)劑 +7	梁茜雯 2026-03-19	7/350	2026-03-23 16:21 by lingjue
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5	困困困困坤坤 2026-03-20	6/300	2026-03-22 17:41 by hxsm
[考研] 318求調(diào)劑 +4	plum李子 2026-03-21	7/350	2026-03-22 14:17 by ColorlessPI
[考研] 286求調(diào)劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO

24小時(shí)熱門(mén)版塊排行榜

[求助] matlab矢量寫(xiě)法，這樣可以嗎？急求啊�。�

» 猜你喜歡

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] matlab矢量寫(xiě)法，這樣可以嗎？急求啊�。�