 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 2429 | 回復(fù): 2 | ||
菲斐新蟲 (初入文壇)
|
[求助]
求matlab中 diff Difference order N must be a positive integer scalar什么意思 已有1人參與
|
|
我是matlab新手，不知道程序中為什么會(huì)出現(xiàn)??? Error using ==> diff Difference order N must be a positive integer scalar. Error in ==> chern at 90 phpx2=diff(H2,'x'); 我只是把第一個(gè)程序中的H變了一下，就有錯(cuò)誤，可是如果把這2段分到2個(gè)m文件里，就木有問題，咋回事呀。大神幫忙呀，程序如下 syms x y a=1.0; t=1; t1=0.1*t; u0=0.6*sqrt(3)*t; u2=0.6*sqrt(3)*t; f=2*cos(x*a/2)*cos(sqrt(3)*a*y/2)-2*i*cos(sqrt(3)*a*y/2)*sin(a*x/2)+cos(x*a)+i*sin(x*a); g=-i*(2*i*sin(sqrt(3)*a*y)-4*i*sin(sqrt(3)*a*y/2+3*a*x/2)); H=[u0/2+t1*g-u2/2,2*t*f;2*t*conj(f),-u0/2-t1*g+u2/2]; phpx=diff(H,'x'); phpy=diff(H,'y'); h=6.63*10^(-34);%此為普朗克常數(shù) hall=0; for kx = -2*pi/(3*a) : 4*pi/(30*a) : 2*pi/(3*a) for ky = -4*pi*sqrt(3)/(9*a) : 8*pi*sqrt(3)/(90*a) : 4*pi*sqrt(3)/(9*a) if -2*pi/(3*a) <= kx <= 0 && (-sqrt(3)*kx/3-4*pi*sqrt(3)/(9*a)) <= ky <=(sqrt(3)*kx/3+4*pi*sqrt(3)/(9*a)) % f=exp(-i*((kx*a/2)+sqrt(3)*ky*a/2))+ exp(-i*((kx*a/2)-sqrt(3)*ky*a/2))+exp(-i*(-kx*a)) ; %g1=exp(i*(ky*sqrt(3)*a))-exp(i*(-kx*3*a/2+ky*sqrt(3)*a/2))+exp(i*(-kx*3*a/2-ky*sqrt(3)*a/2)); %g2=-exp(i*(-ky*sqrt(3)*a))+exp(i*(kx*3*a /2-ky*sqrt(3)*a/2))-exp(i*(kx*3*a/2+ky*sqrt(3)*a/2)); %g=-i*(g1+g2); f=2*cos(kx*a/2)*cos(sqrt(3)*a*ky/2)-2*i*cos(sqrt(3)*a*ky/2)*sin(a*kx/2)+cos(kx*a)+i*sin(kx*a); g=-i*(2*i*sin(sqrt(3)*a*ky)-4*i*sin(sqrt(3)*a*ky/2+3*a*kx/2)); H=[u0/2+t1*g-u2/2,2*t*f;2*t*conj(f),-u0/2-t1*g+u2/2]; [P,E]=eig(H); mk=[P(1,1);P(2,1)]; %用eig求得的本征矢和本征值不是按從小到大的順序排列的 nk=[P(1,2);P(2,2)]; Em=E(1,1); En=E(2,2); vx=subs(phpx,{x,y},{kx,ky}); vy=subs(phpy,{x,y},{kx,ky}); s=i*exp(2)*9*a^2/(4*sqrt(3)*pi*h);%此為系數(shù) if Em<En ha=s*((mk'*vx*nk)*(nk'*vy*mk)-(mk'*vy*nk)*(nk'*vx*mk)) /(Em-En)^2; else ha=s*((nk'*vx*mk)*(mk'*vy*nk)-(nk'*vy*mk)*(mk'*vx*nk)) /(En-Em)^2; end hall=hall+ha; end if 0 <= kx <= 2*pi/(3*a) && (sqrt(3)*kx/3-4*pi*sqrt(3)/(9*a)) <= ky <=(-sqrt(3)*kx/3+4*pi*sqrt(3)/(9*a)) %f=exp(-i*((kx*a/2)+sqrt(3)*ky*a/2))+ exp(-i*((kx*a/2)-sqrt(3)*ky*a/2))+exp(-i*(-kx*a)) ; %g1=exp(i*(ky*sqrt(3)*a))-exp(i*(-kx*3*a/2+ky*sqrt(3)*a/2))+exp(i*(-kx*3*a/2-ky*sqrt(3)*a/2)); %g2=-exp(i*(-ky*sqrt(3)*a))+exp(i*(kx*3*a /2-ky*sqrt(3)*a/2))-exp(i*(kx*3*a/2+ky*sqrt(3)*a/2)); %g=-i*(g1+g2); f=2*cos(kx*a/2)*cos(sqrt(3)*a*ky/2)-2*i*cos(sqrt(3)*a*ky/2)*sin(a*kx/2)+cos(kx*a)+i*sin(kx*a); g=-i*(2*i*sin(sqrt(3)*a*ky)-4*i*sin(sqrt(3)*a*ky/2+3*a*kx/2)); H=[u0/2+t1*g-u2/2,2*t*f;2*t*conj(f),-u0/2-t1*g+u2/2]; [P,E]=eig(H); mk=[P(1,1);P(2,1)]; nk=[P(1,2);P(2,2)]; Em=E(1,1); En=E(2,2); vx=subs(phpx,{x,y},{kx,ky}); vy=subs(phpy,{x,y},{kx,ky}); s=i*exp(2)*9*a^2/(4*sqrt(3)*pi*h); if Em<En ha=s*((mk'*vx*nk)*(nk'*vy*mk)-(mk'*vy*nk)*(nk'*vx*mk)) /(Em-En)^2; else ha=s*((nk'*vx*mk)*(mk'*vy*nk)-(nk'*vy*mk)*(mk'*vx*nk)) /(En-Em)^2; end hall=hall+ha; end end end c1=hall*h/(2*pi*exp(2));%第一段 H2=[u0/2-t1*g+u2/2,2*t*f;2*t*conj(f),-u0/2+t1*g-u2/2]; phpx2=diff(H2,'x'); phpy2=diff(H2,'y'); hall2=0; for kx = -2*pi/(3*a) : 4*pi/(30*a) : 2*pi/(3*a) for ky = -4*pi*sqrt(3)/(9*a) : 8*pi*sqrt(3)/(90*a) : 4*pi*sqrt(3)/(9*a) if -2*pi/(3*a) <= kx <= 0 && (-sqrt(3)*kx/3-4*pi*sqrt(3)/(9*a)) <= ky <=(sqrt(3)*kx/3+4*pi*sqrt(3)/(9*a)) % f=exp(-i*((kx*a/2)+sqrt(3)*ky*a/2))+ exp(-i*((kx*a/2)-sqrt(3)*ky*a/2))+exp(-i*(-kx*a)) ; %g1=exp(i*(ky*sqrt(3)*a))-exp(i*(-kx*3*a/2+ky*sqrt(3)*a/2))+exp(i*(-kx*3*a/2-ky*sqrt(3)*a/2)); %g2=-exp(i*(-ky*sqrt(3)*a))+exp(i*(kx*3*a /2-ky*sqrt(3)*a/2))-exp(i*(kx*3*a/2+ky*sqrt(3)*a/2)); %g=-i*(g1+g2); f=2*cos(kx*a/2)*cos(sqrt(3)*a*ky/2)-2*i*cos(sqrt(3)*a*ky/2)*sin(a*kx/2)+cos(kx*a)+i*sin(kx*a); g=-i*(2*i*sin(sqrt(3)*a*ky)-4*i*sin(sqrt(3)*a*ky/2+3*a*kx/2)); H2=[u0/2-t1*g+u2/2,2*t*f;2*t*conj(f),-u0/2+t1*g-u2/2]; [P,E]=eig(H2); mk=[P(1,1);P(2,1)]; %用eig求得的本征矢和本征值不是按從小到大的順序排列的 nk=[P(1,2);P(2,2)]; Em=E(1,1); En=E(2,2); vx2=subs(phpx2,{x,y},{kx,ky}); vy2=subs(phpy2,{x,y},{kx,ky}); s=i*exp(2)*9*a^2/(4*sqrt(3)*pi*h);%此為系數(shù) if Em<En ha=s*((mk'*vxb*nk)*(nk'*vyb*mk)-(mk'*vyb*nk)*(nk'*vxb*mk)) /(Em-En)^2; else ha=s*((nk'*vxb*mk)*(mk'*vyb*nk)-(nk'*vyb*mk)*(mk'*vxb*nk)) /(En-Em)^2; end hall2=hallb+ha; end if 0 <= kx <= 2*pi/(3*a) && (sqrt(3)*kx/3-4*pi*sqrt(3)/(9*a)) <= ky <=(-sqrt(3)*kx/3+4*pi*sqrt(3)/(9*a)) %f=exp(-i*((kx*a/2)+sqrt(3)*ky*a/2))+ exp(-i*((kx*a/2)-sqrt(3)*ky*a/2))+exp(-i*(-kx*a)) ; %g1=exp(i*(ky*sqrt(3)*a))-exp(i*(-kx*3*a/2+ky*sqrt(3)*a/2))+exp(i*(-kx*3*a/2-ky*sqrt(3)*a/2)); %g2=-exp(i*(-ky*sqrt(3)*a))+exp(i*(kx*3*a /2-ky*sqrt(3)*a/2))-exp(i*(kx*3*a/2+ky*sqrt(3)*a/2)); %g=-i*(g1+g2); f=2*cos(kx*a/2)*cos(sqrt(3)*a*ky/2)-2*i*cos(sqrt(3)*a*ky/2)*sin(a*kx/2)+cos(kx*a)+i*sin(kx*a); g=-i*(2*i*sin(sqrt(3)*a*ky)-4*i*sin(sqrt(3)*a*ky/2+3*a*kx/2)); H2=[u0/2-t1*g+u2/2,2*t*f;2*t*conj(f),-u0/2+t1*g-u2/2]; [P,E]=eig(H2); mk=[P(1,1);P(2,1)]; nk=[P(1,2);P(2,2)]; Em=E(1,1); En=E(2,2); vx2=subs(phpx2,{x,y},{kx,ky}); vy2=subs(phpy2,{x,y},{kx,ky}); s=i*exp(2)*9*a^2/(4*sqrt(3)*pi*h); if Em<En ha=s*((mk'*vxb*nk)*(nk'*vyb*mk)-(mk'*vyb*nk)*(nk'*vxb*mk)) /(Em-En)^2; else ha=s*((nk'*vxb*mk)*(mk'*vyb*nk)-(nk'*vyb*mk)*(mk'*vxb*nk)) /(En-Em)^2; end hall2=hall2+ha; end end end c2=hallb*h/(2*pi*exp(2)); c=c1+c2;%第二段 |
回復(fù)此樓» 猜你喜歡
材料與化工專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
086500 325 求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料080500調(diào)劑求收留
已經(jīng)有7人回復(fù)
-
一志愿南昌大學(xué)，327分，材料與化工085600
已經(jīng)有6人回復(fù)
-
環(huán)境工程調(diào)劑
已經(jīng)有9人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有9人回復(fù)
-
287求調(diào)劑
已經(jīng)有5人回復(fù)
-
22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù)
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有3人回復(fù)
-
317求調(diào)劑
已經(jīng)有10人回復(fù)
zhchh008
金蟲 (正式寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 1602.2
- 散金: 871
- 紅花: 13
- 帖子: 835
- 在線: 202.7小時(shí)
- 蟲號(hào): 99456
- 注冊(cè): 2005-11-11
- 性別: GG
- 專業(yè): 資源化工
2樓2014-10-30 19:53:58
菲斐
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 21
- 帖子: 2
- 在線: 1.5小時(shí)
- 蟲號(hào): 3247358
- 注冊(cè): 2014-06-01
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2014-11-04 16:03:16
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 菲斐 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +3 | solanXXX 2026-03-20 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +8 | 梁梁校校 2026-03-17 | 8/400 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-18 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +5 | 曼殊2266 2026-03-18 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級(jí)已過，有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|