matlab求解一元四次方程区间值,急,f(x)=-x^4+32x^3-208x^2-1436x+122860
来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/24 00:02:40
matlab求解一元四次方程区间值,急,f(x)=-x^4+32x^3-208x^2-1436x+122860
matlab求解一元四次方程区间值,急,
f(x)=-x^4+32x^3-208x^2-1436x+12286
0
matlab求解一元四次方程区间值,急,f(x)=-x^4+32x^3-208x^2-1436x+122860
y=@(x)-((16-x).^2*(576-(24-x).^2)./(24-x).^2);
[x,fval] = fminunc(y,8)
Optimization terminated:relative infinity-norm of gradient less than options.TolFun.
x =
7.3593
fval =
-80.6403
在x=7.359时取得最大值80.64
原式问题的回答
该方程可以使用下面的两种方法求出方程所有的根
可以看出多项式在[0 24]之间没有实根
%by dynamic
%
>> syms x;
>> y=x^4+32*x^3-208*x^2-1436*x+12286;
>> p=sym2poly(y)
p =
1 32 -208 -1436 12286
>> roots(p)
ans =
-36.3775
-7.6853
6.0314 + 2.7511i
6.0314 - 2.7511i
>> root=solve(y)
>> vpa(root)
ans =
6.031361313677396281550706226790+2.7511228762923327275995084867862*i
6.031361313677396281550706226790-2.7511228762923327275995084867862*i
-7.685265104236527554399686795075
-36.377457523118265008701725658505
可用符号方程求解的命令
solve('(16-x)^2*(576-(24-x)^2)/(24-x)^2=0','x')
得到
ans =
[ 16]
[ 16]
[ 0]
[ 48]
可知在0-24区间内的饥解为16
r =
-36.3775
-7.6853
6.0314 + 2.7511i
6.0314 - 2.7511i
>>
所以,x再次范围内无解
x=0:.1:24;
y=-x.^4+32*x.^3-208*x.^2-1436*x+12286;
plot(x,y)
如果想看具体的值,前面两句去掉后面的分号即可.