منبع : زلزله

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

در ادامه مطلب

برنامه رسم تابع چگالی احتمال بزبان متلب (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;

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  که تعداد ستونها نیز مهم نمیباشد.

نمونه مشابه

برنامه ای برای ترسیم تابع چگالی طیفی و تابع خودهمبستگی زلزله تحت matla