应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 多元多项式重构问题请教
101/1返回列表
查看: 1933  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lovvv002569

铁虫 (小有名气)

[求助] 多元多项式重构问题请教 已有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？
谢谢！
回复此楼

» 猜你喜欢
1楼 2017-12-10 17:22:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zaq123321

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lovvv002569

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by zaq123321 at 2017-12-11 10:52:43
比如 ∂f/∂x =a10, ∂f/∂y=a01...，∂f/∂x ≈/2.0/dx,  Is this correct?

If your model is fixed, like here you choose multiple variable polynomial form, you just need ...

你好，非常感谢你的回复。
(1) 这个表达式的确有问题，正切的写法是: a10=∂f/∂x|i,j, ∂f/∂x|i,j ≈[f(i+1)-f(i-1)]/2.0/dx.
(2) 您说的linear square method是指linear square fiting吗? 貌似linear square 最多只有二阶精度?
赞一下 回复此楼
3楼2017-12-11 13:49:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zaq123321

专家顾问 (著名写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by lovvv002569 at 2017-12-11 13:49:35
你好，非常感谢你的回复。
(1) 这个表达式的确有问题，正切的写法是: a10=∂f/∂x|i,j, ∂f/∂x|i,j ≈/2.0/dx.
(2) 您说的linear square method是指linear square fiting吗? 貌似linear sq ...

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lovvv002569

铁虫 (小有名气)
引用回帖:
4楼: Originally posted by zaq123321 at 2017-12-13 10:25:48
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 ord ...

前辈，感谢您的回复，第一个方法我略去了一些条件，这边就不再详细叙述了；但我已经在重构算例中测出了相应的精度。
关于二阶精度是我搞错了linear的概念，我以为linear指的是linear polynomial。惭愧。。。
我这两天认真研究了下least square，对于polynomial，用least square确定系数的确是一个非常不错的方法。
感谢您的帮助！
赞一下 回复此楼
5楼2017-12-13 13:07:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zaq123321

专家顾问 (著名写手)
The (1) you wrote is not correct. . It's different with yours

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
小木虫给我温暖，给我希望，爱就要爱小木虫。
6楼2017-12-13 17:50:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lovvv002569

铁虫 (小有名气)
引用回帖:
6楼: Originally posted by zaq123321 at 2017-12-13 17:50:07
The (1) you wrote is not correct.  \partial f/\partial x=a10+2*a20*x+a11*y+.... It's different with yours

前辈,
这个重构的基本思想是对变量f在点(i,j)处进行泰勒展开，f(x,y)=f(i,j)+∂f/∂x(i,j)*(x-xi)+∂f/∂y(i,j)*(y-yi)+ ...
我把上述式子改写成 f=f(i,j)+a10*X+a01*Y+...
其中a10=∂f/∂x(i,j), X=x-xi
您写的公式里，后面几项在离散时都可以约去。
赞一下 回复此楼
7楼2017-12-13 18:20:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zaq123321

专家顾问 (著名写手)
Where is f(i,j) in your first formulation? Do not mix Taylor expansion and fitting function. They are different.

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
小木虫给我温暖，给我希望，爱就要爱小木虫。
8楼2017-12-13 20:32:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lovvv002569

铁虫 (小有名气)
引用回帖:
8楼: Originally posted by zaq123321 at 2017-12-13 20:32:48
Where is f(i,j) in your first formulation? Do not mix Taylor expansion and fitting function. They are different.

前辈，不好意思，第一个公式写的时候少写了这一项。
对网格(i,j)内的光滑变量f重构，f(x,y)=f(i,j)+∂f/∂x(i,j)*(x-xi)+∂f/∂y(i,j)*(y-yi)+ ...

我不太明白您说的多项式重构和泰勒展开的具体区别，可以再解释下吗？
赞一下 回复此楼
9楼2017-12-15 08:20:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zaq123321

专家顾问 (著名写手)
I didn't not say that.

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
小木虫给我温暖，给我希望，爱就要爱小木虫。
10楼2017-12-15 09:07:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lovvv002569 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:52 by o3EMxHUAQDtS
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:40 by o3EMxHUAQDtS
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉？ +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