2013数学建模

2013年数学建模B题思路

2013年数学建模B题程序求助高手解答 ....

2013高教社杯全国大学生数学建模竞赛B题
评阅要点[说明]本要点仅供参考，各赛区评阅组应根据对题目的理解及学生的解答，自主地进行评阅。
本题要求对数据提取合适的特征、建立合理有效的碎纸片拼接复原模型。可以考虑的特征有邻边灰度向量的匹配、按行或按列对灰度求和、行距等。关于算法模型，必须有具体的算法过程（如流程图、算法描述、伪代码等）及设计原理。虽然正确的复原结果是唯一的，但不能仅从学生提供的复原效果来评定学生解答的好坏，而应根据所建的数学模型、求解方法和计算结果（如复原率）三方面的内容做出评判。另一方面，评判中还需要考虑人工干预的多少和干预时间节点的合理性。问题1.仅有纵切文本的复原问题由于“仅有纵切”，碎纸片较大，所以信息特征较明显。一种比较直观的建模方法是：按照某种特征定义两条碎片间的（非对称）距离，采用最优Hamilton路或最优Hamilton圈（即TSP）的思想建立优化模型。关于TSP的求解方法有很多，学生在求解过程中需要注意到非对称距离矩阵或者是有向图等特点。还可能有种种优化模型与算法，只要模型合理，复原效果好，都应当认可。本问题相对简单，复原过程可以不需要人工干预，复原率可以接近或达到100%。问题2. 有横、纵切文本的复原问题一种较直观的建模方法是：首先利用文本文件的行信息特征，建立同一行碎片的聚类模型。在得到行聚类结果后，再利用类似于问题1中的方法完成每行碎片的排序工作。最后对排序后的行，再作纵向排序。本问题的解法也是多种多样的，应视模型和方法的合理性、创新性及有效性进行评分。例如，考虑四邻近距离图，碎片逐步增长，也是一种较为自然的想法。问题3.正反两面文本的复原问题这个问题是问题2的继续，基本解决方法与问题2方法相同。但不同的是：这里需要充分利用双面文本的特征信息。该特征信息利用得好，可以提升复原率。在阅卷过程中，可以考虑学生对问题的扩展。例如，在模型的检验中，如果学生能够自行构造碎片，用以检验与评价本队提出的拼接复原模型的复原效果，可考虑适当加分。阅卷时应有程序，程序的运行结果应和论文给出的结果一致。

clear %释放空间
clc %清屏
%图片数据读取
left_col = [];
right_col = [];
for fp = 0 : 208
str = int2str(fp);
if fp < 10
name = ['0' '0' str '.bmp'];
elseif fp >= 10 & fp < 100
name = ['0' str '.bmp'];
else
name = [str '.bmp'];
end
a = imread(name);
[m,n] = size(a);
left_col = [left_col a(:,1)];
right_col = [right_col a(:,n)];
end
%读取完毕
left_col = double(left_col);%类型转换
right_col = double(right_col);
% 找纸片最左边(left_col)像素全为255(空白)的所有列
row = 1;
for bi=1:209;
number=length(find(left_col(:,bi)==255));
if number == 180
S(row,1)=bi;%保存第一列像素为空(灰度值:255)的放在数组S第一列
row = row + 1;
end
end
S = [S(:,1) zeros(row-1,18)];%矩阵初始化
O = [ones(row-1,19)]; %初始化一个单位矩阵
sign = 1;
w = 0;
for r=1:row-1;%行
for p=1:18;%列p+1
num = 10000000000;%使num足够大
for j=1:209;
count = 0;
count = length(find(S==j));%除去重复
if count ~= 0
continue;
else
blank = length(find(right_col(:,S(r,p)) == 255));%如果碎纸片右边界全为255(即空白)，则跳出，终止此行后面拼接
if blank == 180
sign = 0;
break;%跳出本循环，进入p循环
else
ri=right_col(:,S(r,p));%计算左右拼接精确度
le=left_col(:,j);
c=ri-le;
c = c.^2;
error=sum(c(:));
end
if num >= error %找出差值最小的，精确度最高
num = error;
w = j;
end
end
end
if sign == 0
sign = 1;
break; %跳出p循环，进入r循环
else
S(r,p+1)=w; %二维数组储存每个碎纸片拼接位置
end
end
end
S = S - O; %数据整理，图片从000.bmp开始，数组下标从1开始

