برنامه رسم تابع چگالی احتمال بزبان متلب (matlab)

clc;

F=load('path file');

N=0;

[R C]=size(F);

for i=1:R

for j=1:C

Acc(1,(i-1)*C+j)=F(i,j);

N=N+1;

end

end

dt=0.02;

T=(N-1)*dt;

A(1)=min(Acc);

Amax=max(Acc);

k=1;n=0;dx=0.0001;

sumpr(1)=0;

while A(k)<=Amax

for i=1:N

if Acc(1,i)>=A(k) & Acc(1,i)<(A(k)+dx)

n=n+1;

end

end

proden(k)=(n*dt)/T;

sumpr(k+1)=sumpr(k)+proden(k);

A(k+1)=A(k)+dx;

k=k+1;

n=0;

end

A(k)=[];sumpr(k)=[];

subplot(1,2,1);

plot(A,sumpr);

title('Distribution function');

subplot(1,2,2);

plot(A,proden);

title('Probability density function');

لازم بذکر است که فایل ورودی از نوع شتاب زلزله با پسوند txt.* میباشد شبیه فایل دانلود شده از سایت peer  که تعداد ستونها نیز مهم نمیباشد.

برگرفته از وبلاگ مهندس فراهانی