答案详解实验指导书Matlab软件应用与开发副本.docx
- 文档编号:13106805
- 上传时间:2023-06-11
- 格式:DOCX
- 页数:15
- 大小:497.68KB
答案详解实验指导书Matlab软件应用与开发副本.docx
《答案详解实验指导书Matlab软件应用与开发副本.docx》由会员分享,可在线阅读,更多相关《答案详解实验指导书Matlab软件应用与开发副本.docx(15页珍藏版)》请在冰点文库上搜索。
答案详解实验指导书Matlab软件应用与开发副本
《Matlab软件应用与开发》实验指导书
实验1Matlab操作基础及矩阵运算
(一)、实验类型:
验证型
(二)、实验类别:
基础实验
(三)、每组人数:
1
(四)、实验要求:
选修
(五)、实验学时:
3个学时
(六)、实验目的:
(1)熟悉MATLAB软件中关于矩阵运算的各种命令;
(2)学会运用MATLAB软件自定义函数,并求出函数值;
(3)学会在MATLAB环境下编写函数。
(七)、预备知识:
线性代数中的矩阵运算;高等数学中微积分知识。
本实验所用MATLAB命令
●矩阵输入格式:
A=[a11,a12;a21a22];b=初始值:
步长:
终值
●求A的转置:
A’
●求A加B:
A+B
●求A减B:
A-B
●求A乘B:
A*B
●求A的行列式:
det(A)
●求A的逆:
inv(A)
●求A的特征向量和特征值:
[KL]=eig(A)
●求A的秩:
rank(A)
●符号变量定义:
symsxy
●求函数的极限limit(.)
●求函数的导数diff(.)
●求函数的积分int(.)
●求代数方程的解solve(.)
●求微分方程的解dsolve(.)
●functiony=f()
(八)、内容与要求:
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。
可建立以下线性规划模型:
改写为:
minf*x
S.t.
源程序为:
f=[1391011128];
A=[0.41.11000;
0000.51.21.3];
b=[800;900];
aeq=[100100;010010;001001];
beq=[400;600;500];
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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 答案详解 实验指导书Matlab软件应用与开发 副本 答案 详解 实验 指导书 Matlab 软件 应用 开发