MATLAB第二次上机实验报告.doc

电子科技大学 电子工程 学院 标 准 实 验 报 告 （实验）课程名称 MATLAB与数值分析 学生姓名： 学 号： 指导教师： 实验名称 实验二 线性方程组求解和函数的数值逼近 实验目的 通过上机实验，使学生对病态问题、线性方程组求解和函数的数值逼近方法有一个初步的理解。 实验涉及的核心知识点病态方程求解、矩阵分解和方程组求解、Lagrange插值。 实验重点与难点算法设计和MATLAB编程 1. 对高阶多项式 编程求下面方程的解 并绘图演示方程的解与扰动量的关系。 2. 对，生成对应的Hilbert矩阵，计算矩阵的条件数；通过先确定解获得常向量b的方法，确定方程组 最后，用矩阵分解方法求解方程组，并分析计算结果。 3. 对函数 的Chebyshev点 ， 编程进行Lagrange插值，并分析插值结果。 实验数据及结果分析 对高阶多项式 编程求下面方程的解 并绘图演示方程的解与扰动量的关系。 p=[1,-1]; for i=2:20 n=[1,-i]; p=conv(p,n); % 求多项式乘积 end m=zeros(1,21); % m的最高次幂为20，有21项 hold on x=1:20; d=[-1,0,0.1,0.5,1]; for i=1:5 delt=d(i); m(2)=delt; y=(roots(p+m)); % 求多项式的根 plot(x,y,-o,color,[i/5,i/20,i/10]); end title(方程p(x)=0的解与扰动量delt的关系) legend(delt=-1,delt=0,delt=0.1,delt=0.5,delt=1) 2.对，生成对应的Hilbert矩阵，计算矩阵的条件数；通过先确定解获得常向量b的方法，确定方程组 最后，用矩阵分解方法求解方程组，并分析计算结果。 for n=2:20 h=hilb(n); fprintf(\n\nn=%-10dcond(Hn)∞=%d,n,cond(h,inf))%cond：求矩阵范数 X=1:n; b=h*X; [l u]=lu(h); %lu分解 x=u\(l\b); %利用lu分解求线性其次方程组的根 x=x; fprintf(\n X’=) for i=1:n fprintf(%-8.2f,X(i)) end fprintf(\n x’=) for i=1:n fprintf(%-8.2f,x(i)) end end 输出结果如下： n=2 cond(Hn)∞=2.700000e+001 X’=1.00 2.00 x’=1.00 2.00 n=3 cond(Hn)∞=7.480000e+002 X’=1.00 2.00 3.00 x’=1.00 2.00 3.00 n=4 cond(Hn)∞=2.837500e+004 X’=1.00 2.00 3.00 4.00 x’=1.00 2.00 3.00 4.00 n=5 cond(Hn)∞=9.436560e+005 X’=1.00 2.00 3.00 4.00 5.00 x’=1.00 2.00 3.00 4.00 5.00 n=6 cond(Hn)∞=2.907028e+007 X’=1.00 2.00 3.00 4.00 5.00 6.00 x’=1.00 2.00 3.00 4.00 5.00 6.00 n=7 cond(Hn)∞=9.851949e+008 X’=1.00 2.00 3.00 4.00 5.00 6.00 7.00 x’=1.0