声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 4248|回复: 26

[编程技巧] [求助]请教一个曲面拟合问题

[复制链接]
发表于 2006-3-11 10:17 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
数据如下,那个例子我参考过了,也得到了拟合函数,不过用插值法得到的数据不对,与原数据差别太大,请各位高手再指点小弟一下。
x=[6 13 14 26 29 32 50 59 69 74 79 82];

y=[0.1 0.15 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 4.6 4.8 5.0 5.5 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10 11 12 13 14 15 16 17 18 19 20 22 24 26 28 30 32 34 36 38 40 45 50 55 60 65 70 75 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280 290 300 ];


z=[792 598 578 412 384 360 261 230 202 191 181 176 ...
624 485 470 344 322 303 223 197 175 165 157 152 ...
526 416 404 301 282 266 199 176 156 148 141 137 ...
505 311 303 257 242 227 157 134 117 112 109 109 ...
417 273 265 218 206 194 136 117 103 98.9 96.8 96.4 ...
360 248 239 192 181 171 122 105 93.3 89.8 88 87.7 ...
319 228 220 174 164 155 111 96.7 86.1 83 81.4 81.2 ...
288 213 205 159 150 142 103 90 80.5 77.7 76.2 76 ...
264 197 189 148 140 132 96.3 84.2 75.4 72.8 71.5 71.3 ...
245 183 177 139 131 124 90.6 79.4 71.2 68.8 67.5 67.3 ...
229 172 166 131 124 117 85.9 75.3 67.6 65.4 64.2 64 ...
202 152 147 117 111 105 77.6 68.3 61.6 59.6 68.6 58.4 ...
181 138 133 107 101 95.8 71.2 63 56.9 55.2 54.3 54.1 ...
165 126 122 98.2 93.2 88.6 66.2 58.7 53.2 51.6 50.8 50.6 ...
152 116 113 91.4 86.8 82.5 62 55.1 50.1 48.6 47.9 47.7 ...
141 109 105 85.7 81.4 77.5 58.5 52.1 47.5 46.1 45.4 45.3 ...
132 102 98.8 80.8 76.9 73.2 55.5 49.5 45.2 44 43.3 43.2 ...
124 96.2 93.3 76.6 72.9 69.4 52.9 47.3 43.3 42.1 41.5 41.4 ...
117 91.3 88.6 72.9 69.5 66.2 50.6 45.3 41.5 40.4 39.8 39.7 ...
111 86.9 84.3 69.7 66.4 63.3 48.5 43.6 40 38.9 38.4 38.3 ...
106 83 80.6 66.8 63.7 60.8 46.7 42 38.6 37.6 37.1 37 ...
101 79.2 76.9 63.9 60.9 58.1 44.7 40.2 37 36 35.5 35.4 ...
96.2 75.8 73.7 61.2 58.4 55.7 42.9 38.6 35.5 34.6 34.1 34 ...
92 72.8 70.7 58.8 56.1 53.5 41.3 37.2 34.2 33.3 32.8 32.8 ...
88.3 70 68 56.6 54 51.6 39.8 35.8 32.9 32.1 31.7 31.6 ...
84.8 67.4 65.5 54.7 52.1 49.8 38.4 34.6 31.8 31 30.6 30.5 ...
81.7 65.1 63.3 52.8 50.4 48.1 37.2 33.5 30.8 30 29.6 29.5 ...
78.8 62.9 61.2 51.1 48.8 46.6 36 32.4 29.9 29.1 28.7 28.6 ...
76.1 60.9 59.3 49.6 47.3 45.2 34.9 31.5 29 28.2 27.9 27.8 ...
73.7 59.1 57.5 48.1 45.9 43.9 33.9 30.6 28.2 27.5 27.1 27 ...
71.4 57.4 55.8 46.8 44.6 42.6 33 29.8 27.4 26.7 26.4 26.3 ...
66.3 53.6 52.1 43.8 41.8 39.9 30.9 27.9 25.7 25.1 24.7 24.7 ...
62 50.3 49 41.2 39.3 37.6 29.1 26.3 24.2 23.6 23.3 23.3 ...
58.2 47.5 46.2 39 37.2 35.5 27.6 24.9 23 22.4 22.1 22.1 ...
55 45 43.8 37 35.3 33.8 26.2 23.7 21.9 21.3 21 21 ...
52 42.6 41.5 35.2 33.7 32.2 25.2 22.8 21.1 20.6 20.4 20.3 ...
49.3 40.5 39.5 33.6 32.2 30.8 24.3 22.1 20.4 20 19.7 19.7 ...
46.9 38.6 37.7 32.2 30.9 29.6 23.4 21.3 19.8 19.4 19.2 19.1 ...
44.7 36.9 36 30.9 29.7 28.4 22.6 20.7 19.3 18.9 18.7 18.6 ...
42.8 35.4 34.6 29.8 28.6 27.4 21.9 20.1 18.8 18.4 18.2 18.1 ...
41 34 33.2 28.7 27.6 26.5 21.3 19.6 18.3 17.9 17.7 17.7 ...
38 31.6 30.8 26.7 25.6 24.7 19.9 18.3 17.1 16.8 16.6 16.6 ...
35.4 29.5 28.8 25 24 23.1 18.7 17.2 16.1 15.8 15.7 15.6 ...
33.2 27.7 27.1 23.5 22.6 21.8 17.1 16.3 15.3 15 14.8 14.8 ...
31.3 26.1 25.6 22.2 21.4 20.6 16.8 15.5 14.5 14.3 14.1 14.1 ...
29.6 24.8 24.2 21.1 20.3 19.6 15.9 14.7 13.9 13.6 13.5 13.4 ...
28.1 23.5 23 20.1 19.4 18.6 15.2 14.1 13.3 13 12.9 12.9 ...
26.7 22.5 22 19.2 18.5 17.8 14.6 13.5 12.7 12.5 12.4 12.3 ...
25.5 21.5 21 18.4 17.7 17.1 14 13 12.2 12 11.9 11.9 ...
24.4 20.6 20.1 17.7 17 16.4 13.5 12.5 11.8 11.6 11.5 11.5 ...
23.4 19.8 19.4 17 16.4 15.8 13 12.1 11.4 11.2 11.1 11.1 ...
21.7 18.4 18 15.8 15.2 14.7 12.1 11.3 10.7 10.5 10.4 10.4 ...
20.2 17.1 16.8 14.8 14.3 13.8 11.4 10.6 10 9.88 9.79 9.77 ...
19 16.1 15.8 13.9 13.4 13 10.8 10 9.51 9.35 9.27 9.26 ...
17.9 15.2 14.9 13.2 12.7 12.3 10.2 9.54 9.04 8.89 8.82 8.8 ...
16.9 14.4 14.1 12.5 12.1 11.7 9.74 9.09 8.62 8.48 8.41 8.4 ...
16.1 13.7 13.4 11.9 11.5 11.1 9.29 8.68 8.23 8.1 8.03 8.02 ...
15.3 13.1 12.8 11.4 11 10.6 8.89 8.31 7.88 7.76 7.7 7.69 ...
14.6 12.5 12.3 10.9 10.5 10.2 8.53 7.97 7.57 7.45 7.39 7.38 ...
14 12 11.8 10.5 10.1 9.78 8.2 7.67 7.28 7.17 7.11 7.1 ...
13.4 11.5 11.3 10.1 9.73 9.41 7.9 7.39 7.02 6.91 6.86 6.85 ...
12.2 10.5 10.3 9.21 8.9 8.62 7.25 6.79 6.45 6.36 6.31 6.3 ...
11.2 9.72 9.54 8.51 8.23 7.97 6.71 6.29 5.99 5.9 5.85 5.84 ...
10.4 9.03 8.87 7.92 7.66 7.42 6.26 5.87 5.59 5.51 5.47 5.46 ...
9.72 8.45 8.29 7.41 7.18 6.95 5.88 5.52 5.26 5.18 5.14 5.13 ...
9.12 7.94 7.8 6.98 6.76 6.55 5.54 5.21 4.96 4.89 4.85 4.85 ...
8.6 7.5 7.37 6.6 6.39 6.19 5.25 4.94 4.71 4.64 4.61 4.6 ...
8.17 7.11 6.99 6.27 6.08 5.89 5 4.71 4.49 4.43 4.39 4.39 ...
7.79 6.77 6.65 5.98 5.8 5.62 4.78 4.5 4.29 4.23 4.2 4.2 ...
7.14 6.18 6.08 5.49 5.32 5.16 4.4 4.14 3.96 3.9 3.87 3.87 ...
6.6 5.7 5.61 5.08 4.93 4.78 4.08 3.85 3.68 3.63 3.6 3.6 ...
6.2 5.36 5.27 4.78 4.63 4.5 3.84 3.62 3.46 3.42 3.39 3.39 ...
5.85 5.07 4.98 4.52 4.38 4.25 3.64 3.43 3.28 3.24 3.21 3.21 ...
5.55 4.81 4.73 4.29 4.16 4.04 3.46 3.26 3.12 3.08 3.05 3.05 ...
5.29 4.58 4.51 4.09 3.97 3.85 3.3 3.11 2.98 2.94 2.91 2.91 ...
5.05 4.38 4.32 3.92 3.8 3.69 3.16 2.98 2.85 2.81 2.79 2.79 ...
4.84 4.21 4.14 3.76 3.64 3.54 3.03 2.86 2.73 2.7 2.68 2.68 ...
4.65 4.04 3.98 3.61 3.51 3.4 2.91 2.75 2.63 2.6 2.58 2.58 ...
4.48 3.9 3.84 3.48 3.38 3.28 2.81 2.65 2.54 2.5 2.49 2.48 ...
4.32 3.76 3.7 3.36 3.26 3.17 2.72 2.56 2.45 2.42 2.4 2.4 ...
4.18 3.64 3.58 3.26 3.16 3.07 2.63 2.48 2.37 2.34 2.33 2.32 ...
4.05 3.53 3.47 3.16 3.06 2.97 2.55 2.41 2.3 2.27 2.26 2.25 ...
3.93 3.42 3.37 3.06 2.97 2.89 2.47 2.34 2.24 2.21 2.19 2.19 ...
3.81 3.33 3.27 2.98 2.89 2.8 2.41 2.27 2.17 2.15 2.13 2.13 ...
3.71 3.23 3.18 2.9 2.81 2.73 2.34 2.21 2.12 2.09 2.07 2.07 ...
3.61 3.15 3.1 2.82 2.74 2.66 2.28 2.15 2.06 2.04 2.02 2.02 ...
3.52 3.07 3.02 2.75 2.67 2.59 2.23 2.1 2.01 1.99 1.97 1.97 ...
3.43 3 2.95 2.69 2.61 2.53 2.17 2.05 1.96 1.94 1.93 1.92 ...
3.35 2.93 2.88 2.62 2.55 2.47 2.12 2.01 1.92 1.9 1.88 1.88 ...
3.27 2.86 2.82 2.57 2.49 2.42 2.08 1.96 1.88 1.85 1.84 1.84 ...
3.2 2.8 2.76 2.51 2.44 2.37 2.03 1.92 1.84 1.81 1.8 1.8]
回复
分享到:

