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

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

11楼: Originally posted by dbb627 at 2015-01-03 23:33:10
你想看的是否是P关于W;K;V的等值面？...

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)

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);
              ...

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);
              ...

15楼: Originally posted by wzp0098 at 2015-01-04 16:46:47
问一下，这P(n, m, v)的四维图吗...

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

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

18楼: Originally posted by dbb627 at 2015-01-04 18:29:15
这P(n, m, v)的四维图是第一个帖子里面画的图，后面这个是等值面
他们的关系就像三维曲面和等值线之间的关系...