//////////////////////////////////////////////////////////////////////////
第三题碎纸片特征分类代码：
clear %释放空间
clc %清屏
%图片数据读取
char namea = (209,7);
char nameb = (209,7);
for fpa = 0 : 208
str = int2str(fpa);
if fpa < 10
fpa = fpa + 1;
namea(fpa,:) = ['0' '0' str 'a.bmp'];
elseif fpa >= 10 & fpa < 100
fpa = fpa + 1;
namea(fpa,:) = ['0' str 'a.bmp'];
else
fpa = fpa + 1;
namea(fpa,:) = [str 'a.bmp'];
end
end
for afp = 1:209
a= imread(namea(afp,:));
fdataa(:,:,afp) = a;
end
%%%读取反面b的数据
for fpb = 0 : 208
str = int2str(fpb);
if fpb < 10
fpb = fpb + 1;
nameb(fpb,:) = ['0' '0' str 'b.bmp'];
elseif fpb >= 10 & fpb < 100
fpb = fpb + 1;
nameb(fpb,:) = ['0' str 'b.bmp'];
else
fpb = fpb + 1;
nameb(fpb,:) = [str 'b.bmp'];
end
end
for bfp = 1:209
b= imread(nameb(bfp,:));
fdatab(:,:,bfp) = b;
end
%读取完毕
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
qfdataa = ~fdataa; %取反
qfdatab = ~fdatab; %取反
for lj = 1:209 %行累加求和
Ldataa(:,lj) = sum(qfdataa(:,:,lj),2); %正面(a)累加求和
Ldatab(:,lj) = sum(qfdatab(:,:,lj),2); %反面(b)累加求和
end
Ldataa(Ldataa>0)=1; %正面归一化
Ldatab(Ldatab>0)=1; %反面归一化
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 数据分类横向正面(a)分类
for flta = 1:209;
for pflta = 1:209
numa = 0;
for flha = 1:180;
if Ldataa(flha,flta) == Ldataa(flha,pflta)
numa = numa + 1;
end
end
tsavea(flta,pflta) = numa; %保存每两张图片之间的匹配度
end
end
% 数据分类横向反面(b)分类
for fltb = 1:209;
for pfltb = 1:209
numb = 0;
for flhb = 1:180;
if Ldatab(flhb,fltb) == Ldatab(flhb,pfltb)
numb = numb + 1;
end
end
tsaveb(fltb,pfltb) = numb; %保存每两张图片之间的匹配度
end
end
%%%%%%%%%%%%%%%%%%%%%%%
%总匹配度
%tsave = (tsavea + tsaveb)/2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 找纸片最左边(left_col)像素全为255(空白)的所有列
% fdataa = double(fdataa);
% fdatab = double(fdatab);
row = 1;
for bi=1:209;
number=length(find(fdataa(:,1,bi)~=0 & fdatab(:,72,bi)~=0));
if number == 180
S(row,1)=bi;%保存第一列像素为空(灰度值:255)的放在数组S第一列
row = row + 1;
end
end
S = [S(:,1) zeros(row-1,18)];%矩阵初始化
O = [ones(row-1,19)]; %初始化一个单位矩阵
%%%%%%%%%%%%%%%%%%

数学建模美赛ab离散题和连续题有什么区别

数学建模中的连续和离散题型一般指的是美赛中的AB题。

