%%?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%???FILENAME????????ellipseLS.m
%???FUNCTIPN????????Least-squares?fit?of?ellipse?to?2D?points
%???DATE????????????2015-12-09
%???AUTHOR??????????liuzongming
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
clc;
clear?all;
close?all;

format?long?g
%%??生成?帶噪聲的橢圓
%?參數初始化
g_NumOfPoints?=?500;???%?點數
g_NormDistrVar?=?3;????%?標準偏差
a=100;b=200;?????????????%長軸短軸
angle=?50;??????????????%傾斜角
%%?橢圓生成
beta?=?angle?*?（pi?/?180）;
alpha?=?linspace（0?360?g_NumOfPoints）?.*?（pi?/?180）;?
x?=?（a?*?cos（alpha）?*?cos（beta）-?b?*?sin（alpha）?*?sin（beta）?）+wgn（1length（alpha）g_NormDistrVar^2‘linear‘）;????
y?=?（a?*?cos（alpha）?*?sin（beta）+?b?*?sin（alpha）?*?cos（beta）?）+wgn（1length（alpha）g_NormDistrVar^2‘linear‘）;
plot（x?y?‘.‘?‘Tag‘?‘DATA‘）;