代数方程和数值计算及复杂性理论简介.pptVIP

顾 小 丰 Email：guxf@uestc.edu.cn * §1.1　代数方程的不动点迭代算法 我们知道，形如 ax2＋bx＋c=0，a≠0 的一元二次方程，它的两个根可以按照 一元三次方程的求根公式 一元四次方程的求根公式 x4＋ax3＋bx2＋cx＋d=0 可以配方成 高 次 方 程 上述这种把代数方程的根用方程系数经有限次加、减、乘、除和开方运算表示出来的方法，叫做代数方程的代数解法（或公式解法）。 但是，数学家已经证明，五次和更高次方程，就找不到普遍适用的代数解法，这就是说，不会有用方程系数经有限次加、减、乘、除和开方运算把方程的根表示出来的公式。这种“无公式解”的本性是和五次以下的方程不同的。由于这个原因，以后我们只把五次和高于五次的代数方程叫做高次方程。 高次方程虽然没有普遍适用的代数解法，但是却有一些非代数的或者说非公式的解法。下面先介绍高次方程的不动点迭代解法。 不动点迭代法 代数方程都可以表示成 f(x)=a0xn＋a1xn-1＋a2xn-2＋···＋an-1x＋an=0，a0≠0 这里f(x)是一个n次多项式。如果能够把方程 f(x)=0 改写成 x=φ(x) 的形式，并能够找到一个x*，使得 x*=φ(x*) 那么，x*就是原代数方程的一个解。 不动点迭代法 把方程f(x)=0改写成x=φ(x)的形式，非常容易，也有许多方式。例如，可以写成 x=-a0xn－a1xn-1－a2xn-2－···－(an-1－1)x－an， 也可以写成 不动点迭代法 从上面可以看出，把代数方程改写成不动点方程是容易的，难的是怎样得到不动点x*。为此，我们采用迭代方法：找一个点，记作x0，代入函数φ，得到φ(x0)，记作x1，再代入函数φ，得到φ(x1)，记作x2，···，如此一直做下去，可以得到一个序列 x0, x1, x2, ···, xn, ··· 其迭代关系可以表示成 xn+1=φ(xn)，n=0,1,2,··· 有趣的是，这个迭代序列有时候可以帮助我们找到所要的不动点，这就是不动点迭代方法。 不动点迭代法例1 考虑5次方程 x5－17x＋2=0 首先把它变成不动点方程 不动点迭代法例1 刚才，我们选x0=0开始迭代，获得成功，这是不是巧合？是不是接受了什么暗示？提出怀疑是完全合理的，应当多做几次试验。下面分别从x0=1，x0=-1，x0=2，x0=-2，开始迭代，4个迭代序列如下： 不动点迭代法例1 到目前为止，5个迭代都是成功的，一共找到2个解。下面，再扩大范围试试，从x0=3和x0=-3开始迭代，数据如下： 例1的图示 不动点迭代方法非常简单，但不一定能够保证成功。迭代是否成功，怎样使迭代成功，我们等以会儿再讨论。为了进一步把上述迭代的情况研究清楚，可以画一张迭代图帮助我们分析。 例1的图启示 从-3开始的迭代发散，从-2开始的迭代收敛。-3和-2之间，应当有一个分界点，分界点在哪里？从分界点开始的迭代，究竟怎样进行？ 从1开始的迭代收敛到0.1176483这个不动点，从2开始的迭代驻守在2这个不动点，它们之间也应当有一个分界点，也有同样的问题。 2和3之间也有同样的问题。 两个结论 点2是个孤立的、很不稳定的不动点，迭代出发点x0与点2差一点点，迭代的结果就完全不同 问题(1)的分界点在-2.0590与-2.0589之间 二分法 我们要解的是方程f(x)=0，如果我们已经知道两点a＜b，使得f(a)·f(b)＜0，即f(a)和f(b)符号相反，那么在(a,b)中取中点c，计算f(c)。如果f(c)=0，则解已求得，不必再迭代下去。否则，如果f(a)·f(c)＜0，知f(c)与f(b)同号，就扔掉原来的b，把c作为新b，仍如上法迭代；如果f(b)·f(c)＜0，知f(c)与f(a)同号，就扔掉原来的a，把c作为新a，仍如上法迭代。这就是同符号顶替的原则。这样，每迭代一次，区间(a,b)的长度就缩小一半，而在区间的两端，函数f(x)的符号总是相反的。于是，根据f(x)的连续性，这些每次缩短一半的区间最后套住一个点，这个点一定是f(x)=0的解。 二分法求解 现在就来试一试。前已知道，-2.0590与-2.0589之间应当有一个分界点，我们就拿a=-2.0590和b=-2.0589开始。f(a)= -3.817×10-3＜0，f(b)=3.469×10-3＞0，正好符合f(a)·f(b)＜0，取(a,b)中取中点c=-2.05895，计算f(c)，这时不必记录具体结果，只要知道正负就可以了，算得f(c)＜0，与f(a)同号，所以按照同号顶替得原则，取c作为新a，下次迭代就取(a,b