how to change this code to plot the frequency response of the desired interval?

Question

Michael Sugiarto 2022-3-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1678254-how-to-change-this-code-to-plot-the-frequency-response-of-the-desired-interval

评论： Star Strider 2022-3-23

Hi! I am using the R2021b Matlab and I need this code to plot the frequency response of the audio in the specified time interval. Could anyone help me please?

The code is displayed below:

recObj = audiorecorder;
Fs=8000;
filename = sprintf('myAudioData.wav');
disp('Start speaking.')
recordblocking(recObj, 10);
disp('End of Recording.');
doubleArray = getaudiodata(recObj);
audiowrite(filename,doubleArray,Fs);
[x,Fs] = audioread('myAudioData.wav');
ts = 1 / Fs;
n = round(0.02 / ts); 
nexttile       
t1 = 1.000; 
t2 = 1.020; 
% start and end index in audio data of segment of interest
idx1 = round(t1/ts);
idx2 = round(t2/ts);
% plot it!
plot(x(idx1:idx2));       % plot 1st 20 ms of audio file
nexttile
t3 = 3.000; 
t4 = 3.020; 
% start and end index in audio data of segment of interest
idx3 = round(t3/ts);
idx4 = round(t4/ts);
% plot it!
plot(x(idx3:idx4));

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2022-3-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1678254-how-to-change-this-code-to-plot-the-frequency-response-of-the-desired-interval#answer_924664

That depends on what you want as the result.

Two possibilities are fft and pspectrum with the 'spectrogram' option.

7 个评论
显示 5更早的评论隐藏 5更早的评论

Star Strider 2022-3-23

在 MATLAB Online 中打开

Try this —

load handel.mat

filename = 'handel.wav';

audiowrite(filename,y,Fs);

clear y Fs

tiledlayout(2,3)

[x,Fs] = audioread('handel.wav');

t = linspace(0, size(x,1)-1, size(x,1))/Fs;

L = numel(t);

ts = 1 / Fs;

nexttile

t1 = 1.000;

t2 = 1.020;

% start and end index in audio data of segment of interest

idx1 = round(t1/ts);

idx2 = round(t2/ts);

% plot it!

x1 = x(idx1:idx2);

plot(x(idx1:idx2)); % plot 1st 20 ms of audio file

title('x_1')

L = numel(x1);

Fn = Fs/2;

NFFT = 2^nextpow2(L);

FTx1 = fft(x1 - mean(x1), NFFT) / L;

Fv1 = linspace(0, 1, NFFT/2+1)*Fn;

Iv1 = 1:numel(Fv1);

nexttile

plot(Fv1, abs(FTx1(Iv1))*2)

grid

xlim([min(Fv1) max(Fv1)])

title('FFT')

nexttile

pspectrum(x1,Fs, 'spectrogram')

colormap(turbo)

nexttile

t3 = 3.000;

t4 = 3.020;

% start and end index in audio data of segment of interest

idx3 = round(t3/ts);

idx4 = round(t4/ts);

% plot it!

x2 = x(idx3:idx4);

plot(x(idx3:idx4));

title('x_2')

L = numel(x2);

Fn = Fs/2;

NFFT = 2^nextpow2(L);

FTx2 = fft(x2 - mean(x2), NFFT) / L;

Fv2 = linspace(0, 1, NFFT/2+1)*Fn;

Iv2 = 1:numel(Fv2);

nexttile

plot(Fv2, abs(FTx2(Iv2))*2)

grid

xlim([min(Fv2) max(Fv2)])

title('FFT')

nexttile

pspectrum(x2,Fs, 'spectrogram')

colormap(turbo)

.

Michael Sugiarto 2022-3-23

Thankyou so much I will try it now!!

Star Strider 2022-3-23

As always, my pleasure!

请先登录，再进行评论。

how to change this code to plot the frequency response of the desired interval?

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

采纳的回答

7 个评论
显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

how to change this code to plot the frequency response of the desired interval?

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

采纳的回答

7 个评论 显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

WeChat

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

7 个评论
显示 5更早的评论隐藏 5更早的评论