當(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人參與

程序如下，可是沒有辦法運(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ù)此樓

» 猜你喜歡

329求調(diào)劑已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有11人回復(fù)
一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一）已經(jīng)有3人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有14人回復(fù)
320求調(diào)劑0856 已經(jīng)有3人回復(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ù)

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

請(qǐng)問三維金相立體圖如何畫出來已經(jīng)有17人回復(fù)
如何用matlab導(dǎo)入導(dǎo)入三維數(shù)據(jù)~ 已經(jīng)有5人回復(fù)
請(qǐng)問采用AFM掃描的文件如何用matlab畫出三維圖？已經(jīng)有7人回復(fù)
求matlab功率譜的畫法，我這樣對(duì)么？已經(jīng)有6人回復(fù)
【已完結(jié)】matlab仿真圖形（交流問題即有獎(jiǎng)，呵呵）已經(jīng)有6人回復(fù)
菜鳥求助用matlab如何畫三維圖已經(jīng)有3人回復(fù)
急用�。。�！求畫三維拓?fù)鋱D 已經(jīng)有10人回復(fù)
求助已知三維坐標(biāo)位置，如何畫出空間圖像已經(jīng)有5人回復(fù)
在matlab中用mesh畫三維圖存在的問題已經(jīng)有11人回復(fù)
MATLAB中畫三維雙對(duì)數(shù)圖形的問題已經(jīng)有10人回復(fù)
MATLAB將三維圖和其中一個(gè)面的等高線畫出來的問題，謝謝了。已經(jīng)有7人回復(fù)
如何用matlab畫二維網(wǎng)格已經(jīng)有5人回復(fù)
求助MATLAB做三維圖問題已經(jīng)有4人回復(fù)
關(guān)于Matlab畫三維曲面圖的問題，求助高手，感激感激已經(jīng)有3人回復(fù)
請(qǐng)問如何用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畫出來的這個(gè)三維圖更漂亮？已經(jīng)有3人回復(fù)
彈性模量各向異性的三維圖。。。求高手解決！！已經(jīng)有10人回復(fù)
【求助】菜鳥求助 matlab 繪圖問題繪制光譜圖已經(jīng)有3人回復(fù)

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

已閱回復(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é)

引用回帖:

6樓: Originally posted by iloveyou1991 at 2014-12-30 15:48:42
你畫的是三維的嗎？
你是不是想畫P=P(x,y,z)，這是四維的

是的，是四維的，請(qǐng)指教！

贊一下

回復(fù)此樓

7樓2014-12-30 16:57:27

已閱回復(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)行不起來，請(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)行不起來，請(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ā)表

[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3	葵梓衛(wèi)隊(duì) 2026-03-18	5/250	2026-03-19 19:35 by 給你你注意休息
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +8	zlingli 2026-03-13	8/400	2026-03-19 16:32 by 輕松不少隨
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 332求調(diào)劑 +3	ydfyh 2026-03-17	3/150	2026-03-19 10:14 by 功夫瘋狂
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike