试试吧。不过我还是建议根据你的问题换一种插值的方法，我这个比较蛮干，而且也不是非常准确。
</p>
<p>a=read.csv("DTT1.csv",head=F)<br />
a=as.matrix(a)<br />
a0=data.frame(x=rep(a[1,-1],each=nrow(a)-1),y=rep(a[-1,1],ncol(a)-1),z=as.numeric(a[-1,-1]))<br />
a0[is.na(a0)]=1</p>
<p>eps0=0.000001;sigma0=2;C0=10000<br />
library(kernlab)<br />
md=ksvm(data=a0,z~x+y,C=C0,epsilon=eps0,kpar=list(sigma=sigma0))</p>
<p>x=seq(0.2,0.5,by=0.005)<br />
y=seq(11800,15600,by=20)<br />
zz=expand.grid(x=x,y=y)<br />
zz$z1=predict(md,zz)</p>
<p>library(lattice)</p>
<p>png("3d-svm.png",height=600,width=600)<br />
par(mar=rep(0,4),ann=F)<br />
wireframe(z1~x*y,data=zz,outer=T,shade=T,zlab="z1")<br />
dev.off()<br />
</p>
