mel-scale.zip_Modified_mel scale

% clear all; % close all; % Q=24; % for f=1:1:4000; % fmel(f)=2146.1-2595*log10(1+(4000-f)/700);%9.635202131004828e+02 % end % k=fmel; % % fmelmax=2146.1-2595*log10(1+(4000-4000)/700); % dfmel = fmelmax/(Q+1) % kb=[]; % % for(i=0:1:Q+1) % kb(i+1)=dfmel*i; % end % % figure; % for m=1:1:4000 % x(m)=128*m/4000; % end % % for(i=2:1:Q+1) % y=[]; % for (j=1:4000) % if k(j)<kb(i-1) % y(j)=0; % end % if k(j)<=kb(i) && k(j)>=kb(i-1) % y(j)=(k(j)-kb(i-1))/(kb(i)-kb(i-1)); % end % if k(j)<=kb(i+1) && k(j)>=kb(i) % y(j)=(kb(i+1)-k(j))/(kb(i+1)-kb(i)); % end % if k(j)>kb(i+1) % y(j)=0; % end % end % hold on; % plot(x,y) % end % clear all; % Q=24; % for f=1:1:4000; % fmel(f)=2595*log10(1+f/700); % end % k=fmel; % % fmelmax=2595*log10(1+4000/700); % dfmel=fmelmax/(Q+1) % kb=[]; % for(i=0:1:25) % kb(i+1)=dfmel*i; % end % figure; % for m=1:1:4000 % x(m)=128*m/4000; % end % % for(i=2:1:25) % y=[]; % for (j=1:4000) % if k(j)<kb(i-1) % y(j)=0; % end % if k(j)<=kb(i) && k(j)>=kb(i-1) % y(j)=(k(j)-kb(i-1))/(kb(i)-kb(i-1)); % end % if k(j)<=kb(i+1) && k(j)>=kb(i) % y(j)=(kb(i+1)-k(j))/(kb(i+1)-kb(i)); % end % if k(j)>kb(i+1) % y(j)=0; % end % end % hold on; % plot(x,y) % end % % ************************* % % ************************* clear all; close all; Q=12; for f=1:1:1300; fmel(f)=2146.1-10.635202131004828e+02-2595*log10(1+(-f+1300)/700);%9.635202131004828e+02 1.082579786899517e+03 end k=fmel; fmelmax=2146.1-10.635202131004828e+02-2595*log10(1+(1300-1300)/700); dfmel = fmelmax/(Q+1); kb=[]; for(i=0:1:Q+1) kb(i+1)=dfmel*i; end figure; for m=1:1:1300 x(m)=41*m/1300; end for(i=2:1:Q+1) y=[]; for (j=1:1300) if k(j)<kb(i-1) y(j)=0; end if k(j)<=kb(i) && k(j)>=kb(i-1) y(j)=(k(j)-kb(i-1))/(kb(i)-kb(i-1)); end if k(j)<=kb(i+1) && k(j)>=kb(i) y(j)=(kb(i+1)-k(j))/(kb(i+1)-kb(i)); end if k(j)>kb(i+1) y(j)=0; end end hold on; plot(x,y) end % ***************** % ***************** % clear all; % Q=12; % start=1231; % endd=2530; % for f=start:1:endd % fmel(f-start+1)=1082.6+2595*log10(1+(f-start)/700); % end % k=fmel; % % % fmelmax=2595*log10(1+(endd-start)/700); % fmelmax=2595*log10(1+(endd-start)/700)-2595*log10(1+(start-start)/700); % dfmel=fmelmax/(Q+1); % kb=[]; % for(i=0:1:Q+1) % kb(i+1)=dfmel*i+1082.6; % end % % figure; % for m= start:1:endd % x1(m-start+1)=128*(m)/4000; % end % % for(i=2:1:Q+1) % y=[]; % for (j=1:1:endd-start+1) % if k(j)<kb(i-1) % y(j)=0; % end % if k(j)<=kb(i) && k(j)>=kb(i-1) % y(j)=(k(j)-kb(i-1))/(kb(i)-kb(i-1)); % end % if k(j)<=kb(i+1) && k(j)>=kb(i) % y(j)=(kb(i+1)-k(j))/(kb(i+1)-kb(i)); % end % if k(j)>kb(i+1) % y(j)=0; % end % end % hold on; % plot(x1,y) % end