Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- clear;
- clc;
- %DANE DO ZADANIA DLA STUDENT O NUMERZE n = 2
- f1 = 20;
- A = 17;
- phi = 0.2 * pi;
- T = 1/f1;
- %SYGNAŁ ANALOGOWY
- y = A * square(2*pi*f1*t+phi,gamma);
- figure;
- subplot(3,1,1);
- plot(t,y)
- title("y(t)")
- xlabel("t")
- ylabel("y")
- %SYGNAŁ DYSKRETNY y1
- N1 = 20;
- k1=N/20;
- t1=t(1:k1:N); %na podstawie funkcji alias.m
- y1=y(1:k1:N);
- subplot(3,1,2)
- stem(k1,y1)
- title("y1[k]")
- xlabel("k")
- ylabel("y1")
- %SYGNAŁ DYSKRETNY y2
- N2 = 100
- k2=N/100;
- t2=t(1:k2:N); %próbkowanie czasu 100 razy na okres
- y2=y(1:k2:N); %próbkowanie wartości y 100 razy na okres
- subplot(3,1,3)
- stem(k2,y2)
- title("y2[k]")
- xlabel("k")
- ylabel("y2")
- %WIDMO AMPLITUDOWE I FAZOWE y1
- n1 = length(y1)
- cm1 = (1/n1)*fft(y1)
- c_abs1 = abs(cm1);
- phi1 = angle(cm1);
- m1=1:1:n1;
- figure;
- subplot(2,1,1);
- stem(m1,c_abs1)
- title('Widmo amplitudowe dla y1')
- xlabel('m [-]')
- ylabel('c_m [-]')
- subplot(2,1,2);
- stem(m1,phi1)
- title('Widmo fazowe dla y1')
- xlabel('degree')
- ylabel('phi')
- %WIDMO AMPLITUDOWE I FAZOWE y2
- n2 = length(y2);
- cm2 = (1/n2)*fft(y2)
- c_abs2 = abs(cm2);
- phi2 = angle(cm2);
- m2=1:1:n2;
- figure;
- subplot(2,1,1);
- stem(m2,c_abs2)
- title('Widmo amplitudowe dla y2')
- xlabel('m [-]')
- ylabel('c_m [-]')
- subplot(2,1,2);
- stem(m2,phi2)
- title('Widmo fazowe dla y2')
- xlabel('degree')
- ylabel('phi')
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement