数字信号处理实验四
第一题结果：
（1）没有增加过渡点
源码如下：
N = 15;
H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小
%H(3,13) = 0.75;H(5,11) = 0.25; %设置过渡点
k = 0:N-1;
A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小
HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n)
freqz(hn,1,256); %画出幅频相频曲线figure(2);
stem(real(hn),'.'); %绘制单位冲激响应的实部
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');
单位脉冲响应曲线
幅频和相频特性曲线
（2）增加过渡点
源码如下：
N = 15;
H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小
H(3) = 0.75;H(13) = 0.75;H(5) = 0.25;H(11) = 0.25; %设置过渡点
k = 0:N-1;
A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小
HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2);
stem(real(hn),'.'); %绘制单位冲激响应的实部
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');
单位脉冲响应曲线
幅频和相频特性曲线
第二题结果：
源码如下：
N=35;
M = N-1;L = M/2;
F = [0:1/L:1]; %设置抽样点的频率，抽样频率必须含0和1。

A = [ones(1,8),zeros(1,L-7)]; %设置抽样点相应的幅值
B = fir2(M,F,A);
[H,W]=freqz(B); %画出滤波器幅频相频曲线
subplot(2,1,1),plot(W/pi,10*log10(abs(H)));grid on;
subplot(2,1,2),plot(W/pi,unwrap(angle(H)));grid on;
% figure,plot(F,A,W/pi,abs(H));grid on;
figure(2);
stem(real(B),'.'); %绘制单位冲激响应的实部% line([0,25],[0,0]);
xlabel('n');ylabel('Real(h(n))');
单位脉冲响应曲线
幅频和相频特性曲线。