查看: 5273|回复: 2

2020全国大学生数学建模A题思路讲解与核心代码

字体大小: 正常放大

杨利霞

5273 主题	82 听众	17万积分

TA的每日心情

	开心 2021-8-11 17:59

签到天数: 17 天

[LV.4]偶尔看看III

网络挑战赛参赛者

自我介绍: 本人女，毕业于内蒙古科技大学，担任文职专业，毕业专业英语。

群组: 2018美赛大象算法课程

群组: 2018美赛护航培训课程

群组: 2019年数学中国站长建

群组: 2019年数据分析师课程

群组: 2018年大象老师国赛优

电梯直达

1^#

发表于 2021-7-16 16:32 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

2020全国大学生数学建模A题思路讲解与核心代码

2020全国大学生数学建模A题思路讲解与核心代码
题目
核心方法：
问题一
问题二
问题三和问题四
答案如下：
题目

核心方法：
热传导
有限差分法
遍历法

问题一
建立焊接区域中心温度变化规律模型，推出焊接区域中心温度与其厚度和PCB 板所走过的时间的关系。查阅相关资料可知，由于自动焊接过程中热量传递复杂，因此对模型进行简化，只考虑一维方向的热量传导，即单侧单方向小温区对PCB 板的热量传导。利用能量守恒定律和 Fourier 热传导定律推出热传导方程，再利用附数据件求出方程中的参数，进而建立了焊接区域中心温度变化规律型，即炉温曲线变化模型。依据建立出的炉温曲线变化模型，根据问题一中所给出的各温区的温度参数T1， T2， T3， T4 及过炉速度v，需要求出过炉曲线，即焊接区中心的温度变化

对于热传导方程的求解，需要先确定热传导方程中的参数—热扩散率，这
可以通过附件提供的炉温曲线数据进行参数估计。热传导方程的求解可以利用差分法进行。

// lamda的计算的部分代码
array=zeros(76,length(x1));
array(1,

=y;
array(:,1)=z(:,1);
for k=1:31
for j=1

(1)-1
      for i=2:75
         array(i,j+1)=array(i,j)+u(k)*(array(i-1,j)-2*array(i,j)+array(i+1,j));
      end
      array(76,j+1)=array(74,j+1);
end
e1=1

(1);
e2=time(1:5,

;
[C,ia,ib]=intersect(e1,e2*100);
for i=1:5
b(i)=array(75,ia(i));
end
for i=1:5
c(i)=(temperature(i)-b(i))^2;
end
rss(k)=sum(c(

);
end
result=[u;rss];
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
有限差分的核心代码：

//有限差分的核心代码
array=zeros(76,length(x1));
array(1,

=y;
array(:,1)=z(:,1);
for j=1

(1)
for i=2:75
array(i,j+1)=array(i,j)+u(1)*(array(i-1,j)-2*array(i,j)+array(i+1,j));
end
array(76,j+1)=array(74,j+1);
end
z(:,2)=array(:,2143);
for k=1:9
for j=L(k)

(k+1)
      for i=2:75
         array(i,j+1)=array(i,j)+u(k+1)*(array(i-1,j)-2*array(i,j)+array(i+1,j));
      end
         array(76,j+1)=array(74,j+1);
end
end
array(:,length(array))=[];

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
得模拟数据和真实数据对比得炉温曲线：

问题二
问题二中，基于问题一中所建立的炉温曲线模型，在四个温度参数给定的条件求取传送带的最大过炉速度为优化问题。此问题可以看做是问题一所建立模型的反问题，即在温度分布
已知的条件下，要求通过该分布计算最大过炉速度v。在具体求解该反问题时，可以利用遍历法对过炉速度进行遍历搜索，这样就将反问题转化为了正问题的求解，从而问题一中模型方法都可以继续使用。

问题三和问题四
问题三和问题四仍然和问题二类似，也是对过炉曲线提出了不同的要求，进而在这些要求之下确定影响炉温曲线的 5 个参数 T1， T2， T3， T4， v ，求解也可以采用与问题二相同的遍历法进行，但由于此时遍历的变量个数增多，如果遍历步长较小，必然会使得计算量增大，因而必要情况下，可采用分阶段的遍历，即：大范围，大步长，小范围，小步长。需要考虑的就是对于面积和对称性的数学描述，面积可以采用积分的离散化表示，对称性可以采用以最大峰值温度两侧取对称点，使对称点的温度差值尽可能小来实现。

答案如下：

注：用以上方法算出的结果均在最优解范围内，详细解读请待下次，困了该sleep了
————————————————
版权声明：本文为CSDN博主「盈博简」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/Sandm_Hou/article/details/112726635

zan