Question

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

比如高斯钟方程zi = A*exp( -0.5* ((X-center_x).^2/(sigma^2) + (Y-center_y).^2/(sigma^2))) ;
想把zi以矩阵形式表示，X,Y表示行列，zi表示深度。即zi（X,Y）=？

Answer 1

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' );