bezier曲面程序matlab

⑴ matlab画Bezier曲线求助

原理可参考网络的【贝塞尔曲线】，当然，英语过得去的话建议参考维基网络的【Bézier curve】（没办法，中文版被墙了）。

参考代码：

%控制点
p=[111;123;234;346;545;589;689];

%n为控制点的数量，t为0~1的参变量，b为平滑后的结果
n=size(p,1);
t=linspace(0,1)';
b=0;
fork=0:n-1
tmp=nchoosek(n-1,k)*t.^k.*(1-t).^(n-1-k);
b=b+tmp*p(k+1,:);
end
plot3(p(:,1),p(:,2),p(:,3),'.:',b(:,1),b(:,2),b(:,3))
gridon

⑵ 急找绘制bezier曲面代码（matlab）

/*
产生三次方贝塞尔曲线的程序码
*/

typedef struct
{
float x;
float y;
}
Point2D;

/*
cp 在此是四个元素的阵列:
cp[0] 为起始点，或上图中的 P0
cp[1] 为第一个控制点，或上图中的 P1
cp[2] 为第二个控制点，或上图中的 P2
cp[3] 为结束点，或上图中的 P3
t 为参数值，0 <= t <= 1
*/

Point2D PointOnCubicBezier( Point2D* cp, float t )
{
float ax, bx, cx;
float ay, by, cy;
float tSquared, tCubed;
Point2D result;

/* 计算多项式系数 */

cx = 3.0 * (cp[1].x - cp[0].x);
bx = 3.0 * (cp[2].x - cp[1].x) - cx;
ax = cp[3].x - cp[0].x - cx - bx;

cy = 3.0 * (cp[1].y - cp[0].y);
by = 3.0 * (cp[2].y - cp[1].y) - cy;
ay = cp[3].y - cp[0].y - cy - by;

/* 计算位于参数值 t 的曲线点 */

tSquared = t * t;
tCubed = tSquared * t;

result.x = (ax * tCubed) + (bx * tSquared) + (cx * t) + cp[0].x;
result.y = (ay * tCubed) + (by * tSquared) + (cy * t) + cp[0].y;

return result;
}

/*
ComputeBezier 以控制点 cp 所产生的曲线点，填入 Point2D 结构的阵列。
呼叫者必须分配足够的内存以供输出结果，其为 <sizeof(Point2D) numberOfPoints>
*/

void ComputeBezier( Point2D* cp, int numberOfPoints, Point2D* curve )
{
float dt;
int i;

dt = 1.0 / ( numberOfPoints - 1 );

for( i = 0; i < numberOfPoints; i++)
curve[i] = PointOnCubicBezier( cp, i*dt );
}