kechen4.m

这个是我数字信号处理这门课的课程设计,文件包括用matlab编写的代码和设计报告.这个设计包括卷积演示程序、采样定理演示程序、模拟滤波器设计演示程序、设计切比雪夫I型低通滤波器、切比雪夫I型低通滤波器

fp=100;fs=120;ap=2;as=60;
wp=fp*2*pi;ws=fs*2*pi;
vp=1;vs=fs/fp;
k1=sqrt((10^(0.1*as)-1)/(10^(0.1*ap)-1));
N1=log(k1+sqrt(k1^2-1));
N2=log(vs+sqrt(vs^2-1));
N=ceil(N1/N2);
E=sqrt(10^(0.1*ap)-1);
w=0:300*2*pi;
Ha=1+(E^2)*Cn(w/wp,N).^2;
Ha=1./Ha;
Ha=sqrt(Ha);
subplot(2,1,1)
plot(w/(2*pi),Ha)
title('幅度特性')
hold on
y=1/sqrt(1+E^2);
plot([fp,fp],[0,y],'r')
text(fp-2,0.03,'fp')
plot([0,300],[y,y],'-.g')
text(1,y-0.02,'1/sqrt(1+E^2)')
grid on
subplot(2,1,2)
Ha=20*log10(abs(Ha));
plot(w/(2*pi),Ha)
hold on
title('幅度衰减特性')
plot([100,100],[-200,0],'r');hold on
grid on