【数学建模】实验设计方法

1.1 发展阶段

  • 早期、传统实验设计阶段
    费歇尔在农村进行田间试验的过程中,对高产小麦品种的遗传进行研究。为减少偶然因素对试验的影响,他对各种试验因素的每一水平组进行了试验,并通过方差分析评价指标的优劣(用于排除偶然因素的影响),使小麦大幅度增产。
    方差分析——分析因素重要性的方法
  • 中期:正交实验设计、回归实验设计为代表
  • 现代:均匀试验设计

1.2 正交实验设计

该方法使用一种规范化的表格进行实验设计,可以用较少的试验次数,取得较为准确、可靠的优选结论。
正交实验设计主要可以完成:

  1. 确定出各因素对试验指标的影响规律,得知哪些因素的影响是主要的、哪些因素的影响是次要的、哪些因素之间存在相互影响;
  2. 选出各因素的一个水平组合来确定最佳生产条件
  3. 正交实验的基础是正交表

正交实验法优点

  1. 试验点代表性强,试验次数少
  2. 不需做重复试验就可以估计试验误差
  3. 可以分清因素的主次
  4. 可以使用数理统计的方法处理试验结果,提出展望好条件

正交表

在这里插入图片描述
在这里插入图片描述

正交表特点

  1. 机会均等:各列中每个数字出现的次数相等
  2. 搭配均衡:任意两列中水平的不同搭配方式重复数相同(这两点称为正交性)
  3. 正交:对正交表进行初等变换,不改变其正交性
  4. 均衡分散:试验点在试验范围内排列规律整齐
  5. 整齐可比:试验点在试验范围内散布均匀
  6. 任意列中各水平重复出现的次数相等
  7. 任意两列所构成的水平对是完全有序数字对,各水平对重复出现的次数相等

基本流程

确定试验因素及水平数
选用合适的正交表
列出试验方案及试验结果
对试验结果进行极差分析和方差分析
确定最优或较优因素水平组合

极差分析

在完成试验收集完数据后,将要进行的是极差分析(也称方差分析)。
极差分析就是在考虑A因素时,认为其它因素对结果的影响是均衡的,从而认为,A因素各水平的差异是由于A因素本身引起的。

用极差法分析正交试验结果应引出以下几个结论:
①在试验范围内,各列对试验指标的影响从大到小的排队。
某列的极差最大,表示该列的数值在试验范围内变化时,使试验指标数值的变化最大。所以各列对试验指标的影响从大到小的排队,就是各列极差D的数值从大到小的排队。
②试验指标随各因素的变化趋势。
③使试验指标最好的适宜的操作条件(适宜的因素水平搭配)。
④对所得结论和进一步研究方向的讨论。

1.3 考虑交互作用的正交设计

在前面讨论中我们都是基于一个假设展开的,即在所有被考虑的对试验结果有影响的各因素之间对试验结果的影响是相互独立的,但是工程实践告诉我们这种情况很少出现,因此正交实验设计过程中考虑各因素的相互作用显得十分必要。

综合考虑实验目的、专业知识、以往经验及现有实验条件等多方面情况进行交互作用选择,一般原则是:
①忽略高级交互作用
②有选择地考察一级交互作用。通常只考察那些作用效果明显的,或试验要求必须考察的。
③试验允许的条件下,试验因素尽量取2水平

1.4 多指标试验

单指标试验:衡量试验效果的指标只有一个
多指标试验:衡量试验效果的指标有多个
多个指标之间又可能存在一定矛盾,一项指标好了,另一项指标却差了,这时需要兼顾各个指标,寻找使得每个指标都尽可能好的生产条件

在实际问题中,需要考虑的指标往往不止一个,有时是两三个甚至是多个,这都是多指标问题。解决多指标实验问题可采用两种方法:综合平衡法和综合评分法

1.5 正交实验的方差分析法

在这里插入图片描述
正交实验设计对于水平数较少时可以安排多因素
当水平数较多、因素个数少时使用均匀设计

2 均匀试验设计

