我一直在使用MATLAB对UTM数据(X和Y坐标)执行核密度估计(KDE)。我遇到了一个我似乎不理解的问题。
我用45分的样本执行KDE。一切正常,我用轮廓生成图表。
[bandwidth,density,X,Y]=kde2d(data)
函数kde2d是Zdravko Botev的代码。我是从他在MathWorks上的文件交换中获得的。变量'data'是我的数据的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的任何代码。我希望错误消息可以自己解释。如果没有,我可以提供更多。非常感谢你。