zong1.m

音乐信号处理 实现时域延时和混响的处理 实现频域均衡器的设计

close all;
clear all;
clc;
[x,fs]=wavread('come.wav');
i=0;
while(1)
i=i+1;
temp=input('1.时域处理   2.频域处理    请输入：')
if(temp==1)
    disp('时域处理')
    in_key=input('1.延时（FIR梳妆滤波器）  2.混响（IIR滤波器）  3.混响（改进型）  请输入：')
    switch in_key
       case 1
       y=filter([1,zeros(1,6615),0.8],1,[x;zeros(5000,1)]);
       case 2
       y=filter(1,[1,zeros(1,6615),-0.8],[x;zeros(50000,1)]);
       case 3
       y=filter([0.8,zeros(1,6615),1],[1,zeros(1,6615),0.8],[x;zeros(5000,1)]);
       otherwise
            i=0;
       end
elseif(temp==2)
    disp('频域处理')
    in_key=input('1.一阶均衡器  2.二阶均衡器  请输入：')
    switch in_key
    case 1
    K=2.5,alpha=0.7;
    b1=[alpha,-1];a1=[1,-alpha];
    b=0.5*(1+K)*a1+0.5*(1-K)*b1;
    y=filter(b,a1,x);    
    case 2
    K=2.5,alpha=0.7;beta=0.4;
    b1=[alpha,-beta*(1+alpha),1];a1=[1,-beta*(1+alpha),alpha];
    b=0.5*(1+K)*a1+0.5*(1-K)*b1;
    y=filter(b,a1,x);
    otherwise
        i=0;
    end
else 
    break;
end
sound(y,fs);    
end