row = size(I,1);
col = size(I,2);
for i = 1:20:row
for j = 1:20:col
x0 = j-1;
y0 = 1-i;
[x,y] = meshgrid(-5:5,-5:5);
ezplot('(x.*cos_phi+y.*sin_phi-x0).^2)./a.^2+((x.*sin_phi-
y.*cos_phi-y0).^2)./b.^2 -1');
hold on
end
end
此代码从3个矩阵中调用a
,b
,phi
个参数。