程序设计数值分析.doc

程序设计： 本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x); y=-x*x-sin(x) 写成如上形式即可，下面给出主程序。 二分法源程序： %%%给定求解区间 b=1.5; a=0; %%%误差 R=1; k=0;%迭代次数初值 while (R5e-6) ; c=(a+b)/2; if f12(a)*f12(c)0; a=c; else b=c; end R=b-a;%求出误差 k=k+1; end x=c%给出解 Newton法及改进的Newton法源程序： clear %%%% 输入函数 f=input(请输入需要求解函数,s) %%%求解f(x)的导数 df=diff(f); %%%改进常数或重根数 miu=2; %%%初始值x0 x0=input(input initial value x0); k=0;%迭代次数 max=100;%最大迭代次数 R=eval(subs(f,x0,x));%求解f(x0)，以确定初值x0时否就是解 while (abs(R)1e-8) x1=x0-miu*eval(subs(f,x0,x))/eval(subs(df,x0,x)); R=x1-x0; x0=x1; k=k+1; if (eval(subs(f,x0,x))1e-10); break end if kmax;%如果迭代次数大于给定值，认为迭代不收敛，重新输入初值 ss=input(maybe result is error,choose a new x0,y/n?,s); if strcmp(ss,y) x0=input(input initial value x0); k=0; else break end end end k;%给出迭代次数 x=x0；%给出解 程序设计： 本实验采用Matlab的M文件编写。 Gauss消去法源程序： clear a=input(输入系数阵：\n) b=input(输入列阵b：\n) n=length(b); A=[a b] x=zeros(n,1); %%%函数主体 for k=1:n-1; %%%是否进行主元选取 if abs(A(k,k))yipusilong;%事先给定的认为有必要选主元的小数 yzhuyuan=1; else yzhuyuan=0; end if yzhuyuan; %%%%选主元 t=A(k,k); for r=k+1:n; if abs(A(r,k))abs(t) p=r; else p=k; end end %%%交换元素 if p~=k; for q=k:n+1; s=A(k,q); A(k,q)=A(p,q); A(p,q)=s; end end end %%%判断系数矩阵是否奇异或病态非常严重 if abs(A(k,k)) yipusilong disp(‘矩阵奇异，解可能不正确’) end %%%%计算消元,得三角阵 for r=k+1:n; m=A(r,k)/A(k,k); for q=k:n+1; A(r,q)=A(r,q)-A(k,q)*m; end end end %%%%求解x x(n)=A(n,n+1)/A(n,n); for k=n-1:-1:1; s=0; for r=k+1:n; s=s+A(k,r)*x(r); end t=(A(k,n+1)-s) x(k)=(A(k,n+1)-s)/A(k,k) end 本实验采用Matlab的M文件编写。其中待插值的方程写成function的方式，如下 function y=f(x); y=1/（1+25*x*x）; 写成如上