数学实验分解.doc

数学实验报告 院 系 计算机学院 班 级 2004级二学位 指导教师 张 兴 永 姓 名 田 琳 实 验 目 录 1．计算 3 2. 求的fourier、Laplace、ztrans变换 3 3．矩阵高斯消去法 4 4．画出心型线、星型线、双纽线和四叶玫瑰线的图形 5 5．求fibonacci数 6 6．数据三次拟合曲线 6 7．递推公式的稳定性 7 实验程序: 7 实验结果: 8 8．迭代法的收敛性与收敛速度的比较 9 实验程序 9 实验结果 11 9．雅可比迭代法与高斯塞德尔法的收敛性与收敛速度 11 实验结果 13 10．龙格现象的发生、防止和插值效果的比较 13 结算结果 14 龙格现象分析 19 11．[1]P129页拟合应用题 19 实验结果 20 12．[1]课题1 数值积分的方法的使用、比较 21 实验结果 23 13．数据样条插值 24 14．统计应用 25 15．油罐标尺刻度设计 27 课程体会 致谢 31 参考书 31 图表目录 错误！未定义书签。 1．计算 程序内容：源程序：\s.m； 程序清单： function s(x) s=1; for(i=1:8); s=s*cos(2^i*x); end s 实验结果：S(1)= 3.4160e-004;s(2)= 8.1827e-004;s(pi)=1; s(pi/2)=-1;s(pi/3)= 0.0039 2. 求的fourier、Laplace、ztrans变换 程序内容： fourier、Laplace、ztrans变换 源程序：\bianhuan.m； 程序清单： syms t; f=5*sin(2*t)-3*cos(2*t) fourier_f=fourier(f) laplace_f=laplace(f) ztrans_f=ztrans(f) 实验结果： f =5*sin(2*t)-3*cos(2*t) fourier_f =pi*(5*i*Dirac(w+2)-3*Dirac(w+2)-5*i*Dirac(w-2)-3*Dirac(w-2)) laplace_f =10/(s^2+4)-3*s/(s^2+4) ztrans_f = 10*z*cos(1)*sin(1)/(-4*z*cos(1)^2+z^2+2*z+1)-3*(z+1-2*cos(1)^2)*z/(-4*z*cos(1)^2+z^2+2*z+1) 3．矩阵高斯消去法 对矩阵进行高斯消去法变换 求三次初等变换的总的等价乘子,并用MATLAB求原来A的行列式、秩和迹 源程序 ：\gaosi.m 程序清单： A=[1 0 7;4 1 5;2 -1 9]; A0=A; %输入A，并保留一个备份 A(2,: ) = -4*A(1,: )+A(2,: ) A1=A, B1=A1/A0 % 消去A(2,1),求B1 A(3,: ) = -2*A(1,: )+A(3,: ) A2=A, B2=A2/A1 % 消去A(3,1) A(3,: ) = -A(3,2)/A(2,2)*A(2,: )+A(3,: ) A3=A, B3=A3/A2 % 消去A(3,2) B0 = A3/A0 % 求三次初等变换的总的等价乘子 det_A=det(A0), rank_a=rank(A0), tr_A=trace(A0), %求原来A的行列式、秩和迹 实验结果： ， det_A = -28 rank_a =3 tr_A =11 4．画出心型线、星型线、双纽线和四叶玫瑰线的图形 心形线绘制： Clear t=0:0.001:2*pi; subplot(2,2,1); polar(a, 1+cos(t)) subplot(2,2,2); plot(cos(t).^3,sin(t).^3) subplot(2,2,3); polar(t,abs(sin(t).*cos(t))) subplot(2,2,4); polar(t,(cos(2*t)).^0.5) 图像： 图表 1四曲线 5．求fibonacci数 图表 2 fib 曲线 实验结果：100以内的fibonacci 数：f = 1 1 2 3 5 8 13 21 34 55 89 实验结果：1000以内的fibonacci数：f= 1 1 2 3 5 8 13 21 34 55 89 144 233 377 610