 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 3203 | 回復(fù): 2 | ||
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)行的時(shí)候一直有問題，如果把最前面第一個(gè)for去掉的話只計(jì)算一個(gè)矩陣的本征值和本證矢量，這段程序是沒有問題的。但是加上第一個(gè)for語句就出問題了。第一個(gè)for我的用意是將循環(huán)里的所有矩陣的本征值和本證矢量求出，最后分別整合到兩個(gè)大矩陣中，一個(gè)矩陣裝本征值，一個(gè)矩陣裝本證矢量。如果將第一個(gè)for語句去掉單獨(dú)算一個(gè)矩陣，沒有問題，而且本征值和本證矢量都可以原封不動(dòng)的轉(zhuǎn)移到大矩陣中，位置也是我想要的位置沒有問題。但是加上第一個(gè)for語句就不行了。 我覺得應(yīng)該是 [evrp,evlp]=eig(hsym) 這里有問題，但是我不知道如何解決。菜鳥請(qǐng)大神們幫忙解決一下。含淚感謝。 |
回復(fù)此樓» 猜你喜歡
081700 調(diào)劑 267分
已經(jīng)有8人回復(fù)
-
請(qǐng)教下大家 2026年國家基金申請(qǐng)是雙盲審嗎？
已經(jīng)有5人回復(fù)
-
276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷
已經(jīng)有10人回復(fù)
-
一志愿南航材料專317分求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)學(xué)碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助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ù)
Mokeer
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 20.5
- 帖子: 19
- 在線: 12.4小時(shí)
- 蟲號(hào): 3143218
- 注冊(cè): 2014-04-17
- 專業(yè): 機(jī)器人學(xué)及機(jī)器人技術(shù)
2樓2014-04-24 08:56:38
bdx1989
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 72
- 帖子: 19
- 在線: 11.8小時(shí)
- 蟲號(hào): 1554183
- 注冊(cè): 2011-12-28
- 專業(yè): 原子和分子物理
3樓2014-04-24 10:41:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 bdx1989 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +10 | Zhao anqi 2026-03-22 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研] 材料調(diào)劑 +3 | 匹克i 2026-03-23 | 3/150 |
|
|
|
[考研] 材料專碩找調(diào)劑 +3 | 哈哈哈吼吼吼哈 2026-03-23 | 3/150 |
|
|
|
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +3 | SYA！ 2026-03-23 | 3/150 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛生物生物愛?/a> 2026-03-17 | 5/250 |
|
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|