代码备份：按模式划分椭圆轨迹渐变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