我一直在使用 MATLAB 对 UTM 数据(X 和 Y 坐标)执行核密度估计 (KDE)。我遇到了一个我似乎不理解的问题。
我用 45 个点的样本执行 KDE。一切正常,我制作了带有轮廓的图形。
[bandwidth,density,X,Y]=kde2d(data)
函数 kde2d 是 Zdravko Botev 的代码。我是从他在 MathWorks 上的文件交换中获得的。变量“数据”是我的数据的 45x2 数组。第一列包含 X 坐标,第二列包含 Y。
当我尝试在这 45 个点的子集上执行相同的代码行时,问题就出现了。我收到一个反复出现的错误:
Error using fzero (line 274)
The function values at the interval endpoints must differ in sign.
Error in kde2d (line 101)
t_star=fzero(@(t)(t-evolve(t)),[0,0.1]);
对于一堆不同的 45 点集上的一堆子集,我得到了同样的错误。
完整的集合有以下 45 个值:
1594436.281 572258.1272
1594418.48 572357.5859
1594471.362 572385.5186
1594516.726 572266.8206
1594415.313 572369.2754
1594519.701 572272.7153
1594415.377 572363.4139
1594468.365 572381.5779
1594518.139 572276.6059
1594425.496 572271.6874
1594524.259 572272.7651
1594502.555 572172.8749
1594516.747 572264.867
1594485.314 572360.2689
1594476.027 572375.7997
1594556.087 572419.6609
1594522.718 572274.7021
1594472.775 572395.3039
1594554.568 572419.6443
1594527.255 572276.7054
1594474.315 572393.3669
1594522.697 572276.6557
1594471.319 572389.4262
1594460.854 572373.6799
1594546.022 572228.0609
1594460.79 572379.5414
1594468.323 572385.4855
1594466.953 572371.7926
1594519.722 572270.7614
1594396.76 572398.3826
1594468.131 572403.0693
1594418.288 572375.1697
1594396.377 572433.5499
1594448.287 572271.9361
1594510.541 572276.523
1594424.466 572226.7345
1594413.773 572371.2124
1594511.848 572296.0774
1594513.367 572296.094
1594424.488 572224.7805
1594468.152 572401.1153
1594421.37 572371.2953
1594446.768 572271.9195
1594468.152 572401.1153
1594448.799 572225.0457
我尝试使用的子集之一是:
1594436.281 572258.1272
1594418.48 572357.5859
1594471.362 572385.5186
1594516.726 572266.8206
1594415.313 572369.2754
1594519.701 572272.7153
1594415.377 572363.4139
1594468.365 572381.5779
1594518.139 572276.6059
1594425.496 572271.6874
我不确定是否应该包含任何 Botev 的代码。我希望可以自行解释错误消息。如果没有,我可以提供更多。非常感谢。