[求助]符号计算的菜鸟问题，请教！

shenhongyu · 发表于 2005-8-21 11:26

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

小弟想求一个一元高次方程，但是用符号方程直接解有错误提示，用roots求解误差太大，所以我想用符号做系数，做解析解，再带入数值，这样应该会准确！ 我的问题是一次如何定义很多符号，然后如何让它转化为方程： 例：已知方程的系数， 1.定义符号：a1,a2,a3,........a70； 2.生成方程：equ=a70*x^70+a69*x^69+........+a2*x^2+a1*x 解方程：solve(equ)； 3.之后带回系数：a1,a2,a3,........a70的值求出解析解。 小弟不知怎么实现，不能让我一个一个的输入吧，有什么命令实现上边的三个步骤？？

AaronSpark · 发表于 2005-8-21 12:07

定义矩阵就行了，比如 a=sym('[a1,a2,a3,.....,a70]'); syms x equ=sym('a1*x'); for i=2:70 equ=equ+a(i)*x^i; end solve(equ)

[此贴子已经被作者于2005-8-21 12:19:48编辑过]

hyacinth · 发表于 2005-8-22 14:54

<DIV class=quote>以下是引用shenhongyu在2005-8-21 11:26:11的发言： 小弟想求一个一元高次方程，但是用符号方程直接解有错误提示，用roots求解误差太大，所以我想用符号做系数，做解析解，再带入数值，这样应该会准确！ 我的问题是一次如何定义很多符号，然后如何让它转化为方程： 例：已知方程的系数， 1.定义符号：a1,a2,a3,........a70； 2.生成方程：equ=a70*x^70+a69*x^69+........+a2*x^2+a1*x 解方程：solve(equ)； 3.之后带回系数：a1,a2,a3,........a70的值求出解析解。 小弟不知怎么实现，不能让我一个一个的输入吧，有什么命令实现上边的三个步骤？？</DIV>
除非a1...a70是特殊的表达式，要不基本没有太大意义 我算了一个四次方的，其结果就是NN长了，何况是70次的

zjliu · 发表于 2005-9-21 15:59

应该直接说问题,举例子如果不恰当会累人的

账号		自动登录	找回密码
密码			我要加入

[编程技巧] [求助]符号计算的菜鸟问题，请教！

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（shenhongyu）[求助]符号计算的菜鸟问题，请教...