科学网—标签

相关帖子	版块	作者	回复/查看	最后发表

GrandFT 2019-11-3 21:37

题目：神经网络求解方程主讲人：赵宇霖时间：2019年11月4号上午10：00 地点：天津大学新校区32教302 提纲：1.神经网络基础-多层感知机 2.梯度下降与函数拟合 3.神经网络的应用于求解方程参考索引：周池春师兄提供，mlp示例python程序代码，与相关pdf简介多层感知机(MLP)网页教程“ http://www.weizn.net/m/?post=253 “ 维基百科词条“Gradientdescent” NeuralNetworkMethodforSolvingPartialDifferentialequation，LucieP.aarts

个人分类: 专题讨论班|2303 次阅读|0 个评论

[转载]【MATLAB】解方程相关命令

JerryYe 2017-5-25 20:59

1、求方程X2+5X-10=0的解。 X=solve('X^2+5*X-10') X=double(X) X = -6.53112887414927 1.53112887414927 2、求方程X2+5X+10=0的解。 X=solve('X^2+5*X+10') X=double(X) X = -2.5 - 1.93649167310371i -2.5 + 1.93649167310371i 3.求解微分方程 https://jingyan.baidu.com/article/e52e36154448e940c60c51aa.html

个人分类: Matlab|4697 次阅读|0 个评论

解方程

可变系时空多线矢主人 2016-1-9 17:33

