matlab中如何求有理分式的商多项式和余多项式

Posted 2023-05-04

1、首先使用root()根函数可以查找出多项式的根，编写好的MATLAB求根程序代码为，将此代码写入到MATLAB命令窗口或M文件中，如下图所示。

2、同理，运行程序代码，可以得出求根结果为如下图所示。

3、如果需要对上述解进行进一步计算，比如求倒数，可以使用poly()函数，MATLAB代码如下。

4、再次运行程序，可以得出如下根处理结果，如下图所示就完成了。

参考技术A

使用多项式除法函数即可：[Q, R]=deconv(Y,X) 

其中，Q是商多项式，R是余数多项式，Y是被除数多项式，X是除数多项式。

示例如下：

人工智能数学基础---不定积分4：有理函数求积分的方法

一、引言

在《人工智能数学基础–不定积分2：利用换元法求不定积分》、《人工智能数学基础—不定积分3：分部积分法》分别介绍了换元积分法和分步积分法。但有些函数表达式很复杂，如果直接用换元积分法和分步积分法不好计算积分，这时需要先对函数进行化简。本文介绍的有理函数求积分就是一种化繁为简的不定积分计算方法。

二、有理函数的概念

2.1、定义

两个多项式的商P(x)/Q(x)称为有理函数，又称为有理分式。

2.2、补充说明

1. 上述定义假设P(x)、Q(x)之间没有公因式，因为有就可以相约去除；
2. 当分子多项式P(x)的次数（变量的最高幂次）小于分母多项式Q(x)的次数，则称有理函数为真分式，否则称为假分式；
3. 利用多项式的除法，总可以将一个假分式化为一个多项式和一个真分式的和；
4. 对于真分式P(x)/Q(x)，如果分母Q(x)可以分解为两个多项式的乘积Q1(x)Q2(x)，即Q(x)=Q1(x)Q2(x)，且Q1(x)、Q2(x)没有公因式，则真分式P(x)/Q(x)可以拆分成两个真分式之和，即：P(x)/Q(x) = P1(x)/Q1(x)+P2(x)/Q2(x)。这个过程称为把真分式化成部分分式之和。如果Q1(x)或Q2(x)还可以分解成两个公因式的多项式的乘积，那么可以再拆分成更简单的部分分式。
5. 经过真分式化成部分分式之和的处理后，最后有理函数的分解式中只会出现多项式、P1(x)/(x-a)^k、P2(x)/(x²+px+q)^l等三类函数（这里p² - 4q<0，P1(x)为小于k次的多项式，P2(x)为小于2l次的多项式）。多项式的积分利用不定积分的加法运算即可以化成每个项的积分，就可以比较容易求出。

思考：学这里时，老猿在想为什么有理函数就可以化为这三种类型的函数呢？仔细想了下，应该是如下理由：

1. 多项式以及后面两个真分式的分子就不用说了，主要是后面两个真分式的分母为什么是那样的；
2. 任何一个变量的n次多项式，其构成起决定作用的是其最高次数n，而任何一个整数n，都可以表示成一个整数和一个偶数的和，因此才有这两种形式的分母。

三、真分式的求解积分方法

先看2个例题：

从以上两个例题可以看到，老猿总结真分式求积分的步骤如下：

1. 首先将分母化解为n个没有公因式的多项式的乘积；
2. 将化解后的有理函数拆分成n个真分式的和，每个真分式的分母为上述n个多项式中的一个，分子为比分母低一次的带未知系数的完整多项式；
3. 根据被积函数与上述n个真分式相等的关系，得到变量的每次系数之间的关系，通过这个关系求出所有未知系数，将其带入到步骤2中的n个真分式中，即将原有理函数求积分化为了n个真分式的和求积分；
4. 利用积分加法运算针对每个真分式求积分，得到最后的积分结果。

四、部分非有理函数化为有理函数求积分

有理数函数求积分的方法，不只是能应用于有理函数，也可以应用于用换元法等方法可化为有理函数的的部分函数。

