1、 1 “拍照赚钱”的任务定价 摘要 “拍照赚钱”作为移动互联网下的一种自助式服务模式，用户在 APP 上领取拍照任务并执行，从而获得相应报酬。本文针对任务定价问题，用统计特性分析定价规律，得到距离价格比等模型，并进一步得到新项目的定价计划。 针对问题一，探究任务定价规律和任务执行情况，采用描述性统计量对已分成完成任务和未完成任务这两类的数据 进行 分析，同时对分类后数据做显著性差异分析，得到定价低的任务对应的 任务 完成率也低。用 SPSS 软件绘制两类位置与标价 的 三维立体图，用 Matlab 软件绘制标价与位置范围图，最终定价规律按位置范围可分为四类：北纬约 23至 23.08，东经约
2、113.1至 113.2；北纬约 23.1至 23.2，东经约 113.21至 113.5；北纬约 113.8至 114.1，东经约 22.5至 22.8；北纬约 22.8至 23.9，东经约 113.5至 113.8。这四个 范围分别对应佛山市、清远市、深圳市和东莞市。用二元 Logistic 回归模型得出任务未完成主要与位置有关，纬度越高，任务未完成可能性比完成可能性越大。此外，店铺拒访等原因也会造成任务未完成。 针对问题二，考虑到任务定价与位置 和执行情况有关，故采用聚类分析，按任务与领取该任务的会员间距离 将任务位置主要分为四类，建立距离价格比模型（ DPP 模型），求得 835 个任
3、务的具体定价 ；按任务完成率和定价之间关系，利用 0-1 整数线性规划，建立最小总定价 模型（ TRM 模型）， 同样 得到 每一个任务的具体定价 。最后得到原计划、按距离制定的计划和按完成率制定的计划三者对应的 APP 开发商需支 付的最小总定价分别为 36446 元、 60225 元和 33650 元。最后，结合具体内容分析可得 两个计划均比原计划合理。 针对问题三，由于 会员间会 竞争同一打包分布的任务，则这些会员之间存在非合作博弈关系，用 K-均值 聚 类得到 135 组新的打包分布后的任务，利用 序贯算法 改进问题二建立的 DPP 模型，得到打包任务定价模型（ PTP 模型）。 13
4、5 个任务集合所包含的任务数最多有 21 个，最少仅 2 个，故选择任务数为 10 个的第30 个任务集作为 PTP 模型的建立模板，求得 其 作为新任务的定价为 810 元，不打包情况下该任务集合的总定价为 815.72 元， 两者之间 误差小于 0.1%，说明 PTP模型合理。 针对问题四， 由新项目任务的位置分布图可得，新项目任务主要分布在已结束项目中未完成任务所在区域。为使 APP 开发商可以用最少的成本（任务总定价）得到最多的商检和信息数据（最高的任务完成率），建立双优化定价模型（ DOP模型）。用 Matlab 软件编程 求解得 任务 完成率为 80.1%， APP 开发商应给出的
5、最低总定价为 75827 元， 同时得到 2066 个任务的具体定价 。最后对该模型进行 10次模拟仿真，每次 模拟仿真所得值与实际值的误差都小于 5%，说明该定价计划的实施效果很好。 关键词： 显著性差异分析； K-均值 聚类； 0-1 整数规划；序贯算法 ； 双优化 2 一、 问题重述 当前盛行一种基于移动互联网自助式劳务众包平台 的赚钱方法 “拍照赚钱”。用户下载相应的 APP 并注册成为 APP 会员后，即可在 APP 上通过拍照获得标有具体定价的拍照任务，从而取得报酬。上述平台与传统市场调查方式相比，在为公司提供各类商务考察和信息搜索时，可大幅度减少调查成本。此外，该平台有效保证调查
6、数据的真实性，缩短调查周期。故 APP 变成该平台运行关键，且APP 的任务定价是其核心元素。若价格不合理，一些任务将被忽略，导致商品检查失败。 附件一是一组已结束项目 的 任务数据，包括 各项任务的位置、定价和完成情况（“ 0”为未完成，“ 1”为完成）；附件二是会员信息数据，包括 其位置、信誉值、根据其信誉给定预订任务限额及其开始时间，原则上信誉越高，会员越优先选择任务，配额越高（任务按照预订限额所占比例分配）；附件三是一组新的项目任务检验数据，仅包含任务的位置信息。请根据以上信息解决下述问题： 1. 根据附件一所给的项目任务定价，探究其规律性，分析任务未完成原因。 2. 针对附件一的项目
7、制定新的任务定价计划，并与原计划进行对比。 3. 多个任务可能由于位置较集中，在实际情况下会使得会员之间产生竞争。考虑将这些任务联合在一起打包发布时，该如何修改之前的定价模型，又 会对最终的任务完成情况有何影响？ 4. 针对附件三所给新项目建立任务定价计划，并对该计划的实施效果做出评价。 二、 问题分析 针对问题一，探究项目任务定价规律，属于寻找数据统计特性问题，一般选择用 SPSS 软件对所给数据进行一系列具体分析来解决 此类问题 。首先，按照附件一中任务执行情况一栏将任务标价分为两类，得到完成任务和未完成任务所对应的两组任务标价。分别求解两组任务标价所对应的描述性统计量。利用 M 估计和
8、K-S 检验对两组任务标价分别做探索性分析，判断其是否服从正态分布。其次，由于问题一还需考虑任务完成情况，故需比较两类 任务标价数据对任务完成情况有无显著性差异，则利用方差分析和独立样本 T 检验对分类后两组数据进行显著性差异分析。最后，通过对附件一中任务位置和任务标价进行对比分析，分别画出两类任务位置与标价的三维立体图来初步观察其关系，再进行聚类分析，得到任务位置与标价之间的具体对应几种分类关系。对于任务未完成 的 原因，执行情况取值只有 0和 1，故利用回归分析，建立任务位置与执行情况 间 的二元 Logistic回归模型，判断任务位置是否会造成任务未完成，还可借此模型预测其他位置的任务执
9、行情况。此外，查询一些具体的 “拍照赚钱” APP， 研究 他们对完成任务的具体要求，得出附件一无法给出的其他造成任务未完成的原因。 针对问题二，需制定新的任务定价计划，属于优化问题，解决该问题需找到每个任务的最优定价。由附件一可知任务标价与位置和执行情况都有关，故可按距离关系和任务完成率分别制定一组定价计划。按距离制定的计划，关键在于距离会员近的任务定价低些，而距离会员远的任务定价则按一定比例高些。由问题一得到定价按位置大致分为几类，将任务按这几类划分区域，分别算出各区域内每个任务的最优定价。按 任务 完成率制定的计划，关键在于新计划定价与任务执行情况之间的关系， 可用 0-1 整数线性规划
10、建立相应模型求解。最后，分别计算3 在原计划 和两组新的定价计划下，该平台一组项目需支付的总定价，比较其值大小。三者中总定价 少且任务完成率高的计划 为 最优 定价计划 。 针对问题三，只需考虑任务位置与定价之间的关系，故在问题二按距离关系所建立的定价计划 的 基础上做出改进即可。先将原来的 835 个任务按距离进行聚类分析，利用可打包任务间的距离范围确定聚类个数。对于仍是未打包的任务（单个任务）而言 ， 定价不变；对于打包在一起的多个任务 ， 可整体看成一个任务，聚类的中心 即 这个新任务的位置，即从问题二中的点与点之间距离变成点与集合之间的距离。题目中提到会员之间对任务有竞争关系，则此时的
11、距离不再是任务与最近会员间的距离， 此距离 还与时间有关，可以基于 序贯算法 （优先级的先后次序） 来改进定价模型。首先，对打包后的任务集合以 会员能接受的最远距离 为半径画圆，得到可能会竞争这个任务集合的会员集合。其次，每个会员都有其任务预定限额，若该任务集中任务个数超过某个会员的限额数，则该会员失去竞争力，从而缩小会员集合。每个会员的任务开始预定时间也不相同，挑选预定时间最早的会员得到进一步缩小的会员集。 最后，在上述会员集中按问题二的定价模型，找到与该任务集距离最短的会员，则这个任务集就被这个会员所领取了。按上述算法思想来修改问题三种按距离关系建立的定价计划。另一方面，打包后会员可选择的
12、总任务数减少，之前由于距离太远未被选择的任务可能会因此被没有抢到任务的会员选择，导致任务完成率增大； 而打包 被领取的任务的完成情况不受打包的影响 。 因此，整体任务 完成率增大。 针对问题四， 首先用 Matlab 软件画出附件三中 2066 个任务的位置分布 图，由此初步判断这些任务 的 可能执行情况。对于 APP 开发商而言，希望 在给出最少总 定价的同时满足最多的任务被会员领取，故问题四属于双目标优化 问题， 可用最优策略解决，建立双优化定价模型对新项目给出任务定价计划。对建立的模型进行模拟仿真，从而评价该计划的实施效果。 三、基本假设 1. 假设一个任务只能由一个会员领取，即不能被不
13、同会员重复领取。 2. 假设按一定周期发布一组新的项目任务。 3. 假设会员完全完成任务才能拿到标定的定价，否则得不到任何定价 。 4. 假设不考虑任务的难易程度对任务定价的影响。 5. 假设任务与任务之间、任务与会员之间的距离都为直线距离。 6. 假设打包发布的同一任 务集 合 中的每个任务的定价相同。 7. 假设不同项目的任务定价范围相同。 四 、 符号说明 符号 说明 P 未完成概率 p 每单位距离应增加的价格比例 d 任务位置与领取该任务的会员之间的距离 4 , 1, 2, ,835idi 第 i 个任务与领取该任务的会员之间的最短距离 , 1, 2, ,835ibi 已结束项目中 第
14、 i 个任务的定价 , 1, 2, ,835ixi 已结束项目中 第 i 个任务的执行情况 z 已结束项目的任务 总 定价 * , 1, 2, , 2066iyi 新项目中 第 i 个任务的定价 * , 1, 2, , 2066ixi 新项目中 第 i 个任务的执行情况 y 新项目的任务 总 定价 cp 新项目任务完成率 五、模型的建立与求解 5.1 问题一的模型建立与求解 5.1.1 数据预处理 附件一中罗列了一个已结束项目中 835 个任务的经度和纬度位置、任务标价以及任务的执行情况。在地图上找到具体任务位置主要分布在广东省的广州市、深圳市、佛山市、东莞市、惠州市和清远市，这些都是较发达城
15、市，故“拍照赚钱”任务发布较多，所给的酬金也比其他地区高。 按任务执行情况我们可以将这 835 个任务分为完成和未完成两类，其中完成的任务共 522 个，未完成任务共 313 个，分类后的数 据见附录 3.1。 问题一要探究项目的任务定价规律，可由分类后的任务标价的数字特征入手，故需对相应的统计量进行分析处理。 首先，描述性统计分析可 得出最直观的数据规律。均值、中位数和总和可描述任务定价的数据集中趋势，方差、标准差、极差可描述定价的离散程度，而偏度和峰度则可用来描述完成任务和未完成任务的总体分布形态，从而直观的观察其 是否服从正态分布。利用 SPSS 软件，我们得到两类任务标价所对应的各统计
16、量结果（见表 1），两者的频率直方图（见图 1），以及两类任务标价频率分布表（见附录 3.2）。 表 1 已完成和未完成任务的各统计量值表 N 平均值 标准平均 值误差 中位数 标准 偏差 方差 偏度 峰度 极差 最小值 最大值 平均值 95% 置信区间 下限值 上限 已完成 522 69.82 0.2109 68.833a 4.818 23 1.518 2.3 20 65 85 67.52 68.34 未完成 313 67.928 0.2071 66.282a 3.6647 13 1.965 4.4 20 65 85 69.41 70.23 由表 1 可得，已完 成的任务标价均值、中位数均高
17、于未完成任务的，说明标价高的任务其 完成率也 相对 较高。 5 图 1 完成和未完成任务的频率直方图 由图 1 可 看出两组数据均存在极端值，而 M 估计稳定性高，故用 M 估计值代替均值可得更精确的结果， 表 2 为 M 估计量 结果 表 。 表 2 M 估计量 任务执行情况 休伯 M 估计量 a Tukey 双权 b 汉佩尔 M 估计量 c 安德鲁波 d 任务标价 0 66.19 65.87 65.93 65.87 1 68.91 68.75 69.02 68.75 通过 K-S 检验得到常态性检验表（如表 3 所示）来检验两类任务标价的正态性，从中看出两类任务的显著性水平均为 0.000
18、，小于 0.05，故认为两类任务的数据均不服从正态分布，不满足方差分析的基本假设 1，则不能用方差分析来比较两类任务标价数据对任务完成情况有无显著性差异。 表 3 常态性检验表 任务执行情况 Kolmogorov-Smirnov(K)a Shapiro-Wilk 统计 df 显著性 统计 df 显著性 任务标价 0 .268 313 .000 .737 313 .000 1 .161 522 .000 .825 522 .000 5.1.2 任务定价规律探究 附件一给出了每个任务所对应的位置、标价和执行情况，故可通过分析任务位置与标价的关系、任务执行情况与标价的关系来探究任务定价规律。 首先，
19、观察附件一中数据，发现不论是完成的任务还是未完成的任务，相同的标价下的任务位置都相聚较近，故用 SPSS 软件分别绘制出完成和未完成任务的标价与位置关系图，如图 2 和图 3 所示。 图 2 完成任务的标价与位置关系图 图 3 未完成任务的标价与位置关系图 由图 2 和图 3 可初步得到任务标价与位置间的关系：标价在 65 元（最低价）6 附近的未完成任务较多，且按位置大致分为两类；对已完成的任务而言，任务位置与标价的分布较均匀。接下来对具体某一任务标价和位置进行聚类分析，得到与位置相关的任务定价规律。 由图 1 可得，最低的 4 个任务标价的未完成任务数也最多，用 Matlab 软件（代码见
20、附录 1.1）分别画出这四个标价所对应 的 任务位置范围图，从左往右分别对应于标价 65 元、 65.5 元、 66 元和 66.5 元，如图 4 所示。可得经纬度范围为北纬约 23至 23.08，东经约 113.1至 113.2；北纬约 23.1至 23.2，东经约 113.21至 113.5；北纬约 113.8至 114.1，东经约 22.5至 22.8的这三个地区的未完成任务定价最少，即按这三个地区分为三类。 图 4 四个最低标价对应任务位置范围图 另一方面，图 1 显示任务完成数最多的 7 个标价分别为 65 元、 65.5 元、 66元、 66.5 元、 70 元、 72 元和 75
21、 元，他们所对应的位置范围图见附录 2.1，分布较平均，说明对于已完成任务的定价规律而言， 其值与位置之间的关联 不大。故完成任务对应的位置范围约为北纬 22.8至 23.9，东经 113.5至 113.8。 图 5 任务定价按位置分类图 综合上述分析，任务定价规律按位置可 大致 分为四类，其中已完成任务的定1 1 3 1 1 3 . 2 1 1 3 . 4 1 1 3 . 6 1 1 3 . 8 1 1 4 1 1 4 . 2 1 1 4 . 42 2 . 52 2 . 62 2 . 72 2 . 82 2 . 9232 3 . 12 3 . 22 3 . 32 3 . 4经度纬度未完成任务
22、的标价为 6 5 . 5 元的任务位置散点图1 1 3 1 1 3 . 5 1 1 4 1 1 4 . 52 2 . 42 2 . 52 2 . 62 2 . 72 2 . 82 2 . 9232 3 . 12 3 . 22 3 . 3经度纬度未完成任务的标价为 66 元的任务位置散点图1 1 3 1 1 3 . 2 1 1 3 . 4 1 1 3 . 6 1 1 3 . 8 1 1 4 1 1 4 . 2 1 1 4 . 42 2 . 62 2 . 72 2 . 82 2 . 9232 3 . 12 3 . 22 3 . 3经度纬度未完成任务的标价为 65 元的任务位置散点图1 1 3 1
23、1 3 . 2 1 1 3 . 4 1 1 3 . 6 1 1 3 . 8 1 1 4 1 1 4 . 2 1 1 4 . 42 2 . 52 2 . 62 2 . 72 2 . 82 2 . 9232 3 . 12 3 . 22 3 . 32 3 . 4经度纬度未完成任务的标价为 6 6 . 5 元的任务位置散点图清远 佛山 东莞 深圳 1 1 3 . 5 1 1 4 1 1 4 . 5 2 2 . 5 1 1 3 2 2 . 9 2 3 . 3 经度 纬度 7 价为一类，未完成任务的定价按位置不同分为三类，大致分类情况如图 5 所示 。 其次，为比较两类任务执行情况对定价的影响规律，利用独
24、立样本 T 检验 2做方差齐性检验后，对其作出显著性差异分析，结果如表 4 所示。 表 4 独立样本检验表 列文方差相等性检验 平均值相等性的 t 检验 F 显著性 t 自由度 显著性 （双尾） 平均差 标准误差差值 差值的 95% 置信区间 下限 上限 任务标价 已假设方差齐性 15.511 .000 -5.985 833 .000 -1.892 .316 -2.512 -1.271 未假设方差齐性 -6.400 787.376 .000 -1.892 .296 -2.472 -1.312 根据表 4中的显著性为 0.000，小于显著性水平 0.05，故两类任务方差不等，在此情况下，独立样本
25、 T 检验的结果应看“未假设方差齐性”一栏，第 5 列为相应的双尾检测概率为 0.000，在显著性水平为 0.05 的情况下， T 统计量的概率值小于 0.05，即两类任务执行情况对定价有显著性差异。这与图 1 所示结果相符：虽然标价相同，但是不同任务执行情况下所对应的任务数量不同，且绝大多数情况下，完成任务的数量多于未完成任务的数量。然而，仅由附件一无法得出任务执行情况与定价之间更进一步的关系。 5.1.3 任务未完成原因 在 5.1.2 中，我 们已知仅由附件一无法得到任务执行情况与标价间的具体联系，只能初步得到相同标价下，未完成任务数量大多都少于完成任务数量。另一方面，发达地区的任务会相
26、对多些，则未完成任务也会随之增加；较偏远地区由于 APP 的注册会员较少，可能会导致该区域无人领取任务，从而无法完成任务，故可考虑任务的执行情况与其位置之间的关系。 5.1.3.1 二元 Logistic 回归模型的建立 任务的发布位置和某一任务的执行情况都是无法确定的，属于两组独立的变量，且两变量之间的关系是非确定性的，则利用回归分析来得到两者之间的数学表达式，从而定量地描述任务位置与执行情况间的相关关系。 由于位置是由经度和纬度两个量共同决定的，相当于有两个自变量，无法进行线性回归分析或曲线估计，但可为两自变量产生回归系数，通过这些系数来讨论任务位置与执行情况间关系。另外，我们要考察的是任
27、务未完成的原因，故所建模型的因变量为任务执行情况，其取值仅为 0 或 1。综合上述因素，建立任务执行情况与位置的二元 Logistic 回归模型 。 设 P 为任务未完成概率，其取值范围为 0 到 1，则 1P 为任务完成的概率，可用 Logistic 函数计算这个概率 P ，其表达式为 3 0 1 1 2 20 1 1 2 2e xp ,1 e xp yyP yy （ 1） 其中： 1y 和 2y 分别为一个任务的经度位置和纬度位置， 0 ， 1 和 2 为回归系数。 由于 Logistic 函数是协变的非线性函数，为求得 P 值，首先要求式（ 1）中未知的回归系数，故对式（ 1）进行 lo
28、git 变换，得到 0 1 1 2 2ln .1 P yyP（ 2） 8 表 5 方程式中的变量表 B S.E. Wald 自由度 显著性 Exp(B) 步骤 1a weidu 1.715 .309 30.763 1 .000 5.554 常量 -38.871 7.097 29.997 1 .000 .000 表 6 方程式中没有的变量表 得分 自由度 显著性 步骤 1 变量 jingdu 3.590 1 .058 整体统计信息 3.590 1 .058 利用 SPSS 软件求得所建模型的参数拟合表如表 5 和表 6 所示， 1y 没有出现在方程式中说明 其对应 系数为 0， 2y 的系数为
29、1.715，常数项为 -38.871， 将结果 带入式（ 1），故 为 探究任务未完成原因所建立的任务执行情况与位置的二元Logistic 回归模型为 22e xp 38.871 1.715 .1 e xp 38.871 1.715 yP y （ 3） 2y 的系数为正数，取相应指数后的值为 5.554，大于 1，表示纬度位 置越高，任务未完成的可能性 比任务完成的可能性 就 越大，与 5.1.2 的结果相符。 5.1.3.2 二元 Logistic 回归模型的求解 结合式（ 3），根据附件一中已有的位置与任务执行情况，可预测下一个任务的执行情况，如表 7 所示。可得预测未完成任务的正确率仅为
30、 25.9%，但对完成任务的预测正确率有 94.4%，总体预测正确率为 68.7%。说明任务位置的确会较大程度影响任务执行情况，且任务纬度位置越高，任务未完成性越大。 表 7 预测下一个任务执行情况表 观测值 预测值 任务执行情况 百分比正确 0 1 步骤 1 任务执行情况 0 81 232 25.9 1 29 493 94.4 总体百分比 68.7 从当前模型中除去经度位置这一变量后的改进模型，预测正确率不变，说明任务经度位置对任务执行情况的影响分布较均匀，即完成和未完成的可能性相近。 5.1.3.3 二元 Logistic 回归模型的检验 利用 SPSS 软件中的 Binary Logis
31、tic 模块对模型系数进行检验，具体系数检验结果见表 8。其中常数项系数的相伴概率为 1.668，若以 5%为置信，则常数项不显著； 1y 和 2y 的相伴概率分别为 0.000 和 0.031，均小于 0.05，通过检验，即两个位置变量对任务执行情况影响是显著的。 9 表 8 系数检验结果表 B S.E. Wald 自由度 显著性 Exp(B) 步骤 0 常量 .511 .071 51.187 1 .000 1.668 采用步与 步间、块与块间和模型间的相对似然比对模型全局进行检验，结果见表 9，三种检验方法结果一致，说明模型有显著的统计意义。 表 9 模型系数的 Omnibus 检验表 卡
32、方 自由度 显著性 步骤 1 步长 (T) 32.509 1 .000 块 32.509 1 .000 模型 32.509 1 .000 表 10 为模型拟合情况摘要表，主要给出 -2 对数似然值的两个决定系数，分别为 0.038 和 0.052，都很接近 0，说明模型拟合度较高。 表 10 模型拟合情况摘要表 步长 (T) -2 对数似然 Cox %绘制散点图 xlabel(经度 ); %作出 x 轴的标签 ylabel(纬度 ); %作出 y 轴额标签 title(未完成任务的标价为 65.5 元的任务位置散点图 ); %作出标题 grid on %加网格 clc clear scatte
33、r(x,y); %绘制散点图 xlabel(经度 ); ylabel(纬度 ); title(未完成任务的标价为 66 元的任务位置散点图 ) grid on clc clear scatter(x,y); %绘制散点图 xlabel(经度 ); ylabel(纬度 ); title(未完成任务的标价为 65 元的任务位置散点图 ) grid on clc clear scatter(x,y); %绘制散点图 xlabel(经度 ); ylabel(纬度 ); title(未完成任务的标价为 66.5 元的任务位置散点图 ) grid on 附录 1.2 选取 5 公里内的任务（ Matlab
34、 软件） ；对应支撑材料中的 i1.m clc clear scatter(x2,y2,+) %绘制会员位置的散点图 hold on %保留作图 Cxy=x2 y2 r=0.045*ones(1846,1); %确立半径长度 21 for i=1:length(r) %1846 次 for 循环 theta=0:360 %360 度圆 Circle1=Cxy(i,1)+r(i)*cos(theta); %以每一个会员位置为中心，半径 0.045 画圆 Circle2=Cxy(i,2)+r(i)*sin(theta); plot(Circle1,Circle2) hold on %保留作图 en
35、d scatter(x1,y1,d) %绘制任务位置的散点图 xlabel(经度 ); ylabel(纬度 ); legend(会员位置 ,任务位置 ); title(会员位置和没有执行的任务位置散点图 ) count=0; for i=1:835 %全部任务位置循环 835 次 d=sqrt(x1(i)-x2).2+(y1(i)-y2).2); %计算任务位置和会员位置的距离 q=min(d); %求解每一个任务位置和会员位置距离的最小值 if q0.045; count=count+1; %如果满足距离小于 0.045，则计数相加 e(i)=x1(i); %输出 x 坐标 f(i)=y1(
36、i); %输出 y 坐标 end end 附录 1.3 计算各任务与最近的会员 间的距离（ Matlab 软件） ；对应支撑材料中的 g1.m clc clear for i=1:815 d=sqrt(x(i)-x1).2+(y(i)-y1).2); %计算任务位置和会员位置的距离 q(i)=min(d) 附录 1.4 任务的定价和对应的执行情况代码（ Matlab 软件） ；对应支撑材料中的 f1.m clc clear scatter(x2,y2,+) %绘制会员位置的散点图 hold on Cxy=x2 y2 r=0.045*ones(1846,1); %确立半径长度 for i=1:l
37、ength(r) %1846 次 for 循环 theta=0:360 %360 度圆 Circle1=Cxy(i,1)+r(i)*cos(theta); %以会员位置为圆心，画圆 Circle2=Cxy(i,2)+r(i)*sin(theta); 22 plot(Circle1,Circle2) hold on %保留作图 end scatter(x1,y1,d) %绘制任务位 置的散点图 xlabel(经度 ); %作出 x 宙的标签 ylabel(纬度 ); %作出 y 轴的标签 legend(会员位置 ,任务位置 ); %两个图列 title(会员位置和没有执行的任务位置散点图 ) %
38、作出图像的标题 end 附录 1.5 线性规划计算定价最小（ Matlab 软件） ；对应支撑材料中的 j1.m clc clear for i=1:2066 y=randi(65,85,1,2066);%产生 2066 个在 65 元到 85 元随机的标价 x=randi(0,1,1,2066); %产生 2066 个执行情况 Y=y*x; %计算总价 X=sum(x); %计算执行的个数 XX(i)=min(X); %求解目标 YY(i)=max(Y); %求解目标 end 23 附录 2（图附录） 附录 2.1 任务完成数最多的 7 个标价所对应的位置范围图 1 1 3 1 1 3 .
39、5 1 1 4 1 1 4 . 52 2 . 52 2 . 62 2 . 72 2 . 82 2 . 9232 3 . 12 3 . 22 3 . 32 3 . 4经度纬度已完成任务的标价为 6 5 . 5 元的任务位置散点图1 1 2 . 8 1 1 3 1 1 3 . 2 1 1 3 . 4 1 1 3 . 6 1 1 3 . 8 1 1 4 1 1 4 . 2 1 1 4 . 42 2 . 5232 3 . 524经度纬度已完成任务的标价为 70 元的任务位置散点图1 1 2 . 6 1 1 2 . 8 1 1 3 1 1 3 . 2 1 1 3 . 4 1 1 3 . 6 1 1 3
40、. 8 1 1 4 1 1 4 . 22 2 . 62 2 . 8232 3 . 22 3 . 42 3 . 62 3 . 824经度纬度已完成任务的标价为 75 元的任务位置散点图1 1 3 1 1 3 . 2 1 1 3 . 4 1 1 3 . 6 1 1 3 . 8 1 1 4 1 1 4 . 2 1 1 4 . 42 2 . 42 2 . 52 2 . 62 2 . 72 2 . 82 2 . 9232 3 . 12 3 . 22 3 . 3经度纬度已完成任务的标价为 66 元的任务位置散点图1 1 3 1 1 3 . 2 1 1 3 . 4 1 1 3 . 6 1 1 3 . 8 1
41、 1 4 1 1 4 . 2 1 1 4 . 42 2 . 62 2 . 72 2 . 82 2 . 9232 3 . 12 3 . 22 3 . 32 3 . 42 3 . 5经度纬度已完成任务的标价为 72 元的任务位置散点图1 1 3 1 1 3 . 5 1 1 4 1 1 4 . 52 2 . 52 2 . 62 2 . 72 2 . 82 2 . 9232 3 . 12 3 . 22 3 . 32 3 . 4经度纬度已完成任务的标价为 65 元的任务位置散点图1 1 3 1 1 3 . 5 1 1 4 1 1 4 . 52 2 . 52 2 . 62 2 . 72 2 . 82 2
42、. 9232 3 . 12 3 . 22 3 . 32 3 . 4经度纬度已完成任务的标价为 6 6 . 5 元的任务位置散点图24 附录 3(表附录 ) 附录 3.1 任务表 附录 3.1.1 表 1 未完成任务表（截取部分） 任务号码 任务 gps 纬度 任务 gps 经度 任务标价 任务执行情况 A0001 22.56614225 113.9808368 66 0 A0002 22.68620526 113.9405252 65.5 0 A0004 22.56484081 114.2445711 75 0 A0005 22.55888775 113.9507227 65.5 0 A000
43、6 22.55899906 114.2413174 75 0 A0008 22.56277351 113.9565735 65.5 0 A0009 22.50001192 113.8956606 66 0 A0011 22.52486369 113.9308596 65.5 0 A0012 22.519087 113.9358436 65.5 0 A0018 22.4981901 113.8984817 66 0 A0023 22.51578568 113.9420557 65.5 0 A0024 22.54220963 114.0196493 66.5 0 A0025 22.54458961
44、 113.9986601 67.5 0 A0029 22.56470734 113.9820093 66 0 A0030 22.58193448 114.1563953 68 0 A0031 22.54293611 114.1314072 67 0 附录 3.1.2 表 2 完成任务表（截取部分） 任务号码 任务 gps 纬度 任务 gps 经度 任务标价 任务执行情况 A0003 22.57651183 113.957198 65.5 1 A0007 22.54900371 113.9722597 65.5 1 A0010 22.5437861 113.9239778 66 1 A0013
45、22.54797243 113.977909 65.5 1 A0014 22.50616871 113.9314284 66 1 A0015 22.49962566 113.9365145 66 1 A0016 22.54032142 113.9236456 66 1 A0017 22.52455419 113.9247319 65.5 1 A0019 22.54603946 113.9749684 65.5 1 A0020 22.49772892 113.9373377 66 1 A0021 22.49416247 113.927139 66 1 A0022 22.51592012 113.
46、9356769 65 1 A0026 22.49308313 113.9358391 66 1 A0027 22.52488949 113.9191429 65.5 1 A0028 22.54808582 113.9453119 66.5 1 A0033 22.55101305 113.9567446 65.5 1 A0034 22.5829933 114.1471227 66 1 25 附录 3.2 两类任务定价频率表 附录 3.2.1 表 3 未完成任务定价频率表 频率 百分比 有效百分比 累积百分比 有效 65.0 30 9.6 9.6 9.6 65.5 74 23.6 23.6 33.
47、2 66.0 57 18.2 18.2 51.4 66.5 28 8.9 8.9 60.4 67.0 20 6.4 6.4 66.8 67.5 6 1.9 1.9 68.7 68.0 5 1.6 1.6 70.3 68.5 5 1.6 1.6 71.9 69.0 7 2.2 2.2 74.1 69.5 2 .6 .6 74.8 70.0 19 6.1 6.1 80.8 70.5 2 .6 .6 81.5 71.0 1 .3 .3 81.8 71.5 1 .3 .3 82.1 72.0 18 5.8 5.8 87.9 72.5 2 .6 .6 88.5 73.0 3 1.0 1.0 89.5 7
48、3.5 2 .6 .6 90.1 74.0 4 1.3 1.3 91.4 74.5 1 .3 .3 91.7 75.0 19 6.1 6.1 97.8 80.0 4 1.3 1.3 99.0 85.0 3 1.0 1.0 100.0 总计 313 100.0 100.0 附录 3.2.2 表 4 已完成任务定价频率表 频率 百分比 有效百分比 累积百分比 有效 65.0 35 6.7 6.7 6.7 65.5 76 14.6 14.6 21.3 66.0 46 8.8 8.8 30.1 66.5 35 6.7 6.7 36.8 67.0 18 3.4 3.4 40.2 67.5 17 3.3
49、3.3 43.5 68.0 25 4.8 4.8 48.3 68.5 6 1.1 1.1 49.4 69.0 12 2.3 2.3 51.7 69.5 6 1.1 1.1 52.9 70.0 77 14.8 14.8 67.6 26 70.5 9 1.7 1.7 69.3 71.0 3 .6 .6 69.9 71.5 4 .8 .8 70.7 72.0 42 8.0 8.0 78.7 72.5 7 1.3 1.3 80.1 73.0 7 1.3 1.3 81.4 73.5 3 .6 .6 82.0 74.0 1 .2 .2 82.2 74.5 1 .2 .2 82.4 75.0 59 11.3 11.3 93.7 80.0 9 1.7 1.7 95.4
