空气中PM2.5问题的研究（matlab代码）

%%%%%%%%%空气中PM2.5问题的研究%%%%%%%%%%%%
%%%%%%%%缺失值检测%%%%
clc, clear
a=readmatrix('F:\数学建模\培训模拟\研究生组题目六2013-D题：空气中PM2.5问题的研究\附件1. 数据1.(武汉市一个监测点数据：2013.01.01-2013.08.26).xls','Range','B2:G239');
A = ismissing(a) %显示1为缺失值
[m,n]=size(A)
sumj=sum(A)
j=sumj./m %缺失值所占的比列
%发现可吸入颗粒物有两个缺失值，在表中的D120和D159
% scatter(1:size(a,1),a(:,3))
a([119,158],3)=0
meanValue=sum(a(:,3))/(size(a,1)-2)
a([119,158],3)=79

%%%%%%%%异常值检测%%%%
boxplot(a,'Labels',{'二氧化硫','二氧化氮','颗粒物','一氧化碳','臭氧','PM2.5'})
xlabel("AQI监测指标")
ylabel("AQI指数")
title("AQI监测指标异常值检测")

%%%%%%%%第一题第一问%%%%
%%%%正态性检验%%%%
pd=fitdist(a(:,1),'Normal'),qqplot(a(:,1),pd) %直接画QQ图，检验正态性
%%%%总体相关性分析%%%%
cor=corr(a,'Type','Spearman')
figure(1)
X = categorical({'SO2','NO2','PM10','CO','O3','PM2.5'});
X = reordercats(X,{'SO2','NO2','PM10','CO','O3','PM2.5'});
bar(X,cor(:,1))
title("SO2与AQI监测指标的相关性")
grid
figure(2)
X = categorical({'SO2','NO2','PM10','CO','O3','PM2.5'});
X = reordercats(X,{'SO2','NO2','PM10','CO','O3','PM2.5'});
bar(X,cor(:,2))
title("NO2与AQI监测指标的相关性")
grid
figure(3)
X = categorical({'SO2','NO2','PM10','CO','O3','PM2.5'});
X = reordercats(X,{'SO2','NO2','PM10','CO','O3','PM2.5'});
bar(X,cor(:,3))
title("PM10与AQI监测指标的相关性")
grid
figure(4)
X = categorical({'SO2','NO2','PM10','CO','O3','PM2.5'});
X = reordercats(X,{'SO2','NO2','PM10','CO','O3','PM2.5'});
bar(X,cor(:,4))
title("CO与AQI监测指标的相关性")
grid
figure(5)
X = categorical({'SO2','NO2','PM10','CO','O3','PM2.5'});
X = reordercats(X,{'SO2','NO2','PM10','CO','O3','PM2.5'});
bar(X,cor(:,5))
title("O3与AQI监测指标的相关性")
grid
figure(6)
X = categorical({'SO2','NO2','PM10','CO','O3',