% 绘制3D的托卡马克位形图， (V 0.1 by Jiale Chan for Y. H. Huang)

% Dee Formula

% 特征参数

    rzero = 2.0;

    rmax = 0.75;

    eshape = 2.0;

    xshape = 0.4;

% 极向角与环向角

    theta = linspace(0, 2*pi, 50);

    phi = linspace(0, 1.3*pi, 50);

% 将数据矩阵化

    [phi_g, theta_g ] = meshgrid(phi, theta);

% 参考TOQ网页给的简单公式生成磁面数据：https://fusion.gat.com/THEORY/toq/geometry.html

    z = eshape*rmax*sin(theta_g);

    R = rzero + rmax*( cos(theta_g)-xshape*sin(theta_g).^2 ) ;

    x = R.*cos(phi_g);

    y = R.*sin(phi_g);

% 绘制3D磁面

    surf(x,y,z); shading interp; %或者 surf(x,y,z, 'EdgeColor','none');

    hold on;

%随便画的一条线

    R_line = rzero + rmax*( cos(theta)-xshape*sin(theta).^2 ) ;

    x_line = R_line.*cos(phi*1.5);

    y_line = R_line.*sin(phi*1.5);

    z_line = eshape*rmax*sin(theta);

    plot3(x_line, y_line, z_line, 'k'); 

 

% 设定图的一些其他参数

    colorbar; %绘制色条

    xlabel('x');   ylabel('y');   zlabel('z');%坐标轴的命名

    axis equal; %将坐标轴比例设为“相等”

    view(2,40); %改变视角

% 效果图如下：