连续：问题中的变量因素是连续变化的，比如说今年的美赛A题中浴缸温度随时间连续变化
离散：问题中的变量因素是离散变化的，比如说今年美赛B题中的太空碎片的处理，每个碎片的处理都是离散的
连续和离散指的是问题本身的变量，对解题的影响不大，主要用来区分题型。

谁有2004年数学建模美赛的中文翻译题目？

A题题目是研究人的指纹相同导致确认身份时产生错误的可能性和因为DNA相同导致产生错误的可能性。这道题与生物有关。2004年，美国提出了“科技展望——生物盾国家纳米技术”，针对微分子进行研究，与人体相关的DNA也就成为焦点。指纹是区分人的重要参数，它由DNA决定，由于现在自然环境改变，人体的DNA也有可能发生变异，这样引起相同指纹出现也是有可能的。研究这种变异的可能性有利于刑事案件的侦破，提高政府的防御能力等。本题目可以采用机理分析法来建模。可以采用计算机模拟人体染色体中的基因排列进而检验、统计排列结果与指纹形成情况
B题题目：现在的快通系统在收费站、娱乐公园和其他的地方，正在越来越频繁的使用，来减少人们排队等候的时间，现在我们考虑为一个娱乐公园所设计的快通系统，在一次测试中，这个公园在几个游客比较多的景点旁都设置的快通系统，这个系统的设计创意是对于那些比较热门的景点，可以到旁边的一个机器，将其门票插入后出来一张纸条，上面写着具体的你可以回来时间段，比如说你把你的门票在1：15插到机子里，系统就会告诉你，你可以在3:30-4:30回来，这个时候队伍就比较短，你可以凭你的纸条加入这个队伍，很快就可以进入景点，为了防止游客同时在几个景点使用这个系统。系统的机器只允许你一次在一个景点排队等候。现在你是几个被公园雇佣的互相竞争的一个，你的职责是改善快通系统的运行。很多游客都在抱怨测试期间系统的异常现象，比如说有一次系统提供的回到景点的时间是4小时以后，但是才过一会儿，在相同景点系统提供的时间只有1小时。在另外一些时候根据快通系统组织起来的游客的等候队伍，就和普通的队伍一样长一样慢。现在的问题是要提出并且测试一个模型，这个模型能让快通系统的等候纸条的发放能增加人们在公园的乐趣的目的。问题的一部分就是首先决定衡量不同模型的标准，在你提交的报告里还要附带一份技术性的总结，以便公园的领导在不同的顾问所提出的模型中选择。分析：随着物质生活质量提高，人们开始追求精神生活了。周末或者空闲时间都喜欢去休闲娱乐场所，特别是一些贴近自然的景点，更是门庭若市。景点有限，特别是节假日，公园、动物园门口的排队现象严重影响人们欣赏自然景点的心情，所以构建一个快通系统来缩短排队时间，有利于人们利用自己有限的空闲时间放松自己、享受自然。动态规划

是这个吗？

电子设计大赛一等奖和数学建模一等奖哪个更重要我正

不存在哪个奖更重要的问题，关键要看个人兴趣与个人发展方向:
1、电子设计大赛更注重一个人的动手操作技能，而数学建模更注重一个人的对实际问题的分析能力。参加电子设计大赛一般需要比较强的实际操作技能;而参加数学建模的人一般有比较强的数学基本功，另外有很强的编程能力。你看看自己哪方面比较突出，可以适当进行选择。
2、另外参加电子设计大赛的人一般会在大学四年之后选择工作，而参加数学建模的人一般会往硕士甚至博士这样的级别进行攻读，因此在选择时可以与自己将来的意向结合起来。
3、参加电子设计大赛的人一般会话比较长的时间进行比赛，而相对来说，参加数模的时间周期会比较短一点，你可以根据自己的时间情况进行选择，但是可以肯定的是，在参加这些比赛时可能会花去你放假的时间，希望你有心理准备。
最后希望你不要抱着太功利的心态去参加这些比赛，大学更重要的是锻炼自己的能力，拿不拿到奖是次要的，只要能从这个过程中学到东西，你迟早能成为一个将来很厉害的人物。祝你好运!

