201 年 201 年度 第 学期
机时间
专业年级:数151
姓 名
王行舟 08
实验项目:数学建模算法应第37页题3
31某工厂户提供发动机合规定交货数量日期:第季度末交40台第二季度末交60台第三季度末交80台工厂生产力季100台季生产费f (x)50 x+02 x2(元)处x该季度生产发动机台数工厂生产余发动机移季户交货样工厂需支付存储费台发动机季存储费4元问该厂季应生产少台发动机满足交货合工厂花费费少(假定第季度开始时发动机存货)?
实验目:
1学会利数学建模思想实际问题建立成数学模型
2学会利MATLAB求解二次规划问题
3学会结果进行分析
实验求:(问题分析符号说明基假设)
问题分析:
二符号说明:
实验容:(数学模型理推导1建立程求解包括程序结果分析)
三分类模型求解
四编码程分析
()编写代码分析程
1处理关数
纯文文档cancerdata2中数包含容第列患者病例号病例号分类数关数进行更准确分析利程序语言数第列删
2找分类标签
数进行述1中处理数第列中分找出数值1代表良性肿瘤序号数值1代表恶性肿瘤序号数值1序号gind中数值1序号bind中
3贴标签
观察数现五百行第列11两肿瘤良性恶性进行分类数值利两种数值进行分类包含数值1gind分组样包含数值1bind分组
测数序号gind中输出1样序号bind中输出1判断69名做穿刺采样患者乳腺癌肿瘤良性恶性
4提出已知样测样
(1)测数进行分类首先已知样数调出进行分析前500行第二列开始列已知样细胞数training数提出
(2)样方法xa0提出测样数
5标准化
数分析前通观察发现患者项数中数差距容易结果产生误差消种误差利mapstd数进行标准化
6利支持量机进行分类
(二)编写代码
原始数cancerdatatxt网载数中B换成1M换成1X换成2删分割符*换数命名成cancerdata2txt
clcclear
aload('cancerdata2txt')
a(1)[] 删第列病例号
gindfind(a(1)1) 读出良性肿瘤序号
bind find(a(1)1) 读出恶性肿瘤序号traininga([1500][2end]) 提出已知样点数
traininga([1500][2end])提出已知样点数
training training' 进行数标准化里进行转置
[trainps]mapstd(training) 已分类数标准化
group(gind)1group(bind)1 已分类数标准化
group group' 转化成列量
xa0a([501569][2end]) 提出分类数
xa xa0'xamapstd('apply'xaps) 分类数标准化
ssvmtrain(train'group'Method''SMO''Kernel_Function''quadratic') 序列化方法训练支持量机分类器果二次规划方法训练支持量机法求解
sv_indexsSupportVectorIndices'返回支持量标号
betasAlpha'返回分类函数权系数
bsBias返回分类函数常数项
mean_and_std_transsScaleData第1行返回已知样点均值量相反数第2行返回标准差量倒数
cherksvmclassify(strain') 验证已知样点
err_rate1sum(groupcherk)length(group) 计算错判率
solutionsvmclassify(sxa') 进行判样点分类
solutionsolution'
sgfind(solution1) 求测样点良性编号
sbfind(solution1) 求测样点恶性编号
五结果分析
()结果输出
s
包含字段 struct
SupportVectors [73×30 double]
Alpha [73×1 double]
Bias 05826
KernelFunction @quadratic_kernel
KernelFunctionArgs {}
GroupNames [500×1 double]
SupportVectorIndices [73×1 double]
ScaleData [1×1 struct]
FigureHandles []
sv_index
1 25 列
4 5 11 30 32 41 50 63 64 68 74 81 82 87 90 92 98 100 101 103 107 110 112 120 129
26 50 列
136 137 142 147 148 152 168 173 177 192 193 195 197 198 205 206 216 226 240 256 262 264 266 276 291
51 73 列
292 298 300 330 341 364 368 373 378 397 411 414 422 448 457 458 463 469 482 485 491 492 497
beta
1 15 列
00026 00034 00038 00025 00068 04826 02543 00133 00020 00119 01618 00218 01449 00118 00523
16 30 列
00383 00025 01160 00046 00168 00422 00264 00261 00022 00489 04242 00049 00076 00014 00046
31 45 列
00080 01310 00191 00008 00277 00007 00192 00075 00167 00758 01061 01951 00063 00026 04825
46 60 列
00830 00823 00013 00153 00011 02943 00698 00003 00004 02350 01152 00743 00002 00840 00259
61 73 列
00656 00960 00628 01239 00156 00129 00283 00001 02418 00124 00214 00064 01885
b
05826
mean_and_std_trans
包含字段 struct
shift [1×30 double]
scaleFactor [1×30 double]
err_rate
0
solution
1 25 列
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
26 50 列
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
51 69 列
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
sg
1 25 列
1 3 6 7 8 9 11 12 14 16 20 21 23 24 25 26 28 30 31 32 33 35 38 39 40
26 44 列
41 42 43 44 45 46 47 48 49 50 51 52 57 58 59 60 61 62 69
sb
2 4 5 10 13 15 17 18 19 22 27 29 34 36 37 53 54 55 56 63 64 65 66 67 68
(二)分析结果
输出结果中知数误差分析0输入程序直观反映出数正确分类
n
皇问题算法设计分析回溯算法递更深理解提高
指导老师评语:
程序设计分析理解问题力加强动手力增加编程验
通
n
皇问题算法设计分析回溯算法递更深理解提高
程序设计分析理解问题力加强动手力增加编程验
文香网httpwwwxiangdangnet
《香当网》用户分享的内容,不代表《香当网》观点或立场,请自行判断内容的真实性和可靠性!
该内容是文档的文本内容,更好的格式请下载文档