《2022年重庆大学__数学实验评测报告___实验一基础实验》由会员分享,可在线阅读,更多相关《2022年重庆大学__数学实验评测报告___实验一基础实验(11页珍藏版)》请在金锄头文库上搜索。
1、个人资料整理仅限学习使用1 / 11 重 庆 大 学学 生 实 验 报 告实验课程名称数学实验开课实验室 DS1401 学院 土木工程学院年级2018 专业班土木1 班学 生 姓 名 李潇潇学号 20185453 开 课 时 间 至学年第学期总 成 绩教师签名数学 与 统 计学院制开课学院、实验室:实验时间:年月日课程名称数学实验实验工程名称基础练习实验工程类型验证演示综合设计其他指导教师成绩精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 11 页个人资料整理仅限学习使用2 / 11 实验目的1 熟悉 MATLAB 软件的用户环境;2
2、了解 MATLAB 软件的一般目的命令;3 掌握 MATLAB 数组操作与运算函数;4 掌握 MATLAB 软件的基本绘图命令;5 掌握 MATLAB 语言的几种循环、条件和开关选择结构。基础实验一、实验内容1MATLAB 软件的数组操作及运算练习; 2直接使用MATLAB 软件进行作图练习; 3用 MATLAB 语言编写命令M-文件和函数M-文件。二、实验过程 一般应包括实验原理或问题分析,算法设计、程序、计算、图表等,实验结果及分析)1. 在 E盘建立一个自己的文件夹。2开启软件平台MATLAB ,将你建立的文件夹加入到MATLAB 的搜索路径中。3利用帮助了解函数max, min, su
3、m, mean, sort, length,rand, size和 diag 的功能和用法。4开启 MATLAB 编辑窗口 , 键入你编写的M文件 命令文件或函数文件);5保存文件 注意将文件存入你自己的文件夹)并运行;6若出现错误,修改、运行直到输出正确结果;7写出实验报告,并浅谈学习心得体会。应用实验 。R=rand(3,2 。O=zeros(2,3 。S=diag(rand(1,2。A=E R。O S B=E R+R*S。O S*S H=A*A-B 结果:A = 1.0000 0 0 0.9501 0.4860 0 1.0000 0 0.2311 0.8913 0 0 1.0000 0.
4、6068 0.7621 0 0 0 0.4565 0 0 0 0 0 0.0185 B = 1.0000 0 0 1.3838 0.4950 0 1.0000 0 0.3366 0.9078 0 0 1.0000 0.8838 0.7762 0 0 0 0.2084 0 0 0 0 0 0.0003 H = 0 0 0 0 0 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共 11 页个人资料整理仅限学习使用3 / 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2某零售店有9 种商品的单件进价元)、
5、售价 。A(3,:-A(2,:。str,I=sort(B(2,:.*A(4,:。mlirun=min(str Mlirun=max(str C=I/0.0001。str zshrun=sum(A(3,:.*A(4,: zlirun=sum(B(2,:.*A(4,: 结果:mlirun = 1.2719e+003 Mlirun = 1.3087e+004 C = 1.0e+004 * 5.0000 3.0000 1.0000 4.0000 9.0000 8.0000 7.0000 2.0000 6.0000 0.1272 0.2108 0.2244 0.3451 0.4303 0.5378 0.
6、6075 0.8134 1.3087 zshrun = 1.4294e+005 zlirun = 4.6052e+004 3x2,y4= 1+x+(1/2x2+(1/6x3这四条曲线的图形,要求在图上加各种标注,观察到什么现象?发现有什么规律?程序:x=linspace(0,2,50。y=exp(x 。1+x。1+x+1/2*x.2。1+x+1/2*x.2+1/6*x.3。plot(x,y。grid 。xlabel(x。ylabel(y。text(1.6,exp(1.6,leftarrowy1=ex。text(1.6,1+1.6,leftarrowy2=1+x。text(1.6,1+1.6+1
7、/2*1.62,leftarrowy3=1+x+(1/2x2。text(1.6,1+1.6+1/2*1.62+1/6*1.63,leftarrowy4= 1+x+(1/2x2+(1/6x3。结果:精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 11 页个人资料整理仅限学习使用4 / 11 4用 subplot分别在不同的坐标系下作出下列四条曲线,为每幅图形加上标题,1)概率曲线2xey;2)四叶玫瑰线=sin2;3)叶形线;13,13323ttyttx4)曳物线22111lnyyyx。程序:x1=linspace(-4,4,50。y1=
8、exp(-x1.2。x2=linspace(-pi,pi,50。y2=sin(2*x2。t=linspace(0,10,50。x3=3*t./(1+t.3。y3=3*t.2./(1+t.3。y4=linspace(-1,1,50。x4=log(1+sqrt(1-y4.2./y4-sqrt(1-y4.2。log(1-sqrt(1-y4.2./y4+sqrt(1-y4.2。subplot(2,2,1,plot(x1,y1,title(概率曲线 y=e-x2。subplot(2,2,2,polar(x2,y2,title(四叶玫瑰线=sin2 。subplot(2,2,3,plot(x3,y3,ti
9、tle(叶形线 x=3t/(1+t3,y=3t2/(1+t3。subplot(2,2,4,plot(x4,y4,title(曳物线 x=ln(1(1-y2/y (1-y2。结果:精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 11 页个人资料整理仅限学习使用5 / 11 5作出下列曲面的3 维图形,。Z=sin(pi*sqrt(X.2+Y.2。mesh(X,Y,Z 结果:。x=(1+cos(U.*cos(V。y=(1+cos(U.*sin(V。z=sin(U 。mesh(x,y,z 。精选学习资料 - - - - - - - - - 名
10、师归纳总结 - - - - - - -第 5 页,共 11 页个人资料整理仅限学习使用6 / 11 结果:3)分别作出单位球面在参数为两种不同取值范围的图形, 注意坐标轴的单位长度要相等。提示:附加命令rotate3d可实现 3 维图形旋转。 cos sin ,sinsin ,cos ,xuvyuvzv(0,1.6)(0,)uv程序:u1=0:0.01:1.6*pi。v1= 0:0.01:pi。U1,V1=meshgrid(u1,v1。x1=cos(U1.*sin(V1。y1=sin(U1.*sin(V1。z1=cos(V1 。mesh(x1,y1,z1 rotate3d 结果:精选学习资料
11、- - - - - - - - - 名师归纳总结 - - - - - - -第 6 页,共 11 页个人资料整理仅限学习使用7 / 11 cos sin ,sinsin ,cos ,xuvyuvzv(0,2 )(0.5,)uv程序:u2= 0:0.01:2*pi。v2=0.5*pi:0.01:pi。U2,V2=meshgrid(u2,v2。x2=cos(U2.*sin(V2。y2=sin(U2.*sin(V2。z2=cos(V2 。mesh(x2,y2,z2 rotate3d 结果:精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 11
12、页个人资料整理仅限学习使用8 / 11 。r=X.2+Y.2。z=r 。mesh(X,Y,z。结果: y = -2z,0x。z=linspace(-1,1,500。X,Z=meshgrid(x,z。Y=-Z.2。mesh(X,Y,Z rotate3d 结果:6建立一个命令M-文件:求所有的“水仙花数”,所谓“水仙花数”是指一个三位数,其各位数字的立方和等于该数本身。例如, 153 是一个水仙花数,因为153=13+53+33。程序:function f=shuixianhua_number( for i=1:9 for j=0:9 for k=0:9 if i*100+j*10+k=i3+j3
13、+k3 f=i*100+j*10+k end end end end 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页,共 11 页个人资料整理仅限学习使用10 / 11 结果:f = 153 f = 370 f = 371 f = 407 ans = 407 7编写函数M-文件 sq.m:用迭代法求ax的值。求平方根的迭代公式为)a(211nnnxxx程序:a=2。 x0=a 。 x1=0.5*(x0+a./x0。 while abs(x1-x010.-5 x0=x1。 x1=0.5*(x0+a./x0。 end x1 结果 : x1 =
14、1.4142 8. 求函数的极限、导数或积分:1)xxx1)3lim(当 x时;程序如下所示:syms x limit(x+3x(1/x,x,inf 结果: ans = 3 3)221( ),sin1xxxf xex求( )fx;程序如下所示:syms x diff(x2+2*x-1/(exp(-x*sin(x+1,x 结果: ans = (2*x+2/(exp(-x*sin(x+1-(x2+2*x-1/(exp(-x*sin(x+12*(-exp(-x*sin(x+exp(-x*cos(x 5)已知22lnyxxyarctg,求dydx;程序如下所示:syms x y 精选学习资料 - -
15、- - - - - - - 名师归纳总结 - - - - - - -第 10 页,共 11 页个人资料整理仅限学习使用11 / 11 z=atan(y/x-log(sqrt(x2+y2。zx=diff(z,x。zy=diff(z,y。a=zy/zx 结果:a = (1/x/(1+y2/x2-1/(x2+y2*y/(-y/x2/(1+y2/x2-1/(x2+y2*x 7)dxeexx22。程序如下所示 :syms x int(exp(2*x/(exp(x+2 结果:ans = exp(x-2*log(exp(x+2 总结与体会通过这次实验,我基本掌握了运用MATLAB解决现实生活中问题的思想。学会了怎么样运用solve函数解线性方程组,但是如何将现实生活中的问题转换成数学模型还不熟悉,要多加练习。教师签名年 月 日精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页,共 11 页