Question

我不是MATLAB专业人员，因此，我需要一些帮助来制作迭代定义函数的 ： R ^ 2-0 - ＆gt; ;在 [ - 1,1] 和

中 x，y 值下面定义的 R （伪代码）
每个 I = 1,2,3 和 A =（0.5,0.5）

function f(v in R^2-0) { a=b=0; (a,b in R) for (i=0; i<I; i=i+1) { v = |v| / ||v||^2 - A; a = a + | ||v||-b |; b = ||v||; } return a; }

（| v |表示分量矢量绝对值）

（如果你想，你可以在这里查看我在数学交换问题上由函数生成的分形：

https://math.stackexchange.com/questions/1457733/a-question-about-a-fractal-like-iteratively-defined-function ）

这样做的MATLAB代码将不胜感激。

非常感谢。

Answer 1

保存这个主程序：

clear
clc
close all

% I = 1;
% A = [ 0.5 0.5 ];

I = 10;
A = [ 0.5 0.5 0.5 ];

xmin = -1;
xmax = 1;

ymin = -1;
ymax = 1;

nx = 101;
ny = 101;

dx = (xmax - xmin) / (nx - 1);
dy = (ymax - ymin) / (ny - 1);

x = xmin: dx: xmax;
y = ymin: dy: ymax;

for ix = 1: nx
    for iy = 1: ny
        if (length(A) == 2)
            z(iy, ix) = f([x(ix) y(iy)], A, I);
        elseif (length(A) == 3)
            z(iy, ix) = f([x(ix) y(iy) 0], A, I);
        end
    end
end

pcolor(x, y, z)
shading interp

然后将此功能保存在与主程序相同的目录中f.m：

function result = f(v, A, I)

a = 0;
b = 0;

for i = 1: I
    v = abs(v) / dot(v, v) - A;
    a = a + abs(norm(v) - b);
    b = norm(v);
end
result = a;

end

在MATLAB中绘制迭代定义的函数

1 个答案: