程序：
function f=shuixianhua_number(>
for i=1:9
for j=0:9
for k=0:9
if i*100+j*10+k==i^3+j^3+k^3
f=i*100+j*10+k
end
end
end
end
结果：
f =
153
f =
370
f =
371
f =
407
ans=
407
7．编写函数M-文件sq.m：用迭代法求 的值。求平方根的迭代公式为
单件售价
11.1015.00 6.0016.259.9018.2520.8024.1515.50
销量
5681205 75358039521041538 810 694
程序：
A=[1 2 3 4 5 6 78 9
7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30
程序：
E=eye(3>。
R=rand(3,2>。
O=zeros(2,3>。
S=diag(rand(1,2>>。
A=[E R。O S]
B=[E R+R*S。O S*S]
H=A*A-B
结果：
A =
1.0000 0 0 0.9501 0.4860
0 1.00000 0.2311 0.8913
0 0 1.0000 0.6068 0.7621
1）概率曲线 ；
2）四叶玫瑰线=sin2；
3）叶形线
4）曳物线 。
程序：
x1=linspace(-4,4,50>。
y1=exp(-x1.^2>。
x2=linspace(-pi,pi,50>。
y2=sin(2*x2>。
t=linspace(0,10,50>。
x3=3*t./(1+t.^3>。
y3=3*t.^2./(1+t.^3>。
结果：ans =
3
3） 求 ；
程序如下所示：syms x
diff((x^2+2*x-1>/(exp(-x>*sin(x>+1>,x>
结果：ans =
(2*x+2>/(exp(-x>*sin(x>+1>-(x^2+2*x-1>/(exp(-x>*sin(x>+1>^2*(-exp(-x>*sin(x>+exp(-x>*cos(x>>
重庆大学
学生实验报告
实验课程名称数学实验
开课实验室DS1401
学院土木项目学院年级2018专业班土木1班
学生姓名李潇潇学号20185453
开课时间至学年第学期
总成绩
教师签名
数学与统计学院制开课学院源自实验室：实验时间：年月日课程
名称
数学实验
实验项目
名称
基础练习
实验项目类型
验证
演示
综合
设计
其他
指导
教师
成绩
[U2,V2]=meshgrid(u2,v2>。
x2=cos(U2>.*sin(V2>。
y2=sin(U2>.*sin(V2>。
z2=cos(V2>。
mesh(x2,y2,z2>
rotate3d
结果：
<4）z =y2绕z轴的旋转面图形
程序：
x=-1:0.05:1。
y=x。
[X,Y]=meshgrid(x,y>。
x=(1+cos(U>>.*cos(V>。
y=(1+cos(U>>.*sin(V>。
z=sin(U>。
mesh(x,y,z>。
结果：
<3）分别作出单位球面在参数为两种不同取值范围的图形,注意坐标轴的单位长度要相等。提示：附加命令rotate3d可实现3维图形旋转。
<a>
程序：
u1=0:0.01:1.6*pi。
v1= 0:0.01:pi。
[U1,V1]=meshgrid(u1,v1>。
x1=cos(U1>.*sin(V1>。
y1=sin(U1>.*sin(V1>。
z1=cos(V1>。
mesh(x1,y1,z1>
rotate3d
结果：
<b>
程序：
u2= 0:0.01:2*pi。
v2=0.5*pi:0.01:pi。
教师签名
年月日
程序：
x=linspace(0,2,50>。
y=[exp(x>。1+x。1+x+1/2*x.^2。1+x+1/2*x.^2+1/6*x.^3]。
plot(x,y>。
grid。
xlabel('x'>。
ylabel('y'>。
text(1.6,exp(1.6>,'\leftarrowy1=e^x'>。
text(1.6,1+1.6,'\leftarrowy2=1+x'>。
y4=linspace(-1,1,50>。
x4=[log((1+sqrt(1-y4.^2>>./y4>-sqrt(1-y4.^2>。log((1-sqrt(1-y4.^2>>./y4>+sqrt(1-y4.^2>]。
subplot(2,2,1>,plot(x1,y1>,title('概率曲线y=e^-x^2'>。
[Mlirun]=max(str>
C=[I/0.0001。str]
zshrun=sum(A(3,:>.*A(4,:>>
zlirun=sum(B(2,:>.*A(4,:>>
结果：
mlirun =
1.2719e+003
Mlirun =
1.3087e+004
C =
1.0e+004 *
5.0000 3.0000 1.0000 4.0000 9.0000 8.0000 7.0000 2.0000 6.0000
0 0 0 0.4565 0
0 0 0 0 0.0185
B =
1.0000 0 0 1.3838 0.4950
0 1.0000 0 0.3366 0.9078
0 0 1.0000 0.8838 0.7762
0 0 0 0.2084 0
0 0 0 0 0.0003
H =
0 0 0 0 0
0 0 0 0 0
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 694]。
B=[A(1,:>。A(3,:>-A(2,:>]。
[str,I]=sort(B(2,:>.*A(4,:>>。
[mlirun]=min(str>
实验目的
[1]熟悉MATLAB软件的用户环境；
[2]了解MATLAB软件的一般目的命令；
[3]掌握MATLAB数组操作与运算函数；
[4]掌握MATLAB软件的基本绘图命令；
[5]掌握MATLAB语言的几种循环、条件和开关选择结构。
基础实验
一、实验内容
1．MATLAB软件的数组操作及运算练习；
2．直接使用MATLAB软件进行作图练习；
text(1.6,1+1.6+1/2*1.6^2,'\leftarrowy3=1+x+(1/2>x^2'>。
text(1.6,1+1.6+1/2*1.6^2+1/6*1.6^3,'\leftarrowy4= 1+x+(1/2>x^2+(1/6>x^3'>。
结果：
4．用subplot分别在不同的坐标系下作出下列四条曲线，为每幅图形加上标题，
5）已知 ，求 ；
程序如下所示：
syms x y
z=atan(y/x>-log(sqrt(x^2+y^2>>。
zx=diff(z,x>。
zy=diff(z,y>。
a=zy/zx
结果：
a =
(1/x/(1+y^2/x^2>-1/(x^2+y^2>*y>/(-y/x^2/(1+y^2/x^2>-1/(x^2+y^2>*x>
4．开启MATLAB编辑窗口,键入你编写的M文件<命令文件或函数文件）；
5．保存文件<注意将文件存入你自己的文件夹）并运行；
6．若出现错误，修改、运行直到输出正确结果；
7．写出实验报告，并浅谈学习心得体会。
应用实验<或综合实验）
1．设有分块矩阵 ，其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证 。
程序：
a=2。
x0=a。
x1=0.5*(x0+a./x0>。
while abs(x1-x0>>10.^-5
x0=x1。
x1=0.5*(x0+a./x0>。
end
x1
结果:
x1 =
1.4142
8.求函数的极限、导数或积分：
1） 当x 时；
程序如下所示：
syms x
limit((x+3^x>^(1/x>,x,inf>
7） 。
程序如下所示:syms x
int(exp(2*x>/(exp(x>+2>>
结果：
ans =
exp(x>-2*log(exp(x>+2>
总结与体会
通过这次实验，我基本掌握了运用MATLAB解决现实生活中问题的思想。学会了怎么样运用solve函数解线性方程组，但是如何将现实生活中的问题转换成数学模型还不熟悉，要多加练习。
r=X.^2+Y.^2。
z=r。
mesh(X,Y,z>。
结果：
<5>y = - ,0<x<5柱面图形
程序：
x=linspace(0,1,500>。
z=linspace(-1,1,500>。
[X,Z]=meshgrid(x,z>。
Y=-Z.^2。
mesh(X,Y,Z>
rotate3d
结果：
6．建立一个命令M-文件：求所有的“水仙花数”，所谓“水仙花数”是指一个三位数，其各位数字的立方和等于该数本身。例如，153是一个水仙花数，因为153=13+53+33。
subplot(2,2,2>,polar(x2,y2>,title('四叶玫瑰线ρ=sin2θ'>。
subplot(2,2,3>,plot(x3,y3>,title('叶形线x=3t/(1+t^3>,y=3t^2/(1+t^3>'>。
subplot(2,2,4>,plot(x4,y4>,title('曳物线x=ln((1±√(1-y^2>>/y>±√(1-y^2>>'>。
3．用MATLAB语言编写命令M-文件和函数M-文件。
二、实验过程<一般应包括实验原理或问题分析，算法设计、程序、计算、图表等，实验结果及分析）
1.在E盘建立一个自己的文件夹。
2．开启软件平台——MATLAB，将你建立的文件夹加入到MATLAB的搜索路径中。
3．利用帮助了解函数max, min, sum, mean, sort, length，rand, size和diag的功能和用法。
结果：
5．作出下列曲面的3维图形，
<1） ；
程序：
x=-2:0.01:2。
y=-2:0.01:2。
[X,Y]=meshgrid(x,y>。
Z=sin(pi*sqrt(X.^2+Y.^2>>。
mesh(X,Y,Z>
结果：
<2）环面： 。
程序：
u=0:0.01:2*pi。
v=0:0.01:2*pi。
[U,V]=meshgrid(u,v>。
0.1272 0.2108 0.2244 0.3451 0.4303 0.5378 0.6075 0.8134 1.3087
zshrun =
1.4294e+005
zlirun =
4.6052e+004
3．<5）在同一个坐标下作出y1=ex,y2=1+x,y3=1+x+(1/2>x2,y4=1+x+(1/2>x2+(1/6>x3这四条曲线的图形，要求在图上加各种标注，观察到什么现象？发现有什么规律？
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
2．某零售店有9种商品的单件进价<元）、售价<元）及一周的销量如表1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该10种商品的总收入和总利润。
表1.1
货号
12 345 67 89
单件进价
7.158.253.2010.30 6.6812.0316.8517.519.30
