cai6.m

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

a=7.865;
ap=0.5;
as=50;
wp=0.6*pi;
ws=0.5*pi;

for n=0:63
    b=a*sqrt(1-((2*n/63-1))^2);
    x=1;
    for k=1:20
        s=1;
        for i=1:k
            s=s*i;
        end
        x=((b/2)^(k)/s)^(2)+x;
    end
    y=1;
    for k=1:20
        s=1;
        for i=1:k
            s=s*i;
        end
        y=((a/2)^(k)/s)^(2)+y;
    end
    w(n+1)=x/y;
end
N=64;
alpha=(N-1)/2;
n=[0:(N-1)];
m=n-alpha+eps;
hd=sin((wp-ws)*m)./(pi*m).*exp(j.*ws.*n);
h=hd.*w;
n=1:1:64;
figure(1)
subplot(2,1,1);
stem(n,w);
title('窗函数的幅度|Wn(jw)|特性')

wc=(wp+ws)/2;
N=-20*log10(min(1-10^(-ap/20),10^(-as/20)));
N=N+rem(N,2);
hn=fir1(N,wc/pi,'high',kaiser(N+1,b));%设计的单位取样响应(实际)
[h0,w0]=freqz(hn,1);
subplot(2,1,2); 
plot(w0/pi,20*log10(abs(h0))); 
grid;
title('滤波器的幅度|H(jw)|特性')
figure(2)
n=1:1:64
subplot(2,1,1);
stem(n,hd);
title('理想冲激响应h(n)');
subplot(2,1,2);
stem(n,h,'r');
title('实际冲激响应h(n)');