Radial Distribution function 2D

22 次查看（过去 30 天）

Alakesh Upadhyaya 2021-8-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/892002-radial-distribution-function-2d

评论： Alakesh Upadhyaya 2021-8-4

I have 792 particle cordiantes (x and y) inside a 2D sphere i.e a circle. I want to find the radial distribution function by taking each of the 792 particles as a tagged center and the average them all in the end. But in doing so I need to bin the particles within certain radius. I know the algorithm but I do not know the syntex. This is my code so far:

data=load('data');
box=100;                % size of box
radius=6;           % radius of the circle inside the box where the particles are located
npart=500;        % number of particles
xcord=100.*data(:,3);            % x cordinates
ycord=100.*data(:,4);            % ycordinates
xc=100.*0.66297;                   % x-center             since the center is not exactly in 0,0 as it is not perfect circle.
yc=100.*0.32219;                    % y-center
dr=0.3;                            % width of the bin
nbins=12./drh;                    % number of bins =40 in my case
%% Binning the particle
number=zeros(length(xs),1);
for i=1:length(xs)
          num=sqrt((xc-xs(i)).^2+((yc-ys(i)).^2));                    % This gives the interpartice distance
           number(i)=num';                                                
          idx=number./dr;                                                        %Indexing the distance.
            
        
    end
    
end