Фазовый угол
Создайте комплексное число и вычислите его величину и фазу.
z = 2*exp(i*0.5)
z = 1.7552 + 0.9589i
r = abs(z)
r = 2
theta = angle(z)
theta = 0.5000
Создайте сигнал, состоящий из двух синусоид частот 15 Гц и 40 Гц. Первая синусоида имеет фазу , а вторая - δ/2. Выполните выборку сигнала с частотой 100 Гц в течение одной секунды.
fs = 100; t = 0:1/fs:1-1/fs; x = cos(2*pi*15*t - pi/4) - sin(2*pi*40*t);
Вычислите преобразование Фурье сигнала. Постройте график величины преобразования как функции частоты.
y = fft(x); z = fftshift(y); ly = length(y); f = (-ly/2:ly/2-1)/ly*fs; stem(f,abs(z)) xlabel 'Frequency (Hz)' ylabel '|y|' grid
Вычислите фазу преобразования, удалив значения преобразования малой величины. Постройте график фазы как функции частоты.
tol = 1e-6; z(abs(z) < tol) = 0; theta = angle(z); stem(f,theta/pi) xlabel 'Frequency (Hz)' ylabel 'Phase / \pi' grid
angle принимает комплексное число z = x + iy и использует atan2 функция для вычисления угла между положительной осью x и лучом от начала координат до точки (x, y) в плоскости xy.