读取附件sd.xlsx,计算数据如下统计量:
资源链接:https://download.csdn.net/download/qq_21484461/85801116?spm=1001.2014.3001.5503
>> [num1]=xlsread('F:\sd.xlsx',1,'A1:B191') %读取前1、2列
>>
>> [num2]=xlsread('F:\sd.xlsx',1,'D1:E191') %读取前4、5列
>>
>> [num3]=xlsread('F:\sd.xlsx',1,'G1:H40') %读取前7、8列
>> mean(num1(:,[1]))
ans =
2.9634
>> var(num1(:,[1]))
ans =
0.9337
>> std(num1(:,[1]))
ans =
0.9663
>> max(num1(:,[1]))
ans =
6.4000
>> min(num1(:,[1]))
ans =
0.5000
>> max(num1(:,[1]))-min(num1(:,[1]))
ans =
5.9000
>> median(num1(:,[1]))
ans =
2.8000
>> prctile(num1(:,[1]),40) %40分位数
ans =
2.5000
>> mode(num1(:,[1]))
ans =
2.2000
>> std(num1(:,[1]))/mean(num1(:,[1]))
ans =
0.3261
>> sum(num1(:,[1]).^2)/length(num1(:,[1])) %2阶原点矩
ans =
9.7103
或
>> mean(num1(:,[1]).^2)
ans =
9.7103
>> for i = 1:length(num1(:,[1]))
x1(i) = num1(i) - mean(num1(:,[1]))
end
>>
>> mean(x1.^2)
ans =
0.9288
>> skewness(num(:,[1]))
ans =
0.3339
>> kurtosis(num(:,[1]))
ans =
2.9240
>> cov(num1(:,[1]))
ans =
0.9337
>> cov(num1(:,[2]))
ans =
5.3527
>> cov(num1(:,[1]),num1(:,[2]))
ans =
0.9337 0.3703
0.3703 5.3527
>> corr(num1(:,[1]),num1(:,[2]))
ans =
0.1656