matlab求多项式的实数根_matlab中怎么求多项式的根
2017-10-09 回答
用matlab解方程的三个实例
1、对于多项式p(x)=x3-6x2-72x-27,求多项式p(x)=0的根,可用多项式求根函数roots(p),其中p为多项式系数向量,即
>>p =
p =
1.00 -6.00 -72.00 -27.00
p是多项式的matlab描述方法,我们可用poly2str(p,'x')函数 ,来显示多项式的形式:
>>px=poly2str(p,'x')
px =x^3 - 6 x^2 - 72 x - 27
多项式的根解法如下:
>> format rat %以有理数显示
>> r=roots(p)
r =
2170/179
-648/113
-769/1980
2、在matlab中,求解用符号表达式表示的代数方程可由函数solve实现,其调用格式为:solve(s,v):求解符号表达式s的代数方程,求解变量为v。
例如,求方程(x+2)x=2的解,解法如下:
>> x=solve('(x+2)^x=2','x')
x =
.69829942170241042826920133106081
得到符号解,具有缺省精度。如果需要指定精度的解,则:
>> x=vpa(x,3)
x =
.698
3、使用fzero或fsolve函数 ,可以求解指定位置(如x0)的一个根,格式为:x=fzero(fun ,x0)或x=fsolve(fun,x0)。例如,求方程0.8x+atan(x)-=0在x0=2附近一个根,解法如下:
>> fu=@(x)0.8*x+atan(x)-pi;
>> x=fzero(fu,2)
x =
2.4482
或
>> x=fsolve('0.8*x+atan(x)-pi',2)
x =
2.4482
________________________________________
当然了,对于该方程也可以用第二种方法求解:
>> x=solve('0.8*x+atan(x)-pi','x')
x =
2.4482183943587910343011460497668
对于第一个例子,也可以用第三种方法求解:
>> f=@(x)x^3-6*x^2-72*x-27
f =
@(x)x^3-6*x^2-72*x-27
>> x=fzero(f,10)
x =
12.1229
对于第二个例子,也可以用第三种方法:
>> fun=@(x)(x+2)^x-2
fun =
@(x)(x+2)^x-2
>> x=fzero(fun,1)
x =
0.6983
最近有多人问如何用matlab解方程组的问题,其实在matlab中解方程组还是很方便的,例如,对于代数方程组ax=b(a为系数矩阵 ,非奇异)的求解,matlab中有两种方法:
(1)x=inv(a)*b — 采用求逆运算解方程组;
(2)x=a\b — 采用左除运算解方程组。
例:
x1+2x2=8
2x1+3x2=13
>>a=;b=;
>>x=inv(a)*b
x =
2.00
3.00
>>x=a\b
x =
2.00
3.00;
即二元一次方程组的解x1和x2分别是2和3。
对于同学问到的用matlab 解多次的方程组,有符号解法,方法是:先解出符号解,然后用vpa(f,n)求出n位有效数字的数值解.具体步骤如下:
第一步:定义变量syms x y z ...;
第二步:求解=solve('eqn1','eqn2',...,'eqnn','var1','var2',...'varn');
第三步:求出n位有效数字的数值解x=vpa(x,n);y=vpa(y,n);z=vpa(z,n);...。
如:解二(多)元二(高)次方程组:
x^2+3*y+1=0
y^2+4*x+1=0
解法如下:
>>syms x y;
>>=solve('x^2+3*y+1=0','y^2+4*x+1=0');
>>x=vpa(x,4);
>>y=vpa(y,4);
结果是:
x =
1.635+3.029*i
1.635-3.029*i
-.283
-2.987
y =
1.834-3.301*i
1.834+3.301*i
-.3600
-3.307。
matlab求多项式的实数根_matlab中怎么求多项式的根相关推荐
- matlab计算一组数的函数值_matlab中怎么求一组数据概率最大的数
好象直接使用就可以了,第一次出现就包含了定义 顺便把些命令弄给你看下 matlab命令 matlab commands and functions list A a abs 绝对值.模.字符的ASCI ...
- cms 移动根节点到根节点中的子节点,根节点消失解决办法
cms 移动根节点到根节点中的子节点,根节点消失解决办法 cms使用对应有个NetCMS数据库,找到其中NT_News_Class表,这个表就是保存栏目名称的表 select * from dbo.N ...
- matlab多项式的拟合与插值例题_Matlab中数据处理和多项式插值与曲线拟合
一. 基本统计处理 1.查取最大值 MAX函数的命令格式有: [Y,I]= max (X):将max(X)返回矩阵X的各列中的最大元素值及其该元素的位置赋予行向量Y与I:当X为向量时,则Y与I为单变 ...
- matlab 邻近度 离群点_Matlab 学习记录帖 —— 多项式、插值和数据拟合
关于多项式Matlab命令 一个多项式的幂级数形式可表示为: 也可表示为嵌套形式: 或因子形式 幂系数:在matlab里,多项式用行向量表示,其元素未多项式的系数,并从左至右按降幂排列. Roots: ...
- matlab如何求传递函数的幅值_MATLAB中求开环传递函数的幅值裕度、相位裕度、截止频率的margin()函数用法...
一.幅值裕度.相位裕度 相位裕度γ就是Bode图中幅值为0dB时对应的相位加上180°,如图中0dB时的相位是-99.1°,所以该系统的相位裕度就是 -99.1°+180°=80.9°.所谓相位裕度的 ...
- matlab非同秩矩阵相乘_MATLAB中的矩阵与向量运算
4.1 数组运算和矩阵运算 从外观形状和数据结构来看 , 二维数组和数学中的矩阵没有区别 . 但是 , 矩阵作为一种变换或 映射算符的体现 , 矩阵运算有着明确而严格的数学规则 . 而数组运算是 MA ...
- matlab如何求传递函数的幅值_MATLAB通过传递函数求截止频率
展开全部 可以用MATLAB画伯德图. 有两个函数可以画伯德图 一个是bode函数 格式是[mag,phase,w]=bode(G); G是构建好的系统,62616964757a686964616 ...
- matlab分块矩阵取某一块_matlab中怎么取一个矩阵的一部分
展开全部 matlab中要取一个62616964757a686964616fe4b893e5b19e31333366306562矩阵的一部分,可以参考以下内容: 比如说有一个5×5的矩阵,想取它右下脚 ...
- 求两个数的公约数java_java中怎样求两个数的最大公约数?
方法一:(辗转相除法) 设用户输入的两个整数为n1和n2且n1>n2,余数=n1%n2.当余数不为0时,把除数赋给n1做被除数,把余数赋给n2做除数再求得新余数,若还不为0再重复知道余数为0,此 ...
最新文章
- 2019计算与系统神经科学大会Cosyne 前沿研究汇总
- 错误:”未能加载文件或程序集“System.Web.Mvc, Version=2.0.0.0” 解决方法
- echarts实现双Y轴之散点和折线图
- 刚刚,“国民”APP微信崩了!官方致歉:已经逐步恢复
- 自动化测试验证码代码常用的四种方式
- html5文档加载前调用函数,html调用javascript外部文件显示函数未定义
- Aspose.Cells 使用UnMerge()取消合并单元格(取消合并单行,取消合并单列,取消多行多列)
- 基于 HTML5 + WebGL 的太阳系 3D 展示系统
- SAP中常用到的会计知识
- Android Studio之Activity简介
- 带你快速入门AXI4总线--汇总篇
- X280 工业4g路由器
- c语言校车管理系统,交运智慧驾培管理系统
- python 裁剪图片
- PMP考试六大管理学定律
- 开源为改进服务器、网络、存储及加速技术提供了行之有效的新途径
- 必测的支付漏洞(三)服务器未返回支付结果前进行干预
- 【lua】-- lua编译
- 【AI 简报20201018期】英伟达开源「Imaginaire」、智能手表争夺战是如何打响的?
- 吞食天地2刘蜀传攻略