高分求一段简单的mathematica程序,计算一个高阶多项式的解
来源:学生作业帮 编辑:神马作文网作业帮 分类:数学作业 时间:2024/11/10 13:36:48
高分求一段简单的mathematica程序,计算一个高阶多项式的解
多项式形式为ax^10+bx^9+cx^8+dx^7+ex^6+fx^5+gx^4+hx^3+ix^2+jx^1+k
其中a,b,c,d,f,等为可输入的参数,对这个方程求二阶导数,并让这个二阶导数等于0,计算x的值,.毕业论文计算数据要用.导师说这是个很简单的问题,分数不是问题.
多项式形式为ax^10+bx^9+cx^8+dx^7+ex^6+fx^5+gx^4+hx^3+ix^2+jx^1+k
其中a,b,c,d,f,等为可输入的参数,对这个方程求二阶导数,并让这个二阶导数等于0,计算x的值,.毕业论文计算数据要用.导师说这是个很简单的问题,分数不是问题.
Mathematica代码如下:
co ={a, b, c, d, e, f, g, h, i, j, k};
expr = co.x^Range[10, 0, -1]
NSolve[D[expr, {x, 2}] == 0, x]
因为方程次数高于5,所以一般只有数值解.
第一句是给系数赋初值,对应x次数由高到低(一定要先赋值,否则没法计算).
第二句是产生表达式,这里用了矢量的标量积来缩减表达式.
第三句是求导,求数值解.
作为例子,比如:
co = Range[11];
expr = co.x^Range[10, 0, -1]
NSolve[D[expr, {x, 2}] == 0, x]
co ={a, b, c, d, e, f, g, h, i, j, k};
expr = co.x^Range[10, 0, -1]
NSolve[D[expr, {x, 2}] == 0, x]
因为方程次数高于5,所以一般只有数值解.
第一句是给系数赋初值,对应x次数由高到低(一定要先赋值,否则没法计算).
第二句是产生表达式,这里用了矢量的标量积来缩减表达式.
第三句是求导,求数值解.
作为例子,比如:
co = Range[11];
expr = co.x^Range[10, 0, -1]
NSolve[D[expr, {x, 2}] == 0, x]