使用道具 举报

发表于 2006-3-11 15:41 | 显示全部楼层
看了一下你的数据好像不太合适用
http://forum.vibunion.com/thread-1005-1-1.html
提供的函数进行拟合
你可以找一下NAG Foundation Toolbox,这个工具箱也提供了几个曲面拟合的函数,不过现在我手头没有这个工具箱
 楼主| 发表于 2006-3-12 09:55 | 显示全部楼层
好的,问题已经解决,谢谢,非常谢谢。忘记说明了,上面是12*91=1092的矩阵,现在已经纠正过来了。
现在曲面拟合函数已经得到了,但怎么才能得到最好的P和q值,有没有求误差函数的方法。
谢谢  happy
发表于 2006-3-12 11:41 | 显示全部楼层
p和q需要试的,太大了容易出现数据震荡现象,太小了又拟合的不太好
 楼主| 发表于 2006-3-12 14:33 | 显示全部楼层
谢谢,两维拟合可以返回一个误差函数,请问三维拟合可以返回吗?
如何画出经拟合后的函数图像。
发表于 2006-3-12 19:29 | 显示全部楼层
谢谢,两维拟合可以返回一个误差函数,请问三维拟合可以返回吗?
如何画出经拟合后的函数图像。

好像这个工具包没有提供返回误差函数,如果需要那就要自己变成了
画图用mesh
 楼主| 发表于 2006-3-12 22:09 | 显示全部楼层
