一种画椭圆的方法——MATLAB实现（附源码）

一种画椭圆的方法——MATLAB实现（附源码）

​ 在看 1B3B 视频的时候，对椭圆的实现比较感兴趣，然后就尝试用matlab上实现。代码很简单，直接看就能看懂，其中搞直线旋转搞了一会儿，写了一个小函数。

程序代码：

clear;clc;
r = 10;%半径
n=60;%等分点数
theta = 0:2*pi/n:2*pi; %角度[0,2*pi] 
x = r*cos(theta);
y = r*sin(theta);
xy=[x;y]';
f=[8,0];%给定点
xy2=zeros(2,n+1)';
for i=1:n+1
    xy2(i,:)=(xy(i,:)+f)/2;
end
plot(x,y,'-');
hold on;
plot(xy2(:,1),xy2(:,2));
hold on;

axis equal
for i=1:n+1
    transform(xy(i,:),f);
    hold on
end

transform函数代码：

%% 构建一个函数，用于旋转两点构成的直线绕其中点旋转90度
function transform(a,b) %a,b分别是第一个点和第二个点的坐标,行向量表示，即(x,y)

c=(a+b)/2; %c点为两点中点
xy=[a;b]; %两点坐标矩阵
% plot(xy(:,1),xy(:,2)); %显示原先直线
hold on;
a=a-c;
b=b-c;
xy=[a;b];
T=[cos(pi/2),sin(pi/2);
    -sin(pi/2),cos(pi/2)];
xy2=xy*T;
xy2(1,:)=xy2(1,:)+c;
xy2(2,:)=xy2(2,:)+c;
plot(xy2(:,1),xy2(:,2));

具体效果：

更改n能够获得更加圆滑的椭圆

n=200:

代码写的比较粗糙，感兴趣的伙伴可以改进以下~