Matlab软件使用实验报告.docVIP

重 庆 大 学 学 生 实 验 报 告 实验课程名称 数学实验 开课实验室 DS1401 学 院 计算机学院 年级 2009级 专业班 计科一班 学 生 姓 名 曹柳 学 号 开 课 时 间 2010 至 2011 学年第 2 学期 总 成 绩 教师签名 数 理 学 院 制 开课学院、实验室：数理学院、DS1402 实验时间 ： 2011年 3 月 22 日 课程 名称 数学实验 实验项目 名 称 MATLAB软件入门 实验项目类型 验证 演示 综合 设计 其他 指导 教师 李东 成 绩 √ 实验目的 [1] 熟悉MATLAB软件的用户环境； [2] 了解MATLAB软件的一般目的命令； [3] 掌握MATLAB数组操作与运算函数； [4] 掌握MATLAB软件的基本绘图命令； [5] 掌握MATLAB语言的几种循环、条件和开关选择结构。 通过该实验的学习，使学生能灵活应用MATLAB软件解决一些简单问题，能借助MATLAB软件的绘图功能，对函数的特性进行探讨，广泛联想，大胆猜想，发现进而证实其中的规律。 基础实验 一、实验内容 1．MATLAB软件的数组操作及运算练习； 2．直接使用MATLAB软件进行作图练习； 3．用MATLAB语言编写命令M-文件和函数M-文件。 二、实验过程（一般应包括实验原理或问题分析，算法设计、程序、计算、图表等， 实验结果及分析） 第一题： 实验分析：第一步：应先将矩阵Ｅ、Ｒ、Ｏ、Ｓ分别列出，之后运用矩阵的拼接求出矩阵Ａ的值，将矩阵Ａ求平方。第二步：同样先求出矩阵Ｒ＋Ｒ＊Ｓ、Ｓ＾２， 拼接，最后验证题设。 程序如下：E=eye(3) R=rand(3,2) O=zeros(2,3) temp=rand(1,2) S=diag(temp) A= [E R; O S]; C=[E R+R*S; zeros(2,3) S*S] D=A^2 if D==C Correct! else Error End 第二题 程序如下：cost=[7.15 8.25 3.2 10.3 6.68 12.03 16.85 17.51 9.2]; price=[11.1 15 6 16.25 9.9 18.25 20.8 24.15 15.5]; quantity_of_sale=[568 1205 753 580 395 2104 1538 810 694]; profit=(price-cost).*quantity_of_sale; [min_of_profit min_of_profit_good]=min(profit); 最小利润是： min_of_profit 最小利润的商品是： min_of_profit_good [Sort_profit Sort_profit_good]=sort(profit) 总收入： sum(price.*quantity_of_sale) 总利润 sum(profit) 第三题 程序如下：x=linspace(-2,2,200); y1=exp(x); y2=1+x; y3=1+x+x.^2/2; y4=1+x+x.^2/2+x.^3/6; plot(x,y1) hold on plot(x,y2,r) hold on plot(x,y3,g) hold on plot(x,y4,k) hold off xlabel(x) ylabel(y) title(泰勒公式) text(-1,exp(-1),\leftarrowexp(x)) text(1,2,\leftarrow1+x) text(1.5,1+1.5+1.5^2/2,\leftarrow1+x+x.^2/2) text(2,1+2+2^2/2+2^3/6,1+x+x.^2/2+x.^3/6\rightarrow,HorizontalAlignment,right) 现象：三条曲线的平滑程度和弯曲程度相差不大 结论：在y=1+x后面加上一些x的幂的值对原式不产生很大的影响，可忽略 第四题 程序如下：x=linspace(-pi,pi,200); t=linspace(-pi/2,pi/2,100) t_x=tan(t); y1=exp(-x.^2); ruo=sin(2*x); x3=(3*t_x./(1+t_x.^3)) y3=(3*t_x.^2./(1+t_x.^3)) y4_1=log((1+sqrt(1-x.^2))./x