版塊導(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í)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

田山東

捐助貴賓 (著名寫手)

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

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

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

回復(fù)此樓

» 猜你喜歡

295材料工程專碩求調(diào)劑已經(jīng)有3人回復(fù)
085701環(huán)境工程求調(diào)劑已經(jīng)有5人回復(fù)
083000學(xué)碩274求調(diào)劑已經(jīng)有4人回復(fù)
275求調(diào)劑已經(jīng)有8人回復(fù)
295求調(diào)劑已經(jīng)有4人回復(fù)
一志愿211院校 344分東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑已經(jīng)有8人回復(fù)
085602 307分求調(diào)劑已經(jīng)有5人回復(fù)
322求調(diào)劑已經(jīng)有3人回復(fù)
高溫高壓反應(yīng)求助已經(jīng)有3人回復(fù)
一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑已經(jīng)有4人回復(fù)

everythinghasitsseason.enjoyyourlife.

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

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

田山東

捐助貴賓 (著名寫手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時(shí)
蟲號(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的回帖

田山東

捐助貴賓 (著名寫手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時(shí)
蟲號(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的回帖

田山東

捐助貴賓 (著名寫手)

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

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

贊一下

回復(fù)此樓

everythinghasitsseason.enjoyyourlife.

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

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

田山東

捐助貴賓 (著名寫手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時(shí)
蟲號(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

至尊木蟲 (著名寫手)

驃騎將軍

程序強(qiáng)帖: 40
應(yīng)助: 817 (博后)
金幣: 12914.1
紅花: 64
帖子: 2238
在線: 287.3小時(shí)
蟲號(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又是什么?前面沒有定義啊

贊一下

回復(fù)此樓

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

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

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

田山東

捐助貴賓 (著名寫手)

應(yīng)助: 5 (幼兒園)
金幣: 77
散金: 739
紅花: 8
帖子: 1577
在線: 138.1小時(shí)
蟲號(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又是什么?前面沒有定義啊

謝謝你了。你說(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

木蟲 (職業(yè)作家)

中國(guó)特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時(shí)
蟲號(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í)生活，逃離自己飄忽不定的七情六欲的桎梏。--愛因斯坦

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

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

libralibra

至尊木蟲 (著名寫手)

驃騎將軍

程序強(qiáng)帖: 40
應(yīng)助: 817 (博后)
金幣: 12914.1
紅花: 64
帖子: 2238
在線: 287.3小時(shí)
蟲號(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;
開始循環(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根本沒有用處.
那么修改一下: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寫程序請(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ā)表

[考研] 083000學(xué)碩274求調(diào)劑 +4	Li李魚 2026-03-26	4/200	2026-03-27 19:42 by ACS Nano——
[考研] 一志愿西北大學(xué) 總分282 英語(yǔ)一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 287求調(diào)劑 +10	land xuxu 2026-03-26	10/500	2026-03-27 15:33 by 帕爾馬拉特
[考研] 322求調(diào)劑 +4	我真的很想學(xué)習(xí) 2026-03-23	4/200	2026-03-27 13:51 by 楊楊楊紫
[考研] 材料學(xué)碩333求調(diào)劑 +8	北道巷 2026-03-24	8/400	2026-03-27 10:18 by 我是小康
[考研] 325求調(diào)劑 +5	李嘉圖·S·路 2026-03-23	5/250	2026-03-27 00:42 by wxiongid
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 325求調(diào)劑 +3	Aoyijiang 2026-03-23	3/150	2026-03-26 20:46 by 不吃魚的貓
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 271求調(diào)劑 +6	生如夏花… 2026-03-22	6/300	2026-03-26 16:48 by 張凱十八號(hào)
[考研] 尋找調(diào)劑 +5	倔強(qiáng)芒? 2026-03-21	8/400	2026-03-26 13:25 by 0906ljy
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5	軒小寧—— 2026-03-26	5/250	2026-03-26 12:05 by i_cooler
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機(jī) 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 一志愿北化315 求調(diào)劑 +3	akrrain 2026-03-24	3/150	2026-03-24 19:35 by 了了了了。。
[考研] 307求調(diào)劑 +3	余意卿 2026-03-21	6/300	2026-03-24 15:03 by 余意卿
[考研] 085404電子信息284分求調(diào)劑 +4	13659058978 2026-03-24	4/200	2026-03-24 12:15 by syl20081243
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4	allen-yin 2026-03-23	6/300	2026-03-23 15:04 by 汪��？！
[考研] 291求調(diào)劑 +5	孅華 2026-03-22	5/250	2026-03-23 09:20 by haoshis
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng

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

[求助] matlab矢量寫法，這樣可以嗎？急求啊�。�

» 猜你喜歡

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

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

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