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

6Â¥: Originally posted by iloveyou1991 at 2014-12-30 15:48:42
Äã»­µÄÊÇÈýάµÄÂð£¿
ÄãÊDz»ÊÇÏë»­P=P(x,y,z)£¬ÕâÊÇËÄάµÄ

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])

8Â¥: Originally posted by dbb627 at 2014-12-31 09:37:08
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);
               K ...

8Â¥: Originally posted by dbb627 at 2014-12-31 09:37:08
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);
               K ...