Need explanation on plotting frequency response of a causal discrete-time LTI system implemented using the difference equation.

5 次查看(过去 30 天)
So....
i have a program like this
RH_coeff=[0.1 -0.1176 0.1];% Coefficient of x on RHS of difference equation
LH_coeff=[1 -1.7119 0.81]; % Coefficient of y on LHS of difference equation
%% Frequency response of the filter
[H,w]=freqz(RH_coeff,LH_coeff); % Frequency response
plot(w/(2*pi),20*log10(abs(H)),'LineWidth',1.5); % Plot Frequency response
xlabel('Normalized frequency (f/fs)Hz');
ylabel('Magnitude');
title('Frequency response');
grid;
%% Generating response y[n] due to input x[n]
n=0:100;% X axis
xn=cos(0.1*n*pi);% x[n]
yn=filter(RH_coeff,LH_coeff,xn); %y[n]
figure;
subplot(2,1,1);
stem(n,xn); %Plot input
xlabel('n');
ylabel('x(n)');
title('Input sequence');
subplot(2,1,2);
stem(n,yn);% Plot output
xlabel('n');
ylabel('y(n)');
title('Output sequence');
grid;
  1. Can someone explain to me what [H,w]=freqz(RH_coeff,LH_coeff); means? does it mean H=RH_coeff, w=LH_coeff?
  2. what does this function (plot(w/(2*pi),20*log10(abs(H)),'LineWidth',1.5);) means? why does this function exist? i thought for plotting frequency responce i just stopped at freqz(RH_coeff,LH_coeff);

采纳的回答

Star Strider
Star Strider 2022-5-16
does it mean H=RH_coeff, w=LH_coeff?
No. The first argument are the numerator coefficients of the filter transfer function, and the second argument are the denominator coefficients.
what does this function (plot(w/(2*pi),20*log10(abs(H)),'LineWidth',1.5);) means?
It plots the magnitude section of the Bode plot created by freqz as the frequency in Hz and magnitude in dB.
.
  4 个评论
Evelio Excellenta
Evelio Excellenta 2022-5-18
thankyou~
ok now i know that [H,w] is the output of freqz.
now my question is does [H,w] has name? because im preparing for presentation and i don't know what to call it.
should i just call it by the alphabet H and w? or H means frequency and w is magnitude? or just call it a variable?

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Digital Filter Analysis 的更多信息

产品


版本

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by