数字信号处理报告二.doc

《数字信号处理报告二.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验五functionxt=xtgN=1000;Fs=1000;T=1/Fs;Tp=N*T;t=0:T:(N-1)*T;fc=Fs/10;f0=fc/10;mt=cos(2*pi*f0*t);ct=cos(2*pi*fc*t);xt=mt.*ct;nt=2*rand(1,N)-1;fp=150;fs=200;Rp=0.1;As=70;fb=[fp,fs];m=[0,1];dev=[10^(-As/20),(10^(Rp/20)-1)/(10^(Rp/20)+1)];[n,fo,mo,W]=remez

2、ord(fb,m,dev,Fs);hn=remez(n,fo,mo,W);yt=filter(hn,1,10*nt);xt=xt+yt;fst=fft(xt,N);k=0:N-1;f=k/Tp;subplot(3,1,1);plot(t,xt);grid;xlabel('t/s');ylabel('x(t)');axis([0,Tp/5,min(xt),max(xt)]);title('(a)信号加噪器波形')subplot(3,1,2);plot(f,abs(fst)/max(abs(fst)))

3、;grid;xlabel('f/Hz');ylabel('幅度');axis([0,Fs/2,0,1.2]);title('(b)信号加噪器的频谱')fp=120fs=150tp=0.1ts=60wp=0.24*piws=0.3*piBt=ws-wpN=ceil(6.6*pi/Bt)N=N+mod(N+1,2)wc=(wp+ws)/(2*pi)hn=fir1(N-1,wc,hanning(N))yn=conv(hn,xt)subplot(3,1,3)plot(yn);grid;xlabel('f/H

4、z');ylabel('幅度');title('(c)输出信号yn的频谱')