答案详解实验指导书Matlab软件应用与开发图文精.docx
《答案详解实验指导书Matlab软件应用与开发图文精.docx》由会员分享,可在线阅读,更多相关《答案详解实验指导书Matlab软件应用与开发图文精.docx(12页珍藏版)》请在冰豆网上搜索。
答案详解实验指导书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;
⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=5215722154313521
34153524852421
A,⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎣⎡=232556463515148355372182414512B
[]27252321=b
2、矩阵转置、四则运算。
C1=A’,C2=A+B,C3=A-B,C4=A*B
3、求行列式。
D1=|A|,D2=|B|
4、求矩阵A、B的秩E1,E2
5、求极限11232lim+∞→⎪⎭⎫⎝⎛++xxxx
6、设xxey=,求5(y
7、求dxxx⎰4
1ln
8、求微分方程的通解232++=+'xxyyx
9、自定义函数2122
12211((100,(xxxxxf-+-=,并计算2,1(f。
(九、实验操作
1、>>A=[2124854;2351534;21353154;21721552]
A=
21
248542
35153421
35315421721552
>>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的
一阶导数dxdy或y'、二阶导数xdyd22或y''、三阶导数xdy
d33或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的数量分别为x1、x2、x3,在乙车床上加工工件1、2、3的数量分别为x4、x5、x6。
可建立以下线性规划模型:
源程序为:
f=[1391011128];A=[0.41.11000;0000.51.21.3];b=[800;900];
6543218121110913minxxxxxxz+++++=
⎪⎪⎪⎪
⎩⎪⎪⎪⎪
⎨
⎧=≥≤++≤++=+=+=+6
,2,1,09003.12.15.08001.14.0500600400x..654321635241ixxxxxxxxxxxxtsiS.t.
(X
z812
11
10
9
13
min=
⎪⎪⎭
⎫⎝⎛≤⎪⎪⎭⎫⎝
⎛9008003.12
.15
.00
000011.14.0X⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫⎝
⎛50060040010
1
0010010001001
X,0
6543
21≥⎪⎪⎪⎪⎪⎪⎪
⎪⎭
⎫⎝⎛=xxxxxxX
改写为:
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对下面一组数据作二次多项式拟合。
要求:
Ⅰ.返回多项式的系数;Ⅱ.作出数据点的图形及拟合曲线;
Ⅲ.分析拟合的效果。
(2绘制如下饼状图形标签‘Math',‘English',‘Chinese',‘Music'代表的数值分别为2,4,6,8;
(3用两种方法画出函数sin(x,]
∈
x的图形。
要求:
2,0[π
Ⅰ.利用相关函数加x轴、y轴注解分别为“自变量X”、“函数Y”;Ⅱ.加图形标题“示意图”;Ⅲ.加栅格。
(八实验解答:
(九实验总结