Question

求频率 matlab 信号处理问题 会的同学看 高分！ 我不会求教。 0.25pi 和0.3pi怎么求

求滤波器的通带截至频率 组带截至频率。 最大衰减等 为什么图像是2000多 下降到4000多 fs=22050;x1=wavread('x.wav');t=0:1/22050:(size(x1)-1)/22050;Au=0.03;d=[Au*cos(2*pi*5000*t)]';x2=x1+d;wp=0.25*pi;ws=0.3*pi;Rp=1;Rs=15;Fs=22050;Ts=1/Fs;wp1=2/Ts*tan(wp/2); %将模拟指标转换成数字指标ws1=2/Ts*tan(ws/2); [N,Wn]=buttord(wp1,ws1,Rp,Rs,'s');%选择滤波器的最小阶数[Z,P,K]=buttap(N); %创建butterworth模拟滤波器[Bap,Aap]=zp2tf(Z,P,K);[b,a]=lp2lp(Bap,Aap,Wn); [bz,az]=bilinear(b,a,Fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换[H,W]=freqz(bz,az); %绘制频率响应曲线figure(1)plot(W*Fs/(2*pi),abs(H))gridxlabel('频率／Hz')ylabel('频率响应幅度')title('Butterworth')f1=filter(bz,az,x2);figure(2)subplot(2,1,1)plot(t,x2) %画出滤波前的时域图title('滤波前的时域波形');subplot(2,1,2)plot(t,f1); %画出滤波后的时域图title('滤波后的时域波形');sound(f1,22050); %播放滤波后的信号F0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x2,1024);subplot(2,1,1);plot(f,abs(y2(1:512))); %画出滤波前的频谱图title('滤波前的频谱')xlabel('Hz');ylabel('fuzhi');subplot(2,1,2)F1=plot(f,abs(F0(1:512))); %画出滤波后的频谱图title('滤波后的频谱')xlabel('Hz');ylabel('fuzhi');

Answer 1

代码很难看清，乱，格式对齐一下吧

追问

wp=0.25*pi;ws=0.3*pi;Rp=1;Rs=15; 怎么算出来截至频率是多少HZ 有的资料上写是通带截至 1000 和阻带截至频率是1200 怎么得出的呢

追答

wp1,ws1是归一化后通带截止频率和阻带截止频率
最后输出图像的时候要恢复过来
图中对应的是：
22050/pi * 0.25pi
和 22050/pi *0.3 pi