u=[1620 1639 1262 1309 1445 1627 1860 374 301 739 706 490 99 323];
v=[232 323 228 323 322 616 594 268 364 243 337 351 643 646];
x=[7553 6353 7553 6353 6353 3953 3953 7643 6443 7643 6443 6443 4043 4043];
y=[30931 30451 28531 28531 29251 29251 30451 23059 23059 25459 25459 24259 23059 24259];
p1=y/[v-v+1;v;v.^2;v.^3]%p1=[a1 b1 c1 d1]
p2=x/[u-u+1;u;u.*y]%p2=[a2 b2 c2]
%%验证
V=linspace(min(v),max(v));
U=linspace(min(u),max(u));
plot(v,y,'o');
hold on
Y=p1(1)+p1(2)*V+p1(3)*V.^2+p1(4)*V.^3;
plot(V,p1(1)+p1(2)*V+p1(3)*V.^2+p1(4)*V.^3)
plot(u,x,'o')
plot(U,p2(1)+p2(2)*U+p2(3)*U.*Y)