4.1、三角函数案例

小结：可以看到，通过将三角函数进行u=tan（x/2）的换元，就可以将由三角函数组成的类似有理式的积分计算转换成有理函数的方式来计算，最后结果再将变量x=2arctan u替换回来即可，这种处理方式对三角函数类似有理式的积分都可以应用。

注意：当变量x∈（（2k-1）π，（2k+1）π）时，作变换u=tan（（x-2kπ）/2）=tan(x/2)，x=2kπ+2arctan u，一样可以应用有理分式求积分。

4.2、n次开方案例

案例1：

案例2：

小结：通过上面的案例可以看到，如果被积函数中含有简单根式，可以通过将简单根式设为u，从而将被积函数化为有理分式，用该有理分式求出积分，再用u到x的反变换将结果代换，就可以求得这种带根式的函数积分。

五、小结

本文介绍了有理函数的概念及有理函数求积分的方法，并对于类似有理函数的三角函数形式的被积函数和带根式的被积函数，通过适当地换元变换化为有理函数求积分。

说明：

本文内容是老猿学习同济版高数的总结，有需要原教材电子版以及OpenCV、Python基础知识、、图像处理原理介绍相关电子资料，或对文章内有有疑问咨询的，请扫博客首页左边二维码加微信公号，根据加微信公号后的自动回复操作。

更多人工智能数学基础请参考专栏《人工智能数学基础》。

写博不易，敬请支持：

如果阅读本文于您有所获，敬请点赞、评论、收藏，谢谢大家的支持！

关于老猿的付费专栏

1. 付费专栏《https://blog.csdn.net/laoyuanpython/category_9607725.html 使用PyQt开发图形界面Python应用》专门介绍基于Python的PyQt图形界面开发基础教程，对应文章目录为《 https://blog.csdn.net/LaoYuanPython/article/details/107580932 使用PyQt开发图形界面Python应用专栏目录》；
2. 付费专栏《https://blog.csdn.net/laoyuanpython/category_10232926.html moviepy音视频开发专栏 )详细介绍moviepy音视频剪辑合成处理的类相关方法及使用相关方法进行相关剪辑合成场景的处理，对应文章目录为《https://blog.csdn.net/LaoYuanPython/article/details/107574583 moviepy音视频开发专栏文章目录》；
3. 付费专栏《https://blog.csdn.net/laoyuanpython/category_10581071.html OpenCV-Python初学者疑难问题集》为《https://blog.csdn.net/laoyuanpython/category_9979286.html OpenCV-Python图形图像处理 》的伴生专栏，是笔者对OpenCV-Python图形图像处理学习中遇到的一些问题个人感悟的整合，相关资料基本上都是老猿反复研究的成果，有助于OpenCV-Python初学者比较深入地理解OpenCV，对应文章目录为《https://blog.csdn.net/LaoYuanPython/article/details/109713407 OpenCV-Python初学者疑难问题集专栏目录 》
4. 付费专栏《https://blog.csdn.net/laoyuanpython/category_10762553.html Python爬虫入门 》站在一个互联网前端开发小白的角度介绍爬虫开发应知应会内容，包括爬虫入门的基础知识，以及爬取CSDN文章信息、博主信息、给文章点赞、评论等实战内容。

前两个专栏都适合有一定Python基础但无相关知识的小白读者学习，第三个专栏请大家结合《https://blog.csdn.net/laoyuanpython/category_9979286.html OpenCV-Python图形图像处理 》的学习使用。

对于缺乏Python基础的同仁，可以通过老猿的免费专栏《https://blog.csdn.net/laoyuanpython/category_9831699.html 专栏：Python基础教程目录）从零开始学习Python。

如果有兴趣也愿意支持老猿的读者，欢迎购买付费专栏。

老猿Python，跟老猿学Python！

☞ ░ 前往老猿Python博文目录 https://blog.csdn.net/LaoYuanPython ░