1、重 庆 大 学学 生 实 验 报 告实验课程名称 开课实验室 学 院 年级 专业班 学 生 姓 名 学 号 开 课 时 间 至 学年第 学期总 成 绩教师签名数 学 与 统 计 学 院 制开课学院、实验室： 实验时间 ： 年 月 日实验项目类型课程名称实验项目名 称 验证 演示 综合 设计 其他指导教师成 绩实验目的1 熟悉 MATLAB 软件的用户环境；2 了解 MATLAB 软件的一般目的命令；3 掌握 MATLAB 数组操作与运算函数；4 掌握 MATLAB 软件的基本绘图命令；5 掌握 MATLAB 语言的几种循环、条件和开关选择结构。通过该实验的学习，使学生能灵活应用 MATLAB
2、软件解决一些简单问题，能借助 MATLAB 软件的绘图功能，对函数的特性进行探讨，广泛联想，大胆猜想，发现进而证实其中的规律。基础实验一、实验内容1 熟悉 MATLAB 软件的用户环境；2 了解 MATLAB 软件的一般目的命令；3 掌握 MATLAB 数组操作与运算函数；4 掌握 MATLAB 软件的基本绘图命令；5 掌握 MATLAB 语言的几种循环、条件和开关选择结构。通过该实验的学习，使学生能灵活应用 MATLAB 软件解决一些简单问题，能借助 MATLAB 软件的绘图功能，对函数的特性进行探讨，广泛联想，大胆猜想，发现进而证实其中的规律。二、实验过程（一般应包括实验原理或问题分析，变
3、量说明、程序、调试情况记录、图表等， 实验结果及分析）详见下文四、实验结果及分析1设有分块矩阵 ，其中 E,R,O,S 分别为单位阵、随机阵、零阵和对角阵，232SOREA试通过数值计算验证 。20第一题如上图，经过多次随机数测试，结果都为 0 矩阵可知结论成立。2某零售店有 9 种商品的单件进价（元） 、售价（元）及一周的销量如表 1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该 10 种商品的总收入和总利润。表 1.1货号 1 2 3 4 5 6 7 8 9 单件进价 7.15 8.25 3.20 10.30 6.68 12.03 16.8
4、5 17.51 9.30单件售价 11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50销量 568 1205 753 580 395 2104 1538 810 694M 文件：a=7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30;b=11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50;c=568 1205 753 580 395 2104 1538 810 694;d=(b-a)z=d.*csort(z) %求最大利润sort(b.*c) %对
5、收入进行排序sum(z),sum(b.*c) %对收入和利润求和结果：ans =1.0e+004 *0.1272 0.2108 0.2244 0.3451 0.4303 0.5378 0.6075 0.8134 1.3087ans =1.0e+004 *0.3911 0.4518 0.6305 0.9425 1.0757 1.8075 1.9562 3.1990 3.8398ans =4.6052e+004ans =1.4294e+0051）绘制幂函数 30631,xyxy在区间0,2上的图形。观察图像,列表记录观察现象。以下依次为 图像观察可知以下结论2）比较函数 33)(,)(,)( xh
6、xgxf 在 x0 时函数的性态。观察到什么现象？从观察到的现象，反映了什么结论。M 文件：x=0:0.001:0.1;plot(x,x,+,x,x.3+x,x,x.3,*)legend(y=x,y=x3,y=x+x3)结果：观 察 现 象图像经过的关键点 （0,0） (1,1)函数图形的增减性 同为增函数抛物线的开口方向 开口向上参数 p（指数幂）的影响 p 越大函数增长速度越快结论：y=x，y=x+x3近似相等 y=x3接近于03）比较函数 33)(,)(,)( xhxgxf 在 x时函数的性态。M 文件 x=0:100000000:10000000000;plot(x,x,+,x,x.3
7、+x,x,x.3,*)legend(y=x,y=x3,y=x+x3)结果：结论：y=x3，y=x+x3近似相等 y=x 相对逼近04）在日常生活中我们有这样的经验：与幂函数相比，指数函数是急脾气，对数函数是慢性子。这就是说，当 x时，再小的指数函数也比幂函数变化快，再大的对数函数也比幂函数变化慢。当 x时,比较 10xy与 x.的大小.当 x时,比较 01.xy与 xylg0 的大小.M 文件 x=0:1:100;plot(x,x.10,x,1.1.x,+)legend(y=x10,y=1.1x)结果结论：10xy远大于xy1.M 文件：x=0:100:10000;plot(x,x.0.001
8、,x,1000*log10(x),+)legend(y=x0.001,1000log10(x)结论：1000log10(x)远小于 y=x0.0014用 subplot 分别在不同的坐标系下作出下列四条曲线，为每幅图形加上标题，1）概率曲线 ；2xey2）四叶玫瑰线 =sin2；3）叶形线 ;13,2tytx4）曳物线 。221lnyyxmM 文件 ：subplot(2,2,1)x=0:0.01:2;plot(x,exp(-x.2)subplot(2,2,2)o=0:0.01:2*pi;polar(o,sin(2*o)subplot(2,2,3)t=0:10:100000;ezplot(3*t/(1+t.3),3*t.2/(1+t.3)subplot(2,2,4)x=0:0.00001:0.1plot(log(1+sqrt(1-x.2)./x)-sqrt(1-x.2),x,log(1-sqrt(1-x.2)./x)+sqrt(1-x.2),x)5作出下列曲面的 3 维图形，1） ；)sin(2yxzM 文件：x=-10:0.01:10;y=-10:0.01:10;a,b=meshgrid(x,y);mesh(a,b,sin(pi*sqrt(a.2+b.2)下图为程序运行结果