谢谢,方差函数已经求出。
请问用mesh怎么画图
比如x=[6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82 6 13 14 26 29 32 50 59 69 74 79 82];

y=[0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.6 1.6 1.6 1.6 1.6 1.6 1.6 1.6 1.6 1.6 1.6 1.6 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.2 2.2 2.2 2.2 2.2 2.2 2.2 2.2 2.2 2.2 2.2 2.2 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.6 2.6 2.6 2.6 2.6 2.6 2.6 2.6 2.6 2.6 2.6 2.6 2.8 2.8 2.8 2.8 2.8 2.8 2.8 2.8 2.8 2.8 2.8 2.8 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.4 3.4 3.4 3.4 3.4 3.4 3.4 3.4 3.4 3.4 3.4 3.4 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.8 3.8 3.8 3.8 3.8 3.8 3.8 3.8 3.8 3.8 3.8 3.8 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.2 4.2 4.2 4.2 4.2 4.2 4.2 4.2 4.2 4.2 4.2 4.2 4.4 4.4 4.4 4.4 4.4 4.4 4.4 4.4 4.4 4.4 4.4 4.4 4.6 4.6 4.6 4.6 4.6 4.6 4.6 4.6 4.6 4.6 4.6 4.6 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 5.0 5.0 5.0 5.0 5.0 5.0 5.0 5.0 5.0 5.0 5.0 5.0 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 6.0 6.0 6.0 6.0 6.0 6.0 6.0 6.0 6.0 6.0 6.0 6.0 6.5 6.5 6.5 6.5 6.5 6.5 6.5 6.5 6.5 6.5 6.5 6.5 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.0 7.5 7.5 7.5 7.5 7.5 7.5 7.5 7.5 7.5 7.5 7.5 7.5 8.0 8.0 8.0 8.0 8.0 8.0 8.0 8.0 8.0 8.0 8.0 8.0 8.5 8.5 8.5 8.5 8.5 8.5 8.5 8.5 8.5 8.5 8.5 8.5 9.0 9.0 9.0 9.0 9.0 9.0 9.0 9.0 9.0 9.0 9.0 9.0 9.5 9.5 9.5 9.5 9.5 9.5 9.5 9.5 9.5 9.5 9.5 9.5 10 10 10 10 10 10 10 10 10 10 10 10 11 11 11 11 11 11 11 11 11 11 11 11 12 12 12 12 12 12 12 12 12 12 12 12 13 13 13 13 13 13 13 13 13 13 13 13 14 14 14 14 14 14 14 14 14 14 14 14 15 15 15 15 15 15 15 15 15 15 15 15 16 16 16 16 16 16 16 16 16 16 16 16 17 17 17 17 17 17 17 17 17 17 17 17 18 18 18 18 18 18 18 18 18 18 18 18 19 19 19 19 19 19 19 19 19 19 19 19 20 20 20 20 20 20 20 20 20 20 20 20 22 22 22 22 22 22 22 22 22 22 22 22 24 24 24 24 24 24 24 24 24 24 24 24 26 26 26 26 26 26 26 26 26 26 26 26 28 28 28 28 28 28 28 28 28 28 28 28 30 30 30 30 30 30 30 30 30 30 30 30 32 32 32 32 32 32 32 32 32 32 32 32 34 34 34 34 34 34 34 34 34 34 34 34 36 36 36 36 36 36 36 36 36 36 36 36 38 38 38 38 38 38 38 38 38 38 38 38 40 40 40 40 40 40 40 40 40 40 40 40 45 45 45 45 45 45 45 45 45 45 45 45 50 50 50 50 50 50 50 50 50 50 50 50 55 55 55 55 55 55 55 55 55 55 55 55 60 60 60 60 60 60 60 60 60 60 60 60 65 65 65 65 65 65 65 65 65 65 65 65 70 70 70 70 70 70 70 70 70 70 70 70 75 75 75 75 75 75 75 75 75 75 75 75 80 80 80 80 80 80 80 80 80 80 80 80 90 90 90 90 90 90 90 90 90 90 90 90 100 100 100 100 100 100 100 100 100 100 100 100 110 110 110 110 110 110 110 110 110 110 110 110 120 120 120 120 120 120 120 120 120 120 120 120 130 130 130 130 130 130 130 130 130 130 130 130 140 140 140 140 140 140 140 140 140 140 140 140 150 150 150 150 150 150 150 150 150 150 150 150 160 160 160 160 160 160 160 160 160 160 160 160 170 170 170 170 170 170 170 170 170 170 170 170 180 180 180 180 180 180 180 180 180 180 180 180 190 190 190 190 190 190 190 190 190 190 190 190 200 200 200 200 200 200 200 200 200 200 200 200 210 210 210 210 210 210 210 210 210 210 210 210 220 220 220 220 220 220 220 220 220 220 220 220 230 230 230 230 230 230 230 230 230 230 230 230 240 240 240 240 240 240 240 240 240 240 240 240 250 250 250 250 250 250 250 250 250 250 250 250 260 260 260 260 260 260 260 260 260 260 260 260 270 270 270 270 270 270 270 270 270 270 270 270 280 280 280 280 280 280 280 280 280 280 280 280 290 290 290 290 290 290 290 290 290 290 290 290 300 300 300 300 300 300 300 300 300 300 300 300];
z=108.07-0.55948*y-0.97787*x+0.0051668*x.*y
用mesh(x,y,z )没法画出,该怎么办???
再次感谢。
 楼主| 发表于 2006-3-13 08:58 | 显示全部楼层
