众文网1.简短的幽默故事
【因人而异】 佛雷把刚买的高尔夫球用具拿给朋友看：“我和太太都胖了一点，医生说主要是运动量不足，为此我准备加入高尔夫俱乐部。”
“那你太太呢？”“我已买了一部剪草机给她。” 【晓岚佛前释笑】 一天，纪晓岚（清代大学者）陪同乾隆皇帝游大佛寺。
君臣二人来到天王殿，但见殿内正中一尊大肚弥勒佛，坦胸露腹，正在看着他们憨笑。乾隆问：“此佛为何见朕笑？”纪晓岚从容答道：“此乃佛见佛笑。”
乾隆问：“此话怎讲？” 纪晓岚道：“圣上乃文殊菩萨转世，当今之活佛，今朝又来佛殿礼佛，所以说是佛见佛笑。 ”乾隆暗暗赞许，转身欲走，忽见大肚弥勒佛正对纪晓岚笑，回身又问：“那佛也看卿笑，又是为何？”纪晓岚说：“圣上，佛看臣笑，是笑臣不能成佛。”
乾隆称赞纪晓岚善辩。 【从未失败】 有人问一位科学家：“你试验一种新型电池总是失败，为什么还要继续试验？”科学家回答：“失败？我从来没有失败过，我现在已知道了5万种不能制造这种电池的方法。
” 2005年2月12日 一个富人去拜访一位哲学家，请教他为什么自己有钱后变得越发狭隘自私了。哲学家将他带到窗前问：“向外看，告诉我你看到了什么？”富人说：“我看到了外面世界的很多人。”
哲学家又把他带到一面镜子前问：“现在又看到了什么？”富人答：“我自己。 ”哲学家一笑说：“窗子和镜子都是玻璃做的，区别只在于镜子多了一层薄薄的银子。
但就是因为这一点银子，便叫你只看到自己而看不到世界了。” 。
2.研究线性方程组的jacobi和gauss
①雅克比迭代法：
function [n,x]=jacobi(A,b,X,nm,w)
%用雅克比迭代法求解方程组Ax=b
%输入：A为方程组的系数矩阵，b为方程组右端的列向量，X为迭代初值构成的列向量，nm为最大迭代次数，w为误差精度
%输出：x为求得的方程组的解构成的列向量，n为迭代次数
n=1;
m=length(A);
D=diag(diag(A)）； %令A=D-L-U，计算矩阵D
L=tril(-A)+D； %令A=D-L-U，计算矩阵L
U=triu(-A)+D； %令A=D-L-U，计算矩阵U
M=inv(D)*(L+U)； %计算迭代矩阵
g=inv(D)*b； %计算迭代格式中的常数项
%下面是迭代过程
while n<=nm
x=M*X+g； %用迭代格式进行迭代
if norm(x-X,2)<w
disp（'迭代次数为'）；n
disp（'方程组的解为'）；x
return;
%上面：达到精度要求就结束程序，输出迭代次数和方程组的解
end
X=x;n=n+1;
end
%下面：如果达到最大迭代次数仍不收敛，输出警告语句及迭代的最终结果（并不是方程组的解）
disp（'在最大迭代次数内不收敛！'）；
disp（'最大迭代次数后的结果为'）；x
②高斯赛德尔迭代法：
function [n,x]=gaussseidel(A,b,X,nm,w)
%用高斯-赛德尔迭代法求解方程组Ax=b
%输入：A为方程组的系数矩阵，b为方程组右端的列向量，X为迭代初值构成的列向量，nm为最大迭代次数，w为误差精度
%输出：x为求得的方程组的解构成的列向量，n为迭代次数
n=1;
m=length(A);
I=eye(m)； %生成m*m阶的单位矩阵
D=diag(diag(A)）； %令A=D-L-U，计算矩阵D
L=tril(-A)+D； %令A=D-L-U，计算矩阵L
U=triu(-A)+D； %令A=D-L-U，计算矩阵U
M=inv(D-L)*U； %计算迭代矩阵
g=inv(I-inv(D)*L)*(inv(D)*b）； %计算迭代格式中的常数项
%下面是迭代过程
while n<=nm
x=M*X+g； %用迭代格式进行迭代
if norm(x-X,2)<w
disp（'迭代次数为'）；n
disp（'方程组的解为'）；x
return;
%上面：达到精度要求就结束程序，输出迭代次数和方程组的解
end
X=x;n=n+1;
end
%下面：如果达到最大迭代次数仍不收敛，输出警告语句及迭代的最终结果（并不是方程组的解）
disp（'在最大迭代次数内不收敛！'）；
disp（'最大迭代次数后的结果为'）；x
上面是我上数值分析这门课时用matlab所编的程序，都运行过，没有错误，只要复制运行即可
3.关于高斯精神的议论文
高斯是德国数学家 ，也是科学家，他和牛顿、阿基米德，被誉为有史以来的三大数学家。
高斯是近代数学奠基者之一，在历史上影响之大， 可以和阿基米德、牛顿、欧拉并列，有“数学王子”之称。 他幼年时就表现出超人的数学天才。
1795年进入格丁根大学学习。第二年他就发现正十七边形的尺规作图法。
并给出可用尺规作出的正多边形的条件，解决了欧几里得以来悬而未决的问题。 高斯的数学研究几乎遍及所有领域，在数论、代数学、非欧几何、复变函数和微分几何等方面都做出了开创性的贡献。
他还把数学应用于天文学、大地测量学和磁学的研究，发明了最小二乘法原理。高理的数论研究 总结 在《算术研究》（1801）中，这本书奠定了近代数论的基础，它不仅是数论方面的划时代之作，也是数学史上不可多得的经典著作之一。
高斯对代数学的重要贡献是证明了代数基本定理，他的存在性证明开创了数学研究的新途径。高斯在1816年左右就得到非欧几何的原理。
他还深入研究复变函数，建立了一些基本概念发现了著名的柯西积分定理。他还发现椭圆函数的双周期性，但这些工作在他生前都没发表出来。
1828年高斯出版了《关于曲面的一般研究》，全面系统地阐述了空间曲面的微分几何学，并提出内蕴曲面理论。高斯的曲面理论后来由黎曼发展。
高斯一生共发表155篇论文，他对待学问十分严谨，只是把他自己认为是十分成熟的作品发表出来。其著作还有《地磁概念》和《论与距离平方成反比的引力和斥力的普遍定律》等。
1801年高斯有机会戏剧性地施展他的优势的计算技巧。那年的元旦，有一个后来被证认为小行星并被命名为谷神星的天体被发现当时它好像在向太阳靠近，天文学家虽然有40天的时间可以观察它，但还不能计算出它的轨道。
高斯只作了3次观测就提出了一种计算轨道参数的方法，而且达到的精确度使得天文学家在1801年末和1802年初能够毫无困难地再确定谷神星的位置。高斯在这一计算方法中用到了他大约在1794年创造的最小二乘法（一种可从特定计算得到最小的方差和中求出最佳估值的方法在天文学中这一成就立即得到公认。
他在《天体运动理论》中叙述的方法今天仍在使用，只要稍作修改就能适应现代计算机的要求。高斯在小行星“智神星”方面也获得类似的成功。
由于高斯在数学、天文学、大地测量学和物理学中的杰出研究成果，他被选为许多科学院和学术团体的成员。“数学之王”的称号是对他一生恰如其分的赞颂。
4.高斯赛德尔迭代法matlab编程
去百度文库，查看完整内容>
内容来自用户：w___w000
function[x,k]=GaussSeidel(A,b,x0,eps,M)
%高斯赛德尔迭代法求方程组的解（矩阵公式求解）
%A为方程组的系数矩阵；b为方程组的右端项
%x为线性方程组的解了；x0为迭代初值
%eps为误差限；M为迭代的最大次数
if nargin==3
eps= 1.0e-6；%默认精度
M = 10000；%参数不足时默认后两个条件
elseif nargin ==4
M = 10000；%参数的默认值
elseif nargin<3
error（'参数不足'）；
return
end
[n,m]=size(A);
nb=length(b);
%当方程组行与列的维数不相等时，停止计算，并输出出错信息
if n~=m
error（'矩阵A行数和列数必须相等！'）；
return;
end
%当方程组与右端项的维数不匹配时，停止计算，并输出出错信息
if n~=nb
error（'矩阵A的行数必须和b的长度相等！'）；
return;
end
L =zeros(n,n);
U =zeros(n,n);
D =zeros(n,n);
for i=2:n
for j=1:i-1
L(i,j)=-A(i,j);
end
end
for i=1:n-1
for j=i+1:n
U(i,j)=-A(i,j);
end
end
for i=1:n
D(i,i)=A(i,i);
end
B=inv(D-L)*U; %B为迭代矩阵
g=inv(D-L)*b; %g为右端项
pr=max(abs(eig(B)））； %求迭代矩阵谱半径
if pr>=1
error（'迭代矩阵谱半径大于1迭代法不收敛'）；
return;
end
k=0;
tol=1;
while tol>=eps
x = B*x0+g;
k = k+1； %迭代步数
tol = norm(x-x0)；%前后
转载请注明出处众文网 » 高斯迭代法的毕业论文(简短的幽默故事)