2021-07-26

例3.2-3 根据椭圆方程

$$\begin{array}{c}\left(\begin{array}{cc}x& y\end{array}\right)\left(\begin{array}{cc}3& 1\\ 1& 4\end{array}\right)\left(\begin{array}{c}x\\ y\end{array}\right)\end{array}=5$$

绘制椭圆曲线，并修饰图形。

>> P=[3 1;1 4];
>> r=5;
>> [V,D]=eig(P);   %求特征值，将椭圆化为标准方程
>> a=sqrt(r/D(1));   %椭圆长半轴
>> b=sqrt(r/D(4));  %椭圆短半轴
>> t=linspace(0,2*pi,60);   %等间隔产生一个从0到2pi的包含60个元素的向量
>> xy=V*[a*cos(t);b*sin(t)]; %根据椭圆的极坐标方程计算椭圆上点的坐标
>> plot(xy(1,:),xy(2,:),'k','linewidth',3);    %绘制椭圆曲线，线宽为3，颜色为黑色
>%在当前图形窗口加入带箭头的文本标注框
>> h=annotation('textarrow',[0.606 0.65],[0.55 0.65]);
>%设置文本标注框中显示的字符串，并设字号为15
>> set(h,'string','3x^2+2*x*y+4*y^2=5','fontsize',15);
>%为图形加标题，设字号为18，字体加粗
>> h=title('这是一个椭圆曲线','fontsize',18,'fontweight','bold');
>> set(h,'position',[-0.00345622 1.35769 1.00011]);%设置标题的位置
>> axis([-1.5 1.5 -1.2 1.7]);      %设置坐标轴的显示范围
>> xlabel('X')            %为X轴加标签
>> ylabel('Y')          %为Y轴加标签
>> 

运行结果图：