當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

wzp0098

木蟲 (著名寫手)

應(yīng)助: 2 (幼兒園)
金幣: 2524.4
散金: 39
紅花: 13
帖子: 1218
在線: 226.6小時(shí)
蟲號(hào): 853822
注冊(cè): 2009-09-22
性別: GG
專業(yè): 光學(xué)

[求助] 用matlab 畫三維圖，請(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

回復(fù)此樓

» 猜你喜歡

307求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑已經(jīng)有15人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有14人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有18人回復(fù)
265求調(diào)劑已經(jīng)有3人回復(fù)
296求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有13人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有15人回復(fù)
生物學(xué)調(diào)劑招人�。。� 已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

請(qǐng)問(wèn)三維金相立體圖如何畫出來(lái) 已經(jīng)有17人回復(fù)
如何用matlab導(dǎo)入導(dǎo)入三維數(shù)據(jù)~ 已經(jīng)有5人回復(fù)
請(qǐng)問(wèn)采用AFM掃描的文件如何用matlab畫出三維圖？已經(jīng)有7人回復(fù)
求matlab功率譜的畫法，我這樣對(duì)么？已經(jīng)有6人回復(fù)
【已完結(jié)】matlab仿真圖形（交流問(wèn)題即有獎(jiǎng)，呵呵）已經(jīng)有6人回復(fù)
菜鳥求助用matlab如何畫三維圖已經(jīng)有3人回復(fù)
急用�。。。∏螽嬋S拓?fù)鋱D 已經(jīng)有10人回復(fù)
求助已知三維坐標(biāo)位置，如何畫出空間圖像已經(jīng)有5人回復(fù)
在matlab中用mesh畫三維圖存在的問(wèn)題已經(jīng)有11人回復(fù)
MATLAB中畫三維雙對(duì)數(shù)圖形的問(wèn)題已經(jīng)有10人回復(fù)
MATLAB將三維圖和其中一個(gè)面的等高線畫出來(lái)的問(wèn)題，謝謝了。已經(jīng)有7人回復(fù)
如何用matlab畫二維網(wǎng)格已經(jīng)有5人回復(fù)
求助MATLAB做三維圖問(wèn)題已經(jīng)有4人回復(fù)
關(guān)于Matlab畫三維曲面圖的問(wèn)題，求助高手，感激感激已經(jīng)有3人回復(fù)
請(qǐng)問(wèn)如何用Matlab畫工作空間的輪廓圖已經(jīng)有10人回復(fù)
最近用HB2151表達(dá)遇到以下情況，懇請(qǐng)高手指教已經(jīng)有6人回復(fù)
Matlab三維圖中如何實(shí)現(xiàn)x，y坐標(biāo)名稱與各自的坐標(biāo)軸平行？已經(jīng)有4人回復(fù)
怎么用matlab畫極坐標(biāo)三維圖？各位大俠幫幫忙��！已經(jīng)有5人回復(fù)
如何讓MATLAB畫出來(lái)的這個(gè)三維圖更漂亮？已經(jīng)有3人回復(fù)
彈性模量各向異性的三維圖。。。求高手解決！！已經(jīng)有10人回復(fù)
【求助】菜鳥求助 matlab 繪圖問(wèn)題繪制光譜圖已經(jīng)有3人回復(fù)

1樓 2014-12-26 14:32:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huozhe123

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 698.6
帖子: 182
在線: 44小時(shí)
蟲號(hào): 3556048
注冊(cè): 2014-11-24
專業(yè): 光學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
wzp0098: 金幣+10, ★有幫助, 多謝，能幫忙改改嗎？ 2014-12-30 16:56:21

我查了一下，好像surf是平面圖，你這個(gè)是三維的，所以不能表示
http://bbs.pinggu.org/thread-1171696-1-1.html

贊一下

回復(fù)此樓

5樓2014-12-30 15:18:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 19 個(gè)回答

huozhe123

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 698.6
帖子: 182
在線: 44小時(shí)
蟲號(hào): 3556048
注冊(cè): 2014-11-24
專業(yè): 光學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

K=-3:0.1:3i;這個(gè)什么意思？從-3到0 ？還有，k1你在哪定義的？

贊一下

回復(fù)此樓

2樓2014-12-29 09:09:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wzp0098

木蟲 (著名寫手)

應(yīng)助: 2 (幼兒園)
金幣: 2524.4
散金: 39
紅花: 13
帖子: 1218
在線: 226.6小時(shí)
蟲號(hào): 853822
注冊(cè): 2009-09-22
性別: GG
專業(yè): 光學(xué)

我改了一下，程序還是運(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ù)此樓

3樓2014-12-29 12:45:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wzp0098

木蟲 (著名寫手)

應(yīng)助: 2 (幼兒園)
金幣: 2524.4
散金: 39
紅花: 13
帖子: 1218
在線: 226.6小時(shí)
蟲號(hào): 853822
注冊(cè): 2009-09-22
性別: GG
專業(yè): 光學(xué)

引用回帖:

2樓: Originally posted by huozhe123 at 2014-12-29 09:09:14
K=-3:0.1:3i;這個(gè)什么意思？從-3到0 ？還有，k1你在哪定義的？

我改了一下，程序還是運(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ù)此樓

4樓2014-12-29 12:46:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4	脫穎而出 2026-03-16	12/600	2026-03-19 16:17 by 脫穎而出
[考研] 324分 085600材料化工求調(diào)劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無(wú)際的草原
[考研] 304求調(diào)劑 +6	司空. 2026-03-18	6/300	2026-03-18 23:03 by 星空星月
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無(wú)懈可擊111
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無(wú)際的草原
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo