重庆大学 数学实验报告 机械机自 实验一.docx
- 文档编号:26051328
- 上传时间:2023-06-17
- 格式:DOCX
- 页数:16
- 大小:251.33KB
重庆大学 数学实验报告 机械机自 实验一.docx
《重庆大学 数学实验报告 机械机自 实验一.docx》由会员分享,可在线阅读,更多相关《重庆大学 数学实验报告 机械机自 实验一.docx(16页珍藏版)》请在冰豆网上搜索。
重庆大学数学实验报告机械机自实验一
重庆大学数学实验报告机械2009机自实验一
实验课程名称数学实验开课实验室DS1402
学院2009年级机自专业班6班
学生姓名学号2009
开课时间2010至2011学年第2学期
总成绩
教师签名
数理学院制
开课学院、实验室:
实验时间:
年月日
实验项目类型课程实验项目数学实验数学实验一
验证演示综合设计其他名称名称
指导何仁斌成绩
教师
实验目的
[1]熟悉MATLAB软件的用户环境;
[2]了解MATLAB软件的一般目的命令;
[3]掌握MATLAB数组操作与运算函数;
[4]掌握MATLAB软件的基本绘图命令;
[5]掌握MATLAB语言的几种循环、条件和开关选择结构。
基础实验
一、实验内容
1(MATLAB软件的数组操作及运算练习;
2(直接使用MATLAB软件进行作图练习;
3(用MATLAB语言编写命令M-文件和函数M-文件。
二、实验过程(一般应包括实验原理或问题分析,算法设计、程序、计算、图表等,实
验结果及分析)
ER,,3,33,2A,1(设有分块矩阵,其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵,试通过,,OS2,32,2,,
ER,RS,,2A,数值计算验证。
2,,0S,,
程序:
E=eye(3);R=rand(3,2);
O=zeros(2,3);S=diag(2,1);
A2=[E,R+R*S;O,S*S].^2
结果:
2(某零售店有9种商品的单件进价(元)、售价(元)及一周的销量如表1.1,问哪种商品的利润最大,哪种商品的利润最小;按收入由小到大,列出所有商品及其收入;求这一周该10种商品的总收入和总利润。
表1.1
货号123456789
单件进价7.158.253.2010.306.6812.0316.8517.519.30
单件售价11.1015.006.0016.259.9018.2520.8024.1515.50
销量568120575358039521041538810694
程序:
a=[7.158.253.210.306.6812.0316.8517.519.30];
b=[11.1015.006.0016.259.9018.2520.8024.1515.50];
c=[568120575358039521041538810694];
d=(b-a).*c,
sum=sum(d),
[g,h]=sort(d),
结果:
3.近景图将x的取值范围局限于较小的区间内可以画出函数的近景图,用于显示函数的局部特性。
局部放大在绘图时,把x的范围逐渐缩小,可把函数的细节部分展现的很清楚.特别是观察极限问题时,这种方法比较便利.
远景图函数的远景图,是把x的范围取得比较大,使我们能够在大范围内观察函数图像.当研究x趋向于?
时,这种方法给我们带来方便.
13630y,x,y,x,y,x,y,x1)绘制幂函数在区间[0,2]上的图形。
观察图像,列表记录观察现象。
观察现象
图像经过的关键点
函数图形的增减性
抛物线的开口方向
参数p(指数幂)的影响
程序:
x=linspace(0,2,130);y=x;z=x.^3;m=x.^6;n=x.^30;
subplot(2,2,1),plot(x,y),title('y=x')
subplot(2,2,2),plot(x,z),title('z=x.^3')
subplot(2,2,3),plot(x,m),title('m=x.^6')
subplot(2,2,4),plot(x,n),title('n=x.^30')
结果:
现象:
随着指数的增加,曲线在1之前越来越趋近于零,曲线斜率越来越小,后半段曲线斜率越来越大。
33f(x),x,g(x),x,x,h(x),x2)比较函数在x?
0时函数的性态。
观察到什么现象,从观察到的现象,反映了什么结论。
程序:
x=linspace(0,2,130);p=linspace(0,0.2,130);y=x;z=x.^3+x;m=x.^3;n=p.^3;
subplot(2,2,1),plot(x,y),title('y=x')
subplot(2,2,2),plot(x,z),title('z=x.^3+x')
subplot(2,2,3),plot(x,m),title('m=x.^3')
subplot(2,2,4),plot(p,n),title('n=x.^3趋近于0')
结果:
现象:
越趋近于一,曲线斜率越小,指数越大,曲线斜率越小。
33f(x),x,g(x),x,x,h(x),x3)比较函数在x?
?
时函数的性态。
程序:
x=linspace(100000,100010,1300);
;y=x;z=x.^3+x;m=x.^3;
n=p.^3;
subplot(2,2,1),plot(x,y),title('y=x无穷大')
subplot(2,2,2),plot(x,z),title('z=x.^3+x无穷大')subplot(2,2,3),plot(x,m),title('m=x.^3无穷大')subplot(2,2,4),plot(x,n),title('n=x.^3趋近于无穷大')
结果:
现象:
指数越大,曲线斜率越大
4(用subplot分别在不同的坐标系下作出下列四条曲线,为每幅图形加上标题,
2,xy,e1)概率曲线;
2)四叶玫瑰线,=sin2,;
3t,x,,3,,1,t3)叶形线,23t,y,;3,1,t,
21,1,y2x,ln,1,y4)曳物线。
y
程序:
x=linspace(-pi,pi,130);y1=exp(-x.^2);y2=sin(2*x);
t=linspace(-15,10,100);x3=3.*t./(1+t.^3+eps);y3=3.*t.^2./(1+t.^3+eps);
m=linspace(0,1,100);y4=log((1+sqrt(1-m.^2))./(m+eps))-sqrt(1-m.^2);
y5=log((1+sqrt(1-m.^2))./(m+eps))-sqrt(1-m.^2);
y=[y4;y5];
subplot(2,2,1),plot(x,y1),title('概率曲线')
subplot(2,2,2),polar(x,y2),title('四叶玫瑰线')
subplot(2,2,3),plot(x3,y3),title('叶形线')
subplot(2,2,4),plot(m,y),title('曳物线')
结果:
5(作出下列曲面的3维图形,
22z,sin(,x,y)1);
程序:
x=-7.5:
0.3:
7.5;y=x;
[X,Y]=meshgrid(x,y);
Z=sin(sqrt(X.^2+Y.^2));
mesh(X,Y,Z);
结果:
x,(1,cosu)cosv,,,u,(0,2),2)环面:
y,(1,cosu)sinv,。
v,(0,2,),z,sinu,,
3)分别作出单位球面在参数为两种不同取值范围的图形,注意坐标轴的单位长度要相等。
提示:
附加命令rotate3d可实现3维图形旋转。
xuv,cossin,,u,(0,1.6),,a);yuv,sinsin,,v,(0,),,zv,cos,,
程序:
ezsurfc('(cos(u))*sin(v)','(sin(u))*sin(v)','cos(v)',[0,1.6*pi,0,pi])
axisequal
gridon
xlabel('x轴');ylabel('y轴');zlabel('z轴');
title('a图')
rotate3don
结果:
xuv,cossin,,
yuv,sinsin,u,(0,2),,
zv,cos,v,(0.5,),,,b)
程序:
ezsurfc('(cos(u))*sin(v)','(sin(u))*sin(v)','cos(v)',[0,2*pi,0.5*pi,pi])
axisequal
gridon
xlabel('x轴');ylabel('y轴');zlabel('z轴');
title('a图')
rotate3don
结果:
24)z=y绕z轴的旋转面图形
程序:
x=-10:
0.1:
10;y=x;[X,Y]=meshgrid(x,y);Z=X.^2+Y.^2;mesh(X,Y,Z);
2z5)y=-,0 程序: [X,Z]=meshgrid(0: 0.1: 5); Y=-Z.^2; mesh(X,Y,Z); xlabel('x轴');ylabel('y轴');zlabel('z轴'); 6(建立一个命令M-文件: 求所有的“水仙花数”,所谓“水仙花数”是指一个三位数,其各位数字的立 333方和等于该数本身。 例如,153是一个水仙花数,因为153=1+5+3。 程序: fora=1: 9; forb=0: 9; forc=0: 9; d=a^3+b^3+c^3; ifd==100*a+b*10+c disp(d); end end end end 结果: 7(编写函数M-文件sq.m: 用迭代法求的值。 求平方根的迭代公式为x,a 1ax,x,()n,1nx2n 5,迭代的终止条件为前后两次求出的x的差的绝对值小于10。 程序: a=input('...'); x=1; whileabs(x-0.5*(x+a/x))>10^(-5) x=0.5*(x+a/x); end disp(x) 结果: 8.求函数的极限、导数或积分: 1xxlim(x,3)1)当x时;,, 程序: symsx s=(x+3^x)^(1/x) limit(s,x,inf) 结果: xesinx,x(x,1)lim,x,0;2)3x 程序: symsx s=(exp(x)*sin(x)-x*(x+1))/x^3; limit(s,x,0) 结果: 2xx,,21fx'()fx(),,3)求;,xexsin1, 程序: symsx s=(x^2+2*x-1)/(exp(-x)*sin(x)+1); diff(s) 结果: 2x(n)4)已知(),求;fx,f(0)21,x程序: symsx; f=x^2/(1-x^2); s=diff(f,x,4); x=0 eval(s) dxy225)已知,求;arctg,lnx,ydyx 程序: symsx; symsy; F=log(sqrt(x^2+y^2))-atan(y/x); dy=diff(F,y); dx=diff(F,x); dx/dy 结果: z,zz,xarctgy,求,,6)画函数图;,x,y 程序: symsx; symsy; z=x*atan(y); diff(z,x) diff(z,y) 结果: 2xedx7);x,e,2 程序: symsx; int(exp(2*x)/(exp(x)+2),x) 程序: 439.作出函数y=x-4x+3x+5(x,[0,6])的图形,用小红点标出其在[0,6]之间的最小值点,并在最小值点附近标出该最小值点的坐标值; 程序: x=0: 0.01: 6;y=x.^4-4*x.^3+3*x+5;plot(x,y,'-b') holdon [y,i]=min(y),x=(i-1)*0.01 plot(x,y,'.r') gtext('(2.91,-13.13)') holdoff 总结与体会 通过本次实验,我熟悉了关于MATLAB的几个基本操作,比如基本函数以及绘图命令的运用。 教师签名 年月日备注: 1、同一章的实验作为一个实验项目,每个实验做完后提交电子稿到服务器的“2011全校任选 课数学实验作业提交”文件夹,文件名为“学院学号姓名实验几”,如“机械20073159张 新实验一”。 2、提交的纸质稿要求双面打印,中途提交批改不需要封面,但最后一次需将该课程所有实验 项目内页与封面一起装订成册提交。 3、综合实验要求3人合作完成,请在实验报告上注明合作者的姓名。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 重庆大学 数学实验报告 机械机自 实验一 数学 实验 报告 机械