2013大学生数学建模B题编程

求2013年数学建模B题答案给500分

fileform = 'E:\Desktop\B\附件2\*.bmp';
filepathsrc = 'E:\Desktop\B\附件2\';
file = dir(fileform);
A = cell(1,19);
B=cell(4,19);
for i = 1:19
A{i} = imread([filepathsrc, file(i).name]);
B{1,i}=A{i}(:,1);
B{2,i}=A{i}(:,72);
B{3,i}=A{i}(1,:); % shang
B{4,i}=A{i}(1,:)
end
C=zeros(19,19);
for i=1:19
for j=1:19
if(i~=j)
T=size(find(B{1,i}-B{2,j}==0));
C(i,j)=T(1,1);
end
end
end
C
A=C；
[Y_col,Ind_row]=max(A) %每列的最大值及行号
[Y_row,Ind_col]=max(A') %每行的最大值及列号

2013年数模国赛B题问题2代码，要能用的啊，后续再会追加100

addpath增加一条搜索路径 rmpath删除一条搜索路径
demo运行Matlab演示程序 type列出.M文件
doc装入超文本文档 version显示Matlab的版本号
help启动联机帮助 what列出当前目录下的有关文件
lasterr显示最后一条信息 whatsnew显示Matlab的新特性
lookfor搜索关键词的帮助 which造出函数与文件所在的目录
path设置或查询Matlab路径

figure(1);

%

建立图形窗口
1
[u,v] = meshgrid(-1:0.01:1);
%
生成二维频域网格
F1 = abs(sinc(u.*pi));
F2 = abs(sinc(v.*pi));
F=F1.*F2;
%
计算幅度频谱
F=|
F
(
u
v
)|
surf(u,v,F);
%
显示幅度频谱，如图
2.3(b)
shading interp;
%
平滑三维曲面上的小格
axis off;
关闭坐标系
figure(2);
%
建立图形窗口
2
F1=histeq(F);
%
扩展
F
的对比度以增强视觉效果
imshow(F1);
用图像来显示幅度频谱

高分求2014年全国大学生数学建模大赛E题的相关设计程序

sets:
var1/1..6/:y;
var2/1..2/:x;
endsets
min=x(1)+x(2)+y(1)+y(2)+y(3)+y(4)+y(5)+y(6);
47-12*x(1)>=10*y(1);
47-12*x(1)<=10*(y(1)+1);
8*(y(2)-1)<=47-12*x(1)-10*y(1)+33-6*x(1);
47-12*x(1)-10*y(1)+33-6*x(1)<=8*y(2);
10*y(3)<=31-12*x(2);
31-12*x(2)<=10*(y(3)+1);
8*(y(4)-1)<=42-6*x(2)+31-12*x(2)-10*y(3);
42-6*x(2)+31-12*x(2)-10*y(3)<=8*y(4);
8*(y(5)-1)<=50-(8*y(4)-(42-6*x(2)+31-12*x(2)-10*y(3)));
50-(8*y(4)-(42-6*x(2)+31-12*x(2)-10*y(3)))<=8*y(5);
8*(y(6)-1)<=41-(8*y(2)-(47-12*x(1)-10*y(1)+33-6*x(1)))-(8*y(5)-50+(8*y(4)-42-6*x(2)+31-12*x(2)-10*y(3)));
41-(8*y(2)-(47-12*x(1)-10*y(1)+33-6*x(1)))-(8*y(5)-50+(8*y(4)-42-6*x(2)+31-12*x(2)-10*y(3)))<=8*y(6);
x(1)+x(2)<=0.2*(y(1)+y(2)+y(3)+y(4)+y(5)+y(6));
@gin(x(1));
@gin(x(2));
@gin(y(1));
@gin(y(2));
@gin(y(3));
@gin(y(4));
@gin(y(5));
@gin(y(6));
end