MATLAB软件入门 第一次数学实验

-0.9960 -0.2624 0.7736 0.9541 0.4971 -0.6833
-0.6321 -0.9685 -0.4943 -0.2902 -0.6833 -0.9915
图片：
5．建立一个命令M-文件：求所有的“水仙花数”，所谓“水仙花数”是指一个三位数，其各位数字的立方和等于该数本身。例如，153是一个水仙花数，因为153=13+53+33。
>> R=rand(3,2);
>> O=zeros(2,3);
>> a=[2 3];
>> S=diag(a);
>> A=[E R;O S];
>> B=[E R+R*S;O S^2];
>> A^2-B
结果：
ans =
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
x=linspace(-1,1,20);
y1=x.^2;
y2=x.^3;
y3=x.^4;
y4=x.^5;
plot(x,y1,x,y2,x,y3,x,y4)
gtext('y1=x^2');
gtext('y2=x^3');
gtext('y3=x^4');
gtext('y4=x^5')
xlabel('x')
结果：
>> ht3
ans =
-0.0187 -0.8474 -0.9951 -0.9628 -0.9960 -0.6321
-0.8474 -0.7842 0.0060 0.2656 -0.2624 -0.9685
-0.9951 0.0060 0.9687 0.9684 0.7736 -0.4943
-0.9628 0.2656 0.9684 0.3893 0.9541 -0.2902
text(-0.9,-0.1,'y3=x^4\rightarrow')
text(-0.8,-0.6,'\leftarrowy4=x^5')
xlabel('x')
ylabel('y')
title('curves')
结果：
>> ht0
ans =
Columns 1 through 8
-1.0000 -0.8947 -0.7895 -0.6842 -0.5789 -0.4737 -0.3684 -0.2632
li =
5 3 1 4 9 2 8 7 6
zs =
1.4294e+005
zl =
4.6052e+004
ans =
1.4294e+005
3.用两种方法在同一个坐标下作出y1=x2,y2=x3,y3=x4y4=x5这四条曲线的图形，并要求用两种方法在图上加各种标注。
（1）M文件程序：
function[x,y]=ht0(x)
2．直接使用MATLAB软件进行作图练习；
3．用MATLAB语言编写命令M-文件和函数M-文件。
二、实验过程（一般应包括实验原理或问题分析，算法设计、程序、计算、图表等，实验结果及分析）
1．设有分块矩阵 ，其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证 。
>> E=eye(3);
迭代的终止条件为前后两次求出的x的差的绝对值小于105。
M文件程序：
functionx1=sq(a)
x0=eps;
x1=(1/2)*(x0+a/x0);
while(abs(x1-x0)>=10^(-5))
x0=x1;
x1=(1/2)*(x0+a/x0);
end
结果:
>> sq(4)
ans =
2.0000
注行距：选最小值16磅，每一图应有简短确切的题名，连同图号置于图下。每一表应有简短确切的题名，连同表号置于表上。图表的题名及其中的文字采用小5号宋体。
教师签名
年月日
备注：
1、同一章的实验作为一个实验项目.
2、提交的纸质稿要求双面打印，中途提交批改不需要封面，但最后一次需将该课程所有实验项目内页与封面一起装订成册提交。
x=linspace(-1,1,20);
y1=x.^2;
y2=x.^3;
y3=x.^4;
y4=x.^5;
plot(x,y1)
holdon,plot(x,y2);plot(x,y3);plot(x,y4);holdoff
text(-0.6,0.4,'\leftarrowy1=x^2')
text(-0.9,0.2,'y2=x^3\rightarrow')
zl=sum((b-a).*c)
结果：
ht
m =
6.7500
mi =
2
n =
2.8000
ni =
3
l =
1.0e+004 *
Columns 1 through 8
0.3911 0.4518 0.6305 0.9425 1.0757 1.8075 1.9562 3.1990
Column 9
3.8398
M文件程序：
functionf=fun(t)
fort=100:999
x=fix(t/100);
y=fix((t-x*100)/10);
z=t-x*y*10;
if(x^3+y^3+z^3)==t
A=t
end
end
结果:
>> fun
A =
153
A =
370
A =
371
A =
407
6．编写函数M-文件sq.m：用迭代法求 的值。求平方根的迭代公式为
肖剑
成绩
实验目的
[1]熟悉MATLAB软件的用户环境；
[2]了解MATLAB软件的一般目的命令；
[3]掌握MATLAB数组操作与运算函数；
[4]掌握MATLAB软件的基本绘图命令；
[5]掌握MATLAB语言的几种循环、条件和开关选择结构。
基础实验
一、实验内容
1．MATLAB软件的数组操作及运算练习；
重庆大学
学生实验报告
实验课程名称数学实验
开课实验室
学院年级专业班
学生姓名学号
开课时间2011至2012学年第2学期
总成绩
教师签名
数学与统计学院制
开课学院、实验室：数学与统计DS1422实验时间：2012年3月29日
课程
名称
数学实验
实验项目
名称
MATLAB软件入门
实验项目类型
验证
演示
综合
设计
其他
指导
教师
ylabel('y')
title('curves')
图片：
4．作出下列曲面的3维图形，
；
M文件程序：
function[Z]=ht3(x,y)
x=-1.41:0.5:1.41;y=x;
[X,Y]=meshgrid(x,y);
R=sqrt(X.^2+Y.^2)+eps;
Z=sin(R*pi);
mesh(X,Y,Z);
Columns 9 through 16
-0.1579 -0.0526 0.0526 0.1579 0.2632 0.3684 0.4737 0.5789
Columns 17 through 20
0.6842 0.7895 0.8947 1.0000
图片：
（2）M文件程序：
function[x,y]=ht2(x)
b=[11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50];
c=[568 1205 753 580 395 2104 1538 810 694];
[m,mi]=max(b-a)
[n,ni]=min(b-a)
[l,li]=sort(b.*c)
zs=sum(b.*c)
2．某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如表1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该10种商品的总收入和总利润。
表1.1
货号
1 2 3 4 5 6 7 8 9
单件进价
7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30
单件售价
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
M文件程序：
function[zs,zl]=ht(t)
a=[7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30];
总结与体会
通过这次的实验操作让自己对MATLAB的工作环境有了一个大概了解，并通过一定的练习和看书，是自己对MATLAB的基本函数有了了解，并能进行使用。但在练习与操作过程中与暴露自己的不足，这表明自己在有些方面还是存在缺陷的，自己要注意看书与练习，弥补自己的不足之处，只有这样才能使自己有进一步的提高。加强练习，孰能生巧。