如何将函数绘制到网格上

Question

我是 MATLAB 新用户，我正在尝试绘制一个函数：

function [ uncertainty ] = uncertain(s1, s2, p)
%UNCERTAIN calculates the measurement uncertainty of a triangulation
% provide two coordinates of known stations and a target coordinate 
% of another point, then you get the uncertainty 
 [theta1, dist1] = cart2pol(p(1)-s1(1), p(2)-s1(2));
 [theta2, dist2] = cart2pol(p(1)-s1(1), p(2)-s2(2));
 theta=abs(pi-theta2-theta1);
 uncertainty = dist1*dist2/abs(sin(theta));
end

调用方式：

uncertain([0 0],[8 0],[4 4])

我得到一个结果。 但我想要一个整个表面并调用：

x=-2:.1:10;
y=-2:.1:10;
z = uncertain([0 0],[8 0],[x y]);
mesh(x,y,z)

我收到错误：“Z 必须是矩阵，而不是标量或向量。”

如何修改我的代码以便我的函数绘制表面？

提前致谢。 拉尔夫.

Answer 1

首先，我认为您的函数有一个错误：您的 [theta2, dist2] = cart2pol(p(1)-s1(1), p(2)-s2(2)); 应该有 th第一个 s1 是 s2。

接下来，要获得矢量输入的矢量答案，您必须将 p(i) （选择 p 的第 i 个元素）更改为 p (i,:)，它将选择 p 的第 i 行。

之后，将乘法 (*) 更改为按元素乘法 (.*)。

总之：

function [ uncertainty ] = uncertain(s1, s2, p)
%UNCERTAIN calculates the measurement uncertainty of a triangulation
% provide two coordinates of known stations and a target coordinate 
% of another point, then you get the uncertainty
% target coordinates p are 2xn
% output uncertainty is 1xn
 [theta1, dist1] = cart2pol(p(1,:)-s1(1), p(2,:)-s1(2));
 [theta2, dist2] = cart2pol(p(1,:)-s2(1), p(2,:)-s2(2));
 theta=abs(pi-theta2-theta1);
 uncertainty = dist1.*dist2./abs(sin(theta));
end

唯一的变化是 p(i) -> p(i,:) 和 *->.* 和 /->./.

要获得表面，您可以使用 meshgrid 获取网格中的所有 (x,y) 坐标集，并将它们展平为 2xn 矩阵对于不确定，然后将它们展开回网格以进行绘制。示例：

x=-2:.1:10;  % 121 elements
y=-2:.1:10;  % 121 elements
[xs,ys]=meshgrid(x,y); % xs and ys are each 121 x 121
zs = uncertain([0 0],[8 0],[xs(:) ys(:)]'); %get zs, being 1x(121*121) ie 1x14641
% Reshape zs to be 121x121 in order to plot with mesh
mesh(xs,ys,reshape(zs,size(xs)))

注意：您会得到很多非常大的数字，因为当 theta 为 0 或 pi （或非常接近）时，因为那时您'再除以（几乎）0。