test.m

电磁场实验中

close all;
clear;
load outside.txt;    %调入外部数据文件
f=outside(:,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);          %统计各数据值出现的概率
end
c=sum(F);                                 %统计实验累计概率
for  j=1:length(F)
    d(j)=sum(F(1:j))/c;   
end

figure(1)
    set(1,'Position',[10,100,500,400])     %设定窗口位置及大小                            
figure(2)
    set(2,'Position',[550,100,500,400])    %设定窗口位置及大小 
    
figure(1)
t=linspace(a,b,length(F));
plot(t,d,'*')                                 %画实验概率累积分布曲线
hold on
plot(xx,h,'r')                                %画理论概论累积分布曲线
xlabel('电平值/（ -dbmw）')                %画理论概率密度函数
ylabel('概率累积值')
title('概率累积分布图')

figure(2)
hist(f,(b-a)*10);                          %画频数直方图
hold on
plot(xx,p*100,'r');
xlabel('电平值/（ -dbmw）')                %画理论概率密度函数
ylabel('概率密度值')
title('概率密度分布图')