资源描述：
《数字信号处理实验五答案》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、1・信号产生函数mstg清单functionxt=xtg(N)%实验五信号x(t)产生，并显示信号的幅频特性曲线%xt=xtg产生一个长度为N,冇加性高频噪声的单频调幅信号xt,N二1000.,%采样频率Fs二1000Hz%载;波频率fc=Fs/10=100Hz,调制正弦波频率f0=fc/10=10Hz.N=1000;Fs=l000;T=1/Fs;Tp=N*T;t二0:T:(N-1)*T;fc=Fs/10;f0=fc/10;%载波频率fc=Fs/10,单频调制信号频率为f0=Fc/10;int=cos(2*pi*f0*t);%产生讥频止弦波调制信号mt,频率为f0ct=co
2、s(2*pi*fc*t);%产生载波iF弦波信号ct,频率为fcxt=mt.*ct;%相乘产生单频调制信号xtnt=2*rand(l,N)T;%产生随机噪声nI%=======设计高通滤波器hn,用于滤除噪声nt中的低频成分，牛成高通噪声=======fp=150;fs=200;Rp=0.l;As=70;%滤波器指标fb=[fp,fs];m二[0,1];%计算remezord函数所需参数f,m,devdev=[107-As/20),(10"(Rp/20)-l)/(10"(Rp/20)+l)];[n,fo,mo,W]=remezord(fb,m,dev,Fs);%确定rcin
3、ezP^i数所盂参数hn=rcmcz(n,fo,mo,W);%调用rcme函数进行设计，用J：滤除噪声nt中的低频成分yt二f订ter(hn,1,10*nt);%滤除随机噪声中低频成分，牛成高通噪声yt%========================以下为绘图部分=====================================xt=xt+yt;%噪声加信号fst=fft(xt,N);k=0:N-l;f=k/Tp;subplot(2,1,1);plot(t,xt);grid;xlabel('t/s');ylabel('x(t)');axis([0,Tp/5,min
4、(xt),max(xt)]);titleC(a)信号加噪声波形')subplot(2,1,2);plot(f,abs(fst)/max(abs(fst)));grid;titleC(b)信号加噪声的频谱')axis(f0,Es/2,0,1.2]);xlabelCf/Hz');ylabelCif®度')2•产生FIR数字滤波器设计及软件实现clearall;%==usingthefunction,xtg£「toplotthepictureofxtanditsspectrogram==xt=xtg;N=1000;fp=120;fs=150;Rp=0.2;As=60;Fs=100
5、0;%inputthedefinedparameters%todesignthefilterusingth巳Windowsmethodwc=(fp+fs)/Fs;Ethecut-offfrequenceofthefilterB=2*pi*(fs-fp)/Fs;%thewidthofprototypefilterNb=cei丄(11*pi/B);%thelengthofwindowofblackmanhn=fir1(Nb-1,wc,blackman(Nb));Hw=abs(fft(hn,1024));%thefrequencycharacteristicofthefilter
6、%tofilterxtusingfftfiltywt=fftfi]_t(hnzxt』N);%todrawthefigureusingtheWindowsmethodf=[0:1023]*Fs/1024;~figure(2)subplot(2f1,1);plot(f,20*logl0(Hw/max(Hw)));grid;xlabel(1f/Hz1);ylabel(1Amplitude1);title(1thefrequencycharacteristicofthelow-passfilter1);axis([0,Fs/2z-120,20]);t=[0:N-1]/Fs;Tp=N
7、/Fs;subplot(2f1,2);plot(t,ywt);grid;xlabel(1t/s1);ylabel(1y_w(t)T);title(1thesignalwavesafterfilteringnoises1);axis([0,Tp/2,-1,1]);%todesignthefilterusingthebestuniformapproximationmethodfb=[fpzfs];m=[1,0];%definetheparametersofremezorddev=[(10A(Rp/20)-1)/(10A
