function?tuoyuan=creat_ellipse（ellipsen）
p?=?zeros（n）;

xax?=??（?（0:n-1）-（n-1）/2?）?/?（（n-1）/2）;?
xg?=?repmat（xax?n?1）;???%?x?coordinates?the?y?coordinates?are?rot90（xg）

for?k?=?1:size（ellipse1）????
???asq?=?ellipse（k2）^2;???????%?a^2
???bsq?=?ellipse（k3）^2;???????%?b^2
???phi?=?ellipse（k6）*pi/180;??%?rotation?angle?in?radians
???x0?=?ellipse（k4）;??????????%?x?offset
???y0?=?ellipse（k5）;??????????%?y?offset
???A?=?ellipse（k1）;???????????%?Amplitude?change?for?this?ellipse
???x=xg-x0;????????????????????%?Center?the?ellipse
???y=rot90（xg）-y0;??
???cosp?=?cos（phi）;?
???sinp?=?sin（phi）;
???idx=find（（（x.*cosp?+?y.*sinp）.^2）./asq?+?（（y.*cosp?-?x.*sinp）.^2）./bsq?<=?1）;?
???p（idx）?=?p（idx）?+?A;
???tuoyuan=p;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3722??2012-07-28?13:35??shanshu2.m

?????文件????????745??2012-07-24?17:35??creat_ellipse.m

?????文件????????506??2012-07-25?23:03??gtheta.m

-----------?---------??----------?-----??----

?????????????????4973????????????????????3