Simulating Plane Wave Expansion in Matlab

[Vangjush Komini]

Dear Sir Madam,

I am a new member of this forum and I was trying to understand Plane Wave expansion simulation in Matlab. So far I cannot understand how to simulate Gz vector in K-space consistently to the numbers of plane waves.
Hereby you can have the MATLAB code for Gx and Gy.nx_arr=repmat(n1,1,2*nmax+1);
ny_arr=round(n2/(2*nmax+1));
nz_arr=1:3;

[nx,nxp]=meshgrid(nx_arr,nx_arr);
[ny,nyp]=meshgrid(ny_arr,ny_arr);

nx_arr=repmat(n1,1,2*nmax+1);
ny_arr=round(n2/(2*nmax+1));
nz_arr=1:3;

[nx,nxp]=meshgrid(nx_arr,nx_arr);
[ny,nyp]=meshgrid(ny_arr,ny_arr);

Thank you a lot for your consideration

Regards

 

[Achy47]

,Thank you for your question and welcome to the forum! I understand that you are trying to simulate Gz vector in K-space consistently with the number of plane waves using MATLAB. I would be happy to assist you with this.

First, let me explain what the code you provided is doing. The lines where you define nx_arr and ny_arr are creating arrays of numbers to be used as indices for your nx and ny matrices. These arrays are then used to create meshgrid matrices, which are used for plotting purposes.

In order to simulate Gz vector in K-space, you will need to define the values for n3 or nz_arr in your code. This will depend on the specific simulation you are trying to run and the number of plane waves you want to use. Once you have defined nz_arr, you can use it in a similar way as the nx_arr and ny_arr to create the meshgrid matrix for nz.

I hope this helps and if you have any further questions, please don't hesitate to ask. Best of luck with your simulation!