A..10 meshgrid()

>> nx=5; ny=3;
>> X=linspace(0,1,nx+1)

X =
         0    0.2000    0.4000    0.6000    0.8000    1.0000

>> Y=linspace(2,3,ny+1)

Y =
    2.0000    2.3333    2.6667    3.0000

>> [x,y]=meshgrid(X,Y)

x =
         0    0.2000    0.4000    0.6000    0.8000    1.0000
         0    0.2000    0.4000    0.6000    0.8000    1.0000
         0    0.2000    0.4000    0.6000    0.8000    1.0000
         0    0.2000    0.4000    0.6000    0.8000    1.0000

y =
    2.0000    2.0000    2.0000    2.0000    2.0000    2.0000
    2.3333    2.3333    2.3333    2.3333    2.3333    2.3333
    2.6667    2.6667    2.6667    2.6667    2.6667    2.6667
    3.0000    3.0000    3.0000    3.0000    3.0000    3.0000

>> whos
  Name      Size            Bytes  Class     Attributes

  X         1x6                48  double              
  Y         1x4                32  double              
  nx        1x1                 8  double              
  ny        1x1                 8  double              
  x         4x6               192  double              
  y         4x6               192  double

X $ =[x_1,x_2,\dots,x_{N_x+1}]\in\mathbb{R}^{N_x+1}$, Y $ =[y_1,y_2,\dots,y_{N_y+1}]\in\mathbb{R}^{N_y+1}$とするとき、 [x,y]=meshgrid(X,Y) とすると、 xy $ \mathbb{R}^{(N_y+1)\times(N_x+1)}$ の要素となる。 x(j,i)$ x_i$, y(j,i)$ y_j$ である。

MATLAB の2次元配列は column first で並んでいるので、 例えば x の成分は、メモリー内で、

x(1,1), x(2,1), x(3,1), x(4,1),  
x(1,2), x(2,2), x(3,2), x(4,2),  
$ \cdots$  
x(1,5), x(2,5), x(3,5), x(4,5),  
x(1,6), x(2,6), x(3,6), x(4,6)  
と並んでいる。
>> x(:)
ans =
         0
         0
         0
         0
    0.2000
    0.2000
    0.2000
    0.2000
    0.4000
    0.4000
    0.4000
    0.4000
    0.6000
    0.6000
    0.6000
    0.6000
    0.8000
    0.8000
    0.8000
    0.8000
    1.0000
    1.0000
    1.0000
    1.0000

>> y(:)'

ans =
  1 列から 9 列
    2.0000    2.3333    2.6667    3.0000    2.0000    2.3333    2.6667    3.0000    2.0000

  10 列から 18 列
    2.3333    2.6667    3.0000    2.0000    2.3333    2.6667    3.0000    2.0000    2.3333

  19 列から 24 列
    2.6667    3.0000    2.0000    2.3333    2.6667    3.0000

>>



桂田 祐史