| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1258 | 回復(fù): 18 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
wzp0098木蟲(chóng) (著名寫(xiě)手)
|
[求助]
用matlab 畫(huà)三維圖,請(qǐng)高手指教. 已有2人參與
|
|
|
程序如下,可是沒(méi)有辦法運(yùn)行,請(qǐng)高手指教,重謝。 W=-3:0.1:3;K=-3:0.1:3i;V=-3:0.1:3; P = []; for n=1:length(W); Ky=W(n); for m=1:length(K); Kx=K(m); for v=1:length(V); Kz=V(v); w = 1; K2 = K1+2*w; x = 2*(sin(Kx)+sin(Ky)); r1 = 8; r2 = 0.001*r1; T = 3*sin(Kz); A=[ -r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B = [-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y = A\B; P(n, m, v) = imag(Y(3)); % <---------------------------P向量中的元素賦值 end end end [X,Y,Z]=meshgrid(W,K,V); surf(X,Y,Z,P); shading interp |
木蟲(chóng) (著名寫(xiě)手)
金蟲(chóng) (小有名氣)
木蟲(chóng) (著名寫(xiě)手)
|
我改了一下,程序還是運(yùn)行不起來(lái),請(qǐng)大蝦指教! W=-3:0.1:3;K=-3:0.1:3;V=-3:0.1:3; P = []; for n=1:length(W); Ky=W(n); for m=1:length(K); Kx=K(m); for v=1:length(V); Kz=V(v); w = 1;K1=0; K2 = K1+2*w; x = 2*(sin(Kx)+sin(Ky)); r1 = 8; r2 = 0.001*r1; T = 3*sin(Kz); A=[ -r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B = [-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y = A\B; P(n, m, v) = imag(Y(3)); % <---------------------------P向量中的元素賦值 end end end [X,Y,Z]=meshgrid(W,K,V); surf(X,Y,Z,P); shading interp |
木蟲(chóng) (著名寫(xiě)手)
|
我改了一下,程序還是運(yùn)行不起來(lái),請(qǐng)大蝦指教! W=-3:0.1:3;K=-3:0.1:3;V=-3:0.1:3; P = []; for n=1:length(W); Ky=W(n); for m=1:length(K); Kx=K(m); for v=1:length(V); Kz=V(v); w = 1;K1=0; K2 = K1+2*w; x = 2*(sin(Kx)+sin(Ky)); r1 = 8; r2 = 0.001*r1; T = 3*sin(Kz); A=[ -r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B = [-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y = A\B; P(n, m, v) = imag(Y(3)); % <---------------------------P向量中的元素賦值 end end end [X,Y,Z]=meshgrid(W,K,V); surf(X,Y,Z,P); shading interp |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開(kāi)始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|