解方程中国科学院力学研究所吴中祥提要具体分析“伽罗华理论”、创新解法，解得任意 1 元 n 次不可约代数方程仅引进 2 次根式，就能，也才能，给出认识、运用各种事物客观特性和运动规律所需的各种互不相依方程组，仅由实数、虚数或复数表达的解。关键词：多元互不相依方程组，伽罗华理论， 1 元 n 次不可约代数方程，仅由实数、虚数或复数表达的解 1 ．认识、运用各种事物客观特性和运动规律所需的各种互不相依方程组及其解各种事物相互作用的客观规律形成各种代数或解析的多元互不相依方程组，它们及其解是确切认识、运用各种事物特性和规律的关键。例如，分别由 4 种时空自然力与时空运动力建立的方程，并由它们的初始条件和边界条件，确定其相应的宏观运动轨迹，特别是强力和弱力，都是高于 12 维的，多元、高次，互不相依的，偏微分方程组。各种多元、高次，互不相依的，微分或偏微分方程组都可由相应的变换，成为相应的一元高次不可约代数方程。都须求得它们的公式解。早在公元前 3 世纪，引进 2 次的根式，就已得出 2 次不可约代数方程的根式解。直到公元 16 世纪，才先后引进 2 次和 3 次的根式，得到 3 次和 4 次不可约代数方程的根式解。而 1 元 5 次以上的不可约代数方程却历经数学家 5 百多年的努力，仍未能解得，并主要由“伽罗华理论”，分析各根式群的特点，而给出“代数方程能够求得根式解的判据”之后，阿贝尔，据此，首先提出 n4 的不可约代数方程不能根式求解，学术界就似已公认 n4 的不可约代数方程没有根式解。就只能在具体分析其各“解” 所在数域的基础上，数值地逼近，或引入某些特殊函数求解，不能有公式解。这当然就给许多实际问题和理论工作，因没有相应方程的公式解，而造成许多限制、不便和困难。 2 ．具体分析“伽罗华理论”，创建新的解法，求得 1 元任意 n 次不可约代数方程仅由实数、虚数，或复数表达的公式解。具体分析伽罗华理论，确可证明方程根式解的可解性是相应于将方程各系数作有理运算与逐次添加相应根式的变换群的可解性，而这种变换群的阶数等于其整个求解过程中添加根式的最大指数，而当这种变换群的阶数 4 时，其对称置换群，及其子群，就都是非交换群的单群，就都是不可解的。伽罗华所证明的，实际上，只是“在求解 n 次不可约代数方程的整个过程中，所添加根式的指数， n* ，应是小于 4 ”，并非所解方程的次数， n ，应是小于 4 ，并非方程的次数 n 大于 4 就不能有根式解。显然，其整个求解过程中添加根式的最大指数， n* ，并非所解方程的次数 n ，按伽罗华理论，完全得不出其整个求解过程中添加根式的最大指数， n* ，等于所解方程的次数 n ，或两者有任何关系的根据。阿贝尔也未能给出 n4 的不可约代数方程就没有根式解的任何根据。因而，按伽罗华理论，认为“ n4 的不可约代数方程没有根式解”；本身就是混淆 n* 与 n ，而得出的错误结论。而且：一般而言，引进 j 次根式，就会有任意负实数， -s ，的 j 次根式， (-1)^(1/j)s^(1/j), 而产生以 (-1)^(1/j) 所标志的各种与实数不同的数类。当 j=2 ， (-1)^(1/2) 就被定义为由“ i ”标志的所谓“虚数”，而可以与实数一起，成为复数。当 j 是大于 2 的其它数，则不能确定它们是什么数。这就表明：方程的解，应仅由实数、虚数，或复数，表达，就不能出现 j 大于 2 的 (-1)^(1/j) ，就不仅不能有大于 4 的根式，甚至是不能有大于 2 次的根式。现有 3 次和 4 次不可约代数方程的解，都含有 3 次的根式，就实际上，并非确切的解。本文，创建如下新的解法，仅逐次引进 2 次根式，不仅不违反伽罗华理论，也不致产生 -1 的高于 2 次的根式，而能仅以实数、虚数或复数表达，全面、具体给出 1 元任意 n 次不可约代数方程的确切的公式解。 3 ．任意 1 元 n 次不可约代数方程的创新解法 (3,1) 2 次不可约代数方程的创新解法 2 次不可约代数方程： x^2+a1x+a0=0 ，都可由变换 y=x+a1/2 ， x=y-a1/2 ， x^2=y^2-a1y+a1^2/4, 而使原方程变换为： y^2-a1y+a1^2/4+a1(y-a1/2)+a0=0 ， 1 次项的系数 =0 的如下形式： y^2+b0=0, b0=a0-a1^2/4, 由此，解得： y1=+i(b0)^(1/2), (1’) y2=-i(b0)^(1/2), (2’) x1=-a1/2+((a1/2)^2-a0)^(1/2), (1*) x2=-a1/2-((a1/2)^2-a0)^(1/2) ， (2*) 当 (a1/2)^2-a0=0 ， 2 解为相同的实数： x1= x2=-a1/2, (1*’) 当 (a1/2)^2-a00 ， 2 解为不同的实数： x1=-a1/2+((a1/2)^2-a0)^(1/2), (1*”) x2=-a1/2-((a1/2)^2-a0)^(1/2) ， (2* ” ) 当 (a1/2)^2-a00 ， 2 解为互共轭的复数： x1=-a1/2+i((a1/2)^2-a0)^(1/2), (1*”’) x2=-a1/2-i((a1/2)^2-a0)^(1/2) ， (2* ”’ ) 即得任意 2 次不可约代数方程，仅由实数、虚数或复数表达的公式解。 (3,2) 2m+1 次不可约代数方程的创新解法对于 2m+1 次方程也总可变换变量而成为： x^(2m+1)+a x^(2m-1)+a x^(2m-2)+ … +a x+ a =0 ，有 2m+1 个根： x 、 x 、…、 x ，也还都可以表达为： (x^(2m)+a’ x^(2m-1)+…+a’ x+a’ )(x-x )=0, 由此式与原方程对比，得到的各关系式，即可解得：各 a’ ； j=0,1, … 2m-1 均为 a ； j=0,1, … 2m-1 ，的函数。只要解得： x^(2m)+a’ x^(2m-1)+ … +a’ x+a’ =0 ，都是实数、虚数或复数的解，就得到 2m+1 次方程各解中的 2m 个仅由实数、虚数或复数表达的解。其中，都含有 x ，可利用 x 满足原方程的条件，对它逐次降幂，而消去。方程其它各根仅由各系数表达的方程，又可利用各根满足原方程的条件，对各根逐次降幂，而解得各根。再利用原 2m+1 次不可约代数方程各根之和 =0 ，确定 x ，即得任意 2m+1 次不可约代数方程不含任何根式的 2m+1 个解。 (3,3) 2m 次不可约代数方程的创新解法对于任意 2m 次不可约代数方程都总可表达为： x^(2m)+a x^(2m-2)+ … +a x+a =0 ，当 m1 ，其多项式还都可改写为如下 2 个 m 次不可约代数方程的多项式乘积的形式， (x^m+a’ x^(m-1)+ … +a’ x+a’ ) (x^(m)+a” x^(m-1)+ … +a” x+a” )=0 ，由此都可解得：各 a” 、 a’ ； j=0,1, … m-2 为 a ； j=0,1, … 2m-2 的函数。即可，分别解得 a’ 、 a” ； j=0,1, … m-2 的 2 个 m 次方程仅由实数、虚数或复数表达的， 2m 次不可约代数方程的根式解。由于已经解得任意 2 次不可约代数方程，仅由实数、虚数或复数表达的公式解。按此法，就逐次具体给出了 1 元 3 次、 4 次， 5 次、 6 次，直到 2m 次、 2m+1 次，仅由实数、虚数或复数表达的公式解。 4 ． 1 元 3 次不可约代数方程的创新解法任意 3 次不可约代数方程都可以表达为如下形式： y^3+b1y+b0=0, 有 3 个根 y0 ， y1 ， y2 ，且各根与各系数有如下关系式： -(y0+y1+y2)=0, (1) y0(y1+y2)+y1y2=b1, (2) –y0y1y2=b0, (3) 还可表达为： (y^2+c1y+c0)(y-y0)=0 ，再与原方程对比，有： c1=y0 ， c0=-b0/y0 ，而由 y^2+c1y+c0=0 ，即： y^2+y0y- b0/y0=0 ，解得： y1=-y0/2+i((y0/2)^2+b0/y0)^(1/2), (4) y2=-y0/2-i((y0/2)^2+b0/y0)^(1/2) ， (5) 由 (4) ，有： (y1+y0/2)^2+(y0/2)^2+b0/y0=0 ，即： y0^3+2y1y0^2+2y1^2y0+2b0=0 ， (4’) y0 是原方程的根，有： y0^3+b1y0+b0=0, (6) 由 (4’) 与 (6) ，对 y0 逐次降幂，直到消去 y0 ，即得： 16b0y1^9-16b1^2y1^8-16b1b0y1^7-8b0y1^6+32b1^2b0y1^5 -(4b1^4+20b1b0^2-4b1)y1^4+(2b1^3b0+2b0^3+4b0) y1^3 -12b1^2b0^2y1^2+(2b1^4b0+6b1b0^3-6b1b0)y1 -(b1^3(b0^2-1)+b0^5-b0^2)=0 ， (7) y1 也是原方程的根，也有： y1^3+b1 y1+b0=0, (8) 由 (7) 与 (8) ，对 y1 逐次降幂，直到 1 次，即得： y1=b1((400b1^10b0-b1^8(3300b0^2+20)+b1^6(6696b0^3+72b0) +b1^5(300b0^4+160b0^3-60b0^2) -b1^3(1956b0^5+590b0^4-282b0^3+16b0^2-6b0) -b1(14400b0^6+960b0^5+48b0^4) -76b0^7+32b0^6+126b0^5+48b0^4+18b0^3) /(b0^2(400b1^10-2880b1^8b0+5184b1^6b0^2 +b1^3(2260b0^3+160b0^2-160b0) -b1^2(25484b0^4+2664b0^3+652b0^2+72b0+20) +2296b0^6+1112b0^5+1786b0^4+56b0^3+48b0^2)) ，同样，由 (5) ，有： (y2+y0/2)^2+(y0/2)^2+b0/y0=0 ，即： y0^3+2y2y0^2+2y2^2y0+2b0=0 ， (5’) 仅将 y1 换为 y2 ，就如与 y1 同样地，得到： y2=y1 ，由 (1) ： y0+y1+y2=0, 即得： y0=-(y1+y2)=2 y1, 由此，还可将 y1 简化，即令： (x-2y1)(x-y1)^2=0, 有： x^3-4y1x^2+5y1^2x-2y1^3=0， x^3+a2x^2+a1x+a0=0 ，作变换： , x=y+4y1/3, x^2=y^2+8y1y/3+(4y1)^2/9, x^3=y^3+4y1y^2+(4y1)^2y/3-(4y1)^3/27 x 方程变换为： y^3+b1y+b0=0, 则有： y1^2-3b1/16=0, (a) y1^3+b0/(2-20/3+128/27)=0, (b) (a) 、 (b)2 式对 y1 逐次降幂，即得 y1 ： y1=8b0/(b1(7-64/9)), y2=8b0/(b1(7-64/9)), y0=16b0/(b1(7-64/9)), 于是，解得 1 元 3 次不可约代数方程不含任何根式的 3 个解。类似地，按此法，即可解得 5 次 ( 见： 4 次后的补充说明 ) ，直到 2m+1 次， 1 元不可约代数方程不含任何根式的解。 5 ． 1 元 4 次不可约代数方程的创新解法任意 4 次的不可约代数方程总可表达为： x^4+a2x^2+a1x+a0=0 ，还可将原方程表达为： (x^2+a’1x+a’0)(x^2+a”1x+a”0)=0 ，而有： a’1+a”1=0, a”1 =-a’1, (1) a’0+a”1a’1+a”0=a2, (2) a”0a’1+a”1a’0=a1, (3) a”0a’0=a0, a”0=a0/a’0, (4) (1) 、 (4) 均代入 (2) 、 (3) ： a’0-a’1^2+a0/a’0=a2, a’0^2-(a’1^2+a2)a’0+a0=0, (2’) a’0=(a’1^2+a2)/2+((a’1^2+a2)^2/4-a0)^(1/2), (2”) a0a’1/a’0-a’1a’0=a1, a’0^2+a1a’0/a’1-a0=0, (3’) a’0=(a1/a’1)/2+((a1/a’1)^2/4+a0)^(1/2), (3”) 由 (2”)-(3”) ： a’1^6a0+2a’1^4a2a0+a’1^2(a2^2a0-4a0^2)-a1^2a0=0 ， (5) (2’)-(3’) ，消去 a’0 ： a’1^3+ a’1a2 +a1 =0 ， (6) 由 (2’) 与 (6) 消去 a’1 ： (a’0^2+a1a’0+a0)(a’0^2-a2a’0+a0)=0 ，即： a’0^4-(a2-a1)a’0^3-(a2a1-2a0)a’0^2-(a2a0-a1a0)a’0+a0^2=0 ， (7) 由 (3’) 与 (6) 消去 a’1 ： a’0^6-3a0a’0^4+3a0^2a’0^2-a0^3=0, 即： a’0^2-a0=0 ， (8) 由 (5) 与 (6) ，按 a’1 逐次降幂，即得 : 以方程各系数表达的 a’1 , 即： a’1=-(a2a1-a1)/(a2^2-a2-4a0) ， (1*) (1*) 代入 (1) ，即得 : 以方程各系数表达的： a”1=(a2a1-a1)/(a2^2-a2-4a0) ， (2*) 由 (7) 与 (8) ，按 a’0 逐次降幂，即得 : 以方程各系数表达的 a’0 , 即： 2(a2a0-a1a0)a’0+(a2a1a0-4a0^2)=0 ， (9) 即得： a’0=-(a2a1/2-2a0)/(a2-a1) ， (3*) (3*) 代入 (4) ，即得 : 以方程各系数表达的： a”0=a0/a’0=-a0(a2-a1)/(a2a1/2-2a0) ， (4*) 于是，对于各相应系数的 2 个 2 次方程 x^2+a’1x+a’0=0 ， x^2+a”1+a”0=0 ，即可分别解得各 2 个解： x1=-a’1/2+((a’1/2)^2-a’0)^(1/2), x2=-a’1/2-((a’1/2)^2-a’0)^(1/2), x3=-a”1/2+((a”1/2)^2-a”0)^(1/2), x4=-a”1/2-((a”1/2)^2-a”0)^(1/2), 而得到，没有高于 2 次的根式，仅由实数、虚数或复数表达的， 1 元 4 次不可约代数方程的 4 个根式解。类似地，按此法，即可解得没有高于 2 次的根式，仅由实数、虚数或复数表达的， 6 次，直到 2m 次， 1 元不可约代数方程的根式解。 6 ． 1 元 5 次不可约代数方程的创新解法 ( 对第 4 节的补充说明 ) 任意 5 次不可约代数方程总可表达为： y^5+b3y^3+b2y^2+b1y+b0=0, 有 5 个根： y1 、 y2 、 y3 、 y4 、 y5 ，还可表达为： (y^4+b3’y^3+b2’y^2+b1’y+b0’)(y-y5) =0 ，对比原方程，即得： b3’-y5=0, b2’-b3’y5=b3, b1’-b2’y5=b2, b0’-b1’y5=b1,-b0’y5=b0, 即： b3’=y5 ， b2’=b3+y5^2, b1’=b2+b3y5+y5^3, b0’=b1+b2y5+b3y5^2+y5^4=-b0/y5, y^4+b3’y^3+b2’y^2+b1’y+b0’=0 ，成为： y^4+y5y^3+( b3+y5^2)y^2+( b2+b3y5+y5^3)y-b0/y5=0 ， (1) 将 (1) 变换为 3 次方系数 =0 的形式 (1’) ，即可按第 5 节方法，解得由 y5 和 (1’) 各系数，表达的原 5 次方程的其它 4 个根的方程，再利用这些根都满足原 5 次方程的条件逐次降幂，得到原 5 次方程其它 4 个根的解，再利用原 5 次方程的 5 个根之和 =0 ，即得 y5 之值。于是，解得 1 元 5 次不可约代数方程不含任何根式的 5 个解。由于解法已规范化，完全可以编制程序利用计算机求解，解决较高次方程，逐次降幂求解的繁琐工作。就彻底解决了有关需要的“解方程”问题。 7 ．参考文献：数学百科全书编委（顾问）苏步青等（主任）王元等科学出版社 1994 Basic algebra 1-2 Jacobson, N. Freeman1974-1980 Algebra 1-2 B.I. Van Der WaerdenSpringer-Verleg 1955-1959

个人分类: 数理|2486 次阅读|1 个评论

解方程：转基因+瘦肉精+地沟油+X=20年后中国人体质

热度 4 wangxh 2013-3-24 17:28

该方程有实数解吗？转基因+瘦肉精+地沟油+ X =20年后中国人体质俺打算研究X这种材料，行否？

1737 次阅读|4 个评论

解方程的故事

热度 39 daizhancheng 2013-3-10 14:34

上篇博文 “ 一知半解漫说数 ” 本意也是想介绍代数数与方程可解性（伽罗华理论）。后来还没到这里就已太长了，只好忍痛割爱。这篇就算其续篇吧。 1. 代数方程解代数方程是最古老的数学问题之一，也曾经是对人类智力最富挑战性的一个数学问题。代数方程依赖于两个因素：“ 元 ” 和 “ 次 ”。元是变量个数，次指变量最高幂次。所以最简单的是一元一次方程。当元增加时，就有多元一次方程组；当次增加时就有二次、三次或更高次方程；如果元和次都增加，就有多元高次方程组。多元一次方程组可用矩阵记成 $Ax=b$ ，这里 $A$ 可以是 $m\timesn$ 矩阵，$x=(x_1,\cdots,x_n)^T$ ， $b\inR^m$ 。如果系数阵与增广阵同秩，即$rank =rank $ ，方程就有解，否则就没有解。更详细的结果，可参见任意一本线性代数书，这不是本文的重点，就不多讲了。值得一提的是，我国在公元前一世纪出的《九章算术》中已经会解这类方程了。并且，是将系数排成方阵来解的。这就是矩阵的雏形，也是中国人称其为 “ 方程 ” 的原因。线性方程组解常常不唯一，于是通解可分为两部分：齐次方程的通解 $+$ 非齐次方程的特解。这种思想很重要，寻找线性微分或差分方程的一般解，思路是完全类似的。还有一类有趣的问题，称 “ 不定方程的整数解 ”。中国古代数学家就讨论过许多这种问题，例如东汉《孙子算经》中有：“ 今有物不知数，三三数之余二，五五数之余三，七七数之余二，问物几何？” 还有著名的 “ 百鸡百元 ” 问题：“ 公鸡 5 元 1 只，母鸡 3 元 1 只，小鸡 3 只 1 元，今花百元购得百鸡，问公、母、小鸡各几只？” 此题载于公元 5-6 世纪成书的《张邱建算经》，书中的解答开创了一题多解的先例。多元高次方程组解的讨论是一个极其困难的问题。有一个专门的数学分支研究它，称为 “ 代数几何 ”。本人在华盛顿大学时修过数学系博士生的代数几何课，当时 SSM 系修到此课的只有我一人，可见我对数学的痴迷。自己感觉它是修过的数学课中最艰深的，而本人该科成绩为 “ A+ ”。应行仁兄后来选了许多计算机的课，成了计算机专家，也就与我分道扬镳了。本文主要是谈一元高次方程。初中生都知道一元二次方程的解公式。它是公元八、九世纪阿拉伯数学家花拉子米给出的。其实，《九章算术》中也给出了二次方程的数值解法，但本文只对解的封闭形式 (ClosedFormSolution) 感兴趣。即只用系数的有限次加、减、乘、除、开方表示的解。 2 . 三次、四次方程的解三次方程解曾经是十六世纪最富挑战性的问题，现在通常把三次方程解公式叫做卡当公式。那么，三次方程解是卡当 ( Jerome Cardan , 1501-1576) 发现的吗？也是也不是。此话怎讲？且听我慢慢道来。卡当生于意大利，是一位法官与一位寡妇的私生子。那位法官博学多才，这影响了他。但因是私生子，他自幼受到歧视和虐待，长大后也因此在事业上屡遭排斥。坎坷的经历让他性格怪异。他在数学、哲学、物理学和医学中均有相当成就，但也醉心于占星术与赌博。他的自传《我的生平》，用批评的口吻剖析了自己的一生。回到三次方程，早在 1510 年左右，波伦亚大学教授费罗就发现了缺二次项的三次方程 $x^3+px=q$ 的解法。他死前传给他的学生菲奥尔。 1530 年，另一位意大利数学教师塔尔塔利亚得到缺一次项的三次方程的解法。 1535 年，年轻气盛的菲奥尔听说塔尔塔利亚也会解三次方程，他不信，就向塔尔塔利亚提出挑战。塔尔塔利亚潜心研究，终于在比赛前掌握了缺二次项的三次方程的解法。于是，在比赛中塔尔塔利亚完胜。当时卡当正在写一本代数专著，听到这消息后就邀请塔尔塔利亚前来做客。塔来后，卡当当面再三恳求并发誓对此保密后，塔将解法告诉了他。卡当经研究，给出了三次方程一般解及其证明。卡当后来在《大术》一书中发表了这些结果，这导致他以及他的学生与塔的长期论争。现在看，解三次方程 $z^3+az^2+bz=c$ 其实很容易：设 $z=x-a/3$ ，则得到缺二次项的三次方程 $x^3+px=q$ 。然后，设 $x=u+v$ ，就可得 $3uv=-p$ ， $u^3+v^3=-q$ 。于是可解出 $u^3$ 及 $v^3$ ，也就可得 $u$ 及 $v$ 。从而解出$x=u+v\leftrightarrow z=x-a/3$ 。任何一本数学手册都能找到卡当公式，比如。在使用卡当公式时，会遇到一个实数解却要用负数开平方 $\sqrt{-a}$ 来表示的情况，这时，卡当勇敢地使用了虚数及复数运算，因此，有些书说虚数是解三次方程时引入的，指的就是卡当的工作。随后，在他的一个学生 L . Ferrari 的帮助下，卡当进而得到四次方程解。考虑 $x^4+px^2+qx+r=0$ （如有三次项 $ax^3$ ，设 $x=z-a/4$ 就可将三次项去掉），在 $x^4=-px^2-qx-r$ 两边加 $2zx^2+z^2$ ，则左边成完全平方。右边是 $(2z-p)x^2-qx+(z^2-r)$ 。要让右边也成平方数，令 $2\sqrt{2z-p}\sqrt{z^2-r}=-q$(*) ，则右边也成完全平方，于是就可解出 $x$ （含参数 $z$ ）。 (*) 可化简为一个关于 $z$ 的三次方程，称为预解方程。解出 $z$ ，也就得到真正的 $x$ 了。详细介绍可见。据说，卡当为了证明自己的巫术，在预言自己的死亡之日自杀了。虽然可悲，但他没有白活，因为他的名字已经被历史记住了。 3. 五次方程有公式解吗？中国学人大概都知道，华罗庚先生在 20 岁时（ 1930 年）以一篇论文《苏家驹的五次方程式解法不能成立之理由》崭露头角，被熊庆来破格引进清华，开始了他的数学家生涯。苏家驹不知道，“ 五次方程没有公式解 ” 在他苦思冥想之前的一百年，都已经被证明了。中国至今还有许多人在做规矩三等分任意角，在寻找五次方程公式解（我在海淀图书城还见过号称解决了五次方程这个国际难题的正式出版的书），就因为他们缺少这些常识。前几天还有个 “ 民间数学家 ”，拿了一摞手稿，到我办公室，说他知道怎么三等分一个角。我告诉他这是不可能的。他可怜兮兮地说：“ 我研究了近 20 年，才得到这个结果。你无论如何帮我看看…… ” 搞得我哭笑不得。不知这篇文章会不会有助于减少一、两个痴迷的 “ 民数 ”。现在一提起五次方程没有公式解，就会想到伽罗华。其实，最早证明这个结果的是阿贝尔。阿贝尔与伽罗华一样，也是一个英年早逝的数学天才，他的数学贡献，很晚才被人完全发现。他长期找不到工作， 27 岁时在贫困与疾病的交逼下忧伤而死。他死后两天，柏林大学数学教授的聘书才寄到，令人扼腕。伽罗华得到的是一般的代数方程有解的判定定理。用伽罗华理论，不仅能证明五次方程没有公式解，还可以证明古希腊的三大数学难题同样无解。这些问题是： (1) 规矩三等分任意角； (2) 化圆为方：即做一正方形，使其与己知圆等面积； (3) 倍立方：做一个立方体，其体积为已知立方体两倍。 4. 伽罗华的故事伽罗华父亲是法国资产阶级革命的支持者，在拿破仑 “ 百日政变 ” 时期担任过皇后市自由派市长。其母亲是一位法官的女儿，聪明且受过良好教育，但性格孤僻，行事古怪。伽罗华出生于 1811 年。早年由他母亲自行授课。 11 岁进入路易皇家中学学习，成绩优异。十四岁时碰上一个优秀的数学老师，从此疯狂地爱上了数学。十五岁就开始阅读拉格朗日的原著。偏科使他总成绩下降。学校对他的评价是：“ 奇特、怪异，有原创力又封闭。” 1828 年，与阿贝尔一样，伽罗华一度认为自己找到了五次方程解。但不久，他就找出错误并证明了其不可能性。他的关于群的论文被呈送法国科学院，由柯西主审。这年，他才 17 岁。由于文章过于简略，柯西建议修改。 1831 年，他修改后的论文被送交泊松审查。泊松认真审阅后，结论是 “ 不可理解 ”。他在给科学院的报告中说：“ 我们已经尽了最大努力来研究伽罗华的证明，他的推理显得很不清楚。我们甚至不能弄清他的证明思想。” 他超越了他的时代，而时代也背弃了他。也许受到其父亲因政治迫害而自杀的影响，伽罗华是一名激进的共和主义战士。因为抨击校长，他被学校开除。因为嘲讽国王和率众游行，他两次被捕入狱。 1832 年 3 月，由于霍乱流行，伽罗华被转移到一家私人医院中服刑，在那里，他爱上了医生的女儿。一个月后他获释。医生女儿的前男友提出要和他决斗，他明知此人枪法甚好，但不知为了爱情还是为了荣誉，他接受了挑战。决斗的头一天晚上，他自知必死，彻夜整理自己的手稿，手稿边上多次出现：“ 我没有时间了 ” 的字样。 1832 年 5 月 30 日，他在决斗中中弹致伤，许久后他的兄弟才闻讯赶到，将其送进医院。但终因失血过多，他于第二天与世长辞。一代旷世的数学奇才，死时还未满 21 周岁。伽罗瓦的成就是提出群的概念，并用群论彻底解决了代数方程可解性问题。他的开创性工作，今天被称为伽罗华理论，但当的并未得到同时代人的充分赏识和理解，他的思想太超前了，以至柯西、泊松等同辈数学大家都无法理解他。后来，经刘维尔等人的努力，直到 1870 年，若当出版了《置换和代数方程专论》一书，全面介绍了伽罗华的工作，才使伽罗华理论融入了主流数学。 5 . 域的扩张与古希腊三大难题上篇博文提到群和域，这一节只给还记得这些概念的朋友写的。从有理数域 $Q$ 出发，如果添加一个代数数 $\sqrt{2}$ ，那么，容易证明 $E_1=Q(\sqrt{2})=\{a+b\sqrt{2}|a,b\in Q\}$ 是一个域，这就叫域的添加。如果把 $Q$ 当一维空间， $E_1$ 就是 2 维空间，记相对阶 (degree) 作$ =2$ 。如果添加 $2^{1/3}$ ，那么 $E_2=Q(2^{1/3})=\{a+b 2^{1/3}+c 2^{2/3}|a,b,c\in Q\}$ 也是一个域，并且$ =3$ 。相对阶其实就是添加元满足的最低次数有理系数多项式零点。回到尺规作图问题，在平面上建一直角坐标，如果$(p,q)$ 是有理数，则$(p,q)$ 点可用尺规作图找到，这是熟知的。从有理点出发，尺规作图能走多远呢？实际我们能找到三类点： (1) 直线与直线相交； (2) 直线与圆相交， (3) 圆与圆相交。如果两直线都是过两个有理点画出来的，交点也是有理点，所以不会增加新点；后两种情况，由于圆满足二次方程，容易看出，新的点$(p,q)$ 满足： $p,q\in F,~~~~ =2.$ 即$(p,q)$ 只能在相对阶为 $2$ 的域里。实际上，相对阶为 2 的数域里的数都能做出来。从这些点出发，只有相对阶为 4 的数能做出来。如此等等。于是，容易知道：一个数 $r$ 能由尺规做出来，当且仅当 $ =2^k.$ 回头看古希腊三大难题： (1) 倍立方：设原立方体边长为 1 ，则新立方体边长为 $2^{1/3}$ 。因为$ =3$ ，则显然尺规做不出这个数。 (2) 化圆为方 : 设圆半径为 1 ，则正方形边长为 $\sqrt{\pi}$ 。因为 $\sqrt{\pi}$ 是超越数$ =\infty$ ，显然尺规也做不出这个数。 (3) 三等分任意角 : 只要证明某个特定角不能三等分即可。选 $60^o$ 角。设 $x=\cos(20^o)$ ，由 $\cos(3x)=4\cos^3(x)-3\cos(x)$ 可知 $x$ 满足的最小化零多项式为 $P(x)=x^3-\frac{3}{4}x-\frac{1}{8}$ 因此 $ =3$ 显然尺规也做不出这个数，故 $60^o$ 不能用尺规三等分。 6. 伽罗华理论（这一节是为非数学专业的数学爱好者写的，要求一定数学基础，例如，线性代数教材通常会介绍一点群和正规子群等）设 $F_1$ 与 $F_2$ 为两个域，如果存在一个映射$\pi:F_1\rightarrow F_2$ ，它是一一对应的，并且保持加、乘运算关系（即 $\pi(a*b)=\pi(a)*\pi(b)$ ，这里 $*$ 表加法或乘法，等），则称 $F_1$ 与 $F_2$ 同构； $\pi$ 为 $F_1$ 与 $F_2$ 的一个同构映射。如果 $F_2=F_1$ ，则称 $\pi$ 为 $F_1$ 上的自同构（映射）。容易看出，自同构关于复合运算构成一个群，称为自同构群。给定一个域 $F$ （例如 $F=Q$ ）和 $F$ 上的一个多项式 $a(x)$ ，将 $a(x)$ 的根都添加到 $F$ 上，作扩张，得到一个扩张域 $K$ ，称为 $F$ 的关于 $a(x)$ 的根域。 $K$ 上的所有使 $F$ 不动的自同构，形成 $K$ 的自同构群的一个子群，称为$K$ 在 $F$ 上的伽罗华群，记作$Gal(K:F)$ 。设 $G$ 为一个群， $H$ 为其子群，如果对任何 $g\inG$ ，有 $gH=Hg$ ，则称 $H$ 为 $G$ 的正规子群，记作 $H\vartriangleleftG$ 。显然，如果 $G$ 是阿贝尔群，任一子群都是正规的。一个群 $G$ 称为可解的，如果存在子群 $H_i$ ,$i=1,2,\cdots,m$ ，使得 $e=H_0\vartriangleleftH_1\vartriangleleftH_2\vartriangleleft\cdots\vartriangleleftH_m=G$ 并且，商群 $H_{k+1}/H_k$ ，$k=0,1,\cdots,m-1$ 为阿贝尔群。伽罗华的基本理论是：多项式 $a(x)\inF(x)$ 的根（用根式等）可解，当且仅当，其伽罗华群可解。现在，找一个五次多项式 $p(x)=x^5-5x-2$ ，记 $K=Q(p(x))$ 为其根域，容易证明，$Gal(K,Q)$ 与五阶置换群（ PermutationGroup ） $S_5$ 同构。不难证明 $S_5$ 不可解，故伽罗华群$Gal(K,Q)$ 不可解。从而证明 $x^5-5x-2=0$ 的根不都能用根式等表出。欲知其详，可参见。 7. 后记 2006 年，时任自动化学会教材委员会秘书长、清华大学萧德云教授来访，邀我写一本向工科学生介绍近代数学的书。被他的热情和构想所感动，就写了。其目的在书封底 “ 本书特色 ” 做了说明，摘抄如下： “ 本书的目的是向自动化及相关学科的研究生、青年教师及科研人员介绍近代数学。选材以与系统控制关系密切的分析、代数及几何知识为主线，兼顾数学自身的内在逻辑与严谨性。内容包括：勒贝格测度与积分、泛函分析、点集拓扑，抽象代数、代数拓扑、微分流形、李群与李代数、张量场、黎曼几何、辛几何、代数几何等近代数学的主要分支。作者希望本书能使仅具有工程数学基础的青年学者迅速掌握近代数学的工具，以满足阅读相关文献和从事自动化理论研究的需要。除密切联系自动化学科的需求外，本书也注意为读者展示数学的百花园。例如，证明用尺规三等分任意角是不可能的，证明五次方程没有公式解，介绍什么是选择公理、 Zorn 引理、超限归纳法，连续统假定、罗素悖论、四色问题、费马大定理与椭圆曲线、庞加莱猜想、黎曼猜想…… 读完本书，你或许会发现，你登上了一个新的高度，有一个新的视野看世界…… ” 自己很赞同萧教授及自动化学会教材委员会的想法：中国研究生的一个普遍弱点是基础不足，急于做研究、出文章。要想提高研究生水平，要从加强基础理论的学习和训练做起。这本书是这种努力的一次尝试，也许不是很成功，因为有多位读者反映难懂。但 502 所一位清华老学长陆先生几次来电话鼓励我，说写这本书是 “ 功德无量 ”。现在我在研究工作之余，抽空在写第二版，并增加随机部分，但愿能改得通俗易懂一点。参考文献：叶其孝，沈永欢，《实用数学手册》，第 2 版，科学出版社， 2006 . 李文林，《数学史教程》，高等教育出版社，北京， 2000 . V . J . Katz , AHistoryofMathematics , ReprintedbyChinaMachinePress , 20 04. 吴文俊 ( 主编 ) ，世界著名数学家传记，科学出版社， 1995 . C . C . Pinter , ABookofAbstractAlgebra , 2nd Ed ., McGraw-HillPublishingComp ., 1990 . J . Bewersdorff , GaloisTheoryforBeginners , AMS , RhodeIsland , 2006 . 程代展，《系统与控制中的近代数学基础》，清华大学出版社， 2007 .

个人分类: 其他|34463 次阅读|112 个评论

matlab解方程

liujd 2012-7-30 18:06

1，solve('x^2=4','x') 用单引号把包含字母的方程写出来，后面指定未知数是什么字母。其中可以包含非数值的常量参数，如a,b,c等。solve('2*x=40+y','x') 2，方程组 =solve('cy=f*zy','cx=e*zx^2','z=zx+zy','pxpy=f/(2*zx*e)','pxpy*(b/cy)-(a/cx)=0','pxpy,zx,zy,cx,cy') 用逗号分隔每个方程，最后一个字符串中用逗号分隔未知数。函数返回结果是一个列向量。

个人分类: matlab|10553 次阅读|0 个评论

讨论班（有讲稿）：谱方法（Spectral Method）解方程（王璐）09.12.3

GrandFT 2009-10-22 12:17

周四讨论班：谱方法（Spectral Method）解方程时间：2009.12.3 下午 4：306：10 地点：18教学楼405室主讲：王璐文献： Philippe Grandclment and Jrme Novak , Spectral Methods for Numerical Relativity 的第一、二章。这部分不涉及广义相对论的内容，只需了解一些偏微分方程就可以了。文献下载王璐的讲稿下载

个人分类: 周四讨论班|5518 次阅读|1 个评论

帐号		自动登录	找回密码
密码			注册

关闭安全验证

标签: 解方程

相关帖子

相关日志

关闭 安全验证

标签: 解方程

相关帖子

相关日志

关闭安全验证