自己编的程序是

function h1=SCYTNH(x,h0,H0,x1,H1)

N=length(x); %起始拟合点(x 或y)的个数

e=H0-h0; %求高程异常

M=length(e); %起始高程异常点的个数

NN=length(x1); %被拟合点的个数

if(N~=M) %必须保持一一对应关系

disp('警告：x 与e 必须有相同的维数');

return

end

e1=interp1(x,e,x1,'spline'); %求解未联测的水准点的高程异常

h1=vpa((H1-e1),6); %求解未联测的水准点的高程

但是程序出来以后总是有错误。。希望各位大侠帮忙。。不胜感激、、、被拟合点总共有21个。数据如下表

点号        X坐标        Y坐标        GPS大地高        水准正常高        高程异常

GP01        558490.1980         4670319.9740         79.9960         69.1740         10.8220

GP03        557316.3390         4670413.8920         86.0540         75.2060         10.8480

GP05        554876.1630         4672885.3400         64.4910         53.6630         10.8280

GP08        552548.3900         4674983.1190         63.0690         52.2600         10.8090

GP10        550233.3670         4676991.4200         76.2540         65.4290         10.8250

GP12        549686.1290         4678063.9360         61.4030         50.5570         10.8460

GP14        548788.7790         4679568.7950         67.6640         56.8300         10.8340

GP61        547159.8210         4681020.5790         66.2280         55.3650         10.8630

GP72        546454.4980         4682031.2550         77.2810         66.4210         10.8600

GP74        545641.3090         4683435.2830         75.4560         64.5800         10.8760

GP76        545338.0440         4685431.8100         82.7530         71.8770         10.8760

GP78        544950.4070         4686969.6300         76.6350         65.7530         10.8820

GP80        545136.9750         4689067.2420         90.6610         79.7980         10.8630

GP82        544982.2370         4691142.6270         83.5280         72.6500         10.8780

GP84        543988.8660         4692987.8600         97.9600         87.0590         10.9010

GP89        543967.4810         4693797.4610         101.7140         90.8010         10.9130

GP91        543435.3760         4695301.1550         100.9760         90.0520         10.9240

GP94        541845.3280         4696710.3880         107.9760         97.0080         10.9680

GP96        540154.0270         4697375.0850         124.5190         113.5050         11.0140

GP98        539091.4250         4698570.8010         137.5600         126.5330         11.0270

G100        538146.5870         4700027.6770         160.4450         149.3940         11.0510