如何用matlab解一元三次方程

0=x^3-6*x^2+9*x-9谁能给出详细的过程和结果?Undefinedfunctionormethod'solve'forinputargumentsoftype... 0=x^3-6*x^2+9*x-9
谁能给出详细的过程和结果?
Undefined function or method 'solve' for input arguments of type 'char'
展开
 我来答
心众眠5
推荐于2019-11-07 · TA获得超过7.6万个赞
知道答主
回答量:60
采纳率:100%
帮助的人:4.9万
展开全部

把你方程的系数降幂排列,写成一个向量,利用roots命令即可。

>> p=[1,-6,9,-9];roots(p)

ans =

4.4260e+000

7.8701e-001 +1.1891e+000i

7.8701e-001 -1.1891e+000i

扩展资料

MATLAB线性方程求解

a为m*n的系数矩阵。

Ø    m=n,该方程为“恰定”方程

Ø    m>n,该方程为“超定”方程

Ø    m<n,该方程为“欠定”方程

(1)恰定方程求解

方程ax+b(a为非奇异)

x=a-1 b

两种解:

①x=inv(a)*b 采用求逆运算解方程; 

②x=a\b    采用左除运算解方程。     

(2)超定方程求解

方程 ax=b ,m>n时此时不存在唯一解。

方程解 (a'a)x=a'b  

①x=(a'a)-1 a'b ——求逆法;       

②x=a\b  matlab用最小二乘法找一个准确地基本解。    

(3)欠定方程求解

当方程数少于未知量个数时,即不定情况,有无穷多个解存在。

matlab可求出两个解:

①用除法求的解x是具有最多零元素的解;

②是具有最小长度或范数的解,这个解是基于伪逆pinv求得的。

富港检测技术(东莞)有限公司_
2024-04-02 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大(共振点);正弦振动在任一瞬间只包含一种频率的振动,而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率,所以样品上的共振点会同时激发... 点击进入详情页
本回答由富港检测技术(东莞)有限公司_提供
ilovewenping
2009-09-30 · TA获得超过3.5万个赞
知道大有可为答主
回答量:3380
采纳率:100%
帮助的人:2090万
展开全部
谁能编程用MATLAB解一元三次方程
能帮我用MATLAB编个程序解个一元三次方程么?
方程为想x^3+a*x^2-b=0
我知道用solve很容易解 但是我是想能用一个程序一运行就是输入a,b的值,然后直接出结果。不用每次都自己再输入SOLVE(x^3+a*x^2-b)
我想很简单 就是自己还刚接触还不会 希望大家给帮个忙,
问题补充:这怎么会很难 ,只要你定义一下变量,在用输入语句,在把定义的输入变量赋值给 A,B就可以了 你可别来骗我 ,我就是想让大家来帮忙,
要是你感觉很难 那就不用你帮了 ,我可怕浪费你时间。
]可是我不会 ,因为我学的是C

在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function x=mysolve3(a,b)
%%求解x^3+a*x^2-b=0的解
p=(-a*a)/3;
q=(2*a^3+27*(-b))/27;
D=(p/3)^3+(q/2)^2;
if (D>=0)
temp=-q/2+sqrt(D);
if (temp>0) u=(temp)^(1/3); else u=-(-temp)^(1/3); end
temp=-q/2-sqrt(D);
if (temp>0) v=(temp)^(1/3); else v=-(-temp)^(1/3); end
y(1)=u+v;
y(2)=-(u+v)/2+i*sqrt(3)*(u-v)/2;
y(3)=-(u+v)/2-i*sqrt(3)*(u-v)/2;
else
temp=acos(-q/2*(-p/3)^-1.5);
y(1)=2*sqrt(-p/3)*cos(temp/3);
y(2)=-2*sqrt(-p/3)*cos((temp+pi)/3);
y(3)=-2*sqrt(-p/3)*cos((temp-pi)/3);
end
x=y-a/3;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

在Matlab下面输入:
mysolve3(1,3)
得到三个解

你可以通过vpa(solve('x^3+x^2-3'),5)进行检验

参考资料: http://zhidao.baidu.com/question/59314141.html

已赞过 已踩过<
你对这个回答的评价是?
评论 收起
was_ist_das
推荐于2017-09-03 · TA获得超过6897个赞
知道大有可为答主
回答量:1459
采纳率:0%
帮助的人:732万
展开全部
>> solve('x^3-6*x^2+9*x-9')

ans =

1/(45^(1/2)/2 + 7/2)^(1/3) + (45^(1/2)/2 + 7/2)^(1/3) + 2
2 - (45^(1/2)/2 + 7/2)^(1/3)/2 - (3^(1/2)*i*(1/(45^(1/2)/2 + 7/2)^(1/3) - (1/2*45^(1/2) + 7/2)^(1/3)))/2 - 1/(2*(45^(1/2)/2 + 7/2)^(1/3))
(3^(1/2)*i*(1/(45^(1/2)/2 + 7/2)^(1/3) - (1/2*45^(1/2) + 7/2)^(1/3)))/2 - (45^(1/2)/2 + 7/2)^(1/3)/2 - 1/(2*(45^(1/2)/2 + 7/2)^(1/3)) + 2

>> roots([1 -6 9 -9])

ans =

4.4260
0.7870 + 1.1891i
0.7870 - 1.1891i
本回答被提问者采纳
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
百度网友53799a5ac9
2009-09-30 · TA获得超过3556个赞
知道小有建树答主
回答量:1063
采纳率:0%
帮助的人:1094万
展开全部
把你方程的系数降幂排列,写成一个向量,利用roots命令即可。
>> p=[1,-6,9,-9];roots(p)

ans =

4.4260e+000
7.8701e-001 +1.1891e+000i
7.8701e-001 -1.1891e+000i
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
greatdju
2009-09-30 · TA获得超过2.2万个赞
知道大有可为答主
回答量:2486
采纳率:50%
帮助的人:4162万
展开全部
%数值解
roots([1 -6 9 -9])

%解析解
syms x
solve(x^3-6*x^2+9*x-9)
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
收起 更多回答(4)
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式