I want to plot my results individually under each other like a row and not a continues impulse response . Something isn't correct when i try it, Can someone help please?

6 次查看（过去 30 天）

显示更早的评论

codefanta 2019-1-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/440406-i-want-to-plot-my-results-individually-under-each-other-like-a-row-and-not-a-continues-impulse-respo

评论： Star Strider 2019-2-5

采纳的回答： Star Strider

x = zeros(1, 20000);

x(1 : 4000 : 20000);

xresh = reshapex, 4000, 5);

So every 4000fs i want to plotted graph?

采纳的回答

Star Strider 2019-1-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/440406-i-want-to-plot-my-results-individually-under-each-other-like-a-row-and-not-a-continues-impulse-respo#answer_357001

在 MATLAB Online 中打开

I am not certain what you want.

Try this:

x = linspace(0, 20000, 20000)';                             % Create Column Vector
y  = sin(x*pi/500);                                         % Create Data
xresh = reshape(x, 4000, 5);
yresh = reshape(y, 4000, 5);
yresh = bsxfun(@plus, yresh, (0:2:2*(size(yresh,2)-1)));    % Add Offsets
figure
plot(xresh, yresh)
grid

Experiment to get your desired result.

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

Star Strider 2019-1-18

在 MATLAB Online 中打开

I do not have your sound file, so I cannot respond specifically with respect to it.

Copying some code from the audiorecorder documentation, and appending my code to the end of it:

% Record your voice for 5 seconds.
recObj = audiorecorder;
disp('Start speaking.')
recordblocking(recObj, 2.5);
disp('End of Recording.');
% Play back the recording.
play(recObj);
% Store data in double-precision array.
myRecording = getaudiodata(recObj);
% Plot the waveform.
% plot(myRecording);
Fs = 8000;                                                  % Default Sampling Frequency
x = linspace(0, 20000, 20000)'/Fs;                          % Create Column Vector
y = myRecording;
xresh = reshape(x, 4000, 5);
yresh = reshape(y, 4000, 5);
yresh = bsxfun(@plus, yresh, (0:2:2*(size(yresh,2)-1)));    % Add Offsets
figure
plot(xresh, yresh)
grid
figure
plot(xresh(:,1), yresh)
grid
set(gca, 'YTick',(0:2:2*(size(yresh,2)-1)), 'YTickLabel',(1:5))
xlabel('Time (sec)')
ylabel('Signal Segment')

This runs without error.

‘I want to record a sound and plot the results individually rather than one continouse response.’

I am not certain what you are asking for. My code appears to produce the result you requested, as I understand it. The first plot retains the time offsets, the second plot uses the first column of ‘xresh’ and plots everything as a function of it.

codefanta 2019-1-22

Thank you so much it all make sense.

one last question, if i want to change the recording time do i need to change my reshapes ?

Star Strider 2019-1-22

在 MATLAB Online 中打开

As always, my pleasure.

Yes. You can probably make that dynamic by doing something like:

xresh = reshape(x, fix(x/5), 5);

However, you need to choose a recording tiime and sampling frequency that will produce a value for the length of your recorded sound vector that is an integer multiple of the number of segments you want. The reshape function will fail otherwise.

请先登录，再进行评论。

更多回答（2 个）

codefanta 2019-1-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/440406-i-want-to-plot-my-results-individually-under-each-other-like-a-row-and-not-a-continues-impulse-respo#answer_357682

yresh = bsxfun ( @plus , yresh, ( 0 : 2 : 2*(size (yresh,2) -1)));

