dbb627

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

引用回帖:

10樓: Originally posted by wzp0098 at 2015-01-03 22:40:00
我想看到附件一樣的圖，請指教. 是否能用griddata3？

@RN1%DO4A1R%H5DNB`$PFQI.png
...

你想看的是否是P關(guān)于W;K;V的等值面？

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

11樓2015-01-03 23:33:10

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

wzp0098

木蟲 (著名寫手)

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

引用回帖:

11樓: Originally posted by dbb627 at 2015-01-03 23:33:10
你想看的是否是P關(guān)于W;K;V的等值面？...

是的

回復(fù)此樓

12樓2015-01-04 09:55:24

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
wzp0098: 金幣+60, ★★★很有幫助 2015-01-04 16:47:10

CODE:

function slice_4D

W=-3:0.05:3;K=-3:0.05:3;V=-3:0.05: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);

figure;slice(X,Y,Z,P,[-2.5 0 2.5],0,[-2 2])

result = HistRate(P(:));

[Pm,I]=max(result(:,2));

figure;

result(I,1)

Pv=isosurface(X,Y,Z,P,result(I,1));

p=patch(Pv)

untitled1.jpg

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

13樓2015-01-04 15:36:05

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
wzp0098: 金幣+60, ★★★很有幫助, 多謝了 2015-01-04 22:14:08

引用回帖:

13樓: Originally posted by dbb627 at 2015-01-04 15:36:05
function slice_4D
W=-3:0.05:3;K=-3:0.05:3;V=-3:0.05:3;
P = [];
for n=1:length(W);
      Ky=W(n);
   for m=1:length(K);
      Kx=K(m);
         for v=1:length(V);
            ...

HistRate可換成tabulate 函數(shù)主要找出P中出現(xiàn)最多的那個值，
Pv=isosurface(X,Y,Z,P,result(I,1)); result(I,1)可換成其他可以取到的P值

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

14樓2015-01-04 15:42:21

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

wzp0098

木蟲 (著名寫手)

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

引用回帖:

問一下，這P(n, m, v)的四維圖嗎

贊一下

回復(fù)此樓

15樓2015-01-04 16:46:47

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

引用回帖:

15樓: Originally posted by wzp0098 at 2015-01-04 16:46:47
問一下，這P(n, m, v)的四維圖嗎...

Pv=isosurface(X,Y,Z,P,result(I,1));是畫的三維等值面啊

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

16樓2015-01-04 16:50:40

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

wzp0098

木蟲 (著名寫手)

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

引用回帖:

16樓: Originally posted by dbb627 at 2015-01-04 16:50:40
Pv=isosurface(X,Y,Z,P,result(I,1));是畫的三維等值面啊...

也就P(n, m, v)有關(guān)于X,Y,Z的圖形？

贊一下

回復(fù)此樓

17樓2015-01-04 16:59:07

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

引用回帖:

17樓: Originally posted by wzp0098 at 2015-01-04 16:59:07
也就P(n, m, v)有關(guān)于X,Y,Z的圖形？...

這P(n, m, v)的四維圖是第一個帖子里面畫的圖，后面這個是等值面
他們的關(guān)系就像三維曲面和等值線之間的關(guān)系

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

wzp0098

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

18樓2015-01-04 18:29:15

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

wzp0098

木蟲 (著名寫手)

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

送紅花一朵

引用回帖:

18樓: Originally posted by dbb627 at 2015-01-04 18:29:15
這P(n, m, v)的四維圖是第一個帖子里面畫的圖，后面這個是等值面
他們的關(guān)系就像三維曲面和等值線之間的關(guān)系...

請教高人，我想畫出P(n, m, v)在整個四維空間的分布，將它所有的值都在圖形里面體現(xiàn)出來，怎么修改？
function slice_4D
W=-3:0.05:3;K=-3:0.05:3;V=-3:0.05: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)); %
         end
   end
end
[X,Y,Z]=meshgrid(W,K,V);
figure;slice(X,Y,Z,P,[-2.5 0 2.5],0,[-2 2])
result = HistRate(P(

);
[Pm,I]=max(result(:,2));
figure;
result(I,1)
Pv=isosurface(X,Y,Z,P,result(I,1));
p=patch(Pv)

贊一下

回復(fù)此樓

19樓2015-01-11 21:39:56

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

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

[考研] 271材料工程求調(diào)劑 +6	.6lL 2026-03-18	6/300	2026-03-19 15:41 by haoshis
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-19	3/150	2026-03-19 13:36 by houyaoxu
[考研] 281求調(diào)劑（0805） +9	煙汐憶海 2026-03-16	19/950	2026-03-19 11:42 by laoshidan
[考研] 346求調(diào)劑[0856] +3	WayneLim327 2026-03-16	6/300	2026-03-19 11:21 by WayneLim327
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎覬?/a>
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 288求調(diào)劑 +4	奇點0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky