D2_Space_离子束_相空间_二维相空间_相空间分布_源码

%参照 %-----二维相空间 alpha = -4.52; % beta = 90.3; % vareps = 0.2; %归一化均方根发射度 epslon = 0.2*4; %发射度 gamma = (alpha^2+1)/beta; % x_max = sqrt(beta*epslon); xp_max = sqrt(gamma*epslon); x_valu = -8.5:0.1:8.5; xp_valu = -0.4:0.1:0.4; xp1 = -alpha*x_valu/beta+sqrt(epslon*beta-x_valu.^2)/beta; xp2 = -alpha*x_valu/beta-sqrt(epslon*beta-x_valu.^2)/beta; figure; plot(x_valu,xp1,'b'); hold on; plot(x_valu,xp2,'b'); grid; xlabel('x'); ylabel('px');