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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 多元多项式重构问题请教
101/1返回列表
查看: 1871  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上提前没消息，有中的吗 +9 archvillain 2026-08-02 11/550 2026-08-02 15:44 by alexon
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 8/400 2026-08-02 15:01 by jHcrHMhbKCFV
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 14:58 by jHcrHMhbKCFV
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 18/900 2026-08-02 14:52 by jHcrHMhbKCFV
[找工作] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 12/600 2026-08-02 13:43 by jHcrHMhbKCFV
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +17 低垂的野花 2026-07-31 23/1150 2026-08-02 13:41 by RRRKKK
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 13:34 by jHcrHMhbKCFV
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 13:34 by jHcrHMhbKCFV
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[有机交流] 一个有机合成实验室都需要哪些设备？ 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
信息提示
关闭
请填处理意见
关闭