《matlab软件应用与开发》实验指导书Word下载.docx
- 文档编号:19981602
- 上传时间:2023-01-13
- 格式:DOCX
- 页数:14
- 大小:106.11KB
《matlab软件应用与开发》实验指导书Word下载.docx
《《matlab软件应用与开发》实验指导书Word下载.docx》由会员分享,可在线阅读,更多相关《《matlab软件应用与开发》实验指导书Word下载.docx(14页珍藏版)》请在冰豆网上搜索。
>
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
inv
eig
trace
(四)实验总结
所得及建议。
实验二、MATLAB程序设计
(1)学会用MATLAB软件中M文件的编写程序;
(2)学会运用MATLAB软件解决一个实际问题。
Matlab程序设计,线性规划,动态规划。
(三)实验内容与结果:
(1)定义函数f(x1,x2)=100(x2-x12)2+(1-x1)2并计算f(1,2).
解答:
建立函数文件:
fun.m
functionf=fun(x)
f=100*(x
(2)-x
(1)^2)^2+(1-x
(1))^2;
建立主程序:
main.m
x=[12]
fun(x)
(2)任务分配问题:
某车间有甲、乙两台机床,可用于加工三种工件。
假定这两台车床的可用台时数分别为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];
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
400.0000
500.0000
fval=1.3800e+004
即在甲机床上加工600个工件2,在乙机床上加工400个工件1、500个工件3,可在满足条件的情况下使总加工费最小为13800。
(3)投资问题
引言:
问题描述
设有某种资源(或资金)M个单位(M为整数)欲分配用于N个生产项目,已知第k个生产项目获得u(k)个单位(u(k)为非负整数)这种资源后可创利润为L(u(k),k)。
L(u(k),k)是u(k)的不减函数。
问如何分配这些资源可使所获得总利润最大。
把资源分配过程分为N个阶段。
第K个阶段是向第K个生产项目分配资源。
状态变量x(k)表示分配完第1,2,….k-1个生产项目的资源数量。
x
(1)=M.,决策变量取为u(k),它是分配给第k个生产项目的资源数量。
L(u(k),k)为阶段指标。
最优值函数V(x(k+1),k+1)表示把剩余的资源x(k)分配给第k+1,k+2,…N个生产项目所能获得的最大利润。
根据动态规划方法,利用动态规划基本方程
和状态转移方程
逆向递推可求得最优化决策序列
和总利润的最大值
。
这里
应用示例
某公司新购置了某种设备6台,欲分配给下属4个企业。
已知各企业获得这种设备后年创利润如表1所示,单位为千万元。
问如何分配这些设备能使年创总利润最大,最大利润是多少?
表1各企业获得不同台数设备的年创利润数
甲
乙
丙
丁
1
2
3
6
5
7
8
9
10
程序及使用说明
functionzjsy()
M=6
N=4
L=[0,0,0,0;
4,2,3,4;
6,4,5,5;
7,6,7,6;
7,8,8,6;
7,9,8,6;
7,10,8,6]
[u,s]=ssuu(M,N,L)
function[u,I]=ssuu(M,N,L)
%M为资源数量
%N为生产项目数
%L为(M+1)XN阶利润矩阵
fori=1:
M+1
V(i,N)=L(i,N);
W(i,N)=i-1;
end
fork=N-1:
-1:
fori=1:
V(i,k)=0;
W(i,k)=i-1;
forj=1:
i
ifV(i,k)<
L(j,k)+V(i+1-j,k+1)
V(i,k)=L(j,k)+V(i+1-j,k+1);
W(i,k)=j-1;
end
V(M+1,1)=0;
W(M+1,1)=M;
forj=1:
ifV(M+1,1)<
L(j,1)+V(M+2-j,2)
V(M+1,1)=L(j,1)+V(M+2-j,2)
W(M+1,1)=j-1;
I=V(M+1,1);
u
(1)=W(M+1,1);
s=M-u
(1);
fori=2:
N
u(i)=W(s+1,i);
s=s-u(i);
end
在命令窗口命令运行输入zjsy回车后得到输出
u=
1131
I=
17
这表明,分配给甲、乙和丁企业各1台,分配给丙企业3台可使年创利润最大。
最大利润为1亿7千万元。
实验三、MATLAB绘图与数据模拟
(1)掌握MATLAB软件中的数据拟合与插值;
(2)熟练运用MATLAB软件进行绘图。
数据拟合与插值,Matlab绘图
(三)实验内容与结果
1、对下面一组数据作二次多项式拟合
xi
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
yi
-0.447
1.978
3.28
6.16
7.08
7.34
7.66
9.58
9.48
9.30
11.2
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、绘制如下饼(柄)状图形:
'
math'
'
english'
chinese'
music'
分别是2,4,6,8。
x=[2,4,6,8];
pie(x,{'
})
3在区间[0,2*pi]画sin(x)的图形,并加注图例“自变量X”、“函数Y”、“示意图”,并加格栅.
x=linspace(0,2*pi,30);
y=sin(x);
plot(x,y)
xlabel('
自变量X'
)
ylabel('
函数Y'
title('
示意图'
gridon
4、先利用随机产生的数据x,y,通过函数
生成数据z,然后利用数据x,y,z绘出函数
的图形.
编写程序如下:
首先,生成随机数据x,y,z,它们的空间分布图见图1.1.2.
rand('
seed'
0);
x=rand(50,1)*4-2;
y=rand(50,1)*4-2;
z=x.*y.*exp(-x.^2-y.^2);
figure
(1)
zmax=max(z);
zmin=min(z);
bm=zmin-0.25*(zmax-zmin);
top=zmax+0.25*(zmax-zmin);
plot3(x,y,z,'
.'
markersize'
4*6);
holdon
z0=bm*ones(1,length(z));
plot3([x'
;
x'
],[y'
y'
],[z'
z0],'
k:
holdoff
set(gca,'
Zlim'
[bm,top]);
xlabel('
);
ylabel('
zlabel('
z'
gridoff
再进行二元插值,求出绘图要使用矩形网格点上的值,然后绘图,见图1.1.3.
T=-2:
2;
[Xi,Yi]=meshgrid(T,T);
Zi=griddata(x,y,z,Xi,Yi);
figure
(2)
mesh(Xi,Yi,Zi);
hiddenoff
colormap([0.5,0.5,0.5])
holdoff
实验上机要求
1、遵守实验室一切规章制度,爱护设备;
2、认真完成每次实验任务,并按要求写好实验报告;
3、报告内容:
认真填写报告前面的内容
系:
理学院课程名称:
Matlab软件设计与应用技术日期:
2008年XX月XX日
姓名
XXXX
组号
计算机编号XXXX
学号
实验室
数学专业实验室
专业
班号
老师签名
成绩评定
实验器材
一台计算机
(注:
在实验报告纸里填写内容栏按顺序写如下四个部分)
(三)实验内容与操作结果:
2、操作结果:
(并写出源代码的注释)
3、结果分析:
(有实际模型的写)
(四)实验总结:
考核成绩
该门课程的考核成绩=操作成绩*60%+报告成绩*40%
y=[-0.4471.9783.286.167.087.347.669.589.489.3011.2];
-9.824820.1469-0.0338
x=[2,4,6,8];
pie(x,{'
x=linspace(0,2*pi,30);
自变量x'
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab软件应用与开发 matlab 软件 应用 开发 实验 指导书
![提示](https://static.bdocx.com/images/bang_tan.gif)