Добавил:
Lybov_Vernis
2202 2050 2250 3772 Сб
Песня посвящается героическим защитникам курсовой по ЦСП в апреле 2025 года
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:261 / pfc2612
.m % Define constants
b0 = 1;
b1 = 0;
b2 = -0.84;
T = 1/8000;
% Create frequency vector
f = 0:0.01:8;
% Frequency from 0 to 8 KHz
% Calculate Ф(f)
Phi = atan2((-b1 .* sin(2 * pi * f * T * 1000) - b2 .* sin(4 * pi * f * T * 1000)), (b0 + b1 .* cos(2 * pi * f * T * 1000) + b2 .* cos(4 * pi * f * T * 1000)));
% Plotting
figure; % Create a new figure
plot(f, Phi, 'color', [0 0 0], 'LineWidth',2); % Plot Ф(f) against f
set(gca,'FontName','Helvetica','FontSize',11)
xlabel('Frequency [kHz]', 'color', [0 0 0]); % Label x-axis
ylabel('Phase [rad]', 'color', [0 0 0]); % Label y-axis
title(''); % Title of the plot
grid on; % Add grid for better visualization
b0 = 1;
b1 = 0;
b2 = -0.84;
T = 1/8000;
% Create frequency vector
f = 0:0.01:8;
% Frequency from 0 to 8 KHz
% Calculate Ф(f)
Phi = atan2((-b1 .* sin(2 * pi * f * T * 1000) - b2 .* sin(4 * pi * f * T * 1000)), (b0 + b1 .* cos(2 * pi * f * T * 1000) + b2 .* cos(4 * pi * f * T * 1000)));
% Plotting
figure; % Create a new figure
plot(f, Phi, 'color', [0 0 0], 'LineWidth',2); % Plot Ф(f) against f
set(gca,'FontName','Helvetica','FontSize',11)
xlabel('Frequency [kHz]', 'color', [0 0 0]); % Label x-axis
ylabel('Phase [rad]', 'color', [0 0 0]); % Label y-axis
title(''); % Title of the plot
grid on; % Add grid for better visualization