24小时热门版块排行榜    

查看: 1956  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

KGFan

银虫 (小有名气)

[求助] 想用mathematica求多项式的乘积,求指点

如图,有两个系数待定,阶数为n的x,y的多项式,怎么求出它们的乘积的表达式,要求结果也按照阶数展开。本人是mathematica小白,误鄙视。。。
想用mathematica求多项式的乘积,求指点

Photo Mar 31, 12 26 28 PM.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KGFan

银虫 (小有名气)

引用回帖:
2楼: Originally posted by dollapym at 2015-04-01 10:04:01
可以求出前面几项的所求目标,比如n=2,3,总结n项时规律.下面是n=2时的情况,系数矩阵可以自己写出来,总结n项的规律
A2 = {1, x, y, x^2, x*y, y^2}; B2 = {a1, a2, a3, a4, a5,
  a6}; C2 = {b1, b2, b3, b4,  ...

非常感谢!再请教一下,有没有办法合并同类项,然后还是按照x^3, x^2*y, x*y^2, y^3这个顺序排列?因为我主要是要得到x^i*y^j这些项的系数表达式。
3楼2015-04-01 11:43:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

dollapym

铁虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
KGFan: 金币+15, ★★★很有帮助 2015-04-01 11:43:31
可以求出前面几项的所求目标,比如n=2,3,总结n项时规律.下面是n=2时的情况,系数矩阵可以自己写出来,总结n项的规律
A2 = {1, x, y, x^2, x*y, y^2}; B2 = {a1, a2, a3, a4, a5,
  a6}; C2 = {b1, b2, b3, b4, b5, b6}; Q = A2.B2
P = A2.C2

Out[3239]= a1 + a2 x + a4 x^2 + a3 y + a5 x y + a6 y^2

Out[3240]= b1 + b2 x + b4 x^2 + b3 y + b5 x y + b6 y^2

In[3251]:= Q1[x_, y_] = Q;

In[3247]:= S = Expand[P*Q]

Out[3247]=
a1 b1 + a2 b1 x + a1 b2 x + a4 b1 x^2 + a2 b2 x^2 + a1 b4 x^2 +
a4 b2 x^3 + a2 b4 x^3 + a4 b4 x^4 + a3 b1 y + a1 b3 y + a5 b1 x y +
a3 b2 x y + a2 b3 x y + a1 b5 x y + a5 b2 x^2 y + a4 b3 x^2 y +
a3 b4 x^2 y + a2 b5 x^2 y + a5 b4 x^3 y + a4 b5 x^3 y + a6 b1 y^2 +
a3 b3 y^2 + a1 b6 y^2 + a6 b2 x y^2 + a5 b3 x y^2 + a3 b5 x y^2 +
a2 b6 x y^2 + a6 b4 x^2 y^2 + a5 b5 x^2 y^2 + a4 b6 x^2 y^2 +
a6 b3 y^3 + a3 b6 y^3 + a6 b5 x y^3 + a5 b6 x y^3 + a6 b6 y^4

In[3255]:= S1 = Expand[P*D[Q1[x, y], x, y]]

Out[3255]=
a5 b1 + a5 b2 x + a5 b4 x^2 + a5 b3 y + a5 b5 x y + a5 b6 y^2
2楼2015-04-01 10:04:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
feixiaolin: 金币+20, http://emuch.net/bbs/viewthread.php?tid=8740247&fpage=1 2015-04-01 22:53:13
KGFan(feixiaolin代发): 金币+10 2015-04-01 22:53:21
引用回帖:
3楼: Originally posted by KGFan at 2015-04-01 11:43:09
非常感谢!再请教一下,有没有办法合并同类项,然后还是按照x^3, x^2*y, x*y^2, y^3这个顺序排列?因为我主要是要得到x^i*y^j这些项的系数表达式。...

Collect加模式匹配:
CODE:
Collect[a1 b1 + a2 b1 x + a1 b2 x + a4 b1 x^2 + a2 b2 x^2 +
  a1 b4 x^2 + a4 b2 x^3 + a2 b4 x^3 + a4 b4 x^4 + a3 b1 y + a1 b3 y +
  a5 b1 x y + a3 b2 x y + a2 b3 x y + a1 b5 x y + a5 b2 x^2 y +
  a4 b3 x^2 y + a3 b4 x^2 y + a2 b5 x^2 y + a5 b4 x^3 y +
  a4 b5 x^3 y + a6 b1 y^2 + a3 b3 y^2 + a1 b6 y^2 + a6 b2 x y^2 +
  a5 b3 x y^2 + a3 b5 x y^2 + a2 b6 x y^2 + a6 b4 x^2 y^2 +
  a5 b5 x^2 y^2 + a4 b6 x^2 y^2 + a6 b3 y^3 + a3 b6 y^3 +
  a6 b5 x y^3 + a5 b6 x y^3 + a6 b6 y^4, x^_. y^_.]
(*
a1 b1 + a2 b1 x + a1 b2 x + a4 b1 x^2 + a2 b2 x^2 + a1 b4 x^2 +
a4 b2 x^3 + a2 b4 x^3 + a4 b4 x^4 + a3 b1 y +
a1 b3 y + (a5 b1 + a3 b2 + a2 b3 + a1 b5) x y + (a5 b2 + a4 b3 +
    a3 b4 + a2 b5) x^2 y + (a5 b4 + a4 b5) x^3 y + a6 b1 y^2 +
a3 b3 y^2 +
a1 b6 y^2 + (a6 b2 + a5 b3 + a3 b5 + a2 b6) x y^2 + (a6 b4 + a5 b5 +
    a4 b6) x^2 y^2 + a6 b3 y^3 + a3 b6 y^3 + (a6 b5 + a5 b6) x y^3 +
a6 b6 y^4
*)

小木虫Mathematica版块已毁(当然原本也不咋的),建议大家前往百度贴吧或Stackexchange。
4楼2015-04-01 12:15:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)

【答案】应助回帖

引用回帖:
3楼: Originally posted by KGFan at 2015-04-01 11:43:09
非常感谢!再请教一下,有没有办法合并同类项,然后还是按照x^3, x^2*y, x*y^2, y^3这个顺序排列?因为我主要是要得到x^i*y^j这些项的系数表达式。...

仔细看了一下你的问题,于是,有个函数叫做CoefficientList楼主知不知道。
小木虫Mathematica版块已毁(当然原本也不咋的),建议大家前往百度贴吧或Stackexchange。
5楼2015-04-01 12:17:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:23 by CgyNCDVNhGVg
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:15 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 16:26 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-15 15:27 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:50 by VMDqgnlKimtZ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 14:38 by VMDqgnlKimtZ
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-15 14:02 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见