shujuchuli.rar

%输入基本数据 y=[113.5 110.2 115.3 16.2 115.2 140.2 135.2 103.6 150.7 120.3 135.2 132.5 137.1 120.9 127.5 130.5 136.5 124.8 120.4 134.5]; %光强数据 x1=[19.5 24.7 30.7 29.8 19.1 25.6 31.4 27.9 22.1 25.5 31.1 30.4 18.7 19.7 14.6 29.5 27.7 30.2 22.7 25.2]; %色温数据 x2=[43 49 51 54 42 59 58 52 49 53 56 57 45 42 47 54 55 58 48 51]; %寿命数据 x3=[201 215 210 206 203 208 204 206 203 204 201 203 205 206 205 209 211 204 216 218 ]; %光通量数据 %创建所有测试数据矩阵 A(1,:)=x1; A(2,:)=x2; A(3,:)=x3; A1=mean(A); %样本平均 A2=nanmean(A); %算数平均 A3=geomean(A); %几何平均 A4=harmmean(A); %和谐平均 A5=trimmean(A,1); %调整平均 %绘制均值曲线 plot(A1,'bo') hold on plot(A2,'r-') plot(A3,'c--') plot(A4,'y:') plot(A5,'g-..') hold off title('均值曲线') xlabel('测试数据'),ylabel('白炽灯测量数据') legend('样本平均','算数平均','几何平均', '和谐平均','调整平均') %数据条形图 M(1:4,1)=0.05.*sum(y); M(1:4,2)=0.05.*sum(x1); M(1:4,3)=0.05.*sum(x2); M(1:4,4)=0.05.*sum(x3); bar(M) %样本方差的分析 miu=mean(A); sigma=var(A,1); %协方差计算 B=cov(A); %绘制对比图 subplot(1,2,1),plot(A) title('样本数据') subplot(1,2,2),plot(B) title('协方差结果') %相关系数计算 corrcoef(A) plot(B)