hi everyone , I want to compare the voice of male and female and wants to find difference between them,
for example :frequency on per cycle , any suggestion? what to do ?how to do.
MATLAB: How to get frequency of a wav file for per cycle in matlab
frequencySignal Processing Toolbox
Related Solutions
Try this:
% Initialization steps.
clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
clear; % Erase all existing variables. Or clearvars if you want.
workspace; % Make sure the workspace panel is showing.
format long g;format compact;fontSize = 15;% read file into memory */
[wave, fs]=audioread('1.wav'); % sound(wave, fs); % see what it sounds like */
t=0:1/fs:(length(wave)-1)/fs; % and get sampling frequency */
subplot(3, 1, 1)%plot(t,wave)
plot(t, wave, 'b-');title('Wave File', 'FontSize', fontSize);ylabel('Amplitude', 'FontSize', fontSize);xlabel('Length (in seconds)', 'FontSize', fontSize);grid on;L = length(wave);NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(wave,NFFT)/L;f = fs/2*linspace(0,1,NFFT/2+1);% Plot single-sided amplitude spectrum.
subplot(3, 1, 2);y = 2*abs(Y(1:NFFT/2+1));plot(f, y, 'b-');title('Single-Sided Amplitude Spectrum of y(t)', 'FontSize', fontSize);xlabel('Frequency (Hz)', 'FontSize', fontSize);ylabel('|Y(f)|', 'FontSize', fontSize);grid on;subplot(3, 1, 3);y = 2*abs(Y(1:NFFT/2+1));plot(f, y, 'b-');title('Single-Sided Amplitude Spectrum of y(t)', 'FontSize', fontSize);xlabel('Frequency (Hz)', 'FontSize', fontSize);ylabel('|Y(f)|', 'FontSize', fontSize);grid on;% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0, 0.04, 1, 0.96]);xlim([0, 400]);% Smooth with a Savitzky Golay Filter
smoothy = sgolayfilt(y, 2, 101);hold on;plot(f, smoothy, 'r-', 'LineWidth', 2);% Get the upper envelope
uppery = movmax(y,51);% Smooth it out some
windowWidth = 31;uppery = conv(uppery, ones(windowWidth, 1)/windowWidth, 'same');plot(f, uppery, 'r-', 'LineWidth', 2);
What you've put doesn't really make sense. See if this is what you want. It computes the histograms of x and y and then plots them.
clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
clear; % Erase all existing variables.
workspace; % Make sure the workspace panel is showing.
fontSize = 15;x=rand(10,1); y=rand(10,1);% Plot the original x data.
subplot(2, 2, 1);bar(x);xlabel('Index of X', 'fontSize', fontSize);ylabel('Value of X', 'fontSize', fontSize);title('X as a function of Index', 'fontSize', fontSize);% Plot the original y data.
subplot(2, 2, 2);bar(y);xlabel('Index of Y', 'fontSize', fontSize);ylabel('Value of Y', 'fontSize', fontSize);title('Y as a function of Index', 'fontSize', fontSize);% Get the histograms.
[xCounts xValues] = hist(x)[yCounts yValues] = hist(y)% Plot X histogram.
subplot(2, 2, 3);bar(xValues, xCounts);xlabel('Value of X', 'fontSize', fontSize);ylabel('Frequency of Occurrence for X', 'fontSize', fontSize);grid on;title('Histogram of X', 'fontSize', fontSize);subplot(2, 2, 4);bar(yValues, yCounts);xlabel('Value of Y', 'fontSize', fontSize);ylabel('Frequency of Occurrence forY', 'fontSize', fontSize);grid on;title('Histogram of Y', 'fontSize', fontSize);% Enlarge figure to full screen.
set(gcf, 'Position', get(0,'Screensize'));
Best Answer