Matlab大作业
MATLAB大作业
1.设计一个程序，要求用户用键盘输入10次，输入内容可以是数字、英文、汉字、或者
其他文字或者符号。
如果输入的是数字，则把存入一个数字数组当中；如果输入的是英文的字符，则存入一个英文字符数组当中；如果是其他字符，则存入一个字符串数组中。
当10次内容都输入完毕后，先输出英文数组中的内容，然后输出字符串数组中的内容。
如果数字数组中元素数量大于3个，则由数字数组中的数字组成一个a0 X n+a1X n-1+………an=f（X），其中a0…an是数组中第一到最后一个元素，n为元素的个数，并且画出这个函数在区间[0,20]之间的曲线。
如果数字数组中的元素数量小于等于3个，则由数字数组中的数字组成一个a0 X n+a1Y n-1+anZ=0的方程，并且画出这个函数在[0,10]之间的曲线。
具体要求：
（1）画出程序的结构图
（2）画出程序的流程图
（3）每一行或者每一小段都需要标注
（4）附上完整的程序
2.按如下条件，设计一个完整的程序。
（1）有一片草原，草原上有人、狼、羊、鸡和蝗虫生活；
（2）被吃过的草，三年之后可以修复；
（3）一个人一年吃一只羊，10个人一年杀一只狼，一个人一年吃一只鸡；
（4）一头狼一年吃一只羊，一头狼三年成活一个幼崽；
（5）一只羊一年吃掉100平方米的草原；
（6）一只鸡一年吃掉1000只蝗虫；
（7）一只蝗虫一年吃掉一平方米的草原。
按如上条件，模拟出草原及其上述的生物从现在起500年间的状态。
具体要求：
（1）写出具体的解题思路
（2）列出详细的变量表，写出每一个变量的含义
（3）写出完整的程序，最好有图。
MATLAB大作业
M A T L A B大作业作业要求：（1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。
（2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。
打印稿不要求MATLAB程序，但电子稿要包含MATLAB程序。
（3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。
不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。
（4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。
（5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。
第一类：绘制图形。
（B级）问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。
斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契螺旋线，如图所示。
问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢尔宾斯基在1915年提出,它是一种典型的自相似集。
其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三角形，然后去掉中间的那一个小三角形。
接下来对其余三个小三角形重复上述操作，如图所示。
问题三：其他分形曲线或图形。
分形曲线还有很多，教材介绍了科赫曲线，其他还有皮亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。
这方面的资料很多（如），请分析构图原理并用MATLAB 实现。
问题四：模拟掷骰子游戏：掷1000次骰子，统计骰子各个点出现的次数，将结果以下表的形式显示，并绘制出直方图。
Matlab大作业
Matlab 大作业（组内成员：彭超杰、南彦东、江明伟）一、研究模型（电车）通过控制油门（保持一定角度）来调节电动机能输出稳定的转速，从而控制车速稳定。
数学依据说明如下：由图可知存在以下关系：a d a au w k R i dtdi L =++ (w k e d d =) L M M dtdw J-=a m i k M =La m M i k dtdw J -=k为反电势常数，m k为电动机电磁力矩常数，这里忽略阻尼力矩。
d二、数学模型再看整个研究对象，示意图以课本为依据，不同点是这里将数控的进给运动，转换为汽车行驶所需要的扭矩。
（这里不说明扭矩的具体产生过程，仅仅说明输出车轮旋转的角速度w ）对照课本不同，() s θ变为()s N ，1221z z w w =，1w 为电动机的转速，2w 为轮胎的转速，1z 为电动机的光轴齿轮的齿数，2z 为与轮胎相连光轴的齿轮齿数。
)(*10110w x w k x ==，121z z k =()c a m m d ba m x K K K k s k k JRs JLs K K K k s G i 1231+++=()ca m m d M K K K k s k k JRs JLs R Ls K s G L 1231)(++++-=同理，忽略电枢绕组的电感L ，简化系统传递函数方框图如下()JRK K K k JR s k k s JR K K K k s G c a m m d ba m x i 121++=()JRK K K k JR s k k s K K K K k s k k Rs R K s G c a m m d c a m m d M L 121121++-=++-=三、系统分析1.分析时间响应其传递函数如下：（1）系统时间响应令τ=0、τ=0.0125、τ=0.025，应用impulse函数，可得到系统单位脉冲响应；应用step函数，可得系统单位跃阶响应。
Matlab大作业
Matlab的基本操作一、使用函数实现对下列矩阵的左旋和右旋以及反转已知答案：如图：矩阵如图矩阵的左旋如图一矩阵的右旋如图二矩阵的左右反转如图三矩阵的上下反转如图四二、已知A=[8 9 5 ] B=[-1 3 -2][36 -7 11] [2 0 3][21 -8 5] [-3 1 9 ]计算：1)A+5*B;2)A*B和A.*B3)A^3和A.^34)A/B和B\A5)[A,B]答案：A=[8 9 5；36 -7 11；21 -8 5];B=[-1 3 -2；2 0 3；-3 1 9];A+5*BA*BA.*BA^3A.^3A/BA\B[A,B]程序运行结果如下：第二章、Matlab程序设计一、已知S=1+2+22+23+…+263求S的值答案：代码如下：ClearClcs=0,j=2for i=1:63s=s+j^iends程序运行结果如下：二、第三章、Matlab绘图一、在同一坐标系中画出下列函数的图像：x2，-x2，xsin（x）在[0,2π]上的函数图像程序代码如下：clearclcx=0:pi/100:2*pi;y1=x.^2;y2=-x.^2;y3=x.*sin(x);title('同一坐标下的函数图像')plot(x,y1,':',x,y2,'h',x,y3,'--')程序运行结果如下：二、绘制极坐标图像：程序代码如下：clearclcx=0:0.01:2*pi;y=2*cos(2*(x-pi/8));title('极坐标图像')xlabel('x')ylabel('y')polar(x,y)程序运行结果如下：第四章、Matlab符号运算一、求函数y=e-x2 的傅立叶变换及其逆变换程序代码如下：clearclcsyms x t;y=exp(-x^2);Ft=fourier(y,x,t)fx=ifourier(Ft,t,x)程序运行结果如下：二、求下列极限值：答案：1、程序代码：clearclcsyms x;s1=sin(2*x)/sin(5*x);limit(s1,x,0)s2=(1+1/x)^(2*x);limit(s2,x,inf)程序运行结果如下：第五章、Matlab数值运算一、建立一个3*4阶随机矩阵，求矩阵的最大值、最小值、方差和标准差、极差、协方差，和自相关阵程序代码如下：clearclcA=rand(3,4)B=var(A)C=std(A)D=range(A)E=cov(A)F=corrcoef(A)程序运行结果如下：二、求函数f（x）=x3-2x+1在x=[-1 1]之间的极小值和x=-1附近的零点程序代码如下：clearclc[x,y]=fminbnd('x.^3-2.*x+1',-1,1)[x,y]=fzero('x.^3-2.*x+1',-1)程序运行结果如下：第六章、Matlab图形用户界面一、建立一个具有三个输入框的窗口对话框程序代码如下：clearclcprompt={'姓名','年龄','班级'}；DTitle='注册学生信息';line[1;1;1];def{'乔阳'，'18','100412105'};info=inputdlg(prompt,DTitle,line,def,'ON')程序运行结果如下：二、设计一个表现下载进度的进度条程序代码如下：clearclch = waitbar(0,'正在下载，请稍等……');for i=1:1000waitbar(i/10000,h)endclose(h)程序运行结果如下：第七章、Matlab Simulink 仿真一、仿真信号x（t）=sin(t)sin(10t)的波形仿真步骤：1.建立模型窗口：生成无标题（intitled）的模型窗口；2.添加信号源模块（sine）、输出模块（scope）、数学模块（Dot product）3.设置模块参数：如下图：4.编辑模块即将各个模块连接起来。
matlab大作业
2014年春季学期MATLAB 课程考查题一.必答题(80分)1.如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有何区别？(2)答：方法一：在MATLAB命令窗口中输入editpath或pathtool命令或通过【File】/【SetPath】菜单，进入“设置搜索路径”对话框，通过该对话框编辑搜索路径。
方法二：在命令窗口执行“path(path,…D:\Study ‟)”，然后通过“设置搜索路径”对话查看“D:\Study”是否在搜索路径中。
方法三：在命令窗口执行“addpathD:\Study- end”，将新的目录加到整个搜索路径的末尾。
如果将end改为begin，可以将新的目录加到整个搜索路径的开始。
区别：当前文件目录是正在运行的文件的目录，显示文件及文件夹的详细信息，且只有将文件设置为当前目录才能直接调用。
搜索路径中的文件可以来自多个不同目录，在调用时不用将其都设置为当前目录，为同时调用多个文件提供方便。
2.创建符号变量和符号表达式有哪几种方法？(4)答：定义符号变量：方法一：sym函数，可以定义单个符号变量，调用格式为：符号量名=sym（'符号字符串'）；方法二，syms函数，定义多个符号变量，调用格式：syms 符号变量名1 符号变量名2 符号变量名3 符号变量名4 。
空格隔开。
定义符号表达式：方法一：用单引号定义符号表达式；方法二：用sym函数定义符号表达式；方法三：用已经定义好的符号变量组成符号表达式；3.GUIDE提供哪些常用的控件工具，各有什么功能？（5分）答：一、控件风格和外观（1）BackgroundColor：设置控件背景颜色，使用[R G B]或颜色定义。
（2）CData：在控件上显示的真彩色图像，使用矩阵表示。
（3）ForegroundColor：文本颜色。
（4）String属性：控件上的文本，以及列表框和弹出菜单的选项。
（5）Visible：控件是否可见。
matlab大作业实验报告
matlab大作业学号姓名：年级：专业：1、产生一个10 10的随机矩阵A，要求A中元素均为整数，范围[1,50]。
1)求出A中所有元素之和S，平均值M。
2)找到所有小于平均值，且能被3整除的元素。
3)绘制出A的二维纵向柱状图，横坐标为[8 5 9 1 2 3 4 7 10 13],条形宽度为0.7的“stacked”样式。
代码如下：clc,clear all,close allA=round(rand(10,10)*50);disp(A)S=sum(sum(A));P=mean(mean(A));disp(S)disp(P)disp('所有小于平均数且能被三整除的元素')XPS=H((mod(H,3)==0)&(H<P));disp(XPS')subplot(1,1,1),bar(A,0.7,'stacked'),title('ygh');set(gca,'XTickLabel',{'8','5','9','1','2','3','4','7','10','13'})2、产生一个随机四位密码。
用户用“input”进行输入对比。
猜错提示“WRONG”，正确提示“RIGHT”同时退出程序，最多五次机会。
代码如下：clc,clear all,close alldisp('请输入密码')A=round(8999*rand(1,1))+1000;m=1;while m<=5N=input('请输入一个四位数：');if A==N;disp('RIGHT');breakelsedisp('WRONG');endm=m+1;enddisp('密码是：')disp(A)disp('输入结束')3、按照脚本文件的编程风格，用for和while循环嵌套输出如下的乘法口诀表。
matlab课程设计大作业
matlab课程设计大作业一、教学目标本课程的教学目标是使学生掌握MATLAB基本语法、编程技巧以及MATLAB 在工程计算和数据分析中的应用。
通过本课程的学习，学生将能够熟练使用MATLAB进行简单数学计算、线性方程组求解、函数图像绘制等。
1.掌握MATLAB基本语法和编程结构。
2.了解MATLAB在工程计算和数据分析中的应用。
3.熟悉MATLAB的函数库和工具箱。
4.能够使用MATLAB进行简单数学计算。
5.能够使用MATLAB求解线性方程组。
6.能够使用MATLAB绘制函数图像。
7.能够利用MATLAB进行数据分析和处理。
情感态度价值观目标：1.培养学生对计算机辅助设计的兴趣和认识。
2.培养学生团队合作和自主学习的能力。
二、教学内容本课程的教学内容主要包括MATLAB基本语法、编程技巧以及MATLAB在工程计算和数据分析中的应用。
1.MATLAB基本语法：介绍MATLAB的工作环境、基本数据类型、运算符、编程结构等。
2.MATLAB编程技巧：讲解MATLAB的函数调用、脚本编写、函数文件编写等编程技巧。
3.MATLAB在工程计算中的应用：介绍MATLAB在数值计算、线性方程组求解、图像处理等方面的应用。
4.MATLAB在数据分析中的应用：讲解MATLAB在数据采集、数据分析、数据可视化等方面的应用。
三、教学方法本课程采用讲授法、案例分析法、实验法等多种教学方法相结合的方式进行教学。
1.讲授法：通过讲解MATLAB的基本语法、编程技巧以及应用案例，使学生掌握MATLAB的基本知识和技能。
2.案例分析法：通过分析实际工程案例，使学生了解MATLAB在工程计算和数据分析中的应用。
3.实验法：安排上机实验，使学生在实际操作中巩固所学知识，提高实际编程能力。
四、教学资源本课程的教学资源包括教材、实验设备、多媒体资料等。
1.教材：选用《MATLAB教程》作为主要教材，辅助以相关参考书籍。
2.实验设备：为学生提供计算机实验室，配备有MATLAB软件的计算机。
MATLAB期末大作业
1.龟兔赛跑本题旨在可视化龟兔赛跑的过程。
比赛的跑道由周长为P面积为A的矩形构成。
每单位时间，乌龟沿跑道缓慢前进一步，而兔子信心满满，每次以一个固定的概率决定走或不走。
如果选择走，就从2-10步中等概率选择一个步长。
每个单位时间用一个循环表示。
赛跑从矩形跑道左上点（0，0）开始，并沿顺时针方向进行。
不管是乌龟或兔子，谁先到达终点，比赛就告结束。
要求：编写MATLAB程序可视化上述过程。
程序以P，A以及兔子每次休息或前进的概率为输入参量。
程序必须可视化每个时刻龟兔赛跑的进程，并以红色“*”表示乌龟，蓝色的“—”表示兔子。
测试时可取P=460，A=9000。
通过上述例子，可否从理论和实验角度估计兔子休息或前进的概率，是的兔子和乌龟在概率意义下打平手。
2.黄金分割Fibonacci数列F n通过如下递推格式定义F n=F n-1+F n-2，其中F0=F1=1要求：1.计算前51项Fibonacci数，并存入一个向量2.利用上述向量计算比值F n/F n-13.验证该比值收敛到黄金比例(1+5)/23.图像处理此题旨在熟悉图像处理的基本操作，请各位自己选择一张彩色图像pMATLAB以三维数组读取一张彩图。
该彩图上每个像素位置分别存放一个取值0-255的三维向量，其三个分量分别表示该点的红（R）绿（G）蓝（B）强度信息。
要求：编写MATLAB程序，读入原始彩色图像，并且在一个图形窗口界面下显示六张图像。
这六张图分别是原始RGB彩图，及其5各变形：RBG，BRG，BGR，GBR和GRB。
每张子图要求以其对应变形命名。
最后将图像以a.jpg形式保存并黏贴至报告中。
提示：imread, imshow, cat。
MATLAB大作业
贵州大学实验报告姓名学号实验组实验时间指导教师成绩实验项目名称MATLAB大作业实验通过综合考察，让同学们能对学到的知识综合运用，更好地灵活地解决问题目的实验利用matlab提供的一些函数实现数据分析和数值计算。
原理实验运行MATLAB7.0或以上的计算机一台仪器实验内容实验数据1、试编写名为test01.m的MATLAB函数，用以计算下述的值：⎪⎩⎪⎨⎧-<->=ttntttntf的)4/sin()(si对所有)4/sin(其他情况)sin(的)4/sin()(si对所有)4/sin()(ππππ要绘制t关于函数f(t)的图形，其中t的取值范围为ππ66≤≤-t，间距为10/π。
（提示：注意要产生一系列的点，这里可考虑t的输入是向量形式，可以利用find函数找出所需限定值的元素的位置，对其按需要赋值后，再进行绘图；其次，另外一种思路，也可考虑使用循环的形式来实现）①test01.m代码如下：function y=test01(t)if sin(t)>sin(pi/4)y=sin(pi/4);elseif sin(t)<sin(-pi/4)y=sin(-pi/4);elsey=sin(t);end②调用test01函数计算值③绘制t关于函数f(t)的图形，其中t的取值范围为ππ66≤≤-t，间距为10/π代码如下：t=-6*pi:pi/10:6*pi;y=sin(t);y(find(y>sin(pi/4)))=sin(pi/4);%f(t)=sin(pi/4),y>sin(pi/4) y(find(y<sin(-pi/4)))=sin(-pi/4);%f(t)=sin(pi/4),y>sin(pi/ 4)plot(t,y) %绘制图形title('t from -6{\pi} to 6{\pi}') %添加图形标题xlabel('t') %x轴标签ylabel('f(t)') %y轴标签实验结果：2、编写函数，在同一窗口的4个子图中利用plot等语句绘制y=at2图像，其中a=[1 2 5 10],t范围[-2，5]。
matlab综合大作业(附详细答案)
m a t l a b综合大作业(附详细答案)-标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII《MATLAB语言及应用》期末大作业报告1．数组的创建和访问（20分，每小题2分）：1)利用randn函数生成均值为1，方差为4的5*5矩阵A；实验程序：A=1+sqrt(4)*randn(5)实验结果：A =0.1349 3.3818 0.6266 1.2279 1.5888-2.3312 3.3783 2.4516 3.1335 -1.67241.2507 0.9247 -0.1766 1.11862.42861.5754 1.6546 5.3664 0.8087 4.2471-1.2929 1.3493 0.7272 -0.6647 -0.38362)将矩阵A按列拉长得到矩阵B；实验程序：B=A(:)实验结果：B =0.1349-2.33121.25071.5754-1.29293.38183.37830.92471.65461.34930.62662.4516-0.17665.36640.72721.22793.13351.11860.8087-0.66471.5888-1.67242.42864.2471-0.38363)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的矩阵C；实验程序：C=[A(2,2),A(2,4);A(3,2),A(3,4)]实验结果：C =3.3783 3.13350.9247 1.11864)寻找矩阵A中大于0的元素；]实验程序：G=A(find(A>0))实验结果：G =0.13491.25071.57543.38183.37830.92471.65461.34930.62662.45165.36640.72721.22793.13351.11860.80871.58882.42864.24715)求矩阵A的转置矩阵D；实验程序：D=A'实验结果：D =0.1349 -2.3312 1.2507 1.5754 -1.29293.3818 3.3783 0.9247 1.6546 1.34930.6266 2.4516 -0.1766 5.3664 0.72721.2279 3.1335 1.1186 0.8087 -0.66471.5888 -1.67242.4286 4.2471 -0.38366)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E；实验程序：E=flipud(fliplr(A))实验结果：E =-0.3836 -0.6647 0.7272 1.3493 -1.29294.2471 0.80875.3664 1.6546 1.57542.4286 1.1186 -0.1766 0.9247 1.2507-1.6724 3.1335 2.4516 3.3783 -2.33121.5888 1.2279 0.6266 3.3818 0.13497)删除矩阵A的第2列和第4列得到矩阵F；实验程序：F=A;F(:,[2,4])=[]实验结果：F =0.1349 0.6266 1.5888-2.3312 2.4516 -1.67241.2507 -0.17662.42861.5754 5.3664 4.2471-1.2929 0.7272 -0.38368)求矩阵A的特征值和特征向量；实验程序：[Av,Ad]=eig(A)实验结果：特征向量Av =-0.4777 0.1090 + 0.3829i 0.1090 - 0.3829i -0.7900 -0.2579 -0.5651 -0.5944 -0.5944 -0.3439 -0.1272-0.2862 0.2779 + 0.0196i 0.2779 - 0.0196i -0.0612 -0.5682 -0.6087 0.5042 - 0.2283i 0.5042 + 0.2283i 0.0343 0.6786 0.0080 -0.1028 + 0.3059i -0.1028 - 0.3059i 0.5026 0.3660 特征值Ad =6.0481 0 0 0 00 -0.2877 + 3.4850i 0 0 00 0 -0.2877 - 3.4850i 0 00 0 0 0.5915 00 0 0 0 -2.30249)求矩阵A的每一列的和值；实验程序：lieSUM=sum(A)实验结果：lieSUM =-0.6632 10.6888 8.9951 5.6240 6.208710)求矩阵A的每一列的平均值；实验程序：average=mean(A)实验结果：average =-0.1326 2.1378 1.7990 1.1248 1.24172．符号计算（10分，每小题5分）：1)求方程组20,0++=++=关于,y z的解；uy vz w y z w实验程序：S = solve('u*y^2 + v*z+w=0', 'y+z+w=0','y,z');y= S. y, z=S. z实验结果：y =[ -1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))-w] [ -1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))-w] z =[ 1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))] [ 1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))]2)利用dsolve 求解偏微分方程,dx dyy x dt dt==-的解； 实验程序：[x,y]=dsolve('Dx=y','Dy=-x')实验结果：x =-C1*cos(t)+C2*sin(t)y = C1*sin(t)+C2*cos(t)3．数据和函数的可视化（20分，每小题5分）：1)二维图形绘制：绘制方程2222125x y a a +=-表示的一组椭圆，其中0.5:0.5:4.5a =；实验程序：t=0:0.01*pi:2*pi; for a=0.5:0.5:4.5; x=a*cos(t); y=sqrt(25-a^2)*sin(t); plot(x,y) hold on end实验结果：2) 利用plotyy 指令在同一张图上绘制sin y x =和10x y =在[0,4]x ∈上的曲线；实验程序：x=0:0.1:4; y1=sin(x); y2=10.^x;[ax,h1,h2]=plotyy(x,y1,x,y2); set(h1,'LineStyle','.','color','r'); set(h2,'LineStyle','-','color','g'); legend([h1,h2],{'y=sinx';'y=10^x'});实验结果：3)用曲面图表示函数22z x y =+；实验程序：x=-3:0.1:3; y=-3:0.1:3; [X,Y]=meshgrid(x,y); Z=X.^2+Y.^2; surf(X,Y,Z)实验结果：4)用stem 函数绘制对函数cos 4y t π=的采样序列；实验程序：t=-8:0.1:8;y=cos(pi.*t/4); stem(y)实验结果：4. 设采样频率为Fs = 1000 Hz ，已知原始信号为)150π2sin(2)80π2sin(t t x ⨯+⨯=，由于某一原因，原始信号被白噪声污染，实际获得的信号为))((ˆt size randn x x+=，要求设计出一个FIR 滤波器恢复出原始信号。
MATLAB大作业
M A T L A B大作业作业要求：（1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。
（2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。
打印稿不要求MATLAB程序，但电子稿要包含MATLAB程序。
（3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。
90问题五：利用MATLAB软件绘制一朵鲜花，实现一定的仿真效果。
提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。
第二类：插值与拟合。
（B级）问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。
注意，整个实验过程中从未减速。
在一组时间段50个时间点的速度。
（2）绘制插值图形并标注样本点。
问题二：估算矩形平板各个位置的温度。
已知平板长为5m，宽为3m，平板上3×5栅格点上的温度值为44，25，20，24，30；42，21，20，23，38；25，23，19，27，40。
（1）分别使用最近点插值、线性插值和三次样条插值进行计算。
（2）用杆图标注样本点。
（3）绘制平板温度分布图。
对a,b,c,d的值。
提示：曲线拟合并绘图分析第三类：定积分问题。
（B级）问题一：地球密度随着离中心（r=0）距离的变化而变化，不同半径处的密度如表所示，试估问题二：河道平均流量Q（m3/s）可使用速度和深度的乘积的积分来计算（河道横截面不规则），公式如下。
其中V(x)是离岸x（m）距离处的水速（m/s）,H(x)是离岸x距离处的水深（m）。
根据收集到过5（1（2（3（Q，单位是m（1（2（1（2（3）将节点1的力改为方向向上，计算这种改变对H2和V2的影响。
（4）将节点1的力撤销，而在节点1和2处施加1500N的水平外力，求节点3处垂直反作用力（V3）。
MATLAB大作业题目
一、 程序分析题1、命令窗口中输入如下命令，在%号后写出每行命令的功能。
（不用写出运算结果） A=1; f0=4; phi=pi/3; w0=2*pi*f0;t=0:0.01:1; %__________________________________________ y=A*cos(w0*t+phi) ;axis off %__________________________________________ PT=plot(t,y); %__________________________________________ set(PT,’LineWidth ’,[5]); %__________________________________________ AX=gca; %__________________________________________ set(AX,’FontSize ’,14); %__________________________________________ axis square %__________________________________________ xlabel(‘时间[s]’); %__________________________________________ T=title(‘余弦波’) %__________________________________________ set(T,’FontSize ’,’16’); %__________________________________________ legend(‘cos ’); %__________________________________________ grid on %__________________________________________二、 程序设计题2、已知两个矩阵：3765213100132610A ---⎡⎤⎢⎥---⎢⎥=⎢⎥-⎢⎥-⎣⎦，1468237021570010B ⎡⎤⎢⎥⎢⎥=⎢⎥--⎢⎥-⎣⎦将矩阵A 中所有等于-1的元素改为3，将矩阵B 中等于0的元素值改为新的矩阵A 中相应位置元素的值，最后将新的矩阵A 的第三、四列的元素全改为-1，编写相应的程序。
MATLAB大作业
1. 某控制系统为单位负反馈系统，其开环传递函数为：()()10.5KG s s s =+。
（1） 要求设计一校正装置，使校正后系统的静态速度误差系数20v K >，超调量%15%σ<，调整时间5s t s <，同时使用Simulink 仿真（示波器）验证。
（2） 绘制校正前后系统的根轨迹。
解：1）查看符合条件的zeta 。
zeta=0:0.001:0.99;sigma=exp(-zeta*pi./sqrt(1-zeta.^2))*100; plot(zeta,sigma) xlabel('\zeta'); ylabel('\sigma');title('\sigma%=e^{-\zeta*\pi/sqrt(1-\zeta^2)}*100%', 'fontsize',16) gridz=spline(sigma,zeta,15)运行结果： z=0.5169σ与ζ的关系曲线如图1-1图1-12）求系统期望主导极点sigma=0.15;zeta=((log(1/sigma))^2/((pi)^2+(log(1/sigma))^2))^(1/2);wn=1/zeta;p=[1 2*zeta*wn wn*wn];s=roots(p)运行结果：s = -1.0000 + 1.6560i，-1.0000 - 1.6560i3）查看根轨迹G0=tf([1],[0.5 1 0]);rlocus(G0);sgrid(.5169,[])根轨迹如图1-2和局部放大根轨迹图1-3图1-2图1-3读出增益为1.88，则校正前系统的稳态误差Kv=1.88/1=1.88。
按要求，偶极子的零点和极值比值应为20/1.88=10.6383,取Zc=0.01，Pc=0.01/11=0.0009,因此校正环节为0009.0)01.0(88.1)(++=s s s Gc4)得出校正后的系统，并进行验证。
matlab期末大作业题目及答案
matlab 期末大作业（30分，每题6分）1. 积分运算（第四数值和五章符号）(1)定积分运算：分别采用数值法（quad ，dblquad ）和符号运算（syms, int ）一重定积分π⎰1. 数值法（quad ）a) 运行代码：b) 运行结果：2. 符号运算（syms ）a) 运行代码：b) 运行结果：二重定积分112200()x y dxdy+⎰⎰1.数值法（dblquad）：a)运行代码：b)运行结果：2.符号运算（syms）：a)运行代码：b)运行结果：(2) 不定积分运算sin dxdy ⎰⎰（(x/a)+b/y） i.运行代码：ii.运行结果：2. 用符号法和数值法求解线性代数方程 （第五章和第二章）⎩⎨⎧=+=+12*22x *213*12x *a11y a a y a （1） 用syms 定义待解符号变量x,y 和符号参数a11,a12,a21,a22，用符号solve 求x,y 通解 1. 运行代码：2. 运行结果：（2） 用subs 带入a11=2,a12=4,a21=6,a22=8,求x 和y 特解，用vpa 输出有效数值4位的结果 1. 运行代码：2. 运行结果：（3） 采用左除(\)和逆乘法求解符号参数赋值后的方程 ⎩⎨⎧=+=+12*8x *63*4x *2y y1. 运行代码：2. 运行结果：3.数值法和符号法求解非线性方程组(第四数值和五章符号 )（1）采用数值法（fsolve ）求解初始估计值为x0 = [-5; -5]的数值解1. 运行代码：2. 运行结果：21x 21x 21e x 2x e x x 2--=+-=-（2）符号法（solve ）的符号结果用eval 或double 转化为数值结果.1. 运行代码：2. 运行结果：4. 解二阶微分方程 (第四数值和五章符号 )⎪⎩⎪⎨⎧===++6)0(',0)0(09322y y y dx dy dx y d（1）数值ode 求特解，用plot (x,y) 画t 在[0,10]范围内（x ，y ）数值曲线 1. 运行代码：2. 运行结果：（2）符号运算dsolve求通解，用ezplot画t在[0,10]范围内（x，y）符号曲线1. 运行代码：2. 运行结果：5. 三维绘图（第六章）已知：x和y都在[-8,8]范围内，采用subplot(3,1,x)绘制三个子图，它们分别是用meshgrid和mesh绘制网格图、用c=contour 绘制等位线和用surf 绘制曲面图1.运行代码：2.运行结果：。
MATLAB大作业
实验七图像的形态学运算一、实验目的1.掌握图像的奇异值分解与合成2.能够利用模板操作实现图像的腐蚀运算3.能够利用模板操作实现图像的膨胀运算4.掌握图像的开、闭运算二、实验内容1. 读下面的程序，分析该程序要完成什么样的功能。
A=[1 2 1;2 3 2;1 2 1];[u d v]=svd(A);[nx,ny]=size(A);I=eye(nx,nx);for i=nx:-1:2I(i,i)=0;endB=u*I*d*v';BX=imread('cameraman.bmp');X=double(X);A=X(20:256,20:256);[nx,ny]=size(A);I=eye(nx,nx);for i=nx:-1:21I(i,i)=0;end[u d v]=svd(A);B=u*I*d*v';A=uint8(A);`B=uint8(B);subplot(2,2,1);imshow(A)subplot(2,1,2);imshow(B)该程序的功能是对图像分别进行膨胀和腐蚀的操作！2. 利用模板操作，编写程序实现二值图像的腐蚀运算函数imerode的功能。
程序：a=imread('D:\005.bmp');b=~a;m=[1;1;1];n=[1 1 0];p=[0 0 0;0 1 0;0 0 0];q=[0 0 1;0 0 0;1 0 0];c1=imerode(b,m);c2=imerode(b,n);c3=imerode(b,p);c4=imerode(b,q); subplot(1,5,1);imshow(b)subplot(1,5,2);imshow(c1)subplot(1,5,3);imshow(c2)subplot(1,5,4);imshow(c3)subplot(1,5,5);imshow(c4)3. 利用模板操作，编写程序实现二值图像的膨胀运算函数imdilate的功能。
MATLAB大作业
若系统是完全可控的，则可引入状态反馈调节器，且：
u v Kx
这时，闭环系统的状态空间模型为： (A BK)x Bu x y Cx 系统方框图如下：
此时，相稳定裕量Pm=46.7deg,满足题目中430<γ<480的要求。
（5）计算校正后系统阶跃响应曲线及其性能指标。 程序如下： % BodeJiaoZheng4 % 系统校正后性能指标及阶跃响应曲线 global y t; k0=1000;n1=1;d1=conv([1 0],[1 2]); s1=tf(k0*n1,d1); n2=[0.05443 1];d2=[0.009339 1]; s2=tf(n2,d2); sope=s1*s2; sys=feedback(sope,1); step(sys); %绘制阶跃响应曲线 [y,t]=step(sys); % 求出阶跃响应的函数值及其对应时间 [sigma,tp,ts]=ste(y,t) % 调用函数ste()
1 s(s 2)
2）作原系统的波特图与阶跃响应曲线，检查是否满足要求。
在MATLAB命令窗口输入：（SysExpBode.M) >>k0=1000;n1=1;d1=conv([1 0],[1 2]); [mag,phase w]=bode(k0*n1,d1); s1=tf(k0*n1,d1); figure(1);
根据校正后的结构与参数，调用函数ste()，给出上述例题系统的性 能指标及阶跃响应曲线。程序如下：
%计算系统超调量，峰值时间，调节时间。 global y,t s1=tf(20700,[1 20 75 0]); s2=tf([1 3.841],[1 73.12]); sope=s1*s2;sys=feedback(sope,1); step(sys) [y,t]=step(sys); [sigma,tp,ts]=ste(y,t) 运行该程序，可得到系统的阶跃响应曲线，并有系统性能指标： 超调量：sigma=0.2489<30% 峰值时间：tp=0.2209 调节时间：ts=0.3037<0.5s
matlab大作业例子
Matlab作业题目1:(1)程序部分：>> x=0:pi/50:2*pi;>> y=cos(0.5+((3*sin(x))./(1+x.^2))); >> plot(x,y)>> xlabel('x');>> ylabel('y');>> title('x-y');(2)运行结果截图：（1）程序部分：>> clear>> clc>> x=0:pi/100:4*pi;>> y1=sin(x);>> y2=cos(x);>> plot(x,y1,'r-',x,y2,'g:') >> hold on>> y3=x(find(abs(y1-y2)<0.001)); >> plot(y3,sin(y3),'*')（2）运行结果截图：（1）程序部分：>> t=(0:2*pi/100:2*pi)'; >> y1=sin(t)*[1,-1]; >> y2=sin(t).*sin(9*t); >> subplot(2,1,1);>> plot(t,[y1,y2]);>> subplot(2,1,2);>> plot(t,[y1,y2]) （2）运行结果截图：题目4（1）程序部分：>> t=0:pi/50:4*pi;>> y=exp(-t/3);>> y0=exp(-t/3).*sin(3*t); >> plot(t,y,'r-',t,y0,'b:') >> xlabel('\bf\it t')>> ylabel('\bf\it y')>> grid on;（2）运行结果截图：题目5（1）程序部分：>> n=0;>> sum=0;>> while sum<2000n=n+1;sum=sum+n;end>> n（2）运行结果截图：题目6（1）程序部分：for j=1:3n=input('n='); sum=0;for i=1:na=1/(i.^2);sum=sum+a;endPI=sqrt(6*sum) end（2）运行结果截图：题目7（1）程序部分：n0=0;y0=0;while 3*y0<5n0=n0+1;y0=y0+1/(2*n0-1);y=y0-1/(2*n0-1);n=n0-1;endn（2）运行结果截图：题目8（1）程序部分：for i=1:7x=input('put x:');if x<0&x~-3y=x^2+x-6;elseif x>=0&x<5&x~=3&x~=2 y=x*x-5*x+6;elsey=x*x-x-1;endend（2）运行结果截图：题目9①switch语句（1）程序部分：>> score=input('请输入成绩'); switch fix(score/10)case {9,10}disp('A');case {8}disp('B');case {7}disp('C');case {6}disp('D');case {0,5}disp('E');otherwisedisp('error');end（2）运行结果截图：如右图②if语句（1）程序部分：score=input('请输入成绩:') n=score/10;if n>=9&n<10disp Aelseif n>=8&n<9disp Belseif n>=7&n<8disp Celseif n>=6&n<7disp Delseif n>=0&n<6disp Eelsedisp errorend（2）运行结果截图：题目10（1）程序部分：t=input('员工的工作时间:')if t<60m=84*t-700;elseif t>120m=84*(t-120)*1.15+84*120;elsem=84*t;endm2）运行结果截图：题目11（1）程序部分：a=round(20*rand(5,6))n=input('请输入n的值:') tryb=a(n,:);catchb=a(5,:);endblasterr（2）运行结果截图：。
matlab大作业1
计算机实习班级：姓名：学号：题目一：使用while循环及for循环计算1+2+3+…+100分析：无。
解答：MATLAB程序：%使用for-end语句>> clear>> sum=0;>> for i=1:1:100sum=sum+i;end>> sumsum =5050%使用while-end语句>> clear>> sum=0;i=0;>> while i<100i=i+1;sum=sum+i;end>> sumsum =5050小结：了解matlab中的基本语句，与C语言中的相关内容进行对比，从而更好的掌握及运用语句。
题目二：圆轮在平面上滚动时，圆周上任一点的轨迹为摆线，圆内任一点的轨迹为内摆线，圆外任一点的轨迹为外摆线，建立模型，并作出图。
分析：建立模型：设r为轮的半径，R为点的半径，则由分析得方程x=rt-Rsinty=r-Rcost；（t为参数）MATLAB程序：t=0:0.1:10;r=1;R=input('R=');x=r*t-R*sin(t); y=r-R*cos(t); plot(x,y),axis('equal');hold ont=0:0.1:10;r=1;R=input('R=');x=r*t-R*sin(t); y=r-R*cos(t); plot(x,y),axis('equal');hold ont=0:0.1:10;r=1;R=input('R=');x=r*t-R*sin(t); y=r-R*cos(t); plot(x,y),axis('equal');hold on令R=1, 0.5, 2，作出如下图形小结：逐渐将matlab应用到解决实际问题当中，了解matlab强大的绘图功能。
MATLAB大作业
M A T L A B大作业(总15页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--MATLAB大作业作业要求：（1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。
（2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。
打印稿不要求MATLAB程序，但电子稿要包含MATLAB程序。
（3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。
不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。
（4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。
（5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。
第一类：绘制图形。
（B级）问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。
斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契螺旋线，如图所示。
问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢尔宾斯基在1915年提出,它是一种典型的自相似集。
其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三角形，然后去掉中间的那一个小三角形。
接下来对其余三个小三角形重复上述操作，如图所示。
问题三：其他分形曲线或图形。
分形曲线还有很多，教材介绍了科赫曲线，其他还有皮亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。
MATLAB 大作业
MATLAB 大作业请各位同学，自己完成matlab 大作业的内容。
禁止相互抄袭，如有雷同，零分计算。
大作业的格式按照实验报告的格式书写，务必标明题号，作业完成后，将生成的报告打印出来提交。
正文的字号以宋体五号字，1.5倍行距的格式打印。
请与18周五前将大作业报告由班级负责人统一收齐交给我，过期不收。
1、 试编写名为test01.m 的MATLAB 函数，用以计算下述的值：⎪⎩⎪⎨⎧-<->=t t n t t t n t f 的对所有其他情况的对所有)4/sin()(si )4/sin()sin()4/sin()(si )4/sin()(ππππ要绘制t 关于函数f (t )的图形，其中t 的取值范围为ππ66≤≤-t ，间距为10/π。
（提示：注意要产生一系列的点，这里可考虑t 的输入是向量形式，可以利用find 函数找出所需限定值的元素的位置，对其按需要赋值后，再进行绘图；其次，另外一种思路，也可考虑使用循环的形式来实现）2、 编写函数，在同一窗口的4个子图中利用plot 等语句绘制y=at 2图像，其中a=[1 2 5 10],t 范围[-2，5]。
3、 求函数32)(3-+=x x x f 在区间[-5,5]上的最大值和最小值。
4、 求解函数⎰102dx e x 的数值积分和符号积分，并比较结果。
5、 求解微分方程3|;1|2)1(002='='=''+==x x y y y x y x 的精确解和解析解，并绘制图形。
假设求解区间为[0,10] 。
6、 说说你对MATLAB 及应用这门课程学习后的体会，另外请说明在所学章节中哪一章的内容你最感兴趣，为什么？哪一章的内容你认为是没有必要学习的，为什么？如果可以选择MATLAB 的学习的内容的话，谈谈你所期望学到的知识类别的前三种。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
q ((q-1)p+1)s
((q-1)p+2)s
((q-1)p+3)s
…
(qp)s
…
…
…
…
…
…
3
(2p+1)s
(2p+2)s
(2p+3)s
…
(3p)s
2
(p+1)s
(p+2)s
(p+3)s
…
(2p)s
1
1s
2s
3s
…
ps
层/ 列
1
2
3
…
p
这种对应关系是指将货物代码为1g的货物放在货位编码为1s的货位上，货物
3、货位优化模型的建立
根据“提高作业效率”和“满足限制条件”的原则，进行货位的优化调整， 使各货品处于理想的货位。在货位优化中．满足限制条件是以损失作业效率为代 价换取的。因此，货位优化是一个多目标决策问题。
3、1模型中的变量定义
仓储区货位规划有许多目标，但主要目标有两个: 一定时期内全部货物出入 库搬运距离之和最小;降低货架重心，定性而言，要实现这两个目标，可以将周 转率高的货物安排在接近仓储区出入口的位置，将重货放在低层货架，定量而言， 可以根据货位规划问题的两个主要目标建立货位规划模型，设仓储区的每个货架 有货位p列q层，每个货位的长和高分别为a、b，最靠近出入口的为第1列，最接 近地面的为第1层，仓储区采用叉车作业。 (1) 主目标函数F1
左边的一列作为第1列，将所有的货物按顺序编一个代码1g、2g、3g、4g 、5g、
6g、7g..,,，货物代码不变，然后随机生成一个对应的货位编码，比如1s、3s、
4s、6s、7s、8s、9s、…,将货物的代码和货位的编码对应起来，表2 给出了货
物与货位的对应关系。
表1 p 列q 层货架货位的编码顺序
定存放的位置，而同属一类的不同货品又按一定的法则来分配 货位。也有与定位储放同样的缺点。
(4)分类随机储放。每一类货品有固定存放位置，但在各类的储区内，每个 货位的分配是璃机的。分类随机储放使货品出入库管理及盘点工作的难度较高。
2、当前弊端及解决策略
分析现存的几种存储策略，发现有如下不足之处： 其一，货位相对比较固定，刚性比较大。其二，货位的分配原则中只是针对 生活资料而言，对于体积大、吨位高的生产资料(例如钢材)等大件货物并不适用。 其三，货位的存储策略中，没有一套对应的货位调度方案。针对以上存在的问题， 结合目前物流企业的实际情况，设计出了一套优化的存储策略和货位调度方案， 并从以下几点进行分析。 (一)货位规划的原则。在现有的物流环境下，货位规划要考虑的原则是很多 的。一般应考虑的分配原则大致如下： (1)货架承载均匀。物料分散存放在货场的不同位置，避免因集中存放而造 成货架受力不均匀，致使货架受损。(2)较重的物品存放在下面的货位，较轻的 物品存放在高处的货位，使货架承载稳定。(3)加快周转，先人先出，同种物料 出库时，先入库者先提取出库，避免因物料长期积压产生锈蚀、变形、变质及其 他损坏造成的损失。(4)提高可靠性，分巷道存放。仓库有多个巷道时，同种物 品分散在不同的巷道进行存放，以防止因某巷道堵塞影响某种物料的出摩，造成 生产中断。(5)提高效率，就近入出库。为保证快速响应出库请求，仓库～般将 物料就近放置在出库台附近。(6)产品相关性。商品相关性大的在订购时经常被 同时订购，所以应尽可能存放在相邻或相近的位置。 (二)货位的划分方法。企业一般是根据库区(房)的建筑形式、面积大小、库 房楼层或固定通道的分布和设施设备状况，结合储存货品需要的条件，将储存场 所划分为若干货区(库)，每一货区(库)再划分为若干货位，每一货位相对的固定 存放一类货品，或者几类品种相近的货物。赁区(库)的具体划分，通常以货场(或 库房)为单位，即把每一个独立的货场(或者每一座独立的仓库建筑)划分为一个 或几个货区(库)。 (1)临时存放区。在基本的货位划分号之后，从货场中开辟出一块空的区域， 用来存放已经到库，但还没有安排货位的物品，这个区域我们称为临时存放区。 临时存放区中的货物在安排好货位之后，可以一次到位将其安置在给定的货位，
z 1 q
z 1
xz
z
xz
(5)
1 x p S.t.1 z q
3、求解货位规划模型
货位规划是多目标优化问题，本文拟采用遗传算法求解，遗传算法是一种高
效、并行、随机的全局搜索算法，其中有许多解决多目标优化问题的方法，应用
较为普遍的是权重系数变化法，但是对于货位规划问题，由于两个目标函数的量
1、货位优化简述
货位优化(Slotting Optimization)是用来确定每一货品的恰当储存方式， 在恰当的储存方式下的空间储位分配。货位优化追求不同设备和货架类型特征、 货品分组、货位规划、人工成本内置等因素以实现最佳的货位布局，能有效掌握 商品变化，将成本节约最大化。因此，货位优化的最终目标就是通过提高生产力 和将无用的运动减到最少来降低成本货位优化是一个多目标优化问题，可以采用 多目标决策来解决。目前国内对货位优化系统(Slotting OptimizationSystem) 的研究还十分有限，关于sO的解决方案包括“分层序列法”和“线性加权和法”。 但这几种方法在进行货位优化时有以下不足之处：首先，在进行货位优化时主要 考虑两个因素：货物的周转率和重量。而对于其他因素的考虑很少，例如：货物 的分类存放、货物的存储要求等。特别是没有考虑到货物按照产品族聚集放置的 要求：其次，容易陷入局部优化，全局优化的效果不佳。因此，这些优化方法在 实际的应用中略显不足，针对药品仓库的实际情况，对于货位优化问题我们需要 一种新的算法。
t ax bz
xz vx vz 式中： a---------货位长度(m)；
（2）
b---------货位高度( m)；
x---------货物所在货位的列数；
z---------货物所在货位的层数；
vx--------叉车水平运行平均速度( m/s)；
vz--------叉车垂直升降平均速度( m/s)。
纲不同,采用此法并不合适，这里将遗传算法中的并列选择法、多性别进化法和
条件约束法相结合，求解的具体步骤如下:
(1) 编码
通常采用二进制字符串对个体(染色体)进行编码，但是这种编码方式不能较
好地适用于货位规划这个多目标优化问题，本文采用自然数的编码方式对货物进
行编码，表1给出了p列q层货架的货位编码顺序，最下面的一层作为第1 层，最
保证了货物一次作业的成功率，同时也减少翻堆作业，降低了作业成本。如果在 没有安排货位之前，就需要发货的物品，可以在其发货单上货位选项中选择“临 时存放区”以做区别，不会因为没有安排货位而不能发货(发货单申货位为关键 信息，必须填写)。(2)实货位。根据货位与出库台的距离远近将货位顺序排列编 号，形成一个货位的单链。为了平衡各专用线(或货架)的物料存放及龙门吊(或 堆跺机)的运行，还须对各专用线联合编号，形成了一条货位链。把按与出库台 的距离远近排列形成的货位链分成若干段，就形成了优先级不同的货位，称为实 货位。(3)虚货位。根据物料出入库的频度大小将物料排序，形成物料链。并根 据物料存取的方便程度以及离出货口的远近程度将物料链划分为优先级不同的 若干货位，这样重新编号的货位称为虚货位。它只是一个虚拟的货位，实际并不 存在。(4)虚实货位耦合的柔性货位。通过虚、实货位优先级的匹配，即出库频 度高的物料对应短程出库距离的货位，就可以完成货位的合理分配。实际上每个 实位的大小不是一成不变的，它是在耦合过程中依据虚货位的大小而定的。因此， 当某种物品的存货数量变化或分类变化时，货区的大小也随之变化。仓库进行货 位划分之后，每个货位相对都有其专门的服务对象，会使货位失去通用性，从而 降低了仓库的利用率和灵活性。为了弥补不足，可将相邻两个货位结为“兄弟货 位”。两相邻货位分别按不同的方向进行入库空货位的查找。在两分区之间的边 界地带的货位可以分别为两个分区的物料服务，即在必要的时候允许物料跨越边 界存放。这样，两分区之间的边界是模糊和可变的，这就大大增加了货位的柔性 消除了刚性货位的缺点。
为了缩短出入库搬运距离，提高作业效率，应将周转率高的货物放在靠近出 入口的货位上，使得
pq
f t F1=min x1 z1
xz xz
(1)
f 式中:
______ 第x列、z层上的货物的周转率(台(件)/月)：
xz
t ________第x列、z层上的货物运到出入口所需时间(s) 。 xz
将式(2)代入式(1)得：
f p q
F1=min
x1 z1
xz
ax vx
bz vz
(3)
(2)次目标函数F2
为使货架受力均匀，保证货架稳定，应该遵循下重上轻的原则, 保证货架的
整体重心最低。即:
M q
p
xz
z
F2=min z1
q
x 1
M z1
Matlab 大作业
智能存储过程优化设计
班 级： 10020741 姓 名： 学 号： 1002074139
2013 年 11 月
智能存储过程优化设计
从长远考察，企业竞争优势来源于比竞争对手以更低的成本、更快的速度去 发展自身的能力。在这个思想的指导下，企业通过提高自己的信誉、改善自己的 工作方法提高工作效率，吸引更多的客户。赢得更多的商业活动，实现最大利润。 为了实现更大的利润，企业除了在改善自身的经营方法，吸引更多的进驻客户， 实现最大的吞吐量。获得更多的利润之外，也可以通过降低自己的作业成本来缩 减开支，这就涉及到货位的优化问题。
xz
(4)
式中:Mxz________存储于第x列、z层货位上货物质量。