图像纵向扭曲算法

竖向扭曲，对每个网格顶点赋值时，算法在u方向均匀取值，v方向采用函数变化，变形后网格以及变形函数如下
u(i,j)=i/499;
v(i,j)=v(i,j-1)+tension*s(i,j)*sin(i/499*f)。
v方向增量变化为tension*s(i,j)*sin(i/499*f)，其中正弦函数频率为其中正弦函数频率为f=1/tension*PI，其中tension控制变形程度。
振幅为可变函数，s(i,j)=(1-fabs(u(i,j-1)-0.5)/0.5) * (1-fabs(v(i,j-1)-0.5)/0.5).
求大神指教啊，用MFC怎么做？感觉老师给的公式有问题啊！！