资源描述:
《《高斯滤波器特性》PPT课件》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、——a=0.25,Bb=2.4——a=0.5,Bb=1.2——a=1,Bb=0.6——a=2,Bb=0.3高斯滤波器的传输函数H(f)f=0:0.05:2.5;y1=exp(-(0.25*f).^2);y2=exp(-(0.5*f).^2);y3=exp(-(f.^2));y4=exp(-(2*f).^2);plot(f,y1,'r',f,y2,'g',f,y3,'b',f,y4,'m');axis([02.50.71]);xlabel('f');ylabel('H(f)');title('TransferFunctionofGau
2、ssianFilter');——a=0.25——a=0.5——a=1——a=2t=0:0.01:1.5;y1=sqrt(pi)/0.25*exp(-((pi*t).^2)/0.25.^2);y2=sqrt(pi)/0.5*exp(-((pi*t).^2)/0.5.^2);y3=sqrt(pi)*exp(-(pi*t).^2);y4=sqrt(pi)/2*exp(-((pi*t).^2)/2.^2);plot(t,y1,'r',t,y2,'g',t,y3,'b',t,y4,'m');xlabel('t/T');ylabel('h(t)
3、');title('ImpulseResponseofGaussianFilter');高斯滤波器的冲激响应h(t)高斯脉冲成形滤波器对矩形脉冲b(t)的响应波形g(t)——a=0.25,Bb=2.4——a=0.5,Bb=1.2——a=1,Bb=0.6——a=2,Bb=0.3t1=-1.5:0.01:1.5;t2=-0.5:0.01:0.5;b1=ones(1,length(t2));y1=sqrt(pi)/0.25*exp(-((pi*t1).^2)/0.25.^2);z1=0.5*(erfc(pi/0.25*(t1-0.5))-
4、erfc(pi/0.25*(t1+0.5)));y2=sqrt(pi)/0.5*exp(-((pi*t1).^2)/0.5.^2);z2=0.5*(erfc(pi/0.5*(t1-0.5))-erfc(pi/0.5*(t1+0.5)));y3=sqrt(pi)*exp(-(pi*t1).^2);z3=0.5*(erfc(pi*(t1-0.5))-erfc(pi*(t1+0.5)));y4=sqrt(pi)/2*exp(-((pi*t1).^2)/2.^2);z4=0.5*(erfc(pi/2*(t1-0.5))-erfc(pi/2*
5、(t1+0.5)));subplot(3,1,1),plot(t2,b1);axis([-1.51.501.5]);xlabel('t/T');ylabel('b(t)');title('rectimpulseinputofGaussianpulseshapingfilter');subplot(3,1,2),plot(t1,y1,'r',t1,y2,'g',t1,y3,'b',t1,y4,'m');xlabel('t/T');ylabel('h(t)');title('impulseresponseofGaussianpulses
6、hapingfilter');subplot(3,1,3),plot(t1,z1,'r',t1,z2,'g',t1,z3,'b',t1,z4,'m');xlabel('t/T');ylabel('g(t)');title('output(rectimpulseresponse)ofGaussianpulseshapingfilter');ASK、MSK信号及其功率谱仿真2、Tb=0.1,fc=20时,比较时域波形和功率谱的情况(差异)注意:采样频率fs=1000,数据率Rb=10数字角频率=模拟角频率*Ts,数字频率=模拟角频率/采
7、样频率=100/1000=0.1N为pwelch的谱计算长度,采样数字频率k/N1、Tb=0.1,fc=100时,观察区间T=1功率谱计算方法:①利用MATLAB的FFT函数②利用MATLAB的pwelch函数d(k)={1,0,0,01,0,1,1,0,0},采样间隔Ts=0.001利用MATLAB的FFT函数计算功率谱t=0.001:0.001:1;Tb=0.1;w0=20*pi/Tb;c=cos(w0*t);subplot(7,1,1),plot(t,c);d=[1-1-1-11-111-1-1];w=d*pi/(2*Tb);
8、s=2*[ones(1,100)zeros(1,100)zeros(1,100)zeros(1,100)ones(1,100)zeros(1,100)ones(1,100)ones(1,100)zeros(1,100)zeros(1,