椭圆拟合的matlab程序

1、clear; clc; close all;% 设出圆锥曲线方程F=(p,x)p(1)*x(:,1)A2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2)A2+p(4)*x(:,1)+p(5)*x(:,2)+p(6);% 离散数据点x1=276227542745273527232710269526792661264326232602258025572532250724812454242523962367233623052273224022072173213821032068203219951959192118841846180717691730169116521613157315

2、341494145414151375133612961256121711781139110010611023985947910873836800764729694660626593561529499468439410383356330305281259237216197179162146132119107978982767270697073778492102113127143160179201224249276304335367401438475514555598642687734782832882933986103910931148120312591315137114271484154015

3、9616521707176118151868192119722022207221202166221222552297233923782415245024842516254725752602262626492670268927062722 2735 2746 2756 2764 2770 2774 2776 2777 2776 2773 2769 ; y1=1792 1786 1781 1775 1769 1764 1758 1753 1747 1742 1737 1732 1727 1722 1717 1713 1708 1704 1700 1696 1692 1688 1684 1681 1

4、677 1674 1671 1668 16651663166016581655165316511650164816461645164416431642164116401640164016391639164016401640164116411642164316441646164716491651165216541657165916611664166716701673167616791683168616901694169817021706171117151720172517301735174017451750175617611767177317781784179017961802180818141

5、82118271833183918451852185818641870187718831889189519001906191219171923192819331938194219471951195519591963196619701972197519771979198119831984198519851986198619851985198419831981198019781975197319701967196319601956195219481943193819341929192319181913190719011896189018841878187218661859 1853184718411835182818221816181018041798x=x1,y1 p0=1 1 1 1 1 1;warning off% 拟合系数，最小二乘方法 p=nlinfit(x,zeros(size(x,1),1),F,p0);plot(x(:,1),x(:,2),'ro');hold on;xmin=min(x(:,1); xmax=max(x(:,1);ymin=min(x(:,2); ymax=max(x(:,2);% 作图ezplot(x,y)F(p,x,y),-1+x