实验指导书Matlab软件应用与开发.docx
《实验指导书Matlab软件应用与开发.docx》由会员分享，可在线阅读，更多相关《实验指导书Matlab软件应用与开发.docx（14页珍藏版）》请在冰豆网上搜索。
实验指导书Matlab软件应用与开发
《Matlab软件应用与开发》实验指导书
实验1Matlab操作基础及矩阵运算
（一）、实验类型：
验证型
（二）、实验类别：
基础实验
（三）、每组人数：
1
（四）、实验要求：
选修
（五）、实验学时：
3个学时
（六）、实验目的：
（1）熟悉MATLAB软件中关于矩阵运算的各种命令；
（2）学会运用MATLAB软件自定义函数，并求出函数值；
（3）学会在MATLAB环境下编写函数。
（七）、预备知识：
线性代数中的矩阵运算；高等数学中微积分知识。
本实验所用MATLAB命令
●矩阵输入格式：
A=[a11a12;a21a22];b=初始值:
步长:
终值
●求A的转置：
A’
●求A加B:
A+B
●求A减B:
A-B
●求A乘B:
A*B
●求A的行列式：
det（A）
●求A的逆:
inv（A）
●求A的秩:
rank（A）
●求函数的极限limit（.）
●求函数的导数diff（.）
●求函数的积分int（.）
●求代数方程的解solve（.）
●求微分方程的解dsolve（.）
（八）、内容与要求：
1、输入矩阵A,B,b;
，
2、矩阵转置、四则运算。
C1=A’，C2=A+B,C3=A-B,C4=A*B
3、求行列式。
D1=|A|,D2=|B|
4、求矩阵A、B的秩E1,E2
5、求极限
6、设
，求
7、求
8、求微分方程的通解
9、自定义函数
，并计算
。
（九）、实验操作
1、>>A=[2124854;2351534;21353154;21721552]
A=
2124854
2351534
21353154
21721552
>>B=[1245124;18725335;4811535;46562523]
B=
1245124
18725335
4811535
46562523
>>b=21:
2:
27
b=
21232527
2、>>C1=A',C2=A+B,C3=A-B,C4=A*B
C1=
2122121
24353572
85153115
4345452
C2=
33698628
201076869
69364689
671284075
C3=
9-2184-20
-16-37-38-1
-27341619
-2516-1029
C4=
4948298226684411
2938452929322580
4854652036914056
4660905653624745
3、>>D1=det（A）,D2=det（B）
D1=
2181568
D2=
-3182276
4、>>E1=rank（A）,E2=rank（B）;
E1=
4
E2=
4
5、>>symsx;%定义符号变量x
>>limit（（（2*x+3）/（2*x+1））^（x+1）,x,inf）%求函数（（2*x+3）/（2*x+1））^（x+1）当x->inf时的极限（inf即英文infinity“无穷”的缩写）inf表示正无穷大
ans=
exp
（1）
6、
>>symsyx%定义符号变量
>>y=x*exp（x）;%定义符号函数
>>diff（y,5）%计算符号函数的五阶导数
ans=
5*exp（x）+x*exp（x）
7、
>>symsxs%定义符号变量
>>s=（log（x））/sqrt（x）%定义符号表达式
>>int（s,1,4）%计算符号表达式在区间[1，4]上的定积分
ans=
8*log
（2）-4
8、
>>[y]=dsolve（'x*Dy+y=x^2+3*x+2','x'）%微分或导数的输入是用Dy、D2y、D3y、…来表示y的一阶导数
或
、二阶导数
或
、三阶导数
或
、…。
如果自变量是系统默认的，则自变量输入部分可省略。
y=
2+1/3*x^2+3/2*x+1/x*C1
9、>>
建立函数文件：
fun.m
functionf=fun（x）
f=100*（x
（2）-x
（1）^2）^2+（1-x
（1））^2;
建立主程序：
main.m
x=[12]
fun（x）
。
。
。
。
。
实验二、Matlab程序设计
（一）、实验类型：
综合型
（二）、实验类别：
基础实验
（三）、每组人数：
1
（四）、实验要求：
选修
（五）、实验学时：
3个学时
（六）、实验目的：
学会运用Matlab软件解决线性规划中的实际问题
（七）、实验内容：
任务分配问题：
某车间有甲、乙两台机床，可用于加工三种工件。
假定这两台车床的可用台时数分别为800和900，三种工件的数量分别为400、600和500，且已知用这两种车床加工单位数量不同工件所需的台时数和加工费用如下表。
问怎样分配车床的加工任务，才能既满足加工工件的要求，又使加工费用最低？
工件1
工件2
工件3
甲
加工费：
13
台时数：
0.4
加工费：
9
台时数：
1.1
加工费：
10
台时数：
1
乙
加工费：
11
台时数：
0.5
加工费：
12
台时数：
1.2
加工费：
8
台时数：
1.3
（八）、实验解答：
设在甲车床上加工工件1、2、3的数量分别为x1、x2、x3，在乙车床上加工工件1、2、3的数量分别为x4、x5、x6。
可建立以下线性规划模型：
S.t.
改写为：
源程序为：
f=[1391011128];
A=[0.41.11000；
0000.51.21.3];
b=[800;900];
Aeq=[100100；
010010；
001001];
beq=[400600500];
vlb=zeros（6,1）;
vub=[];
[x,fval]=linprog（f,A,b,Aeq,beq,vlb,vub）
结果:
x=
0.0000
600.0000
0.0000
400.0000
0.0000
500.0000
fval=1.3800e+004
即在甲机床上加工600个工件2,在乙机床上加工400个工件1、500个工件3，可在满足条件的情况下使总加工费最小为13800。
实验三、Matlab绘图与数据模拟
（一）实验类型：
综合型
（二）实验类别：
基础实验
（三）每组人数：
1
（四）实验要求：
选修
（五）实验学时：
3个学时
（六）实验目的：
（1）掌握Matlab软件中的数据拟合
（2）熟练运用Matlab软件进行绘图
（七）实验内容：
（1）对下面一组数据作二次多项式拟合。
要求：
Ⅰ.返回多项式的系数;Ⅱ.作出数据点的图形及拟合曲线;
Ⅲ.分析拟合的效果。
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
-0.447
1.978
3.28
6.16
7.08
7.34
7.66
9.58
9.48
9.30
11.2
（2）绘制如下饼状图形
标签‘Math',‘English',‘Chinese',‘Music'代表的数值分别为2，4，6，8；
（3）用两种方法画出函数sin（x），
的图形。
要求：
Ⅰ.利用相关函数加x轴、y轴注解分别为“自变量X”、“函数Y”;
Ⅱ.加图形标题“示意图”;Ⅲ.加栅格。
（八）实验解答：
1、解答：
输入以下命令：
x=0:
0.1:
1;
y=[-0.4471.9783.286.167.087.347.669.569.489.3011.2];
A=polyfit（x,y,2）
z=polyval（A,x）;
plot（x,y,'k+',x,z,'r'）%作出数据点和拟合曲线的图形
2、解答：
>>x=[2,4,6,8];
>>pie（x,{'math','English','Chinese','music'}）
3、解答：
>>x=linspace（0,2*pi,30）;
>>y=sin（x）;
>>plot（x,y）
>>xlabel（'自变量x'）
>>ylabel（'函数y'）
>>title（'示意图'）
>>gridon
升级会员 