Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- clc
- clear all;
- T = 2; %okres
- U = 1; %amplituda
- N = 100; %liczba próbek
- f= 1/T; %częstotliwość
- w = 2* pi*f; %pulsacja
- t = -T/2:T/N:T/2-T/N+2*T; %przedział czasu
- for k=1:N
- u(k) = abs((2*U*t(k))/T); %funkcja u1(t)
- end
- u1 = [u, u, u];
- subplot(2,2,1); %rozmieszczenie wykresu
- plot(t, u1); %stworzenie wykresu funkcji
- title("Sygnał u1 w postaci funkcji czasu"); %tytuł wykresu
- xlabel("Czas [s]"); %podpis osi X
- ylabel("Amplituda [V]"); %podpis osi Y
- grid on; %włączenie siatki
- a=fft(u1);
- m=abs(a);
- n=0:N/2;
- f_a=n*f/(3*N/2);
- subplot(2,2,2);
- stem(f_a(1:N/2),m(1:N/2),'.');
- u2 = U/2-4*U/pi*pi*(cos(w*t)+cos(3*w*t)/3^2+cos(5*w*t)/5^2+cos(7*w*t)/7^2+cos(9*w*t)/9^2+ cos(11*w*t)/11^2+...
- cos(13*w*t)/13^2+cos(15*w*t)/15^2+cos(17*w*t)/17^2+cos(19*w*t)/19^2);
- subplot(2,2,3); %rozmieszczenie wykresu
- plot(t, u2);
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement