求助：离散点拟合为光滑曲面，并得到方程

海天天

各位高手、大侠们：
       小弟遇到一棘手问题，恳请大家帮忙！
       马上就要写论文了，迫切需要曲面的方程，很着急，如有高手能抽时间看一下我的数据,提点建议，小弟不胜感激。数据如下，按z向坐标，这些数据可分为11层，每一层都是21个点:
    x                    y                 z
13.0044696        2.05963021        13.8
13.0490684        1.89011848        13.8
13.1272623        1.73188664        13.8
13.2357864        1.59223748        13.8
13.3699352        1.47749755        13.8
13.524613        1.39197574        13.8
13.6931075        1.33930283        13.8
13.8688601        1.32157422        13.8
14.1270913        1.32195656        13.8
14.3823859        1.30911225        13.8
14.6344467        1.27260632        13.8
14.883585        1.21522649        13.8
15.1284864        1.14226302        13.8
15.3695332        1.05849668        13.8
15.6078169        0.96620071        13.8
15.8432355        0.86643117        13.8
16.0758525        0.75982344        13.8
16.3054819        0.64662781        13.8
16.5315845        0.52600994        13.8
16.5374905        0.263067978        13.8
16.5395208        0                13.8

13.0067505        2.05987924        12.4
13.0513621        1.89042583        12.4
13.1295751        1.73229582        12.4
13.2380474        1.59278193        12.4
13.3720852        1.47815604        12.4
13.5265852        1.3926979        12.4
13.6948784        1.33998713        12.4
13.8705061        1.32212398        12.4
14.1286512        1.32234243        12.4
14.3839074        1.30936678        12.4
14.6359317        1.27279276        12.4
14.8849984        1.21540125        12.4
15.1298003        1.14243454        12.4
15.3707537        1.0586896        12.4
15.6089549        0.96643961        12.4
15.8443189        0.8667783        12.4
16.0768089        0.76035666        12.4
16.3062224        0.64709045        12.4
16.5324166        0.52612075        12.4
16.538635        0.262980025        12.4
16.540854        0                12.4

13.0092707        2.06008928        11
13.0538611        1.8905324        11
13.1320089        1.73222884        11
13.2405078        1.59247661        11
13.3746555        1.47761043        11
13.5293733        1.39198139        11
13.6979189        1.3392276        11
13.8737268        1.32145854        11
14.1319741        1.32179168        11
14.387195        1.30893896        11
14.6391021        1.27242841        11
14.8880374        1.21503988        11
15.1327095        1.14211092        11
15.373524        1.05837408        11
15.6115806        0.966115        11
15.8467729        0.86635665        11
16.0791814        0.75971046        11
16.3086993        0.64658803        11
16.5346766        0.52619587        11
16.5402717        0.26319896        11
16.5422466        0                11


13.017122        2.06101067        9.6
13.0618067        1.89169509        9.6
13.1401948        1.73374436        9.6
13.2487414        1.59452394        9.6
13.3827673        1.48025177        9.6
13.5371113        1.39505544        9.6
13.7051266        1.34241269        9.6
13.88054        1.32455697        9.6
14.1383068        1.32498784        9.6
14.3932339        1.31208433        9.6
14.6449518        1.27574297        9.6
14.8937354        1.21851236        9.6
15.1382763        1.14563498        9.6
15.3789768        1.06203668        9.6
15.616887        0.9697916        9.6
15.8519687        0.87017098        9.6
16.0841579        0.76376151        9.6
16.3131087        0.64980265        9.6
16.5389971        0.52702362        9.6
16.5458722        0.26336451        9.6
16.5482727        0                9.6


13.0229399        2.06165775        8.2
13.0677196        1.89250526        8.2
13.1457131        1.73463639        8.2
13.2541864        1.59532313        8.2
13.3881714        1.4811081        8.2
13.5422251        1.39548194        8.2
13.7102482        1.34271784        8.2
13.885709        1.32481954        8.2
14.1435012        1.32545383        8.2
14.3983427        1.3127061        8.2
14.6499301        1.2769147        8.2
14.8985625        1.2200172        8.2
15.1428947        1.1473607        8.2
15.3833638        1.0639987        8.2
15.6209777        0.9716217        8.2
15.8558041        0.8720168        8.2
16.0877135        0.76539667        8.2
16.3164015        0.65146788        8.2
16.5416267        0.52764307        8.2
16.5474447        0.26359153        8.2
16.5507692        0                8.2


13.0257705        2.06202153        6.8
13.0701869        1.89288871        6.8
13.1481719        1.73491569        6.8
13.2565235        1.59544587        6.8
13.3905839        1.48109002        6.8
13.5448521        1.39564361        6.8
13.7131571        1.3433951        6.8
13.888474        1.3255625        6.8
14.1462354        1.3261922        6.8
14.4010652        1.3134608        6.8
14.6526237        1.2776042        6.8
14.9011978        1.2206309        6.8
15.1454257        1.1479123        6.8
15.3857371        1.0644159        6.8
15.6231457        0.9719083        6.8
15.8577241        0.8720492        6.8
16.0893927        0.76522234        6.8
16.3178982        0.65109023        6.8
16.5430717        0.52756029        6.8
16.5487627        0.26360465        6.8
16.5516746        0.0002                6.8


