海天天 发表于 2006-8-21 10:15

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

各位高手、大侠们:
       小弟遇到一棘手问题,恳请大家帮忙!
       马上就要写论文了,迫切需要曲面的方程,很着急,如有高手能抽时间看一下我的数据,提点建议,小弟不胜感激。数据如下,按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 发表于 2006-8-21 21:12

将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 发表于 2006-8-22 10:13

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

dingd 发表于 2006-8-22 10:30

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

mhl516 发表于 2006-8-22 10:44

谢谢
试试看

mhl516 发表于 2006-8-22 10:59

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

mhl516 发表于 2006-8-22 11:22

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

dingd 发表于 2006-8-22 11:56

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

NiuSD_2006 发表于 2006-8-24 00:01

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

海天天 发表于 2006-8-24 17:11

谢谢

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

dingd 发表于 2006-8-24 20:43

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

海天天 发表于 2006-8-24 21:02


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

海天天 发表于 2006-8-27 17:10

to dingd

请问你有1stOpt的软件么,在网上找了很久找不到?
可以传一份给我吗,谢谢
我的邮箱:dtyazlm@yahoo.com.cn
页: [1]
查看完整版本: 求助:离散点拟合为光滑曲面,并得到方程