1、首先将原式写为参数方程,以便作图:
x = 1/a + 1/a * cos(t)
y = -1 + 1/a * sin(t)
其中 0≤t≤2π
2、可以通过循环a的方法作图,这里提供使用更为简洁的arrayfun()函数的方法,代码参考:
>> t = 0:0.01*pi:2*pi;
>> A = -5:0.5:5; % a的范围,可以自行修改
>> hold on;axis equal % 图形参数
>> arrayfun(@(a)plot(1./a+1./a*cos(t),-1+1./a*sin(t)), A)
3、结果如下: