maxima-2 多项式和分式计算

LanAlters

已于 2023-07-15 22:40:47 修改

阅读量160

点赞数

分类专栏： maxima 文章标签：学习

于 2023-07-15 22:34:16 首次发布

本文链接：https://blog.csdn.net/LanAlters/article/details/131745083

版权

maxima 专栏收录该内容

6 篇文章 0 订阅

订阅专栏

文章探讨了在概率论背景下，比较三局两胜和五局三胜比赛中，甲赢的概率，使用wxMaxima软件进行多项式合并与计算。通过expand和factor函数简化表达式，得出在p∈(0.5,1)区间内，五局三胜对甲更有利。此外，还展示了如何使用ratexpand和ratsimp函数对分式进行通分合并，以及factor函数对分式进行因式分解的过程。

摘要由CSDN通过智能技术生成

2 多项式和分式计算

多项式合并

考虑如下的一个推导问题.
甲、乙比赛, 每局比赛甲赢的概率 p \in (0.5, 1).

问:
三局两胜还是五局三胜对甲有利?

三局两胜中甲赢的概率为
$\begin{align*} & P(三局两胜中甲赢) \\ = & P(甲赢 2 局) + P(甲赢 3 局) \\ = & C_{3}^{2}p^{2}(1-p) + p^{3} \\ = & 3p^{2} - 2p^{3} \end{align*}$

五局三胜中甲赢的概率为
$\begin{align*} & P(五局三胜中甲赢) \\ = & P(甲赢 3 局) + P(甲赢 4 局) + P(甲赢 5 局) \\ = & C_{5}^{3}p^{3}(1-p)^{2} + C_{5}^{4}p^{4}(1-p) + p^{5} \\ = & 10p^{3} - 15p^{4} + 6p^{5} \end{align*}$

在 wxMaxima 中用如下两个命令定义了这两个概率:

(%i37) P1: 3*p^2*(1-p)+p^3;
                                3              2
(%o37)                         p  + 3 (1 - p) p
(%i38) P2: 10*p^3*(1-p)^2 + 5*p^4*(1-p) + p^5;
                       5              4             2  3
(%o38)                p  + 5 (1 - p) p  + 10 (1 - p)  p

用 expand 函数可以把多项式合并同类项化简:

(%i39) expand(P1);
                                     2      3
(%o39)                            3 p  - 2 p
(%i40) expand(P2);
                                5       4       3
(%o40)                       6 p  - 15 p  + 10 p

如下命令计算两个多项式的差并合并同类项:

(%i41) P3: expand(P2 - P1);
                             5       4       3      2
(%o41)                    6 p  - 15 p  + 12 p  - 3 p

为判断其是否在 p \in (0.5, 1) 总为正值, 可以分解因式.

Maxima 用 factor 函数分解因式.

(%i42) factor(P3);
                                     2  2
(%o42)                      3 (p - 1)  p  (2 p - 1)

在 p \in (0.5, 1) 为正值.

推导结果导出

wxMaxima 的结果显示为数学公式形式, 选中结果, 可以选"编辑—复制为 LaTeX"菜单, 把结果导出到 LATEX 文件中.

也可以把公式复制为纯文本、图片.

分式通分合并

考虑如下分式:

(%i1) R1:(x+1)^2/(x-1) + 1/(x+1);
                                      2
                               (x + 1)      1
(%o1)                          -------- + -----
                                x - 1     x + 1

expand 函数可以分别展开分子分母中的多项式, 但不能通分合并:

(%i2) expand(R1);
                                   2
                           1      x       2 x      1
(%o2)                    ----- + ----- + ----- + -----
                         x + 1   x - 1   x - 1   x - 1

函数 ratexpand 可以把分式的分子分母中多项式合并同类项, 然后通分, 结果表示为分子每个不同幂次单独一个分式的加法.

(%i3) ratexpand(R1);
                              3         2
                             x       3 x      4 x
(%o3)                      ------ + ------ + ------
                            2        2        2
                           x  - 1   x  - 1   x  - 1

函数 ratsimp 把分式之和通分合并变成一个分式.

(%i4) ratsimp(R1);
                                 3      2
                                x  + 3 x  + 4 x
(%o4)                           ---------------
                                     2
                                    x  - 1

factor 函数可以对分式的分子分母分别进行因式分解.

(%i5) factor(R1);
                                   2
                               x (x  + 3 x + 4)
(%o5)                          ----------------
                               (x - 1) (x + 1)