求助啊 用matlab 读取一个音乐，然后再把这段音乐的低频给去除掉 该怎么编程啊

我没对声音处理过不知道有没有特殊的啊，我觉得过程应该是这样

[y,fs,nbits]= wavread(‘wavFile1’);这个wavefile1就是你的声音文件了，如‘awav'

y就是音频信号；

fs是采样频率，比如说16000就是每秒16000次；

nbit是采样精度，比如说16就是指16位精度的采样；

然后你去低频，再设计一个高通滤波器，高通滤波器很容易设计的啦。

要注意的是滤波器有FIR和IIR，前者是线性相位，后者是非线性相位。

如果非实时处理信号那就加上filtfilt函数就行了。

下面是一个基于Kaiser窗的高通滤波器设计：

function y=gaotong(signal) ;signal是信号

wp=215/360; 这个是之前的设计程序了，有些地方你自己要改掉；15是阻带，360是采样频率，自己改。比如，我要去掉100HZ的，你就可以选 98——100，这里的15用98代替，下面的3就用100代替

ws=23/360; 3是通带

fpts = [wp ws]; 后面的不用大改。

mag = [0 1];

dev = [001 01];

[N,Wn,beta,ftype] = kaiserord(fpts,mag,dev);

kw = kaiser(N+1,beta);

b = fir1(N,Wn,'high',kaiser(N+1,beta));

[h,omega] = freqz(b,1,512);

figure;

plot(omega/pi,20log10(abs(h)));grid;

xlabel('\omega/\pi'); ylabel('Gain, dB');

title('FIR高通滤波器')

y=filtfilt(b,1,signal);

end

%%%%%%%%%%%%%%%%%%%%%

如果不用FIR，可以用一个IIR的滤波器

IIR滤波器又有几种，比如椭圆滤波器，切比雪夫型滤波器，贝塞尔滤波器等，其实在matlab上只是函数调用不同而已

%设计低通滤波器

%%%%%%%%%%

%巴特沃斯低通滤波器

f_N=16000; %采样率

%设计要求指标

f_p=2400;%通带截止频率

f_s=5000;%阻带截止频率

R_p=3;%通带最大衰减

R_s=25;%阻带最大衰减

Ws=f_s/(f_N/2);

Wp=f_p/(f_N/2); %计算归一化角频率

[n, Wn]=buttord(Wp,Ws,R_p,R_s); %计算阶数和截止频率

[b,a]=butter(n, Wn); %计算

freqz(b,a, 1000, 8000) %作出H(z)的幅频相频图, freqz(b,a, 计算点数, 采样率)

subplot(2,1,1); axis([0 4000 -30 3])

y=filtfilt(b,a,signal);这里是滤波了

切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。

根据频率响应曲线波动位置的不同，切比雪夫滤波器可以分为以下两种： 在通带（或称“通频带”）上频率响应幅度等波纹波动的滤波器称为“I型切比雪夫滤波器”；

'''n'''阶第一类切比雪夫滤波器的幅度与频率的关系可用下列公式表示 ：

:

其中：

而是滤波器在[[截止频率]]的放大率 (''注意'': 常用的以幅度下降3[[分贝]]的频率点作为截止频率的定义不适用于切比雪夫滤波器!)

f>：> 是 n阶[[切比雪夫多项式]]

其中

或:

'''切比雪夫滤波器'''的阶数等于此滤波器的电子线路内的电抗元件数。

切比雪夫滤波器的幅度波动 =[[分贝]]

当 ,切比雪夫滤波器的幅度波动= 3分贝。

如果需要幅度在在阻频带边上衰减得更陡峭，可允许在复平面的 轴上存在零点。但结果会使通频带内振幅波动较大，而在阻频带内对信号抑制较弱。 这种滤波器叫[[椭圆函数滤波器]]或考尔滤波器。 切比雪夫多项式是与棣莫弗定理有关，以递归方式定义的一系列正交多项式序列。 通常，第一类切比雪夫多项式以符号Tn表示， 第二类切比雪夫多项式用Un表示。切比雪夫多项式Tn或Un代表n阶多项式。

切比雪夫多项式在逼近理论中有重要的应用。这是因为第一类切比雪夫多项式的根（被称为切比雪夫节点）可以用于多项式插值。相应的插值多项式能最大限度地降低龙格现象，并且提供多项式在连续函数的最佳一致逼近。 切比雪夫多项式n切比雪夫多项式112 3 4 5 6 7 8 9 10　　　　　　　　　　　　　　 在阻带（或称“阻频带”）上频率响应幅度等波纹波动的滤波器称为“II型切比雪夫滤波器”。

也称倒数切比雪夫滤波器，较不常用，因为频率截止速度不如I型快，也需要用更多的电子元件。II型切比雪夫滤波器在通频带内没有幅度波动，只在阻频带内有幅度波动。

II型切比雪夫滤波器的转移函数为：

分贝

5分贝衰减度相当于; = 06801; 10分贝衰减度相当于; = 03333。

-3分贝频率fH 和截止频率 fC 有如下关系：