数字信号处理课程设计

figure(3)

subplot(211);plot(x);title('matlab信号波形'); subplot(212);plot(magX);title('matlab信号幅值'); figure(4)

subplot(211);plot(t);title('三路合成信号波形'); subplot(212);plot(magT);title('三路合成信号幅值');

3.3.4 滤波器设计与滤波

低通滤波器的截止频率为6000HZ，为了获取语音信号。

高通滤波器的截止频率为15000HZ，为了获取matlab信号。

%设计低通滤波器的

h(n)______________________________________________________________ for i=1:N

if i==M/2

h1l_p(i)=wc/pi; % 当i==M/2时的极限值 else

h1l_p(i)=sin((i-M/2)*wc)/((i-M/2)*pi); % 低通滤波器的h(n)表达式 end

end %设计高通滤波器的h(n)______________________________________________________________ for i=1:N

if i==M/2

h1h_p(i)=(pi-wh)/pi; % 当i==M/2时的极限值

else

h1h_p(i)=(sin((i-M/2)*pi)-sin((i-M/2)*wh))/((i-M/2)*pi); % 高通滤波器的h(n)表达式 end end

y_hp=conv(h1h_p,t); %高通滤波 for i=1:1000 y_hp(i)=0;

end

for i=65537:75535 y_hp(i)=0;

end

Y2=fft(y_hp);

magY2=abs(Y2);

8

y_lp=conv(h1l_p,t); %低通滤波 for i=1:1000 y_lp(i)=0;

end

for i=65537:75535 y_lp(i)=0;

End 3.3.5

信号的回放与恢复

figure(5);

subplot(211);plot(y_lp);title('恢复语音信号波形'); subplot(212);plot(magY1);title('恢复语音信号幅值'); figure(6);

subplot(211);plot(y_hp);title('matlab声信号波形'); subplot(212);plot(magY2);title('matlab声信号幅值');

fprintf('Finished recording.\\n');

fprintf('Press any key to play the recording...\\n'); pause;

wavplay(y_lp,fs);

4.实验结果与分析 4.1

9

10

4.2

11