MATLAB编程技巧与数据分析实践

55、你的朋友弗雷德试图在频域中对核进行归一化。如果gausX是一维高斯函数的傅里叶频谱，下面的代码是否可行？如果gausX是二维高斯函数的傅里叶频谱呢？代码为：gausX = gausX ./ max(gausX);

对于一维高斯函数的傅里叶频谱，该代码可以通过将频谱中每个值除以最大值来实现最大值归一化，是可行的。

对于二维高斯函数的傅里叶频谱：

• max(gausX) 默认按列计算最大值，可能无法得到预期的归一化结果。
• 需要使用 max(max(gausX)) 来获取二维矩阵中的全局最大值，以实现正确的归一化。

56、下面的代码行有两个错误（变量cr是二维卷积的结果），找出并修复它们。cr = cr(halfK(1)+1:end - halfK(2),halfK(1)+1:end - halfK(2));

以下是给定文本内容转换为 Markdown 格式 的结果：

---

错误一：

在二维卷积结果的索引中，两次使用的行和列索引范围相同，可能不符合预期，应分别使用不同的行和列索引范围，即分别使用 halfK(1) 和 halfK(2) 来确定行和列的裁剪范围。

错误二：

原代码中裁剪行和列时，列裁剪使用了 halfK(1) ，应改为 halfK(2) 。

修复后的代码为：

cr = cr(halfK(1)+1:end - halfK(1), halfK(2)+1:end - halfK(2));

57、使用intFact参数的不同值进行绘图，参数值范围从5到200，以16个线性间隔取值。将结果显示在一个4×4的子图中。较高的intFact值计算时间更长，但会产生更平滑的图。你认为该参数的合适范围是多少？1000比200更好吗？

由于没有相关内容，无法确定 intFact 参数的合适范围，也无法判断 1000 是否比 200 更好。一般来说，若追求更平滑的图且计算资源和时间允许， 1000 可能更好；若考虑计算效率， 200 可能更合适。

58、在MATLAB中，判断代码 [gridX,gridY] = ndgrid(interpX,interpY); 和 [gridX,gridY] = meshgrid(interpX,interpY); 是否会给出相同的输出。若不同，调整代码使两组输出相同。

在 MATLAB 中测试这两行代码的输出是否相同：

meshgrid 仅支持 2D 和 3D 网格，且会交换输出中的列和行；
ndgrid 适用于更一般的情况。

若输出不同，可根据需求调整代码。若用于更通用场景，可考虑统一使用 ndgrid 函数。

59、如果数据大致呈正态分布，你可以使用标准的参数 t 检验进行统计评估。对两组模拟放电率数据进行双样本 t 检验。如果你不知道 MATLAB 中哪个函数可以进行双样本 t 检验，你需要去搜索一下。

在 MATLAB 中，可以使用 'ttest2' 函数进行双样本 t 检验。

假设两组模拟放电率数据分别存储在变量 'data1' 和 'data2' 中，使用以下代码进行双样本 t 检验：

[h, p, ci, stats] = ttest2(data1, data2);

其中：

• 'h' 是假设检验的结果（ 1 表示拒绝原假设， 0 表示不拒绝原假设）
• 'p' 是 p 值
• 'ci' 是均值差的置信区间
• 'stats' 包含其他统计信息

60、在模拟的放电率数据和脑电图时频数据上进行元置换检验，在20次元置换运行中，阈值的变化大吗？

进行元置换检验，需将置换测试代码放在循环中运行多次（如20次），获取多个z值，然后将这些z值平均得到最终结果。可分别在模拟的放电率数据和脑电图时频数据上运行此循环。

要确定20次元置换运行中阈值的变化情况，需记录每次运行得到的阈值，然后分析这些阈值数据，如计算方差或极差等指标来判断变化大小。

61、在下面的代码行中添加两个字符，以绘制一条穿过 y 轴 50% 位置的垂直线。代码为 plot([5 5],get(gca,’ylim’))

可添加的两个字符是 *0.5 ，修改后的代码为 plot([5 5],get(gca,'ylim')*0.5) ，这样就能绘制穿过 y 轴 50% 位置的垂直线。

62、以下代码行将产生错误。在MATLAB中运行代码之前，找出并修复错误。然后在MATLAB中测试原始代码和你修改后的代码。注意错误消息的文本，并确保你的解决方案有效。r = randn(1000,1); r(r>0