实验指导书实验一

《实验指导书实验一》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、重庆大学学生实验报告实验课程名称数学建模与数学实验开课实验室学院弘深学院年级2010专业班理工05班学生姓名陈晨学号20106464开课时间2012至2012学年第二学期总成绩教师签名数学与统计学院制开课学院、实验室：实验时间：年月日课程名称数学建模与数学实验实验项目名称实验一实验项目类型验证演示综合设计其他指导教师龚劬成绩实验目的[1]熟悉MATLAB软件的用户环境；[2]了解MATLAB软件的一般目的命令；[3]掌握MATLAB数组操作与运算函数；[4]掌握MATLAB软件的基本绘图命令；[5]掌握MATLAB语言的几种循环、条件和开关选择结构。基础实验一、实验内容及实验过程（一

2、般应包括实验原理或问题分析，算法设计、程序、计算、图表等，实验结果及分析）1．MATLAB软件的数组操作及运算练习；2．直接使用MATLAB软件进行作图练习；3．用MATLAB语言编写命令M-文件和函数M-文件。1．设有分块矩阵，其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证。程序E=eye(3);R=rand(3,2);O=zeros(2,3);V=rand(1,2);S=diag(V);A=[ER;OS];X=A*A;Y=[ER+R*S;OS*S];...if(X==Y)...disp('p=1');...elsedisp('p=0');...end输出

3、结果：p=1得证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.2010.306.6812.0316.8517.519.30];...b=[

4、11.1015.006.0016.259.9018.2520.8024.1515.50];...c=[568120575358039521041538810694];...d=c.*(b-a);e=b.*c;[m,im]=min(d),[M,iM]=max(d),...[rd,ir]=sort(e),sum1=sum(b.*c),sum2=sum(d),运行结果m=1.2719e+003im=5M=1.3087e+004iM=6rd=1.0e+004*0.39110.45180.63050.94251.07571.80751.95623.19903.8398ir=531492876s

5、um1=1.4294e+005sum2=4.6052e+004所以，第五种商品利润最小，第六种商品利润最大，商品按收入排序如上所示，总收入为1.4294e+005，总利润为4.6052e+0044．用subplot分别在不同的坐标系下作出下列四条曲线，为每幅图形加上标题，1）概率曲线；2）四叶玫瑰线r=sin2q；3）叶形线4）曳物线。程序：x1=linspace(-3,3,100);y1=exp(-x1.*x1);subplot(2,2,1),plot(x1,y1),title('y=exp(-x*x)');...r=-pi:0.01:pi;subplot(2,2,2),polar

6、(r,sin(2*r),'-r'),title('p=sin(2*theda)');...t=linspace(0,9,900);x2=3*t./(1+t.^3);y2=3*t.*t./(1+t.^3);subplot(2,2,3),plot(x2,y2),title('x=3*t/(1+t^3),y=3*t*t(1+t^3)');...y3=linspace(-2*pi,2*pi,900);w=sqrt(1-y3.*y3);x3=log((1+w)./y3)-w;x4=log((1-w)./y3)+w;subplot(2,2,4),plot(x3,y3,x4,y3),...titl

7、e('y=log((1+[-]sprt(1-x*x))/x)-[+]sqrt(1-x*x),');Warning:ImaginarypartsofcomplexXand/orYargumentsignored.Warning:ImaginarypartsofcomplexXand/orYargumentsignored.图形如图4.1所示疑问：1）为什么将t=linspace(0,9,900)改为t=linspace(-2,9,900)后，运行结果变成了图