13.027002        2.06239745        5.4
13.0714913        1.89313565        5.4
13.1495137        1.73518576        5.4
13.2578947        1.59569793        5.4
13.3920581        1.4814023        5.4
13.5464056        1.3961169        5.4
13.7149455        1.3440866        5.4
13.8901549        1.3263191        5.4
14.1477996        1.3267881        5.4
14.4025825        1.3139338        5.4
14.654102        1.2780009        5.4
14.9026268        1.2209295        5.4
15.1467982        1.1480652        5.4
15.3870799        1.064448        5.4
15.6244965        0.9719039        5.4
15.8590973        0.8720912        5.4
16.0907943        0.76530974        5.4
16.3193151        0.65128888        5.4
16.5443916        0.52763744        5.4
16.5500352        0.26360964        5.4
16.5529589        -0.0002                5.4


13.0285877        2.0624736        4
13.0730015        1.89343196        4
13.1508849        1.7355225        4
13.2592228        1.596127        4
13.3932944        1.4818877        4
13.5475116        1.3966804        4
13.7158299        1.3444781        4
13.8910108        1.3266817        4
14.1486613        1.327153        4
14.4034607        1.314206        4
14.6550008        1.2782445        4
14.9035322        1.2212682        4
15.1476999        1.1484895        4
15.3879504        1.0649439        4
15.6253044        0.9724016        4
15.8598285        0.8724991        4
16.0914363        0.7656414        4
16.3198947        0.65148092        4
16.5449764        0.52780514        4
16.5504319        0.26379337        4
16.5530251        -0.0002                4


13.0285629        2.0623368        2.6
13.0730273        1.89305527        2.6
13.1510057        1.7350543        2.6
13.259371        1.5955375        2.6
13.3934918        1.48115589        2.6
13.5478479        1.39569307        2.6
13.7163638        1.3435555        2.6
13.8916291        1.3259279        2.6
14.1492229        1.3267733        2.6
14.4038513        1.3140227        2.6
14.6552444        1.2781472        2.6
14.9037117        1.221295        2.6
15.1478228        1.1485833        2.6
15.3880454        1.0650951        2.6
15.6254036        0.9726298        2.6
15.8599337        0.8728105        2.6
16.0916081        0.7659131        2.6
16.3201596        0.65188658        2.6
16.5453154        0.52780264        2.6
16.5511322        0.26367869        2.6
16.5541204        -0.0001                2.6

13.0302935        2.06261089        1.2
13.0744498        1.8939664        1.2
13.1521239        1.7363318        1.2
13.2601861        1.5970862        1.2
13.393905        1.4830348        1.2
13.5476147        1.3974946        1.2
13.7156029        1.3451225        1.2
13.8905303        1.3269304        1.2
14.1479814        1.3271048        1.2
14.402671        1.3139423        1.2
14.6541062        1.2779696        1.2
14.9025115        1.2207982        1.2
15.1466393        1.147884        1.2
15.3869249        1.0642773        1.2
15.6243891        0.9718813        1.2
15.8590339        0.8722378        1.2
16.0906189        0.7657392        1.2
16.3187832        0.65109706        1.2
16.5438091        0.5277219        1.2
16.5491008        0.26385475        1.2
16.5512044        0                1.2


13.0235402        2.0616922        -0.2
13.0685179        1.89201803        -0.2
13.1469447        1.73380397        -0.2
13.2554988        1.59421547        -0.2
13.3896375        1.47942474        -0.2
13.5443532        1.39382456        -0.2
13.7128674        1.34115185        -0.2
13.888575        1.32357741        -0.2
14.1464199        1.32401236        -0.2
14.4011176        1.31132073        -0.2
14.6524688        1.27483551        -0.2
14.9009228        1.21755708        -0.2
15.1450764        1.14478062        -0.2
15.3853254        1.06107894        -0.2
15.6228031        0.96884012        -0.2
15.8573704        0.86903276        -0.2
16.0891078        0.76223967        -0.2
16.3178654        0.64796543        -0.2
16.5434788        0.52683629        -0.2
16.5488394        0.26342812        -0.2
16.5507845        0                        -0.2

matlab我刚接触，极度菜鸟，麻烦高手在指点时，尽量详细点，实在不好意思。
甚为感谢！！

dingd

将x看作因变量，y、z为自变量：
x = p1+p2*z+p3*Exp(-0.5*((y-p4)/p5)^2);

p1                 16.6462185832468
p2                 -0.00229226168361568
p3                 -3.79497705846567
p4                 1.81086694299024
p5                 0.520466292917289

mhl516

楼上的，能不能具体点，怎么求p1,P2,p3,p4,p5,能把你的程序发上来吗？

dingd

对这类拟合问题，用1stOpt，不用编程，也不用猜初值。

mhl516

谢谢
试试看

mhl516

没找到这个函数，能不能把你的程序直接发上来，谢谢啦

mhl516

这个函数是不是哪个工具箱里的啊，怎么我找不到呢

dingd

是个独立的优化软件，论坛已有些介绍。

NiuSD_2006

1stOpt 优化软件啊，看看
我有组测绘数据，假定方程很麻烦

海天天

to dingd
    很感谢。请问能让拟合的曲面更接近于数据点吗

dingd

如果想100％吻合，可考虑用神经网络来拟合。

海天天

哦
你还在吗，楼上的大侠
很想问你几个问题

海天天

to dingd

请问你有1stOpt的软件么，在网上找了很久找不到？
可以传一份给我吗，谢谢
我的邮箱：dtyazlm@yahoo.com.cn