代码备份：按模式划分椭圆轨迹渐变matlab

     x=sqrt(a^*(-y^/b^));

     r=ceil(+*y);

     x=ceil(x);

     if r==

         I(m/,m/+x)=+*y;

     elseif r==

         I(m/::m/+,m/+x::m/++x)=+*y;

     elseif r==

         I(m/-::m/+,m/+x)=+*y;

         I(m/,m/-+x::m/++x)=+*y;

     elseif r==

         I(m/-::m/+,m/-+x::m/++x)=+*y;

     elseif r==

         I(m/-::m/+,m/-+x::m/++x)=+*y;

         I(m/-::m/+,m/+x)=+*y;

         I(m/,m/-+x::m/++x)=+*y;

     end