基于matlab的数字信号处理的课程设计报告

基于matlab的数字信号处理的课程设计报告内容摘要：

ubplot(2,2,3)。 plot(w/pi,mag,39。 r39。 )。 axis([0 0 ])。 title(39。 幅频响应 39。 )。 subplot(2,2,4)。 plot(w/pi,db,39。 r39。 )。 axis([0 100 10])。 title(39。 幅频响应 (dB)39。 )。 xlabel(39。 以 \pi为单位的频率 39。 )。 ylabel(39。 对数幅度 /dB39。 )。 %任务二。 %切比雪夫滤波器的设计 Rp=。 As=50。 T=1。 Fs=1/T。 omegap=(2/T)*tan(wp/2)。 omegas=(2/T)*tan(ws/2)。 [cs,ds]=cheb(omegap,omegas,Rp,As)。 %调用函数求解拉普拉斯变换的分子与分母多项式的系数数组 [b,a]=bilinear(cs,ds,Fs)。 %做双线性变换 [db,mag,pha,grd,w]=freqz_m(b,a)。 %频域上采样并记录幅频和相频响应 %作图 figure(2)。 subplot(3,1,1)。 plot(w/pi,mag)。 title(39。 幅度响应 39。 )。 grid。 axis([0 0 ])。 ylabel(39。 幅度 39。 )。 xlabel(39。 以 \pi为单位的频率 39。 )。 subplot(3,1,2)。 plot(w/pi,db)。 title(39。 幅频响应 (dB)39。 )。 grid。 axis([0 60 10])。 xlabel(39。 以 \pi为单位的频率 39。 )。 ylabel(39。 对数幅度 /dB39。 )。 subplot(3,1,3)。 plot(w/pi,pha)。 title(39。 相频响应 39。 )。 grid。 axis([0 4 4])。 xlabel(39。 以 \pi为单位的频率 39。 )。 ylabel(39。 相位 39。 )。 %任务三： fs=10000。 f(1)=120。 f(2)=4980。 f(3)=3 f(4)=1150。 Nq=512。 nq=0:Nq1。 xg=2*cos(2*pi*f(1)*nq/fs)+5*sin(2*pi*f(2)*nq/fs)+8*cos(2*pi*f(3)*nq/fs)+*cos(2*pi*f(4)*nq/fs)。 %模拟信号采样的表达式 %作图 figure(3)。 subplot(2,1,1)。 fplot(39。 xat39。 ,[0,511],39。 c39。 )。 title(39。 x(t)39。 )。 axis([0,512,24,24])。 subplot(2,1,2)。 stem(nq,xg,39。 r39。 )。 title(39。 x(n)39。 )。 axis([0,512,25,25])。 %任务四： Nx=fs/5。 nq1=0:Nx1。 xg1=2*cos(2*pi*f(1)*nq1/fs)+5*sin(2*pi*f(2)*nq1/fs)+8*cos(2*pi*f(3)*n。