Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- clear all;
- close all;
- N = 32
- load("h.mat")
- h = h'
- freq1 = 2
- freq2 = 6
- u3 = zeros(32,1);
- for i = 1:1:32
- u3(i) = (sin(2*pi*freq1*i/32) + sin(2*pi*freq2*i/32));
- end
- u3
- U3 = toeplitz(u3)
- y3 = U3 * h
- %Sposob wykladniczy u3
- F1 = zeros(N);
- for n = 1:1:32
- for m = 1:1:32
- F1(n,m) = (1/sqrt(N))*exp(-j*((2*pi)/N)*n*m);
- end
- end
- F1
- X1 = F1 * u3
- new_X1 = 20*log10(abs(X1))
- figure(1)
- plot(new_X1, 'g');
- title("Sygnal w dziedzinie czestotliwosci - exponenta")
- x1_1 = F1' * X1
- x2_1 = inv(F1) * X1
- %Sposob trygonometryczny u3
- F2 = zeros(N);
- for n = 1:1:32
- for m = 1:1:32
- F2(n,m) = (1/sqrt(N))*cos((2*pi/N)*n*m) -j*(1/sqrt(N))*sin((2*pi/N)*n*m);
- end
- end
- F2
- X2 = F2 * u3
- new_X2 = 20*log10(abs(X2))
- figure(2)
- plot(new_X2, 'r');
- title("Sygnal w dziedzinie czestotliwosci - trygonometria")
- x1_2 = F2' * X2
- x2_2 = inv(F2) * X2
- %Analiza sygnalu y
- F1y = zeros(N);
- for n = 1:1:63
- for m = 1:1:63
- F1y(n,m) = (1/sqrt(N))*exp(-j*((2*pi)/N)*n*m);
- end
- end
- F1y
- y = conv(h, u3)
- Y = F1y * y
- figure(3)
- plot(Y)
- title("Analiza dla sygnalu Y")
- %Random signal
- signal = randn(32,1)
- Signal = F1 * signal
- signal1 = F1' * Signal
- signal2 = inv(F1) * Signal
- figure(4)
- plot(Signal)
- title("Random signal")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement