Phase angle
Создайте комплексное число и вычислите его величину и фазу.
z = 2*exp(i*0.5)
z = 1.7552 + 0.9589i
r = abs(z)
r = 2
theta = angle(z)
theta = 0.5000
Создайте сигнал, который состоит из двух синусоидов частот 15 Гц и 40 Гц. Первая синусоида имеет фазу , и второй имеет фазу . Дискретизируйте сигнал со скоростью 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 + <reservedrangesplaceholder2> <reservedrangesplaceholder1> и использует atan2 функция для вычисления угла между положительным x - осью и лучом от источника до точки (x, y) в xy -плане.
