如何在matlab中为自定义函数f(x,y)制作等高线图,其中x和y不能取矢量值？

2019年8月3日 297次阅读

我们知道在Matlab中为函数Z(x,y)制作等高线图的常用方法是

[X,Y] = meshgrid(-2：.2：2,-2：.2：3);

Z = X. * exp(-X.^ 2-Y.^ 2); (例如)

轮廓(X,Y,Z);

但是,这种方式不适用于以下函数f(x,y)：
假设h_ {ij}(x,y)是大(例如,100×100)矩阵,其中每个分量是x和y的(自定义)函数.我们定义另一个函数

F(X,Y)= DET(H_ {IJ}(X,Y))

并且想要绘制函数f(x,y)的等高线图.

f = det(h)中的行列式要求矩阵h的每个分量都是数字.因此,只有当x和y是数字而不是向量时,才能通过Matlab计算f(x,y).如果我们使用[X,Y] = meshgrid(…),则意味着矩阵h的每个分量都是向量,并且不能计算f(X,Y).

有没有办法为上面的函数f(x,y)制作等高线图,其中x和y不能取矢量值？

最佳答案 假设h被预定义为函数矩阵,每个函数采用两个标量参数并输出矩阵(或对det函数的任何有效输入),而下标i和j指的是X和Y中的索引.该函数的参数,类似下面的代码应该工作(X和Y应该与h相同)：

applyh = @(fn, x, y) fn(x, y);
[I, J] = meshgrid(1:m, 1:n);
Z = arrayfun(@(i, j) det(applyh(h(i, j), X(i), Y(j))), I, J);

我认为你误解了meshgrid的作用 – meshgrid的输出可以很容易地输入到上面的函数中.它们不是每个元素中的向量(只是一个二维矩阵).然后你可以照常绘制Z.

点赞