±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1291  |  »Ø¸´: 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
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

dbb627

ÈÙÓþ°æÖ÷ (ÖøÃûдÊÖ)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ...
wzp0098: ½ð±Ò+70, ¡ïÓаïÖú, лл£¡ 2015-01-02 09:19:29
CODE:
function slice_4D
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);
slice(X,Y,Z,P,[-2.5 0 2.5],0,[-2 2])

ÓÃmatlab »­Èýάͼ£¬Çë¸ßÊÖÖ¸½Ì.
untitled.jpg

» ±¾ÌûÒÑ»ñµÃµÄºì»¨£¨×îÐÂ10¶ä£©

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.
8Â¥2014-12-31 09:37:08
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 19 ¸ö»Ø´ð

huozhe123

½ð³æ (СÓÐÃûÆø)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
K=-3:0.1:3i;Õâ¸öʲôÒâ˼£¿´Ó-3µ½0 £¿»¹ÓУ¬k1ÄãÔÚÄ͍ÒåµÄ£¿
2Â¥2014-12-29 09:09:14
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wzp0098

ľ³æ (ÖøÃûдÊÖ)

ÎÒ¸ÄÁËһϣ¬³ÌÐò»¹ÊÇÔËÐв»ÆðÀ´£¬Çë´óϺָ½Ì£¡
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
3Â¥2014-12-29 12:45:43
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wzp0098

ľ³æ (ÖøÃûдÊÖ)

ÒýÓûØÌû:
2Â¥: Originally posted by huozhe123 at 2014-12-29 09:09:14
K=-3:0.1:3i;Õâ¸öʲôÒâ˼£¿´Ó-3µ½0 £¿»¹ÓУ¬k1ÄãÔÚÄ͍ÒåµÄ£¿

ÎÒ¸ÄÁËһϣ¬³ÌÐò»¹ÊÇÔËÐв»ÆðÀ´£¬Çë´óϺָ½Ì£¡
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
4Â¥2014-12-29 12:46:15
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Äܶ¯µ÷¼Á326ר˶ +4 wan112233 2026-04-04 4/200 2026-04-04 22:47 by yu221
[¿¼ÑÐ] ¿É¿çרҵµ÷¼Á +3 ÖÜµÄµÃµØ 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[¿¼ÑÐ] »·¾³¿ÆÑ§Ó빤³Ì334·ÖÇóµ÷¼Á +9 ÍõÒ»Ò»ÒÀÒÀ 2026-03-30 12/600 2026-04-04 20:55 by dongzh2009
[¿¼ÑÐ] 085701Çóµ÷¼Á +7 ¹¨ÓíÃú 2026-04-04 8/400 2026-04-04 13:49 by ССÊ÷2024
[¿¼ÑÐ] 280Çóµ÷¼Á +21 ¹¾ààÏþÏþ 2026-04-02 22/1100 2026-04-04 11:12 by Öí»á·É
[¿¼ÑÐ] 350Ò»Ö¾Ô¸±±¾©º½¿Õº½Ìì´óѧ08500²ÄÁÏ¿ÆÑ§Ó빤³ÌÇóµ÷¼Á +5 kjnasfss 2026-04-03 5/250 2026-04-03 22:29 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] »¯¹¤Çóµ÷¼Á +11 ÀóÏã֥ʿҬÄÌ 2026-04-03 11/550 2026-04-03 22:06 by à£à£à£0119
[¿¼ÑÐ] 295Çóµ÷¼Á +6 FZAC123 2026-04-03 6/300 2026-04-03 21:01 by zhq0425
[¿¼ÑÐ] 289-Çóµ÷¼Á +4 ÕâÀïÊÇ_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
[¿¼ÑÐ] µ÷¼ÁÇóÖú +8 Ïë»»ÊÖ»ú²»Ïë½âÊ 2026-04-02 11/550 2026-04-03 14:16 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] Çóµ÷¼Á +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[¿¼ÑÐ] 283Çóµ÷¼Á +3 jiouuu 2026-04-03 4/200 2026-04-03 13:28 by jiouuu
[¿¼ÑÐ] 326Çóµ÷¼Á +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[¿¼ÑÐ] ÍÁľ304Çóµ÷¼Á +4 ÍÃͻͻͻ£¬ 2026-04-02 5/250 2026-04-02 21:16 by ÍÃͻͻͻ£¬
[¿¼ÑÐ] 327Çóµ÷¼Á +9 ¼³Òàê» 2026-03-29 9/450 2026-04-02 11:44 by ÔËÆøyunqi
[¿¼ÑÐ] 354Çóµ÷¼Á +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[¿¼ÑÐ] 070300»¯Ñ§279Çóµ÷¼Á +15 ¹þ¹þ¹þ^_^ 2026-03-31 17/850 2026-04-01 21:37 by ¸øÄãÄã×¢ÒâÐÝÏ¢
[¿¼ÑÐ] 085601²ÄÁϹ¤³ÌÕÒµ÷¼Á +20 oatmealR 2026-03-29 21/1050 2026-04-01 21:00 by lijunpoly
[¿¼ÑÐ] 311Çóµ÷¼Á +10 ÀîÜÆÐÂ1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[˶²©¼ÒÔ°] ¿¼Ñе÷¼Á +5 ÂæÍÕÄÐÈË 2026-04-01 5/250 2026-04-01 14:28 by syjjj0321
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û