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

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

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

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

<P>定义矩阵就行了，比如<br><br>a=sym('');<br>syms x<br><br><br>equ=sym('a1*x');<br>for i=2:70<br>equ=equ+a(i)*x^i;<br>end<br><br>solve(equ)</P>
[此贴子已经被作者于2005-8-21 12:19:48编辑过]

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

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

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

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

应该直接说问题,举例子如果不恰当会累人的
页: [1]
查看完整版本: [求助]符号计算的菜鸟问题，请教！