y3.m

基于matlab的数字信号处理经典实验代码及所有完整实验报告

N=40;T1=0.5925;T2=0.1099;
ws1=0.2*pi;wp1=0.35*pi;ws2=0.8*pi;wp2=0.65*pi;
alpha=(N-1)/2;
L=0:N-1;w1=(2*pi/N)*L;
hrs=[zeros(1,5),T2,T1,ones(1,7),T1,T2,zeros(1,9),T2,T1,ones(1,7),T1,T2,zeros(1,4)];
k1=[0:floor((N-1)/2)];
k2=floor((N-1)/2)+1:N-1;
angh=[-alpha*(2*pi)/(0.5*N)*k1,alpha*(2*pi)/(0.5*N)*(N-k2)];
H=hrs.*exp(j*angh)
plot(H);