牛顿法求方程的根matlab程序,matlab 中用牛顿法求方程的根

问题描述：

matlab 中用牛顿法求方程的根

在matlab中,用牛顿法求方程的根的时候,如果所求的方程是一个非多项式方程,即方程里面含有一些含有未知数的特殊函数(例如bessel函数)时,方程不能化简为多项式方程.此时,如果要用牛顿法求解方程的根,要不要经过一些特殊处理?还是就是跟求解普通的方程一样,用基本的牛顿法进行计算就可以了?

如果需要进行特殊处理,

唉,我已经用牛顿法试了很久了,但是答案一直不对,所以我才开始怀疑是不是因为我方程里面的函数太特殊了,用牛顿法有问题.

1个回答

分类：

数学

2014-10-12

问题解答：

我来补答

你是在说符号运算么,亲?

再问： 恩，其实我感觉应该不太算是符号运算。。。因为整个运算，在源程序中是符号运算，但是因为我是计算的物体有具体的尺寸限制，而且我是需要得到具体的数值的，所以我会另外进行一系列的赋值，所以运算到最后的未知数只有一个，会得到一个非多项式的方程(这个方程中含有bessel函数，并且这个函数中也含有我上面说的那个未知数)。

再答： 你不给具体问题我不太敢想啊。。。我觉得你最后想得到一个关系式的话，用用看maple吧，matlab的符号运算我不会啊

再问： 最后的方程如下，这个方程是不是用牛顿法的基本的步骤就可以了？还是要因为bessel函数，要特殊处理？要如何处理？ (2.19e-012*(5.45e+015*w^2-3.35e+020)^(1/2)*besseli(2,(5.98e+004*w^2-3.68e+009)/(8.93e+015*w^2-5.50e+020)^(1/2))- 。。。+1.71e-015*i*w^6*besseli(2,8.57e-012*(5.45e+015*w^2-3.35e+020)^(1/2)))/(5.45e+015*w^2-3.35e+020)^(1/2)/besseli(2,(5.98e+004*w^2-3.68e+009)/(8.93e+015*w^2-5.50e+020)^(1/2))

再答： 我感觉直接牛顿法就行了吧。。。符号运算的牛顿法？你先试试吧

展开全文阅读