光时域检测仪（OTDR）

大学数学建模竞赛校赛暨国赛选拔 B 题 光时域检测仪（OTDR）是检测光缆完整性的重要工具。在测试光缆的过程中，仪器 从光缆的一端注入较高功率的激光或光脉冲，并通过同一侧接收反射信号。当光脉冲通过 光缆传输时，部分散射及反射将返回发射端。光时域反射仪(OTDR)只会测量强度较高的反 射回来的光信号，光时域反射仪(OTDR)根据光的后向散射与菲涅耳反向原理制作，利用光 在光纤中传播时产生的后向散射光来获取衰减的信息，从而间接地测量光缆损耗与故障位 置。在使用 OTDR 进行故障位置检测时，通过 OTDR 反射回来的信号，我们需要检测和识 别有用信号在整个观测到的信号序号中的相对位置，从而确定信号发生故障的位置。信号 是由伪随机二进制序列（PRBS）算法生成的一段伪随机数（PN）序列。我们观测到的信 号全程伴有高斯噪声。值得注意的是，真实信号小于观测到的带噪信号序列长度的一半， 并且它只存在于整个序列的一部分之中。由于无法精确控制信号发射的确切起始点，我们 通常会固定信号截取的起始位置和长度。因此，在不同的相对起始位置发射信号，就形成 了我们所观测到的各种带有不同噪声特性的信号序列。信噪比是衡量信号质量的重要参 数，表示有用信号的功率与背景噪声功率的比值，高信噪比表示信号强，噪声弱，在后面 的题目中，我们固定信噪比 snr=10。 信号的起始位置表示在带噪信号序列中第一次出现信号的位置编号，如长度为 10 的带 噪信号中出现长度为 4 的信号，若信号起始位置为 5，则认为 5-8 位置为信号，其余位置 为纯噪声。 问题一：现有 10 条带噪声的信号序列，长度为 1024，已知这 10 条序列中包含的信号 都是同一列伪随机序列（PN），长度为 511，你能否找到该信号在这 10 条带噪信号序列中 的起始位置？（附件 1） 问题二：由于不知道信号的伪随机序列（PN），我们现在将信号随机的在含有高斯噪 声的光纤中传递 10000 次，每次传递信号的起始位置不同，观测到 10000 条带噪信号序 列，带噪信号序列长度为 1024，其中前 8000 条序列通过专业手段测出了信号的起始位 置，你能否识别剩下带噪信号序列中信号的起始位置？（附件 2） 问题三：在问题二基础上，已知信号长 511，你能否推断出信号的伪随机序列？（伪 随机序列的每个位置只有 1 和-1 两种信号） 问题四：在实际的故障检测中，我们观测到的带噪信号序列长度远大于前面给出的序 列长度，导致在求信号起始位置时耗时太长，实际应用中难以接受，附件三中给出了 10 条 长度为 2 20 的带噪信号序列和它们所包含信号的伪随机码，你有什么方法能在较短的时间 内计算该信号在这 10 条带噪信号序列中的起始位置吗？

源代码：代码