Matlab语言程序设计分析.docxVIP

MATLAB 语言程序设计1、应用题（20 分）某工厂生产 A 和 B 两种产品，它们需要经过三种设备的加工，工时如表 1 所列。设备一、二、三每天可使用的时间分别不超过 12、10 和 8 小时。产品 A 和 B 的利润随市场的需求有所波动，如果预测未来某个时期内 A 和 B 的利润分别为 4000 元/吨和 3000 元/吨，问每天应安排产品 A、B 各多少吨，才能使工厂获利最大？表 1 生产产品工时表产品设备一设备二设备三A/（小时/吨）334B/（小时/吨）432设备每天最多可工作时数/小时12108请写出具体解决方案，并编写该解决方案的具体 MATLAB 程序代码。假设每天应安排生产产品A和B分别为和吨，依题意可建数学模型如下：max z=4000x1+ 3000x23x1+ 4x2 ≤123x1+ 3x2 ≤104x1+ 2x2≤8x1≥0, x2≥0首先将目标函数转换为标准形式:min z=-4000-3000M文件如下：f = [-4000;-3000];A = [3 4;3 3;4 2];b = [12;10;8];lb = zeros(2,1);[x,fval] = linprog(f,A,b,[],[],lb)结果如下图1所示：图12、结合 MATLAB 图像的读写知识，请你写出具体 MATLAB 图像读写方法，结合实例，给出图像读写步骤。（20 分）读写方法：imfinfo(FileNAme) %从FileNAme中获取或显示图像文件的特征数据[X,cmap]=imread(FileNAme) %读取变址图像的数据矩阵与伴随色图，并分别送给X和cmapX=imread(FileNAme) %从FileNAme中读取强度图像或真彩图像的数据矩阵并送给Ximwrite(X,cmap,FileName,Parameter,Value) %将变址图像写入文件imwrite（X ,FileName,Parameter,Value）%将强度图像或真彩图像写入文件[例2-1]M文件如下：I=imread(m2.png);imwrite(I,m2.tif);[X,cmap]=imread(m2.tif);image(X);colormap(cmap);axis image off;结果如下图2所示：图23、结合课本第四章所提供的18种绘图方法，请实现其中8种绘图，提供程序代码及所绘制的图形。（30分）（1）面域图面域图特别适合表达一个数据在总数中占据的比例。[例3-1]某公司1998年-2002年的销售额、成本、利润如下表1所示，试用面域图表达之。表1 某公司1998年-2002年的销售额、成本、利润一览表（单位：万元）年度19981999200020012002销售额520.8824.9931.7920.8901.8成本190.5367.9450.3403.2397.5利润330.3457.0481.3517.6504.3M文件如下：sales=[520.8 824.9 931.7 920.8 901.8];expence=[190.5 367.9 450.3 403.2 397.5];x=1998:2002;area(x,sales,FaceColor,[0.5 0.5 0.9],EdgeColor,b,LineWidth,5);hold on;area(x,expence,FaceColor,[0.5 0.9 0.5],EdgeColor,b,LineWidth,3);hold off;set(gca,XTick,[1998:2002]);set(gca,Layer,top);gtext(\leftarrow销售额);gtext(成本);gtext(利润);xlabel(年份,FontSize,18);ylabel(利润=销售额-成本,FontSize,18);结果如下图3所示，用鼠标选择合适的位置后，文字会自动标注：图3（2）饼图饼图便于比较各个元素在一个整体中的比例。[例3-2]某年级学生成绩为：优秀12人，良好87人，中等75人，及格64人，不及格20人。试利用饼图予以表达和分析。M文件如下所示：X=[21 87 75 64 20];subplot(2,2,1);pie(X,[1 0 0 0 0]);legend({优秀,良好,中等,及格,不及格});subplot(2,2,2);pie(X,[0 0 0 0 1]);subplot(2,2,3);pie3(X,[1 0 0 0 0]);subplot(2,2,4);pie3(X,[0 0 0 0 1]);饼图如下图4所示：图4（3）火柴杆图火柴杆图把图中每个数据点画成一根曲线，且数据点处用圆圈表示，看起来很像火柴杆。[例3-3]横坐标在0~2π之