————————————————————————————————————————————————————
现在我有分布在球表面上的有限个应力值(形式T(x,y,z)或T(r,theta,phi)),怎样通过内插得出整个球表面上的应力值,只分布在球表面上?
我自己内插了一下,可是在球面以内都出现应力值,不知道为什么?跪求高手赐教!
我的矩阵格式如下:
% x y z stress
2.864745084375788e-009 1.275466687660724e-009 2.983565686104821e-008 1433.076568587437
2.098295820241338e-009 2.330393597773816e-009 2.983565686104821e-008 1049.663577521429
9.690321463680712e-010 2.982374283870613e-009 2.983565686104821e-008 484.7542180077301
-3.27785988992915e-010 3.11867536226639e-009 2.983565686104821e-008 -163.9735496533228
-1.567926949014802e-009 2.715729138250094e-009 2.983565686104821e-008 -784.3488009263832
-2.536959095382873e-009 1.843208674605667e-009 2.983565686104821e-008 -1269.103018934113
……
内插的设置和绘制如下图:
标题: COMSOL为什么不能精确插值 [打印本页]
作者: sunnywell 时间: 2011-10-29 11:11:15 标题: COMSOL为什么不能精确插值
————————————————————————————————————————————————————
给出一些离散点的值,通过插值函数进行插值,为什么COMSOL插值的结果不能保证原给出值的离散点处的值不变?
作者: vinwim 时间: 2011-10-29 12:02:41
————————————————————————————————————————————————————
不对吧。这是我试验的线性插值结果,
作者: sunnywell 时间: 2011-10-29 13:17:02
————————————————————————————————————————————————————
vinwim 发表于 2011-10-29 12:02
不对吧。这是我试验的线性插值结果,
我的是二维的,中心点是有数据点的,数据的最大值大于7,而插值以后最大值不到6
作者: hit_hx 时间: 2011-10-31 07:43:03
————————————————————————————————————————————————————
我也遇到了这个问题,希望明白人能帮忙解释一下。
作者: sunnywell 时间: 2011-10-31 21:00:26
————————————————————————————————————————————————————
hit_hx 发表于 2011-10-31 07:43
我也遇到了这个问题,希望明白人能帮忙解释一下。
中仿给的答复:
COMSOL绘制的函数图,是通过格点插值取得的,这些格点大小相等,比方说是1e-9的正方形小格子,那么就您给出的表单里的数值,包括值大于7的中心点,都可能落在正方形格子的内部,那么再插值的话就会使最大值不到6.但是,这只是COMSOL显示函数图形上些许不精确,当我们调用该函数求解模型时,通过加密中心处的网格,就能使插值函数的值非常精确,而不用担心得不到最大值7。
作者: sunnywell 时间: 6 天前
————————————————————————————————————————————————————
hit_hx 发表于 2011-10-31 07:43
我也遇到了这个问题,希望明白人能帮忙解释一下。
这种插值方法可控性还是差了点,虽然有几个可以选择的,其中对于大于一维的数据只能使用其中两个,而且差别不大,都不精确。既然是插值,就应该在我们已有数据间插值,肯定要保证我们的原有数据不变。
你有何进展吗?QQ544732871
https://wap.sciencenet.cn/blog-506324-505794.html
上一篇:
COMSOL使用问题2:散射声场【已解决】下一篇:
COMSOL使用问题4:能带结构在T点突变问题