clear;load floor1.txt;    %調入外部數據文件f=floor1(:,1);      %取圖像信號數據faverage=mean (f);  %求均值fstd=std(f);        %求標準差a=min(f);           %最小值b=max(f);           %最大值xx=(a-2):0.05:b+2;  p=normpdf(xx,faverage,fstd);              %理論正態分布函數h=normcdf(xx,faverage,fstd);              %理論累計概率count=0;for i=a:0.1:b  count=count+1;  aa=find(f==i);   F(count)=length(aa)/length(f);          %統計各數據值出現的概率endc=sum(F);                                 %統計實驗累計概率for  j=1:length(F)    d(j)=sum(F(1:j))/c;   endfigure(1)    set(1,'Position',[10,50,300,200])     %設定窗口位置及大小                            figure(2)    set(2,'Position',[350,50,300,200])    %設定窗口位置及大小     figure(1)t=linspace(a,b,length(F));plot(t,d)                                 %畫實驗概率累積分布曲線hold onplot(xx,h，'r')                                %畫理論概論累積分布曲線xlabel('電平值/（ -dbmw）')                %畫理論概率密度函數ylabel('概率累積值')title('概率累積分布圖')figure(2)hist(f,(b-a)*10);                          %畫頻數直方圖hold onplot(xx,p*20,'r');xlabel('電平值/（ -dbmw）')                %畫理論概率密度函數ylabel('概率密度值')title('概率密度分布圖')