资源描述:
《matlab软件应用与开发》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、《Matlab软件应用与开发》实验指导书实验1Matlab操作基础及矩阵运算(一)、实验类型:验证型(二)、实验类别:基础实验(三)、每组人数:1(四)、实验要求:选修(五)、实验学时:3个学时(六)、实验目的:(1)熟悉MATLAB软件中关于矩阵运算的各种命令;(2)学会运用MATLAB软件自定义函数,并求出函数值;(3)学会在MATLAB环境下编写函数。(七)、预备知识:线性代数中的矩阵运算;高等数学中微积分知识。本实验所用MATLAB命令l矩阵输入格式:A=[a11a12;a21a22];b=初始值:步长:终值l求A的转置:A’l求A加B:A+Bl求
2、A减B:A-Bl求A乘B:A*Bl求A的行列式:det(A)l求A的逆:inv(A)l求A的秩:rank(A)l求函数的极限limit(.)l求函数的导数diff(.)l求函数的积分int(.)l求代数方程的解solve(.)l求微分方程的解dsolve(.)(八)、内容与要求:1、输入矩阵A,B,b;,1、矩阵转置、四则运算。C1=A’,C2=A+B,C3=A-B,C4=A*B2、求行列式。D1=
3、A
4、,D2=
5、B
6、3、求矩阵A、B的秩E1,E24、求极限5、设,求6、求7、求微分方程的通解9、自定义函数,并计算。(九)、实验操作1、>>A=[21248
7、54;2351534;21353154;21721552]A=212485423515342135315421721552>>B=[1245124;18725335;4811535;46562523]B=124512418725335481153546562523>>b=21:2:27b=212325272、>>C1=A',C2=A+B,C3=A-B,C4=A*BC1=212212124353572851531154345452C2=3369862820107686969364689671284075C3=9-2184-20-16-37-38-1-2734
8、1619-2516-1029C4=49482982266844112938452929322580485465203691405646609056536247453、>>D1=det(A),D2=det(B)D1=2181568D2=-31822764、>>E1=rank(A),E2=rank(B);E1=4E2=45、>>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表示
9、正无穷大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)-48、>>[y]=dsolve('x*Dy+y=x^2+3*x+2','x')%微分或导数的输入是用Dy、D2y、D3y、…来表示y的一阶导数或、二阶导数或、三阶导数或、…。如果自变量是系统
10、默认的,则自变量输入部分可省略。y=2+1/3*x^2+3/2*x+1/x*C19、>>建立函数文件:fun.mfunctionf=fun(x)f=100*(x(2)-x(1)^2)^2+(1-x(1))^2;建立主程序:main.mx=[12]fun(x)。。。。。实验二、Matlab程序设计(一)、实验类型:综合型(二)、实验类别:基础实验(三)、每组人数:1(四)、实验要求:选修(五)、实验学时:3个学时(六)、实验目的:学会运用Matlab软件解决线性规划中的实际问题(七)、实验内容:任务分配问题:某车间有甲、乙两台机床,可用于加工三种工件。假定这
11、两台车床的可用台时数分别为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。可建立以下线性规划模型:
12、S.t.改写为:源程序为:f=[1391011128];A=[0.