139

Analiza harmonicznych 139

gdzie:

t    - zbiór posortowanych rosnąco wartości współrzędnych czasu od 0 do

9nax»

x    - zbiór wartości pomiarów odpowiadających chwilom t,

tp = [0:7':f_nlax] - zbiór równoodległych o T chwil czasowych, dla których mają być wykonane obliczenia interpolacyjne,

‘metoda_interpolacji’ - parametr tekstowy mogący przybierać wartości: ‘linear’ - interpolacja liniowa, ‘cubic’ - interpolacja wielomianami 3. stopnia, ‘spline’ - interpolacja funkcjami sklejonymi 3. stopnia.

Należy zwrócić uwagę, że metoda ‘cubic’ wymaga, aby odległości między elementami wektora t były stałe, czyli t = [0:d/:r_max].

Zastosowanie funkcji interpl pokazano na przykładzie interpolacji pomiarów prądu. Funkcja ipomiar zawiera pomiary prądu w jednostkach względnych zarejestrowane w nieregularnych odstępach czasu:

function [t,i] = ipomiar

i funkcja zawiera pomiary chwilowych wartości prądu pomiary= [

% t, s	i(t), pu	t, s	i(t), pu	t, s	i(t), pu
0.0000	-0.09311963 ;	0.0002	-0.09772925	0.0003	-0.10219942;
0.0005	-0.10652674;	0.0006	-0.11070815	0.0014	-0.12934544;
0.0022	-0.14409672 ;	0.0030	-0.15501359	0.0038	-0.16238718 ;
0.0057	-0.16872978;	0.0076	-0.16954201	0.0095	-0.17680811;
0.0114	-0.20024868;	0.0135	-0.25102910	0.0156	-0.32221004;
0.0177	-0.40181197;	0.0197	-0.47459991	0.0221	-0.53105173
0.0245	-0.55561562;	0.0269	-0.55937869	0.0293	-0.56411882;
0.0319	-0.59811287;	0.0344	-0.66748733	0.0369	-0.76100443;
0.0395	-0.85367743;	0.0409	-0.89502590	0.0424	-0.92355900;
0.0439	-0.93917719;	0.0454	-0.94430765	0.0469	-0.94334939;
0.0483	-0.94173459;	0.0498	-0.94418700	0.0513	-0.95389955;
0.0528	-0.97297414;	0.0544	-0.99918131	0.0559	-1.02831530;
0.0574	-1 . 05462879;	0.0588	-1.07026858	0.0601	-1.07497238;
0.0615	-1.06647784;	0.0628	-1.04429568	0.0644	-1 . 00326949;
0.0659	-0.95093290;	0.0675	-0.89397349	0.0691	-0.83944675;
0.0705	-0.79716667 ;	0.0719	-0.76550762	0.0734	-0.74531538;
0.0748	-0.73517102;	0.0765	-0.73158508	0.0782	-0.72939476;
0.0799	-0.72128558;	0.0817	-0.70258445	0.0832	-0.67562715;
0.0847	-0.64193881 ;	0.0863	-0.60640322	0.0878	-0.57517305;
0.0892	-0.55610931;	0.0906	-0 . 54987840	0.0920	-0.55849444]

t=pomiary(:,1) ; i =pomiary(:,2) ; return

Należy przeprowadzić interpolację pomiarów prądu tak, aby wartości prądu były równoległe od siebie o przedział czasowy T = 2IN = 2/500 = 0.004 s. Zadanie to realizuje funkcja interpi. Wynik interpolacji prądu pokazano na rysunku 11.1.