自己编的程序是
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