plotting time as legend values with filter as x and amplitude as y

Question

S 2024-4-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2104411-plotting-time-as-legend-values-with-filter-as-x-and-amplitude-as-y

编辑： S 2024-4-9

I am trying to make a plot where filter number is my x axis is my filter number and my y axis is amplitude. I am also trying to make a legend where the the t1 time values are what the lines of the plot are. I am not getting an error but I am getting a blank graph and am unsure what to do. Thank you for your time!

close all

clear all

clc

%Version with 3D, x axis filter

fs = 20e3;

numFilts = 32; %

filter_number = numFilts;

order = 4;

CenterFreqs = logspace(log10(50), log10(8000), numFilts);

signal_freq = 100; % Hz

Nperiods = 15; % we need more than 1 period of signal to reach the steady state output (look a the IR samples)

t = linspace(0,Nperiods/signal_freq,200*Nperiods); %

input = sin(2*pi*signal_freq*t) + 0*rand(size(t));

%Impulse

for ii = 1:filter_number

IR = gammatone(order, CenterFreqs(ii), fs);

[tmp,~] = freqz(IR,1,1024*2,fs);

% scale the IR amplitude so that the max modulus is 0 dB

a = 1/max(abs(tmp));

% % or if you prefer - 3dB

% g = 10^(-3 / 20); % 3 dB down from peak

% a = g/max(abs(tmp));

IR_array{ii} = IR*a; % scale IR and store in cell array afterwards

[h{ii},f] = freqz(IR_array{ii},1,1024*2,fs); % now store h{ii} after amplitude correction

plot(IR_array{ii})

end

figure

hold on

t1=[0.02 0.05 0.1]

for ii = 1:numel(filter_number)

output(ii,:) = filter(IR_array{ii},1,input);

plot(filter_number,output(ii,:))

LEGs{ii} = ['Time ' num2str(t1)]; %assign legend name to each

end

hold off

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function IR = gammatone(order, centerFrequency, fs)

% Design a gammatone filter

earQ = 9.26449;

minBW = 24.7;

% erb = (centerFrequency / earQ) + minBW;

erb = ((centerFrequency/earQ)^order + (minBW)^order)^(1/order);% we use the generalized

% function (depending on order)

B = 1.019 * erb;

f = centerFrequency;

tau = 1. / (2. .* pi .* B);

% gammatone filter IR

tmax = tau + 20/(2*pi*B);

t = (0:1/fs:tmax);

temp = (t - tau) > 0;

IR = ((t - tau).^(order - 1)) .* exp(-2*pi*B*(t - tau)).* cos(2*pi*f*(t - tau)) .* temp;

end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Hassaan 2024-4-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2104411-plotting-time-as-legend-values-with-filter-as-x-and-amplitude-as-y#answer_1438311

在 MATLAB Online 中打开

close all;

clear all;

clc;

% Initialization

fs = 20e3; % Sampling frequency

numFilts = 32; % Number of filters

CenterFreqs = logspace(log10(50), log10(8000), numFilts); % Center frequencies

signal_freq = 100; % Signal frequency in Hz

Nperiods = 15; % Number of periods

t = linspace(0, Nperiods/signal_freq, 200*Nperiods); % Time vector

input = sin(2*pi*signal_freq*t); % Input signal

% Impulse Response Calculation

IR_array = cell(1, numFilts);

h = cell(1, numFilts);

for ii = 1:numFilts

IR = gammatone(4, CenterFreqs(ii), fs);

[tmp, ~] = freqz(IR, 1, 1024*2, fs);

a = 1/max(abs(tmp));

IR_array{ii} = IR*a;

[h{ii}, ~] = freqz(IR_array{ii}, 1, 1024*2, fs);

end

% Filtering and Plotting

figure;

hold on;

t1 = [0.02, 0.05, 0.1]; % Example time values for legend

for ii = 1:numFilts

output = filter(IR_array{ii}, 1, input);

plot(abs(h{ii})); % Plotting the magnitude response

LEGs{ii} = ['Filter ' num2str(ii) ', Time ' num2str(t1(mod(ii-1, numel(t1)) + 1))]; % Unique legend for each line

end

legend(LEGs);

hold off;

xlabel('Frequency (Hz)');

ylabel('Amplitude');

title('Filter Responses with Time Legends');

% Gammatone Filter Function Remains the Same

function IR = gammatone(order, centerFrequency, fs)

% Design a gammatone filter

earQ = 9.26449;

minBW = 24.7;

% erb = (centerFrequency / earQ) + minBW;

erb = ((centerFrequency/earQ)^order + (minBW)^order)^(1/order);% we use the generalized

% function (depending on order)

B = 1.019 * erb;

f = centerFrequency;

tau = 1. / (2. .* pi .* B);

% gammatone filter IR

tmax = tau + 20/(2*pi*B);

t = (0:1/fs:tmax);

temp = (t - tau) > 0;

IR = ((t - tau).^(order - 1)) .* exp(-2*pi*B*(t - tau)).* cos(2*pi*f*(t - tau)) .* temp;

end

-----------------------------------------------------------------------------------------------------------------------------------------------------

If you find the solution helpful and it resolves your issue, it would be greatly appreciated if you could accept the answer. Also, leaving an upvote and a comment are also wonderful ways to provide feedback.

It's important to note that the advice and code are based on limited information and meant for educational purposes. Users should verify and adapt the code to their specific needs, ensuring compatibility and adherence to ethical standards.

Professional Interests

Technical Services and Consulting
Embedded Systems | Firmware Developement | Simulations
Electrical and Electronics Engineering

Feel free to contact me.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Hassaan 2024-4-8

@S You can adjust as per your needs.

-----------------------------------------------------------------------------------------------------------------------------------------------------

If you find the solution helpful and it resolves your issue, it would be greatly appreciated if you could accept the answer. Also, leaving an upvote and a comment are also wonderful ways to provide feedback.

It's important to note that the advice and code are based on limited information and meant for educational purposes. Users should verify and adapt the code to their specific needs, ensuring compatibility and adherence to ethical standards.

Professional Interests

Technical Services and Consulting
Embedded Systems | Firmware Developement | Simulations
Electrical and Electronics Engineering

Feel free to contact me.

S 2024-4-8

编辑：S 2024-4-9

@Hassaan This is not showing what filter 1,2,3,... is at all 3 times? I am trying to make filter number on the x axis and amplitude y axis while there are 3 lines on the graph 1 representing each time in t1

请先登录，再进行评论。