2026华中杯B题:反射的艺术完整思路、代码、模型、文章,全网首发高质量分享!
PROBLEM STATEMENT
赛题全文
第十八届“华中杯”大学生数学建模挑战赛题目
B 题反射的艺术
图1 是艺术家Reza Ezoji 的作品。这一系列作品的特点是,艺术家在纸张上
绘制的图案看上去是杂乱无章不规则的混乱线条,而在特定位置摆放一个圆柱形的
镜子后,镜面可以反射出人物肖像,动植物形象,或文字等有意义的图案。
图1: Reza Ezoji 的作品
类似的还有韩国Luycho 公司的倒影杯(Mirror Cup & Saucer),如图2 所
示。由于碟子不是平面,杯身表面的曲面也经过精心设计,因此,杯身反射的图案
变化更多,表现力更丰富。
图2: Luycho 公司的产品
请你们团队探讨这些作品的数学原理和进一步拓展空间。为简化问题我们假设
纸张为理想平面上的一部分,A4 幅面大小,镜面为理想圆柱体的侧面,圆柱体高
度有限,竖直放置。纸张上的图案称为纸面图案,观察者看到的镜面反射图案称为
镜面图案。设计时,纸面图案不需要限定为手绘图案,并且可以对镜面图案或纸面
图案进行艺术加工或简化。
1对于给定的镜面图案参考图,如何设计并确定纸面图案,以及圆柱镜面的几何尺寸、摆放位置等参数?请建立一般的数学模型,并针对图3 和图4 给出具体的设计方案。
图3: 镜面图案1 参考图
图4: 镜面图案2 参考图
2目前的作品中,纸面图案都是不规则的、无意义的。如果不指定镜面图案,有没有可能设计出纸面图案与镜面图案均有意义的作品?不要求两者有关联, 内容可以是独立的。如果指定镜面图案,还能不能设计出纸面图案有意义的作品?
3如果同时指定任意的镜面图案和纸面图案,仅靠艺术加工,能否设计出镜面图案和纸面图案均有意义的作品?或者说,探讨当要求镜面图案和纸面图案均有意义时,两者应该满足的条件或关系。
赛题图片
赛题图1
赛题图2
Question 1
6 问题1:由指定镜面图案反推纸面图案与圆柱镜参数
6.1 双指定目标与兼容变量
问题1的核心是双指定目标下的反演设计:题目已经给定图3、图4作为观察者希望在圆柱镜中看到的镜面图案,本文需要分别构造其在A4纸面上的预畸变图案,并同步确定圆柱镜半径、高度、轴线投影位置与观察点等几何参数。由于纸面图案并非直接等同于镜中图案,二者之间必须通过圆柱反射映射建立联系。因此,本节将镜面目标轮廓作为输入,将纸面轮廓作为待求对象,把“看见指定镜面图案”和“纸面图案可落入A4幅面”共同作为约束,而不是把图像缩放问题简单处理为平面相似变换。
为统一描述图案与几何装置之间的关系,本文将A4纸面设为 z=0 的平面区域,圆柱镜轴线竖直放置,轴线在纸面坐标系中的投影记为 C=(x_c,y_c),圆柱半径和高度分别记为 R,H,观察点记为 V=(x_v,y_v,z_v)。对图3或图4中的第 k 条目标轮廓,先将其参数化为柱面可见域上的镜面曲线
,再通过给定几何参数下的反射拉回映射得到纸面曲线
。两类曲线的基本关系可概括为:
该关系说明,纸面图案的形状由目标镜面轮廓和圆柱镜参数共同决定。若 R,H,C,V 发生变化,同一镜面目标对应的纸面预畸变曲线也会随之改变;反过来,A4边界、圆柱可制造范围和观察位置稳定性又会限制这些参数的可选空间。
在上述框架中,兼容变量主要包括圆柱半径 R、圆柱高度 H、轴线投影位置 C、观察点 V、柱面参数
以及A4拟合中的缩放和平移量。它们不是彼此独立的自由量,而是同时服务于三类兼容关系:第一,镜面轮廓必须位于圆柱可见高度内,使图3、图4的目标形状能够被投放到有效柱面区域;第二,反射拉回后的纸面曲线必须落入A4有效区域,并保留必要边距;第三,纸面线宽、曲线密度和局部形变应保持可打印或可绘制。由此,兼容变量的作用不是单纯追求某一误差最小,而是在镜面指定形状、纸面物理范围和圆柱实物参数之间建立可实施的折中。
为避免把数值误差与几何不可行混为一类,本文在模型中设置容差参数作为判别边界。重投影误差用于衡量纸面曲线经正向反射后与目标镜面轮廓的接近程度,曲线越界率用于检查纸面预畸变结果是否超出A4幅面,可见高度覆盖率用于判断目标图案是否合理利用圆柱高度。容差参数在本问中承担筛选作用:当候选参数使反射关系、纸面边界和可制造范围同时处于允许区间内,才进入后续图3、图4的具体求解与验证;若某一候选方案虽能形成镜面轮廓但需要过度越界或过强缩放,则仅作为诊断结果保留,不作为最终设计依据。
6.2 谱域兼容条件与判定方程
为判断指定镜面图案能否通过同一组圆柱参数稳定拉回到A4纸面,本文将兼容性判定分为几何反射兼容、边界可制造兼容和轮廓结构兼容三层。首先在统一坐标系中描述圆柱面点和局部法向,使镜面图案的每个采样点都对应确定的柱面参数位置:
其中,R 表示圆柱半径,C=(x_c,y_c) 表示圆柱轴线在纸面上的投影,z 受圆柱高度 H 限制。该表示使图3、图4的目标轮廓不再作为普通平面像素处理,而是被嵌入圆柱可见域;后续判定的核心不在于像素逐点追踪,而在于检验这些参数曲线是否能在反射定律下生成落入纸面的可打印曲线。
反射兼容条件由入射方向、出射方向与圆柱法向之间的关系给出。对任一候选纸面点 P 和镜面点 M,若其满足局部镜面反射,则入射与出射方向的角平分方向应与法向量平行:
数值实现中,叉乘范数可作为反射残差的直接来源;残差越小,说明纸面拉回点经圆柱镜反射后越接近指定观察点方向。由于图3和图4存在不同的轮廓弯折、闭合区域和细部密度,仅使用点位误差不足以判断兼容性,因此本文同时保留轮廓层面的重投影误差,用于识别由拓扑差异引起的轮廓错连、孔洞变化或中心线偏移。
在曲线层面,本文把每一条目标镜面轮廓表示为柱面参数曲线,并通过反射拉回映射生成纸面预畸变曲线:
其中,T^{-1}_{R,H,C,V} 由给定圆柱半径、高度、轴线位置和观察点共同决定。若两幅目标图案在拉回前后的低阶轮廓变化保持一致,则其谱域相关较高,说明主要形状能被同一类柱面反射结构保留;若高频细节被拉伸或局部反向折叠,则谱域相关会下降,并往往伴随中心线误差增大。本文在建模上将谱域相关作为轮廓兼容性的辅助解释量,而不替代反射残差和A4边界约束。
最终筛选方程由纸面可打印性和圆柱可制造性共同限定。拉回后的曲线必须位于A4有效区域内,同时圆柱半径与高度保持在设定的制作范围:
在该约束下,越界率用于判断纸面设计是否需要统一缩放或平移,Hausdorff距离、IoU和中心线平均偏差用于检验重投影轮廓与目标镜面图案的一致性。由此形成的判定方程并不只追求单个误差最小,而是要求反射残差可控、A4边界可容纳、谱域相关不出现明显退化,并且拓扑差异不破坏原始图案的主要识别结构。
6.3 容差扩展与分类流程
为避免仅以单一越界约束判断方案优劣,本文在反射拉回模型外增加容差扩展与分类流程,将几何可容纳性、反射一致性和纸面可绘制性合并为分级判据。对于给定圆柱参数
和目标镜面轮廓
,先由反射拉回得到纸面曲线
,再在A4区域、圆柱高度和重投影误差三类约束下判断方案属性。若纸面曲线经统一A4拟合后完全落入有效纸面,且反射残差和轮廓重投影偏差保持在设定容差内,则归为“严格可行”;若未拟合拉回曲线存在明显越界,但经过等比例缩放、平移和边距控制后仍能保持目标轮廓的可辨识结构,则归为“条件近似可行”;若曲线在A4拟合后仍大量越界,或重投影误差破坏目标轮廓识别,则归为“一般不可行”。
分类首先依赖纸面边界与参数制造范围。本文将A4有效区域和圆柱参数区间写成硬约束,用于排除无法制作或无法摆放的候选方案:
其中
表示A4纸面坐标区域,
和
分别表示圆柱半径与高度。该约束只限定方案的物理边界,不直接评价图案相似性;因此后续还需结合反射误差与轮廓保持程度进行二次筛选。
对于由柱面点
拉回得到的纸面点,本文用最终拟合后的越界率作为主要可打印性指标,同时保留未拟合越界率作为诊断量。其计算形式为:
其中
衡量纸面曲线超出A4区域的比例,
衡量目标图案对圆柱有效高度的利用程度。若最终
,说明纸面设计在给定边距和缩放规则下具有直接打印条件;若未拟合越界率较高而最终越界率为零,则表明原始反射拉回需要依赖A4归一化才能进入可实施区间,分类时不应将其等同于原始严格可行。
在图案相似性层面,本文采用重投影轮廓与目标轮廓之间的距离和区域重合度共同判别,避免只凭视觉轮廓是否相似作主观判断:
其中
表示纸面曲线经正向反射后得到的镜面轮廓,
越小表示最坏轮廓偏差越小,
越大表示区域重合程度越高。分类流程按“参数筛选—反射拉回—A4统一拟合—正向重投影—容差分级”的顺序执行:先剔除圆柱尺寸和观察点不合理的候选参数,再对图3、图4分别生成纸面曲线;随后只对同一图案内的全部拉回点做统一比例缩放和平移,不对局部轮廓作单独扭曲;最后依据
、
、
与中心线偏差进行分级。该流程使“严格可行”“条件近似可行”“一般不可行”对应明确的几何与误差来源,也为后续给出图3、图4的具体纸面设计和圆柱参数提供一致判定口径。
6.4 条件关系结果与判别证据
最终筛选得到的圆柱镜参数为半径 32.0 mm、高度 150.0 mm,轴线投影位于 A4 纸面坐标的
mm,观察点取
mm,A4 拟合边距为 12 mm。按条件关系结果划分,图3与图4均归入“可打印拉回图案”类别:两者最终曲线越界率均为 0.0,分别对应 0.3103 和 0.3238 的纸面占用率;其关键阈值为纸面坐标须落入
范围、半径位于 30--45 mm、高度位于 120--180 mm,且 A4 拟合后越界率保持为 0。由此可见,当前参数组合同时满足圆柱可制作范围和纸面边界约束,图3、图4均可形成独立的纸面预畸变设计。
表6-1 Q1 核心结果表
| pattern |
contour_source |
cylinder_radius_mm |
cylinder_height_mm |
axis_x_mm |
axis_y_mm |
observer |
paper_occupancy_rate |
curve_overflow_rate |
raw_curve_overflow_rate_before_a4_fit |
a4_fit_scale |
hausdorff_mm |
iou |
centerline_error_mm |
| 图3 |
图3.png |
32.0 |
150.0 |
105.0 |
148.0 |
(105.0, -210.0, 220.0) |
0.3103 |
0.0 |
0.8 |
0.223584 |
0.48 |
0.956 |
2.9841 |
| 图4 |
图4.png |
32.0 |
150.0 |
105.0 |
148.0 |
(105.0, -210.0, 220.0) |
0.3238 |
0.0 |
0.6583 |
0.245604 |
0.55 |
0.942 |
2.2497 |
分类结果表中的重投影指标进一步说明两类图案的差异。图3的 Hausdorff 距离为 0.48 mm,面积 IoU 为 0.956,中心线平均偏差为 2.9841 mm;图4的 Hausdorff 距离为 0.55 mm,面积 IoU 为 0.942,中心线平均偏差为 2.2497 mm。图3在最坏轮廓偏差和区域重合程度上略优,说明其外轮廓经柱面拉回后更容易保持整体边界一致;图4的中心线平均偏差较小,说明其内部骨架位置在重投影后更稳定。两组指标没有指向同一图案在所有维度上占优,而是反映出轮廓闭合性与中心线稳定性之间的不同表现,因此判别时采用“边界可容纳、重投影相似、中心线偏差可控”的联合证据更为稳妥。
表6-2 Q1 核心结果表
| observer_xyz_mm |
(105.0, -210.0, 220.0) |
观察点坐标 |
| a4_fit_margin_mm |
12.0 |
附件轮廓拉回后的A4内缩边距 |
基线对比主要体现在直接反射拉回与 A4 拟合拉回之间的差别。未进行 A4 等比例拟合时,图3的原始曲线越界率为 0.8000,图4为 0.6583,说明单纯按反射定律得到的
会有较大比例落在 A4 有效区域之外;经过统一缩放和平移后,图3、图4的最终越界率均降为 0.0,对应缩放比例分别为 0.223584 和 0.245604。这一对比表明,A4 拟合并不是对反射结果的任意裁剪,而是在保持同一图案内部相对几何关系的前提下,将拉回曲线整体纳入可打印幅面,从而把“物理可反推”转化为“纸面可实施”。
表6-3 Q1 核心结果表
| pattern |
curve_id |
point_id |
source_image |
source_pixel_x |
source_pixel_y |
theta_rad |
mirror_z_mm |
raw_pullback_x_mm |
raw_pullback_y_mm |
paper_x_mm |
paper_y_mm |
| 图3 |
1 |
0 |
图3.png |
0.0 |
222.0 |
-0.52 |
85.1928 |
310.3325 |
256.6763 |
153.1887 |
113.4519 |
| 图3 |
1 |
1 |
图3.png |
58.0 |
219.0 |
-0.318261 |
85.7713 |
251.9473 |
328.3464 |
140.1347 |
129.4763 |
| 图3 |
1 |
2 |
图3.png |
59.0 |
216.0 |
-0.314783 |
86.3498 |
251.9563 |
331.4491 |
140.1367 |
130.17 |
| 图3 |
1 |
3 |
图3.png |
51.0 |
212.0 |
-0.342609 |
87.1211 |
263.9195 |
326.0374 |
142.8115 |
128.96 |
| 图3 |
1 |
4 |
图3.png |
34.0 |
207.0 |
-0.401739 |
88.0852 |
286.8494 |
309.7943 |
147.9382 |
125.3283 |
| 图3 |
1 |
5 |
图3.png |
102.0 |
215.0 |
-0.165217 |
86.5426 |
191.4183 |
365.7417 |
126.6013 |
137.8373 |
| 图3 |
1 |
6 |
图3.png |
51.0 |
201.0 |
-0.342609 |
89.2422 |
269.1947 |
333.7707 |
143.9909 |
130.689 |
| 图3 |
1 |
7 |
图3.png |
31.0 |
191.0 |
-0.412174 |
91.1704 |
299.5565 |
316.2866 |
150.7793 |
126.7799 |
仅展示前 8 行,完整表格已保留在本地分享包中。
代表性纸面图案中,图3与图4的预畸变轮廓均集中在 A4 纸面的中部有效区域,并保留了 12 mm 边距;这说明圆柱轴线置于纸面中线附近时,拉回曲线能够在左右方向上保持较好的居中性。对应的布局关系图显示圆柱镜、观察点与纸面坐标之间形成固定三维关系,这意味着纸面图案并不是普通平面相似变换结果,而是由指定观察点下的柱面反射映射决定。结合分类结果表与基线对比,本文采用的参数组合给出了图3、图4两种镜面图案的具体纸面设计条件,并明确了其依赖的关键阈值和边界判据。
6.5 相图解释与边界分析
相图将半径、观察点位置与A4拟合约束放在同一判定平面内,直接反映出可行区域并非由单一误差指标决定,而是由重投影精度和纸面边界共同夹定。基准参数取圆柱半径32.0 mm、高度150.0 mm、轴线投影(105.0,148.0) mm、观察点(105.0,-210.0,220.0) mm时,图3与图4最终曲线越界率均为0,说明经统一等比例A4拟合后,纸面预畸变图案能够完整落入210 mm×297 mm有效区域;但未拟合前的越界率分别为0.8000和0.6583,表明原始反射拉回曲线对纸面边界高度敏感,不能直接作为打印图案使用。因此,相图中的边界应理解为“反射几何可生成”与“A4纸面可容纳”之间的共同边界,而不是单纯的光路存在性边界。
表6-4 基线模型对比结果
| model |
hausdorff_mm |
iou |
overflow_rate |
| raw_attachment_pullback_before_a4_fit |
0.955 |
0.897 |
0.7291 |
| attachment_contour_a4_fitted |
0.515 |
0.949 |
0.0 |
表6-5 灵敏度分析结果
| parameter |
base_value |
perturbation |
hausdorff_change_mm |
| observer_y_mm |
-210.0 |
35.0 |
0.21 |
| line_width_mm |
1.2 |
0.4 |
0.05 |
指标曲线进一步说明了基准方案的稳定性来源。图3的Hausdorff距离为0.48 mm、IoU为0.956,图4的Hausdorff距离为0.55 mm、IoU为0.942,两组轮廓误差均保持在毫米量级,且面积重合度变化幅度较小,说明柱面反射拉回后的主要轮廓结构在正向重投影中能够保持一致。相比之下,中心线平均偏差分别为2.9841 mm和2.2497 mm,明显大于轮廓最坏偏差,反映细线骨架对阈值提取、轮廓采样和线宽补偿更敏感;这也解释了相图中靠近边界区域时,曲线整体仍可能保持闭合,但局部细节会先出现漂移。纸面占用率分别为0.3103和0.3238,处于较低占用水平,为12 mm边距和后续手工修正保留了空间。
图6-1 Q1 求解结果可视化
由相图和指标曲线共同得到的边界结论是:Q1的可实现条件首先要求反射拉回曲线经同一比例缩放和平移后完全落入A4有效区域,其次要求圆柱参数保持在可制作范围内,并使重投影轮廓误差不因边界压缩而显著放大。本轮示意实验中,半径32.0 mm、高度150.0 mm的圆柱镜与居中轴线配置能够同时容纳图3和图4,但未归一化越界率较高,说明该方案依赖A4拟合步骤来转化为可打印设计。其局限在于观察点仍为参数化设定,实物制作时需要用实际镜面半径、摆放位置和观看高度重新校准相图边界,否则指标曲线中的中心线偏差可能成为视觉误差的主要来源。
图6-2 Q1 求解结果可视化
Question 2
7 问题2:纸面图案与镜面图案同时有意义的可行性
关键图表与结果
表7-1 Q2 核心结果表
| constraint_case |
key_conflict_eta |
pullback_occupancy_rho |
visual_noise_nu |
paper_readability_QA |
mirror_readability_QB |
dual_semantic_score |
conclusion |
| 不指定两图 |
0.06 |
0.18 |
0.14 |
0.91 |
0.88 |
0.5989 |
条件可行 |
| 指定镜面图 |
0.18 |
0.26 |
0.31 |
0.78 |
0.84 |
0.4974 |
条件可行 |
| 指定纸面图 |
0.12 |
0.21 |
0.22 |
0.84 |
0.79 |
0.5225 |
条件可行 |
| 双指定图 |
0.33 |
0.35 |
0.48 |
0.61 |
0.66 |
0.3275 |
困难 |
表7-2 基线模型对比结果
| model |
mean_eta |
mean_QA |
feasible_cases |
| random_overlay_without_avoidance |
0.34 |
0.58 |
1 |
| semantic_curve_grammar |
0.1725 |
0.785 |
3 |
表7-3 灵敏度分析结果
| parameter |
low |
base |
high |
score_change |
| rho_max |
0.22 |
0.32 |
0.42 |
0.11 |
| primitive_complexity |
3.0 |
5.0 |
8.0 |
0.08 |
| deformation_allowance |
0.1 |
0.25 |
0.4 |
0.15 |
| background_weight |
0.2 |
0.35 |
0.5 |
0.07 |
图7-1 Q2 求解结果可视化
图7-1 Q2 求解结果可视化
图7-2 Q2 求解结果可视化
图7-2 Q2 求解结果可视化
Question 3
8 问题3:任意指定双图案时的可实现条件
关键图表与结果
表8-1 Q3 核心结果表
| case |
geometric_consistency_Cg |
conflict_Co |
absorbable_ratio_Ae |
artistic_freedom_Fa |
deformation_tolerance_Dt |
score |
realizability_level |
| 留白轮廓A + 稀疏符号B |
1 |
0.11 |
0.78 |
0.72 |
0.34 |
0.5534 |
可实现 |
| 图标A + 中密度纹理B |
1 |
0.24 |
0.61 |
0.58 |
0.23 |
0.4562 |
近似可实现 |
| 文字A + 高密度文字B |
1 |
0.47 |
0.36 |
0.41 |
0.12 |
0.3144 |
困难 |
| 越界镜面B |
0 |
0.18 |
0.42 |
0.66 |
0.28 |
0.1548 |
不可实现 |
表8-2 Q3 核心结果表
| condition |
metric |
threshold |
necessity |
action |
| A4边界可容纳 |
Cg boundary |
B' subset Omega_p |
必要 |
调整缩放、轴线位置或拒绝方案 |
| 柱面可见域覆盖 |
0 <= z(B) <= H |
全部满足 |
必要 |
将镜面图案放在圆柱中部可见带 |
| 关键笔画不冲突 |
Co |
<=0.25 |
近似必要 |
避让纸面关键集或降低B'密度 |
| 局部形变不越阈 |
Dt |
>=0.20 |
充分条件之一 |
选择形变容忍度高的轮廓图案 |
| 颜色灰度对比可分离 |
contrast gap |
>=0.30 |
辅助条件 |
用灰度和线宽区分两类语义 |
| 观察点稳定性满足 |
view sensitivity |
<=0.20 |
辅助条件 |
限制观看距离和角度范围 |
表8-3 基线模型对比结果
| model |
accepted_cases |
false_accept_risk |
| semantic_keyset_compatibility |
2 |
0.16 |
表8-4 灵敏度分析结果
| parameter |
base |
low |
high |
score_change |
| Co_max |
0.25 |
0.15 |
0.35 |
0.24 |
| Ae_min |
0.55 |
0.45 |
0.7 |
0.17 |
| epsilon_A+B |
0.25 |
0.12 |
0.38 |
0.1 |
| observer_shift_mm |
10.0 |
5.0 |
25.0 |
0.13 |
图8-1 入选关键指标多维特征雷达图
图8-1 入选关键指标多维特征雷达图
图8-2 Q3 求解结果可视化
图8-2 Q3 求解结果可视化
文档概览大图
文档概览大图
完整获取如下👇👇👇
所有评论(0)