| 查看: 1934 | 回复: 9 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
多元多项式重构问题请教 已有1人参与
|
||||
|
各位前辈好，最近在工作中遇到一个问题，想请教大家。 我想用高阶的多项式重构一个光滑变量，比如，f(x,y)≈a10*x+a01*y+a20*x*x+a02*y*y+a11*x*y..。 难点在于多项式系数如何确定。 目前我用是对f求偏导，然后用有限差分近似(在均匀网格上)，从而确定polynomial的系数。 比如 ∂f/∂x =a10, ∂f/∂y=a01...，∂f/∂x ≈[f(i+1)-f(i-1)]/2.0/dx, 这个方法的缺点是计算高阶偏导数不好处理，比如我要算二阶混合导数fxy，需要先算出一阶导数fx,fy，然后再对fx,fy求偏导。 上述做法在用高阶有限差分或者计算高阶导数时，会引入很宽的模板。具体应用时，效果不是特别好。 当然，我也可以由f计算出各阶偏导数，相对上述方法更加compact。比如fxy=[f(i+1,j+1)-f(i+1,j-1)-f(i-1,j+1)+f(i-1,j-1)]/4.0/dx/dy. 不过，高阶混合导数有限差分格式推导的工作量实在太大，特别是三维。。。 想请教大家，有没有compact，并且比较简单的方法，比如least square？ 谢谢！ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
9楼2017-12-15 08:20:34
zaq123321
专家顾问 (著名写手)
-
专家经验: +343 - 数学EPI: 6
- 应助: 298 (大学生)
- 贵宾: 0.247
- 金币: 11171.3
- 散金: 545
- 红花: 29
- 帖子: 1223
- 在线: 539.9小时
- 虫号: 405284
- 注册: 2007-06-17
- 性别: MM
- 专业: 生物大分子结构与功能
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lovvv002569: 金币+10, ★★★★★最佳答案 2017-12-13 13:07:23
感谢参与，应助指数 +1
lovvv002569: 金币+10, ★★★★★最佳答案 2017-12-13 13:07:23
|
比如 ∂f/∂x =a10, ∂f/∂y=a01...，∂f/∂x ≈[f(i+1)-f(i-1)]/2.0/dx, Is this correct? If your model is fixed, like here you choose multiple variable polynomial form, you just need build a overdetermined linear system and then use linear square method to get your best coefficients. [ 发自手机版 http://muchong.com/3g ] |
2楼2017-12-11 10:52:43
3楼2017-12-11 13:49:35
zaq123321
专家顾问 (著名写手)
-
专家经验: +343 - 数学EPI: 6
- 应助: 298 (大学生)
- 贵宾: 0.247
- 金币: 11171.3
- 散金: 545
- 红花: 29
- 帖子: 1223
- 在线: 539.9小时
- 虫号: 405284
- 注册: 2007-06-17
- 性别: MM
- 专业: 生物大分子结构与功能
- 管辖: 数学
【答案】应助回帖
|
Are you sure what you have written done (1) is correct? Do you know how to take the derivative of f(x,y) with respect to x according to your formulation? I don't understand what you mean second order of accuracy. It looks to me this is an over determined problem or underdetermined problem [ 发自手机版 http://muchong.com/3g ] |
4楼2017-12-13 10:25:48
