《北京电子科技学院数学建模实验归纳总结报告》由会员分享，可在线阅读，更多相关《北京电子科技学院数学建模实验归纳总结报告（24页珍藏版）》请在金锄头文库上搜索。
1、北京电子科技学院（BESTI）实验报告课程：数学建模与数学实验班级： 姓名： 成绩：指导教师：实验日期及时间： 必修/选修：选修2.6 习题1. 冒泡法排序 n 个数function y = qp(x) m,n=size(x); fori = 1:nfor j = 1:n-i if(x(j)x(j+1)实验一x(j) = x(j+1); x(j+1)=a; endend end y=x;运行结果：a = x(j);2. 求落地后球反弹高度functions,h = sh(n) z=100;s=z;fori = 1:n-1z=z/2; s=s+2*z;end h=z/2运行结果3. 函数func
2、tion y = f(x,y)y = x2 + sin(x*y)+2*y;运行结果：4. 绘图fplot(cos(tan(pi*x),0,1,1e-4) x=0:1e-4:1plot(x,cos(tan(pi*x)运行结果：5. 绘图ezplot(exp(x*y)-sin(x+y),-3,3,-3,3)运行结果：6. 绘图x = -3:0.1:3;y = -3:0.1:3;X,Y = meshgrid(x,y); Z = 2*X2+Y2;surf(X,Y,Z)运行结果：7. 绘图clear a=100;thita=0:0.1:2*pi; rho=a*thita; polar(thita,rho
3、)运行结果：8. 绘图clear a=2;thita=0:0.1:2*pi; rho=a*sin(3*thita); polar(thita,rho)运行结果：实验二1.某鸡场有 1000 只鸡，用动物饲料和谷物混合喂养，每天每只鸡平均食混合饲料 0.5kg，其中动物饲料所占比例不能少于20%，动物饲料每千克 0.30 元，谷物饲料每千克0.18 元，饲料公司每周仅保证供应谷物饲料 6000kg，问饲料怎样混合，才能使成本最低？解：设每天每只鸡食动物饲料为x1 个单位，谷物饲料为x2 个单位（500克/单位），则minf=1000(0.15x+0.09y)x + y = 1x = 0.2s.t
4、. =7 *1000 * y = 0实验代码：MODEL:min=150*x1+90*x2; x1+x2 = 1;x1=0.2;7000*x1 = 0;x2 = 0;end运行结果：2.某工厂用 A1，A2两台机床加工 B1，B2，B3三种不同零件。已知在一个生产周期内 A1 只能工作 8 机时；A2只能工作 100 机时。一个生产周期内计划加工 B1 为 70 件，B2为 50 件，B3 为 20 件，两台机床加工每个零件的成本，分别如下所示：加工时间：A1:1(B1),2(B2),3(B3)A2:1(B1),1(B2),3(B3)加工成本：A1:2(B1),3(B2),5(B3)A2:3(
5、B1),3(B2),6(B3)问怎样安排两台机床一个周期的加工任务，才能使加工成本最低？ 实验代码：model:min = 2*X11+3*X12+5*X13+3*X21+3*X22+6*X23; 1*X11+2*X12+3*X13=80;1*X21+1*X22+3*X23=70; 2*X12+X22=50;2*X13+3*X23=20;gin(x23); End3.model: max=12*x11+12*x21+5*x12+5*x22+4*x13+4*x23; 4*x11+3*x12+x13=180;2*x21+6*x22+3*x23 60;x1 + x2 = 70;x2 + x3 = 6
6、0; x3 + x4 = 50; x4 + x5 = 20; x5 + x6 = 30;end运行结果：5.某工厂生产 A1，A2两种型号的产品都必须经过零件装配和检验两道工序，如果每天可用于零件装配的工时只有100h，可用于检验的工时只有 120h，各型号产品每件需要占用各工序时数和可获得的利润如下所示：产品可用工时工序A1A2装配23100检验42120利润64写出此问题的数学模型，并求出最优化生产方案。实验代码：MODEL:Min= x1 + x2 + x3 + x4 + x5 + x6; x1 60;x1 + x2 = 70; x2 + x3 = 60; x3 + x4 = 50;x4
7、 + x5 = 20; x5 + x6 = 30;end6.model:max=6*x1+4*x2+5*x3; 2*x1+3*x2+4*x3=100;4*x1+2*x2+2*x3=60;x21+x22+x23=100; x11+x21=45; x12+x22=75; x13+x23=40; 4*x1+2*x2+2*x3=120; gin(x11);gin(x12); gin(x13); gin(x21); gin(x22); gin(x23); end运行结果：实验三1. 某公司在六个城市 C1，C2，C3，C4，C，5C6 中都有分公司，从Ci到 Cj 的直达航班票价由下述矩阵的第 i 行、
8、第j 列元素给出（表示无直达航班），该公司想算出一张任意两个城市之间的最廉价路线表，试作出这样的表来。代 码 ： functionD,R=floyd(a) n=size(a,1);D=afor i=1:n for j=1:n R(i,j)=j; endend Rfor k=1:n for i=1:n for j=1:nif D(i,k)+D(k,j)D(i,j)D(i,j)=D(i,k)+D(k,j); 36453R(i,j)=R(i,k);endend end k D REnd运行结果：2. 在一个城市交通系统中取出一段如图所示，其入口为顶点V1,出口为 V8，每条弧段旁的数组表示通过该路段
9、所需时间，每次转弯需要附加时间为 3，球V1 到 V8 的最短时间路径。答：V1 到V8 的最短时间路径为 15 路径为 12478。实验代码：function y=bijiaodaxiao(f1,f2,f3,f4) v12=1;v23=3;v24=2;v35=1;v47=2;v57=2;v56=6;v68=3;v78=4;turn=3;a1=123568;a2=123578;a3=12478;a4=124568;f1=v12+v23+v35+v56+turn+v68; f2=v12+v23+v35+turn+v57+turn+v78; f3=v12+turn+v24+turn+v47+v78
10、; f4=v12+turn+v24+v47+turn+v57+turn+v56+turn+v68; min=f1;a=a1;if f2min min=f2; a=a2; endif f3min min=f3; a=a3; endif f4 polytool(x,Y,1)Xi02468101214161820Yi0.62.04.47.511.817.123.331.239.649.761.72. 某零件上有一段曲线，为了在程序控制机床上加工这一零件，需要求这段曲线的解析表达式，在曲线横坐标Xi 处测得纵坐标 Yi 共 11 对数据如下：求这段曲线的纵坐标 y 关于 x 的二次多项式回归方程实验代码：x=0 2 4 6 8 10 12 14 16 18 20;y=0.6 2.0 4.4 7.5 11.8 17.1 23.3 31.2 39.6 49.7 61.7;p,S=polyfit(x,y,2)实验结果：p =0.14030.19711.0105S =R: 3x3 double df: 8normr: 1.1097y 关于横坐标 x 的二次多项式回归方程：y=0.1403x*x+0.1971x+1.01053.混凝土的抗压强度随养护时间的演唱而增加，现将一批混凝土做成12 个试块，记录了养护日期x（日）以及抗压强度y（kg/cm
