matlab程序切比雪夫I型低通数字滤波器

以下两个滤波器都是切比雪夫 I 型数字滤波器，不是巴特沃尔滤波器，请使用者注意！1。低通滤波器使用说明：将下列代码幅值然后以 m 文件保存,文件名要与函数名相同，这里函数名:lowp.function y=lowp（x,f1，f3，rp，rs，Fs)％低通滤波％使用注意事项：通带或阻带的截止频率的选取范围是不能超过采样率的一半％即，f1,f3 的值都要小于 Fs/2％x:需要带通滤波的序列％ f 1:通带截止频率% f 3：阻带截止频率％rp:边带区衰减 DB 数设置%rs：截止区衰减 DB 数设置%FS：序列 x 的采样频率％ rp=0.1;rs=30;%通带边衰减 DB 值和阻带边衰减 DB 值％ Fs=2000;％采样率%wp=2＊pi*f1/Fs；ws=2＊pi＊f3/Fs;% 设计切比雪夫滤波器;[n,wn]=cheb1ord(wp/pi，ws/pi,rp，rs)；[bz1，az1]=cheby1（n,rp,wp/pi）;%查看设计滤波器的曲线[h,w］=freqz(bz1，az1，256，Fs）；h=20*log10（abs（h）)；figure;plot(w，h)；title（'所设计滤波器的通带曲线'）;grid on;%y=filter(bz1,az1，x）;%对序列 x 滤波后得到的序列 yend-—--—----—-—-—-—---——-—-———--———--—--—低通滤波器使用例子的代码fs=2000；t=（1：fs)/fs;ff1=100；ff2=400;x=sin(2＊pi*ff1*t）+sin(2*pi＊ff2＊t)；figure;subplot(211）；plot（t,x）；subplot(212);hua_fft（x，fs，1)；％低通测试％ y=filter(bz1,az1，x）;y=lowp(x，300，350,0。1,20,fs）;figure；subplot（211)；plot(t，y）；subplot(212）;hua_fft（y,fs,1）;％hua_fft(）函数是画频谱图的函数,代码在下面给出，要保存为 m 文件调用％这段例子还调用了我自己写的专门画频谱图的函数，也给出，不然得不出我的结果％画信号的幅频谱和功率谱％频谱使用 matlab 例子表示function hua_fft（y，fs，style，varargin）％当 style=1,画幅值谱;当 style=2，画功率谱；当 style=其他的，那么花幅值谱和功率谱%当 style=1 时，还可以多输入 2 个可选参数%可选输入参数是用来控制需要查看的频率段的％第一个是需要查看的频率段起点％第二个是需要查看的频率段的终点％其他 style 不具备可选输入参数，假如输入发生位置错误nfft= 2^nextpow2(length(y）)；％找出大于 y 的个数的最大的 2 的指数值（自动进算最佳 FFT 步长nfft）％nfft=1024;％人为设置 FFT 的步长 nfft y=y-mean（y）;％去除直流重量y_ft=fft（y,nfft)；％对 y 信号进行 DFT，得到频率的幅值分布y_p=y_ft.*c...