2026华中杯B题：反射的艺术完整思路、代码、模型、文章，全网首发高质量分享！

weixin_52051317

12人浏览 · 2026-05-15 02:42:56

weixin_52051317 · 2026-05-15 02:42:56 发布

2026华中杯B题：反射的艺术完整思路、代码、模型、文章，全网首发高质量分享！

赛题全文

第十八届“华中杯”大学生数学建模挑战赛题目

B 题反射的艺术

图1 是艺术家Reza Ezoji 的作品。这一系列作品的特点是，艺术家在纸张上

绘制的图案看上去是杂乱无章不规则的混乱线条，而在特定位置摆放一个圆柱形的

镜子后，镜面可以反射出人物肖像，动植物形象，或文字等有意义的图案。

图1: Reza Ezoji 的作品

类似的还有韩国Luycho 公司的倒影杯（Mirror Cup & Saucer），如图2 所

示。由于碟子不是平面，杯身表面的曲面也经过精心设计，因此，杯身反射的图案

变化更多，表现力更丰富。

图2: Luycho 公司的产品

请你们团队探讨这些作品的数学原理和进一步拓展空间。为简化问题我们假设

纸张为理想平面上的一部分，A4 幅面大小，镜面为理想圆柱体的侧面，圆柱体高

度有限，竖直放置。纸张上的图案称为纸面图案，观察者看到的镜面反射图案称为

镜面图案。设计时，纸面图案不需要限定为手绘图案，并且可以对镜面图案或纸面

图案进行艺术加工或简化。

1. 对于给定的镜面图案参考图，如何设计并确定纸面图案，以及圆柱镜面的

几何尺寸、摆放位置等参数？请建立一般的数学模型，并针对图3 和图4 给出具体

的设计方案。

图3: 镜面图案1 参考图

图4: 镜面图案2 参考图

2. 目前的作品中，纸面图案都是不规则的、无意义的。如果不指定镜面图

案，有没有可能设计出纸面图案与镜面图案均有意义的作品？不要求两者有关联，

内容可以是独立的。如果指定镜面图案，还能不能设计出纸面图案有意义的作品？

3. 如果同时指定任意的镜面图案和纸面图案，仅靠艺术加工，能否设计出镜

面图案和纸面图案均有意义的作品？或者说，探讨当要求镜面图案和纸面图案均有

意义时，两者应该满足的条件或关系。

赛题图片

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 条目标轮廓，先将其参数化为柱面可见域上的镜面曲线 $\gamma_m^k(s)$ ，再通过给定几何参数下的反射拉回映射得到纸面曲线 $\gamma_p^k(s)$ 。两类曲线的基本关系可概括为：

$\gamma_p^k(s)=T^{-1}_{R,H,C,V}(\gamma_m^k(s))$

该关系说明，纸面图案的形状由目标镜面轮廓和圆柱镜参数共同决定。若 R,H,C,V 发生变化，同一镜面目标对应的纸面预畸变曲线也会随之改变；反过来，A4边界、圆柱可制造范围和观察位置稳定性又会限制这些参数的可选空间。

在上述框架中，兼容变量主要包括圆柱半径 R、圆柱高度 H、轴线投影位置 C、观察点 V、柱面参数 $\theta,z$ 以及A4拟合中的缩放和平移量。它们不是彼此独立的自由量，而是同时服务于三类兼容关系：第一，镜面轮廓必须位于圆柱可见高度内，使图3、图4的目标形状能够被投放到有效柱面区域；第二，反射拉回后的纸面曲线必须落入A4有效区域，并保留必要边距；第三，纸面线宽、曲线密度和局部形变应保持可打印或可绘制。由此，兼容变量的作用不是单纯追求某一误差最小，而是在镜面指定形状、纸面物理范围和圆柱实物参数之间建立可实施的折中。

为避免把数值误差与几何不可行混为一类，本文在模型中设置容差参数作为判别边界。重投影误差用于衡量纸面曲线经正向反射后与目标镜面轮廓的接近程度，曲线越界率用于检查纸面预畸变结果是否超出A4幅面，可见高度覆盖率用于判断目标图案是否合理利用圆柱高度。容差参数在本问中承担筛选作用：当候选参数使反射关系、纸面边界和可制造范围同时处于允许区间内，才进入后续图3、图4的具体求解与验证；若某一候选方案虽能形成镜面轮廓但需要过度越界或过强缩放，则仅作为诊断结果保留，不作为最终设计依据。

6.2 谱域兼容条件与判定方程

为判断指定镜面图案能否通过同一组圆柱参数稳定拉回到A4纸面，本文将兼容性判定分为几何反射兼容、边界可制造兼容和轮廓结构兼容三层。首先在统一坐标系中描述圆柱面点和局部法向，使镜面图案的每个采样点都对应确定的柱面参数位置：

$M(\theta,z)=(x_c+R\cos\theta,\ y_c+R\sin\theta,\ z)$

$n(\theta)=(\cos\theta,\sin\theta,0)$

