在matlab环境下不用matlab自带函数设计一个FIR数字滤波器和FFT数字滤波器

好的，我可以帮你回答这个问题。

首先，我们来了解一下FIR数字滤波器和FFT数字滤波器的概念。

FIR数字滤波器是一类只有前向传递路径的数字滤波器，它的输出只依赖于当前和以前的输入值。FIR数字滤波器的特点是具有线性相位和稳定性，可以实现各种滤波器特性。

FFT数字滤波器是一种基于FFT算法实现的数字滤波器。它将数字信号转换成频域信号，通过去除频域中不需要的频率分量，再将信号变换回时域，实现数字滤波的效果。

接下来，我将分别介绍如何在matlab环境下设计FIR数字滤波器和FFT数字滤波器。

设计FIR数字滤波器：

步骤1：确定滤波器类型和截止频率

根据需要滤波的信号特性，选择合适的滤波器类型和截止频率。常用的滤波器类型有低通、高通、带通和带阻等。

步骤2：计算滤波器系数

使用fir1函数可以计算FIR数字滤波器的系数。该函数需要输入滤波器阶数和截止频率等参数，输出滤波器系数向量。

例如，下面的代码可以生成一个阶数为50、截止频率为0.4的低通FIR数字滤波器系数：

fs = 1000; % 采样频率
fc = 100; % 截止频率
N = 50; % 阶数
b = fir1(N, fc/(fs/2), 'low');

步骤3：滤波信号

使用filter函数可以对信号进行滤波。该函数需要输入滤波器系数和待滤波的信号，输出滤波后的信号。

例如，下面的代码可以对一个长度为1000的随机信号进行滤波：

x = randn(1,1000); % 生成随机信号
y = filter(b,1,x); % 对信号进行滤波

设计FFT数字滤波器：

步骤1：读取信号

使用wavread函数可以读取wav格式的音频文件。该函数需要输入文件名，输出采样率和音频数据。

例如，下面的代码可以读取一个wav格式的音频文件：

[x, fs] = wavread('audio.wav');

步骤2：进行FFT变换

使用fft函数可以对信号进行FFT变换。该函数需要输入信号，输出频域的复数值。

例如，下面的代码可以对读取的音频文件进行FFT变换：

X = fft(x);

步骤3：滤波信号

根据需要滤波的信号特性，将不需要的频率分量置为0。例如，对于一个频率范围在100Hz到200Hz的带通滤波器，可以将频率在该范围之外的分量置为0。

例如，下面的代码可以实现一个带通FFT数字滤波器：

X(abs(f)<100 | abs(f)>200) = 0; % 将频率在100Hz到200Hz之外的分量置为0

步骤4：进行IFFT变换

使用ifft函数可以对频域信号进行IFFT变换，得到滤波后的信号。该函数需要输入频域信号，输出时域的复数值。

例如，下面的代码可以对滤波后的频域信号进行IFFT变换：

y = ifft(X);

以上就是在matlab环境下设计FIR数字滤波器和FFT数字滤波器的方法。希望对你有所帮助！