(完整word版)同济大学数值分析matlab编程题汇编.pdfVIP

MATLAB 编程题库 ax b 1.下面的数据表近似地满足函数 y 2 ，请适当变换成为线性最小二乘问题，编程求 1 cx 最好的系数 a ,b,c ，并在同一个图上画出所有数据和函数图像 . xi 0.931 0.586 0.362 0.213 0.008 0.544 0.628 0.995 y i 0.356 0.606 0.687 0.802 0.823 0.801 0.718 0.625 解： x=[-0.931 -0.586 -0.362 -0.213 0.008 0.544 0.628 0.995]; y=[0.356 0.606 0.687 0.802 0.823 0.801 0.718 0.625]; A=[x ones(8,1) -x.^2.*y]; z=A\y; a=z(1); b=z(2); c=z(3); xh=-1:0.1:1; yh=(a.*xh+b)./(1+c.*xh.^2); plot(x,y,r+,xh,yh,b*) 1 2.若在 Matlab 工作目录下已经有如下两个函数文件，写一个割线法程序，求出这两个函数 精度为 10 10 的近似根，并写出调用方式： 文件一 文件二 function v = f(x) function z = g(y) v = x .* log(x) - 1; z = y.^5 + y - 1; 解： edit gexianfa.m function [x iter]=gexianfa(f,x0,x1,tol) iter=0; while(norm(x1-x0)tol) iter=iter+1; x=x1-feval(f,x1).*(x1-x0)./(feval(f,x1)-feval(f,x0)); x0=x1;x1=x; end edit f.m function v=f(x) v=x.*log(x)-1; edit g.m function z=g(y) z=y.^5+y-1; [x1 iter1]=gexianfa(f,1,3,1e-10) x1 = 1.7632 iter1 = 6 [x2 iter2]=gexianfa(g,0,1,1e-10) x2 = 0.7549 iter2 = 8 2 3.使用 GS 迭代求解下述线性代数方程组： ì ?5x + 2x + x = - 12 ? 1 2 3 ? ? - x + 4x + 2x = 10 í 1 2