matlab中如何将连续函数离散化

Xa = -1 ; % X的下限
Xb = 1 ; % X的上限
Ya = -1 ; % Y的下限
Yb = 1 ; % Y的上限
Xn = 30; % 在X的下限与上限中取的点数
Yn = 30; % 在Y的下限与上限中取的点数
% 输入参数
center_x = 0;
center_y = 0;
sigma = 0.3;
A = 100;
Xtspan = linspace(Xa,Xb,Xn);
Ytspan = linspace(Ya,Yb,Yn);
[X,Y] = meshgrid(Xtspan,Ytspan);
% zi以矩阵形式表示
Z = A*exp(-0.5*((X-center_x).^2/(sigma^2) + (Y-center_y).^2/(sigma^2)))
surf(X,Y,Z)
title('高斯钟')
xlabel( 'x' );
ylabel( 'y' );
zlabel( 'z' );