有同学指出这个代码在每一年取值范围不同的情况下图形有误，所以提供了一种解决方法，可以自行判断是否正确，我在公众号中将这些代码均放出来，大家可以去我公众号获取：阿奇宏观与计量课堂。同学修改的代码如下：

clc

clear

%--------------------------------------------

cd('E:\计量指导\核密度分析')  %放置数据的文件夹位置

X=xlsread('核密度分析数据西部实体.xlsx','sheet1'); %数据名称和后缀

【m,n】=size(X);

min=min(min(X));

max=max(max(X));

pts=linspace(min,max,100);

%--------------------------------------------

for i=1:10 %年份

  【f,xi】=ksdensity(X(:,i),pts);

  F(:,i)=f;

end

F=F';

%--------------------------------------------

y=2011:2020; %更改为具体的时间  

【pts,y】=meshgrid(x,y);

%--------------------------------------------

mesh(x,y,f)

ylabel('年份');xlabel('年夜间灯光亮度'); %做出相应更改，其他不变