| 5 | 1/1 | 返回列表 |
| 查看: 1264 | 回復: 18 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
wzp0098木蟲 (著名寫手)
|
[求助]
用matlab 畫三維圖,請高手指教. 已有2人參與
|
||
|
程序如下,可是沒有辦法運行,請高手指教,重謝。 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 |
榮譽版主 (著名寫手)
![]() |
專家經驗: +4 |

木蟲 (著名寫手)
|
我改了一下,程序還是運行不起來,請大蝦指教! 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 |
木蟲 (著名寫手)
|
我改了一下,程序還是運行不起來,請大蝦指教! 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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料求調劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0703化學 305求調劑 +5 | FY_yy 2026-03-14 | 5/250 |
|
|
[考研] 材料工程專碩 348分求調劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 0703化學調劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 306求0703調劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
[考研] 346求調劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 265求調劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 一志愿華南師大 070300(化學)304分求調劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 一志愿南昌大學,327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 324求調劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
[考研] 一志愿蘇州大學材料求調劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 0856調劑,是學校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[論文投稿]
申請回稿延期一個月,編輯同意了。但系統(tǒng)上的時間沒變,給編輯又寫郵件了,沒回復
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 收復試調劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 312求調劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 283求調劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 東南大學364求調劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|