set (gca, 'ytick', (0 : 2 : 2*( size( yresh,2) -1) , 'yticklabel', (1:5))

Would you be able to explain this 2 bits of code from your code as I dont 100% understand what they do?

Thanks

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

Star Strider 2019-1-23

在 MATLAB Online 中打开

My apologies for my delay in responding. It took a few minuites to get the Fourier transform code to work correctly with ‘yresh’.

‘Would it be a simple case of just continuing with "yresh" ??’

That depends on what you want to do. If you want to do a Fourier transform of your entire vector before you reshape it, then using ‘yresh’ is not the best approach.

If you want to do the Fourier transform on the segments of the vector created in ‘yresh’, then the frequency vector would have to be changed to accommodate it.

Example:

First, re-name the assignment of the bsxfun call to:

yresh2 = bsxfun(@plus, yresh, (0:(size(yresh,2)-1))*2); % Add Offsets

and change this section of code to conform to that change:

figure
plot(xresh, yresh2)
grid
figure
plot(xresh(:,1), yresh2)
grid
% set(gca, 'YTick',(0:2:2*(size(yresh,2)-1)), 'YTickLabel',(1:5))
set(gca, 'YTick',((0:(size(yresh2,2)-1))*2), 'YTickLabel',(1:5))
xlabel('Time (sec)')
ylabel('Signal Segment')

then do the Fourier transform of each row of the original ‘yresh’ matrix, and plot the results:

Fn = Fs/2;                                      % Nyquist Frequency
N = size(yresh,1);                              % Vector Lengths
NFFT = 2^nextpow2(N);                           % FFT Length (Optional)
FTyresh = fft(yresh,NFFT,1)/N;                  % Fourier Transform
Fv = linspace(0, 1, fix(size(FTyresh,1)/2))*Fn; % Frequency Vector
Iv = 1:numel(Fv);                               % Index Vector
figure
hold all
for k1 = 1:size(FTyresh,2)
    subplot(size(FTyresh,2),1,k1)
    plot(Fv, abs(FTyresh(Iv,k1))*2)
    title(sprintf('Row %d', k1))
    ylabel('Amplitude')
end
hold off
xlabel('Frequency (Hz)')

That should get you started with the frequency domain analysis of your signal.

The rest of the code is unchanged.

Experiment to get the result you want.

codefanta 2019-2-5

Again thank you very much

Star Strider 2019-2-5

As always, my pleasure.

请先登录，再进行评论。

codefanta 2019-1-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/440406-i-want-to-plot-my-results-individually-under-each-other-like-a-row-and-not-a-continues-impulse-respo#answer_357740

hello

is it not possible to get my frequency domain and have my results plotted one under each other like figure 1 and 2 before the new code?

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

codefanta 2019-1-27

编辑：codefanta 2019-1-27

Hello

Sorry I didnt reply to your last comment, It was very helpful though!

I was hoping you could talk me through how this funcion works:

"yresh = bsxfun(@plus, yresh, ( 0 : size(yresh,2 ) -1) ) * 2 ); "

So what does the (yresh, 2), the -1 and the *2 do in the function?

thanks

Star Strider 2019-1-27

As always, my pleasure.

That bsxfun (link) call multiplies each column of ‘yresh’ by the corresponding element of the ‘( 0 : size(yresh,2 ) -1) ) * 2’ vector. That vector is equivalent to ‘[0, 2, 4, 6, 8]’. The purpose is to separate the columns by multiples of 2. Since the sound vectors are limited to a magnitude of ±1, this guarantees that they will be separated, and not overlap in the plot.

请先登录，再进行评论。

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

类别

Signal Processing DSP System Toolbox Filter Design and Analysis Frequency Transformations

在 Help Center 和 File Exchange 中查找有关 Frequency Transformations 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by

I want to plot my results individually under each other like a row and not a continues impulse response . Something isn't correct when i try it, Can someone help please?

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

采纳的回答

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

更多回答（2 个）

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

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

另请参阅

类别

标签

Community Treasure Hunt

I want to plot my results individually under each other like a row and not a continues impulse response . Something isn't correct when i try it, Can someone help please?

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

采纳的回答

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

更多回答（2 个）

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

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

另请参阅

类别

标签

Community Treasure Hunt

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

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

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

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