 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3207 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
bdx1989新蟲 (初入文壇)
|
[求助]
matlab編程中關(guān)于for語句中[V,D]=eig(A)的用法 已有1人參與
|
|
|
我在編程中有一段語句如下 for J=0:Jmax hsym=zeros(2*J+1); %%規(guī)定零矩陣大小，其大小隨每次循環(huán)的J改變而改變 for Kp=-J:J m=m+1; n=0; for K=-J:J n=n+1; if Kp==K hsym(n,m)=F*(J*(J+1)-K^2)+G*K^2; elseif Kp==K+2 hsym(n,m)=H*(0.25*(J*(J+1)-K*(K+1))*(J*(J+1)-(K+1)*(K+2)))^(0.5); elseif Kp==K-2 hsym(n,m)=H*(0.25*(J*(J+1)-K*(K-1))*(J*(J+1)-(K-1)*(K-2)))^(0.5); end end end %%以上是給矩陣hsym的某些矩陣元賦值。沒有問題 [evrp,evlp]=eig(hsym); %%我需要矩陣的本證矢量，所以需要用到[V,D]=eig(A)的形式。 for e=1  2*J+1)evl(J^2+e,J^2+e)=0.5*(A+C)*J*(J+1)+0.5*(A-C)*evlp(e,e); for ee=1 2*J+1)evr(J^2+e,ee)=evrp(e,ee); end end %%以上是將hsym所得的本征值和本證矢量轉(zhuǎn)移到新矩陣中 end 但是這里運(yùn)行的時候一直有問題，如果把最前面第一個for去掉的話只計(jì)算一個矩陣的本征值和本證矢量，這段程序是沒有問題的。但是加上第一個for語句就出問題了。第一個for我的用意是將循環(huán)里的所有矩陣的本征值和本證矢量求出，最后分別整合到兩個大矩陣中，一個矩陣裝本征值，一個矩陣裝本證矢量。如果將第一個for語句去掉單獨(dú)算一個矩陣，沒有問題，而且本征值和本證矢量都可以原封不動的轉(zhuǎn)移到大矩陣中，位置也是我想要的位置沒有問題。但是加上第一個for語句就不行了。 我覺得應(yīng)該是 [evrp,evlp]=eig(hsym) 這里有問題，但是我不知道如何解決。菜鳥請大神們幫忙解決一下。含淚感謝。 |
回復(fù)此樓» 猜你喜歡
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有31人回復(fù)
-
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑
已經(jīng)有5人回復(fù)
-
085801 總分275 本科新能源 求調(diào)劑
已經(jīng)有17人回復(fù)
-
324求調(diào)劑
已經(jīng)有6人回復(fù)
-
學(xué)碩化學(xué)工程與技術(shù)，一志愿中國海洋大學(xué)320+求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料調(diào)劑
已經(jīng)有5人回復(fù)
-
初試301，代碼085701環(huán)境工程，本碩一致，四六級已過，有二區(qū)一作，共發(fā)表5篇論文
已經(jīng)有6人回復(fù)
-
化學(xué)308分調(diào)劑
已經(jīng)有15人回復(fù)
-
346求調(diào)劑
已經(jīng)有3人回復(fù)
-
環(huán)境285分，過六級，求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助Matlab批處理for語句循環(huán)
已經(jīng)有5人回復(fù)
- matlab中如何通過for語句得到已有矩陣的部分?jǐn)?shù)據(jù)
已經(jīng)有3人回復(fù)
- 跪求matlab的for循環(huán)語句的編寫
已經(jīng)有20人回復(fù)
bdx1989
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 72
- 帖子: 19
- 在線: 11.8小時
- 蟲號: 1554183
- 注冊: 2011-12-28
- 專業(yè): 原子和分子物理
3樓2014-04-24 10:41:28
Mokeer
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 20.5
- 帖子: 19
- 在線: 12.4小時
- 蟲號: 3143218
- 注冊: 2014-04-17
- 專業(yè): 機(jī)器人學(xué)及機(jī)器人技術(shù)
2樓2014-04-24 08:56:38
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +12 | 駱駝男人 2026-04-02 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 270求調(diào)劑 +8 | 小杰pp 2026-03-31 | 10/500 |
|
|
[考研] 能源動力 調(diào)劑 +3 | 不破不立0 2026-04-02 | 3/150 |
|
|
|
[考研] 求調(diào)劑推薦 +3 | 南山南@ 2026-04-01 | 3/150 |
|
|
|
[考研] 292分，材料與化工，申請調(diào)劑 +19 | 程晴之 2026-04-01 | 23/1150 |
|
|
|
[考研] 270調(diào)劑 +7 | maxjxbsk 2026-04-02 | 7/350 |
|
|
|
[考研] 0710生物學(xué)，325求調(diào)劑 +3 | mkkkkkl 2026-04-01 | 3/150 |
|
|
|
[考研] 302求調(diào)劑一志愿北航070300，本科鄭大化學(xué) +8 | 圣日耳曼條 2026-04-01 | 11/550 |
|
|
|
[考研] 求生物學(xué)調(diào)劑 +8 | 15172915737 2026-04-01 | 8/400 |
|
|
|
[考研] 英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑 +12 | 12.27 2026-03-30 | 14/700 |
|
|
|
[考研] 291求調(diào)劑 +20 | Y-cap 2026-03-29 | 25/1250 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪��？！ 2026-03-31 | 7/350 |
|
|
|
[考研] 建環(huán)，能源，土木老師路過看一看�。�！ +4 | 嘿嘿uu 2026-04-01 | 4/200 |
|
|
|
[考研] 土木304求調(diào)劑 +5 | 頂級擦擦 2026-03-31 | 5/250 |
|
|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
|
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
|
[考研] 312，生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|