其中，R 表示圆柱半径，C=(x_c,y_c) 表示圆柱轴线在纸面上的投影，z 受圆柱高度 H 限制。该表示使图3、图4的目标轮廓不再作为普通平面像素处理，而是被嵌入圆柱可见域；后续判定的核心不在于像素逐点追踪，而在于检验这些参数曲线是否能在反射定律下生成落入纸面的可打印曲线。

反射兼容条件由入射方向、出射方向与圆柱法向之间的关系给出。对任一候选纸面点 P 和镜面点 M，若其满足局部镜面反射，则入射与出射方向的角平分方向应与法向量平行：

$d_{in}=\frac{P-M}{\|P-M\|},\quad d_{out}=\frac{V-M}{\|V-M\|}$

$(d_{in}+d_{out})\times n(\theta)=0$

数值实现中，叉乘范数可作为反射残差的直接来源；残差越小，说明纸面拉回点经圆柱镜反射后越接近指定观察点方向。由于图3和图4存在不同的轮廓弯折、闭合区域和细部密度，仅使用点位误差不足以判断兼容性，因此本文同时保留轮廓层面的重投影误差，用于识别由拓扑差异引起的轮廓错连、孔洞变化或中心线偏移。

在曲线层面，本文把每一条目标镜面轮廓表示为柱面参数曲线，并通过反射拉回映射生成纸面预畸变曲线：

$\gamma_m^k(s)=M(\theta_k(s),z_k(s))$

$\gamma_p^k(s)=T^{-1}_{R,H,C,V}(\gamma_m^k(s))$

其中，T^{-1}_{R,H,C,V} 由给定圆柱半径、高度、轴线位置和观察点共同决定。若两幅目标图案在拉回前后的低阶轮廓变化保持一致，则其谱域相关较高，说明主要形状能被同一类柱面反射结构保留；若高频细节被拉伸或局部反向折叠，则谱域相关会下降，并往往伴随中心线误差增大。本文在建模上将谱域相关作为轮廓兼容性的辅助解释量，而不替代反射残差和A4边界约束。

最终筛选方程由纸面可打印性和圆柱可制造性共同限定。拉回后的曲线必须位于A4有效区域内，同时圆柱半径与高度保持在设定的制作范围：

0\le x_p(s)\le210,\quad 0\le y_p(s)\le297

30\le R\le45,\quad 120\le H\le180

在该约束下，越界率用于判断纸面设计是否需要统一缩放或平移，Hausdorff距离、IoU和中心线平均偏差用于检验重投影轮廓与目标镜面图案的一致性。由此形成的判定方程并不只追求单个误差最小，而是要求反射残差可控、A4边界可容纳、谱域相关不出现明显退化，并且拓扑差异不破坏原始图案的主要识别结构。

6.3 容差扩展与分类流程

为避免仅以单一越界约束判断方案优劣，本文在反射拉回模型外增加容差扩展与分类流程，将几何可容纳性、反射一致性和纸面可绘制性合并为分级判据。对于给定圆柱参数 (R,H,C,V) 和目标镜面轮廓 $\Gamma_m$ ，先由反射拉回得到纸面曲线 $\Gamma_p$ ，再在A4区域、圆柱高度和重投影误差三类约束下判断方案属性。若纸面曲线经统一A4拟合后完全落入有效纸面，且反射残差和轮廓重投影偏差保持在设定容差内，则归为“严格可行”；若未拟合拉回曲线存在明显越界，但经过等比例缩放、平移和边距控制后仍能保持目标轮廓的可辨识结构，则归为“条件近似可行”；若曲线在A4拟合后仍大量越界，或重投影误差破坏目标轮廓识别，则归为“一般不可行”。

分类首先依赖纸面边界与参数制造范围。本文将A4有效区域和圆柱参数区间写成硬约束，用于排除无法制作或无法摆放的候选方案：

$\Omega_p=[0,210]\times[0,297]$

30\le R\le45,\quad 120\le H\le180

其中 $\Omega_p$ 表示A4纸面坐标区域，和分别表示圆柱半径与高度。该约束只限定方案的物理边界，不直接评价图案相似性；因此后续还需结合反射误差与轮廓保持程度进行二次筛选。

对于由柱面点 $M(\theta,z)$ 拉回得到的纸面点，本文用最终拟合后的越界率作为主要可打印性指标，同时保留未拟合越界率作为诊断量。其计算形式为：

$\rho_{out}=\frac{\mathrm{len}(\Gamma_p\setminus\Omega_p)}{\mathrm{len}(\Gamma_p)}$

$\rho_H=\frac{z_{max}-z_{min}}{H}$

其中 $\rho_{out}$ 衡量纸面曲线超出A4区域的比例， $\rho_H$ 衡量目标图案对圆柱有效高度的利用程度。若最终 $\rho_{out}=0$ ，说明纸面设计在给定边距和缩放规则下具有直接打印条件；若未拟合越界率较高而最终越界率为零，则表明原始反射拉回需要依赖A4归一化才能进入可实施区间，分类时不应将其等同于原始严格可行。

