3582481275

MIS Ćwiczenie 2: Statystyczna analiza danych - część 1 autor: Jan Kowalski (kowal@.), AGH - EAlilB - MIS III rok. Data modyfikacji: 01.01.2015

Contents

» Generacja sygnałów losowych

■    Obliczanie parametrów sygnału (momentów)

» Wyznaczenie histogramu

» Wyznaczenie estymaty fgp - na podstawie histogramu

•    Porównanie estymaty fgp z rozkładem teoretycznym

•    Naniesienie na rysunki wartości oblicznonych parametrów (średnia, mediana i moda)

■    Obliczenie dystrybuanty jako całki z fgp

elear all; close all; clc

Generacja sygnałów losowych

% parametry próbkowania fpr=1000; dt=l/fpr; t0=1.5;

N=t0/dt;

t=(0:N-l)*dt;

% częstotliwość próbkowania (dobrana tak by było spełnione twierdzenie o próbkowaniu) % czas obserwacji sygnałów % liczba próbek % wektor czasu

df=fpr/N; f=df*(0:N-l); Nwektor częstotliwości (potrzebny przy analizie fft)

% typ=l - rozkład normalny, typ=2 - rozkład Chi-kwadrat typ=2;

świtch(typ) case(l)

% u - wartość oczekiwana, s - odchylenie standardowe % generujemy sygnał losowy o rozkładzie normalnym

% L - liczba stopni swobody

% generujemy sygnał losowy o rozkładzie chi-kwadrat

u=0; s=l; x=s*randn(l,N)+u; case (2)

L=4;

x=chi2rnd(L,l,N); otherwise return;

end

X=abs(fft(x-mean(x)));

figurel = figure('PaperPosition',[0.6345 6.345 20.3 15.23],’PaperSize',[20.98 29.68]); axesl = axes('FontName','MS Sans Serif','FontSize',15,'Parent',figurel) subplot(211); plot(t,x);

grid on; xlabel('Czas [s]');ylabel('Amplituda');title(['Sygnał stochastyczny']);

subplot(212); stem(f(1,(1:N/2)),X(1,(1:N/2)));legend(['FFT']); xlabel('Częstotliwość [Hz]'); ylabel('Amplituda'); grid on;

Sygnał stochastyczny

20,-,-1-_r

co ¹5

0    0.5    1    1.5

Czas [s]