Matlab软件使用实验报告.docx

1、Matlab软件使用实验报告重 庆 大 学学 生 实 验 报 告实验课程名称 数学实验 开课实验室 DS1401 学 院 计算机学院 年级 2009级 专业班 计科一班 学 生 姓 名 曹柳 学 号 20095338 开 课 时 间 2010 至 2011 学年第 2 学期总 成 绩教师签名数 理 学 院 制开课学院、实验室：数理学院、DS1402 实验时间 ： 2011年 3 月 22 日课程名称数学实验实验项目名 称MATLAB软件入门实验项目类型验证演示综合设计其他指导教师李东成 绩实验目的1 熟悉MATLAB软件的用户环境；2 了解MATLAB软件的一般目的命令；3 掌握MATLAB数

2、组操作与运算函数；4 掌握MATLAB软件的基本绘图命令；5 掌握MATLAB语言的几种循环、条件和开关选择结构。通过该实验的学习，使学生能灵活应用MATLAB软件解决一些简单问题，能借助MATLAB软件的绘图功能，对函数的特性进行探讨，广泛联想，大胆猜想，发现进而证实其中的规律。基础实验一、实验内容1MATLAB软件的数组操作及运算练习； 2直接使用MATLAB软件进行作图练习； 3用MATLAB语言编写命令M-文件和函数M-文件。二、实验过程（一般应包括实验原理或问题分析，算法设计、程序、计算、图表等， 实验结果及分析） 第一题：实验分析：第一步：应先将矩阵、分别列出，之后运用矩阵的拼接求

3、出矩阵的值，将矩阵求平方。第二步：同样先求出矩阵、，拼接，最后验证题设。程序如下：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*SD=A2if D=C Correct!else ErrorEnd第二题 程序如下：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 120

4、5 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_goodSort_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

5、+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 offxlabel(x)ylabel(y)title(泰勒公式)text(-1,exp(-1),leftarrowexp(x)text(1,2,leftarrow1+x)text(1.5,1+1.5+1.52/2,leftarrow1+x+x.2/2)text(2,1+2+22/2+23/6,1+x+x.2/2+x.3/6rightarrow,HorizontalAlignment,right)

6、现象：三条曲线的平滑程度和弯曲程度相差不大结论：在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)+sqrt(1-x.2);y4_2=log(1+sqrt(1-x.2)./x)-sqrt(1-x.2);y4_3=log(1-sqrt(1-x.2)./x)+sq

7、rt(1-x.2);y4_4=log(1-sqrt(1-x.2)./x)-sqrt(1-x.2);subplot(2,2,1)plot(x,y1)title(概率曲线)subplot(2,2,2)plot(x,ruo)title(四叶玫瑰线)subplot(2,2,3)plot(x3,y3)axis(-5 5 -5 5)title(叶形线)subplot(2,2,4)plot(x,y4_1)hold onplot(x,y4_2)hold onplot(x,y4_3)hold onplot(x,y4_4)title(曳物线)hold off第五题程序如下： x=-1:0.1:1y=x X,Y=m

8、eshgrid(x,y); R=sqrt(X.2+Y.2); Z=sin(pi*R); surf(X, Y, Z); u=0:0.1*pi:2*pi; v=0:0.1*pi:2*pi; U,V=meshgrid(u,v); x=(1+cos(U).*cos(V); y=(1+cos(U).*sin(V); z=sin(U); figure mesh(x, y, z);第六题程序如下： for i=1:999sum=0;num=i;while num sum=sum+(mod(num,10)3; num=fix(num/10);endif sum=i i endend第七题程序如下： funct

9、ion EXP_sq_1_7(a)x(1)=a;x(2)=(x(1)+a/x(1)/2;i=2;while (abs(x(i)-x(i-1)=(10(-5) i=i+1 x(i)=(x(i-1)+a/x(i-1)/2;endx(i)第八题程序如下： syms x;%1题limit(x+3x)(1/x),x,inf)%2题limit(exp(x)*sin(x)-x*(x+1)/(x3),x,0)%3题f=(x2+2*x-1)/(exp(-x)*sin(x)+1);diff(f)%4题f1=x2/(1-x2);% ff1=diff(f);% ff2=diff(f,2);% ff3=diff(f,3

10、)% eval(ff1)% eval(ff2)% eval(ff3)for i=1:10 ff=diff(f,i); x=0; eval(ff)end%5题syms x y;f=atan(y/x)-log(sqrt(x2+y2)dfx=diff(f,x)dxf=1/dfx;dy=diff(f,y)dxf*dy%6题z=x*atan(y)dzx=diff(z,x)dzy=diff(z,y)% x=1;% y=1;% eval(dzy)X=-10:0.1:10;Y=X;x,y=meshgrid(X,Y);z1=eval(dzx);z2=x./(1+y.2);%计算中有向量的运算，先计算导函数，在求

11、解，按z1的方式处理有错mesh(x,y,z1)figuremesh(x,y,z2)%7题syms x;int(exp(2*x)/(exp(x)+2),x)第九题程序如下： function EXP_9x=0:0.1:6y=fun(x);plot(x,y);y_min,in_min=min(y);hold on;plot(x(in_min),y(in_min),r.)text(x(in_min)+0.1,y(in_min)-12,(,num2str(x(in_min),num2str(y(in_min),)return function y=fun(x)y=x.4-4*x.3+3*x+5;re

12、turn 结果：x = Columns 1 through 14 0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 1.0000 1.1000 1.2000 1.3000 Columns 15 through 28 1.4000 1.5000 1.6000 1.7000 1.8000 1.9000 2.0000 2.1000 2.2000 2.3000 2.4000 2.5000 2.6000 2.7000 Columns 29 through 42 2.8000 2.9000 3.0000 3.1000 3.20

13、00 3.3000 3.4000 3.5000 3.6000 3.7000 3.8000 3.9000 4.0000 4.1000 Columns 43 through 56 4.2000 4.3000 4.4000 4.5000 4.6000 4.7000 4.8000 4.9000 5.0000 5.1000 5.2000 5.3000 5.4000 5.5000 Columns 57 through 61 5.6000 5.7000 5.8000 5.9000 6.0000总结与体会： 通过实验一，我初步掌握了MATLAB软件的用户环境；了解了MATLAB软件的一般目的命令；并掌握MATLAB数组操作与运算函数以及基本的绘图命令，还有几种循环等结构。同时，我深切地体会到，只有多动手，多动脑，才能真正学到知识，程序是经过不断地折腾才提炼出来的，学数学实验要有耐心。附图如下 3题图 4题图 5题图(1) 5题图(2) 8题图(1) 8题图(2) 9题图 教师签名年 月 日