均匀设计是只考虑试验点在试验范围内均匀散布的一种试验设计方法
均匀设计通过配套的均匀设计表和使用表来安排试验;当试验因素变化范围较大,需要取较多水平时,均匀设计可以极大减少试验次数。
对均匀试验所得的数据结果进行分析,可以判定所考察的因素中哪些是主要的,哪些是次要的,从而确定出最好的试验条件,得到最优方案。

优点:

  1. 试验次数大大减少
  2. 自动将各试验因素分类;自动将各试验因素分类为重要与次要,并将因素按重要性排序
  3. 过程数字化;通过电脑对结果与因素条件进行界定与预报(如天气预报),进而控制各因素。
  • 0
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
function ceshi_Excel %利用MATLAB生成Excel文档 % ceshi_Excel % % Copyright 2009 - 2010 xiezhh. % $Revision: 1.0.0.0 $ $Date: 2009/12/14 20:40:00 $ % 设定测试Excel文件名和路径 filespec_user = [pwd '\测试.xls']; % 判断Excel是否已经打开,若已打开,就在打开的Excel中进行操作,否则就打开Excel try % 若Excel服务器已经打开,返回其句柄Excel Excel = actxGetRunningServer('Excel.Application'); catch % 创建一个Microsoft Excel服务器,返回句柄Excel Excel = actxserver('Excel.Application'); end; % 设置Excel服务器为可见状态 Excel.Visible = 1; % set(Excel, 'Visible', 1); % 若测试文件存在,打开该测试文件,否则,新建一个工作簿,并保存,文件名为测试.Excel if exist(filespec_user,'file'); Workbook = Excel.Workbooks.Open(filespec_user); % Workbook = invoke(Excel.Workbooks,'Open',filespec_user); else Workbook = Excel.Workbooks.Add; % Workbook = invoke(Excel.Workbooks, 'Add'); Workbook.SaveAs(filespec_user); end % 返回当前工作表句柄 Sheets = Excel.ActiveWorkbook.Sheets; % Sheets = Workbook.Sheets; Sheet1 = Sheets.Item(1); % 返回第1个表格句柄 Sheet1.Activate; % 激活第1个表格 % 页面设置 Sheet1.PageSetup.TopMargin = 60; % 上边距60磅 Sheet1.PageSetup.BottomMargin = 45; % 下边距45磅 Sheet1.PageSetup.LeftMargin = 45; % 左边距45磅 Sheet1.PageSetup.RightMargin = 45; % 右边距45磅 % 设置行高和列宽 % 定义行高向量RowHeight RowHeight = [26,22,15,29,37,29,29,25,25,36,280,31,40,29,15,24]'; % 设置Range对象(从A1到A16)的行高 Sheet1.Range('A1:A16').RowHeight = RowHeight; % 设置Range对象(从A1到H1)的列宽 Sheet1.Range('A1:H1').ColumnWidth = [9,15,9,9,9,9,9,9]; % 合并单元格 Sheet1.Range('A1:H1').MergeCells = 1; Sheet1.Range('A2:H2').MergeCells = 1; Sheet1.Range('A8:A9').MergeCells = 1; Sheet1.Range('B8:D8').MergeCells = 1; Sheet1.Range('E8:H8').MergeCells = 1; Sheet1.Range('B9:D9').MergeCells = 1; Sheet1.Range('E9:H9').MergeCells = 1; Sheet1.Range('A10:H10').MergeCells = 1; Sheet1.Range('A11:H11').MergeCells = 1; Sheet1.Range('A12:H12').MergeCells = 1; Sheet1.Range('A13:H13').MergeCells = 1; Sheet1.Range('A14:H14').MergeCells = 1; Sheet1.Range('D16:H16').MergeCells = 1; % 设置单元格的边框 Sheet1.Range('A4:H14').Borders.Weight = 3; Sheet1.Range('A10:H12').Borders.Item(3).Linestyle = 0; Sheet1.Range('A10:H12').Borders.Item(4).Linestyle = 0; Sheet1.Range('A13:H13').Borders.Item(4).Linestyle = 0; Sheet1.Range('A14:H14').Borders.Item(3).Linestyle = 0; % 设置单元格对齐方式 Sheet1.Range('A1:H9').HorizontalAlignment = 3; Sheet1.Range('A4:A9').HorizontalAlignment = 6; Sheet1.Range('C4:C7').HorizontalAlignment = 6; Sheet1.Range('E4:E7').HorizontalAlignment = 6; Sheet1.Range('G4:G7').HorizontalAlignment = 6; Sheet1.Range('A10:H10').HorizontalAlignment = 6; Sheet1.Range('A11:H11').HorizontalAlignment = 6; % -4130 Sheet1.Range('A11:H11').VerticalAlignment = 1; Sheet1.Range('A12:H12').HorizontalAlignment = 4; Sheet1.Range('A13:H13').VerticalAlignment = 1; Sheet1.Range('A14:H14').HorizontalAlignment = 4; Sheet1.Range('D16:H16').HorizontalAlignment = 4; % 写入单元格内容 Sheet1.Range('A1').Value = '试 卷 分 析'; Sheet1.Range('A2').Value = '( 2009 — 2010 学年 第一学期)'; Sheet1.Range('A4:H4').Value = {'课程名称','','课程号',... '','任课教师学院','','任课教师',''}; Sheet1.Range('A5:H5').Value = {'授课班级','','考试日期',... '','应考人数','','实考人数',''}; Sheet1.Range('A6:H6').Value = {'出卷方式','','阅卷方式',... '','选用试卷A/B','','考试时间',''}; Sheet1.Range('A7:H7').Value = {'考试方式','','平均分',... '','不及格人数','','及格率',''}; Sheet1.Range('A8').Value = '成绩分布'; Sheet1.Range('B8').Value = '90分以上 人占 %'; Sheet1.Range('E8').Value = '80---89分 人占 %'; Sheet1.Range('B9').Value = '70---79分 人占 %'; Sheet1.Range('E9').Value = '60---69分 人占 %'; Sheet1.Range('A10').Value = ['试卷分析(含是否符合教学大纲、难度、知识覆'... '盖面、班级分数分布分析、学生答题存在的共性问题与知识掌握情况、教学中'... '存在的问题及改进措施等内容)']; Sheet1.Range('A12').Value = '签字 : 年 月 日'; Sheet1.Range('A13').Value = '教研室审阅意见:'; Sheet1.Range('A14').Value = '教研室主任(签字): 年 月 日'; Sheet1.Range('D16').Value = '主管院长签字: 年 月 日'; % 设置字号 Sheet1.Range('A4:H12').Font.size = 10.5; % 设置单元格A4至H12的字号为10.5 Sheet1.Range('A1').Font.size = 16; % 设置单元格A1的字号为16 Sheet1.Range('A1').Font.bold = 2; % 单元格A1的字体加粗 % 插入图片,如果当前工作表中有图形存在,通过循环将图形全部删除 Shapes = Sheet1.Shapes; % 返回第1个工作表的Shapes接口的句柄 if Shapes.Count ~= 0; for i = 1 : Shapes.Count; Shapes.Item(1).Delete; % 删除第1个Shape对象 end; end; % 产生均值为75,标准差为4的正态分布随机数,画直方图,并设置图形属性 zft = figure('units','normalized','position',... [0.280469 0.553385 0.428906 0.251302],'visible','off'); % 新建图形窗口,设为不可见 set(gca,'position',[0.1 0.2 0.85 0.75]); % 设置坐标系的位置和大小 data = normrnd(75,4,1000,1); % 产生均值为75,标准差为4的正态分布随机数 hist(data); % 绘制正态分布随机数的频数直方图 grid on; % 添加参考网格 xlabel('考试成绩'); % 为X轴加标签 ylabel('人数'); % 为Y轴加标签 hgexport(zft, '-clipboard'); % 将直方图复制到剪贴板 % 选中工作表Sheet1的A11单元格,插入由MATLAB命令作出的直方图 Sheet1.Range('A11').Select; Sheet1.Paste % Sheet1.PasteSpecial; delete(zft); % 删除图形句柄 Workbook.Save % 保存文档

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Lydia.na

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值