4特征值线性方程组求解数值方法求定积分拟合和插值..docVIP

《数学实验》报告 实验名称?? 数学实验四 ? 学 院 专业班级??? ??? 姓??? 名?? ?? 学??? 号??? ??? ?? 2013年5月 【实验目的】 掌握矩阵的基本运算、特征值、特征向量和线性方程组的求解；能熟练运用数值方法求定积分。掌握多项式的拟合和插值。 【实验任务】P114 12,14 P115 21 (1),(2) P167 17 (2),18 P130 8,10,12 【实验程序】 P114 12,14 12题： a=20*rand(6,6) A=round(a) A det(A) rank(A) rref(A) 14题： A=[2 1 1;1 2 1;1 1 2] p=poly(A) poly2str(p,x) [V,D]=eig(A) P115 21 (1),(2) 21（1）题： A=[1 1 2 -1;-1 1 3 0;2 -3 4 -1] rref(A) 21（2）题： B=[1 -1 -1 1 0;1 -1 1 -3 1;1 -1 -2 3 -1/2] rref(B) P167 17 (2),18 17（2）题： x=0:0.01:pi; y=x.*sin(x)./(1.+(cos(x)).^2); format long t=length(x); s1=sum(y(1:(t-1)))*0.01 s2=sum(y(2:t))*0.01 s3=trapz(x,y) s4=quad(x.*sin(x)./(1.+(cos(x)).^2),0,pi) 18题： x=0:0.01:pi/4; y=1./(1-sin(x)); format long t=length(x); z1=sum(y(1:(t-1)))*0.01 z2=sum(y(2:t))*0.01 z3=trapz(x,y) z4=quad(1./(1-sin(x)),0,pi/4) format short u1=z1-2^0.5, u2=z2-2^0.5, u3=z3-2^0.5, u4=z4-2^0.5 P130 8,10,12 8题： clf x=[0.10 0.30 0.40 0.55 0.70 0.80 0.95]; y=[15 18 19 21 22.6 23.8 26]; p1=polyfit(x,y,1); p3=polyfit(x,y,3); p5=polyfit(x,y,5); disp(一阶拟合函数),f1=poly2str(p1,x); disp(三阶拟合函数),f3=poly2str(p3,x); disp(五阶拟合函数),f5=poly2str(p5,x); x1=0.1:0.01:0.95; y1=polyval(p1,x1); y3=polyval(p3,x1); y5=polyval(p5,x1); plot(x,y,rp,x1,y1,b--,x1,y3,k-,x1,y5,g-); legend(拟合点,一次拟合,三次拟合,五次拟合) 10题： x=10:5:30; y=[25.2 29.8 31.2 31.7 29.4]; x1=10:1:30; y1=interp1(x,y,x1,*nearest); y2=interp1(x,y,x1,*linear); y3=interp1(x,y,x1,*spline); y4=interp1(x,y,x1,*cubic); subplot(221),plot(x,y,rp,x1,y1),title(最近点插值) subplot(222),plot(x,y,rp,x1,y2),title(线性插值) subplot(223),plot(x,y,rp,x1,y3),title(样条插值) subplot(224),plot(x,y,rp,x1,y4),title(立方插值) 12题： [x,y]=meshgrid(-3:0.5:3); z=x.^2/16-y.^2/9; [x1,y1]=meshgrid(-3:0.1:3); z1=x1.^2/16-y1.^2/9; figure(1) subplot(121),mesh(x,y,z),title(数据点) subplot(122),mesh(x1,y1,z1),title(函数图象) [x2,y2]=meshgrid(-3:0.125:3); z2=interp2(x,y,z,x2,y2,*nearest); z3=interp2(x,y,z,x2,y2,*linear); z4=interp2(x,y,z,x2,y2,*splin