最小误差函数已经找到,但是带到原数据里边得到的数值与原数据差别很大,这说明最小二乘法对本组数据不太适合,请问有没有更好的曲面拟合方法。<BR>谢谢。
发表于 2006-3-13 09:04 | 显示全部楼层
谢谢,方差函数已经求出。
请问用mesh怎么画图
用mesh(x,y,z )没法画出,该怎么办???
再次感谢。

你给的数据太乱了,参考mesh的帮助就能画出来
发表于 2006-3-13 09:06 | 显示全部楼层
最小误差函数已经找到,但是带到原数据里边得到的数值与原数据差别很大,这说明最小二乘法对本组数据不太适合,请问有没有更好的曲面拟合方法。谢谢。

实施用样条曲面拟合
曾经看到过一篇文章《样条曲面拟合及其Matlab实现》,可以去期刊网找找看
 楼主| 发表于 2006-3-13 09:30 | 显示全部楼层
谢谢
 楼主| 发表于 2006-3-17 18:43 | 显示全部楼层
谢谢,那篇文章找到了。可是没法使用。有没有直接进行曲面拟合的程序,数据如下所示。横坐标为x,竖坐标为y,y与x相对应的就是z
100 200 300 400 500
100 636 697 624 478 450
200 698 712 630 478 420
300 680 674 598 412 400
400 662 626 552 334 310
发表于 2006-3-17 23:32 | 显示全部楼层
不规则数据的曲面拟合问题
我现在有如下的测量数据:
[x,y,z]=[1256.89 27.45 23.47;
1278.23 29.69 24.00;
1269.53 28.32 25.00;
1260.23 28.96 29.30;]
怎么样运用拟合的方法,得到一个曲面,从而插值出其他的值,并不是运用多项式拟合的方法!
望楼主指点!
不胜感谢!
发表于 2006-3-18 08:16 | 显示全部楼层
不规则数据的曲面拟合问题
我现在有如下的测量数据:
[x,y,z]=[1256.89 27.45 23.47;
1278.23 29.69 24.00;
1269.53 28.32 25.00;
1260.23 28.96 29.30;]
怎么样运用拟合的方法,得到一个曲面,从而插值出其他的值,并不是运用多项式拟合的方法!
望楼主指点!
不胜感谢!


数据太少了
发表于 2006-3-18 08:18 | 显示全部楼层
谢谢,那篇文章找到了。可是没法使用。有没有直接进行曲面拟合的程序,数据如下所示。横坐标为x,竖坐标为y,y与x相对应的就是z
100 200 300 400 500
100 636 697 624 478 450
200 698 712 630 478 420
300 680 674 598 412 400
400 662 626 552 334 310


你去找一下NAG Foundation Toolbox吧,里边有几个拟合函数,不过我也一直没找到,好运
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-20 07:30 , Processed in 0.064412 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表