在图案相似性层面，本文采用重投影轮廓与目标轮廓之间的距离和区域重合度共同判别，避免只凭视觉轮廓是否相似作主观判断：

$e_H=d_H(T(\Gamma_p),\Gamma_m)$

$IoU=\frac{|A_p'\cap A_m|}{|A_p'\cup A_m|}$

其中 $T(\Gamma_p)$ 表示纸面曲线经正向反射后得到的镜面轮廓， e_H 越小表示最坏轮廓偏差越小， IoU 越大表示区域重合程度越高。分类流程按“参数筛选—反射拉回—A4统一拟合—正向重投影—容差分级”的顺序执行：先剔除圆柱尺寸和观察点不合理的候选参数，再对图3、图4分别生成纸面曲线；随后只对同一图案内的全部拉回点做统一比例缩放和平移，不对局部轮廓作单独扭曲；最后依据 $\rho_{out}$ 、 e_H 、 IoU 与中心线偏差进行分级。该流程使“严格可行”“条件近似可行”“一般不可行”对应明确的几何与误差来源，也为后续给出图3、图4的具体纸面设计和圆柱参数提供一致判定口径。

6.4 条件关系结果与判别证据

最终筛选得到的圆柱镜参数为半径 32.0 mm、高度 150.0 mm，轴线投影位于 A4 纸面坐标的 (105.0,148.0) mm，观察点取 (105.0,-210.0,220.0) mm，A4 拟合边距为 12 mm。按条件关系结果划分，图3与图4均归入“可打印拉回图案”类别：两者最终曲线越界率均为 0.0，分别对应 0.3103 和 0.3238 的纸面占用率；其关键阈值为纸面坐标须落入 $210\mathrm{ mm}\times297\mathrm{ mm}$ 范围、半径位于 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 核心结果表

parameter	value	meaning
A4_width_mm	210.0	纸面短边
A4_length_mm	297.0	纸面长边
R_mm	32.0	圆柱半径
H_mm	150.0	圆柱高度
axis_x_mm	105.0	轴线投影x坐标
axis_y_mm	148.0	轴线投影y坐标
observer_xyz_mm	(105.0, -210.0, 220.0)	观察点坐标
a4_fit_margin_mm	12.0	附件轮廓拉回后的A4内缩边距

基线对比主要体现在直接反射拉回与 A4 拟合拉回之间的差别。未进行 A4 等比例拟合时，图3的原始曲线越界率为 0.8000，图4为 0.6583，说明单纯按反射定律得到的 $P_{raw}$ 会有较大比例落在 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
R_mm	32.0	4.0	0.18
H_mm	150.0	15.0	0.06
axis_x_mm	105.0	8.0	0.11
axis_y_mm	148.0	8.0	0.14
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边距和后续手工修正保留了空间。

由相图和指标曲线共同得到的边界结论是：Q1的可实现条件首先要求反射拉回曲线经同一比例缩放和平移后完全落入A4有效区域，其次要求圆柱参数保持在可制作范围内，并使重投影轮廓误差不因边界压缩而显著放大。本轮示意实验中，半径32.0 mm、高度150.0 mm的圆柱镜与居中轴线配置能够同时容纳图3和图4，但未归一化越界率较高，说明该方案依赖A4拟合步骤来转化为可打印设计。其局限在于观察点仍为参数化设定，实物制作时需要用实际镜面半径、摆放位置和观看高度重新校准相图边界，否则指标曲线中的中心线偏差可能成为视觉误差的主要来源。

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
eta0	0.15	0.25	0.35	0.19
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-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
geometry_only	3	0.42
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
Fa_min	0.5	0.4	0.65	0.15
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-2 Q3 求解结果可视化

AtomGit开源社区

AtomGit 是由开放原子开源基金会联合 CSDN 等生态伙伴共同推出的新一代开源与人工智能协作平台。平台坚持“开放、中立、公益”的理念，把代码托管、模型共享、数据集托管、智能体开发体验和算力服务整合在一起，为开发者提供从开发、训练到部署的一站式体验。

更多推荐

深入理解C/C++混合编程

在一起的，比如，RTP视频传输，live555多媒体播放等都是C++下的，他需要调用JRTPLIB库，再比如，我那邮件发送，我也用C++写的，定义了一个Email对象，包含了。，上面代码的意思是：如果是C++文件（*.cpp）后缀，则使用extern “C”，在C++项目中应用的非常广泛。在工作中，C、C++密不可分，做我们嵌入式方面的，当然更多的是C，但，有时候却少不了C++，而且是C、C++混

AtomGit开源社区

【C++篇】map和set的使用

代码语言：javascriptAI代码解释set的声明如上，T就是set底层关键字(key)的类型。set默认要求T是支持比较大小的，如果不支持或者想按自己的比较方式走，可以传仿函数给第二个模板参数。set底层存储数据的内存是从空间配置器申请的，如果需要可以自己实现内存池，传给第三个参数一般情况下是不需要传后两个参数的set底层是用红黑树实现的，增删查的效率为O(logN),迭代器遍历走的是中序遍