高斯濾波器特性ppt課件-資料下載頁(yè)
2025-04-29 04:15本頁(yè)面
【正文】 (Fn,fftshift(Y))。 bc Rf 10?利用 MATLAB的 pwelch函數(shù)計(jì)算功率譜 利用 MATLAB的 pwelch函數(shù)計(jì)算功率譜 bc Rf ?t=::1。Tb=。w0=10*2*pi/Tb。 c=cos(w0*t)。 d=[1 1 1 1 1 1 1 1 1 1]。 w=d*pi/(2*Tb)。 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,100)]1。 y0=c.*s。 t1=::Tb y1=cos((w0+w(1)).*t1)。 for i=1:9 z(1)=0 z(i+1)=z(i)+(d(i)d(i+1))*i*pi/2。 end y2=cos((w0+w(2)).*(t1+Tb)+z(2))。 y3=cos((w0+w(3)).*(t1+2*Tb)+z(3))。 y4=cos((w0+w(4)).*(t1+3*Tb)+z(4))。 y5=cos((w0+w(5)).*(t1+4*Tb)+z(5))。 y6=cos((w0+w(6)).*(t1+5*Tb)+z(6))。 y7=cos((w0+w(7)).*(t1+6*Tb)+z(7))。 y8=cos((w0+w(8)).*(t1+7*Tb)+z(8))。 y9=cos((w0+w(9)).*(t1+8*Tb)+z(9))。 y10=cos((w0+w(10)).*(t1+9*Tb)+z(10))。 y=[y1 y2 y3 y4 y5 y6 y7 y8 y9 y10]。 N=4096。 Fn=:1/N:。 Y=pwelch(y,hamming(33),16,N,39。whole39。)。 Y=10*log10(Y/max(Y)+)。 S=pwelch(s,hamming(33),16,N,39。whole39。)。 S=10*log10(S/max(S)+)。 Y0=pwelch(y0,hamming(33),16,N,39。whole39。)。 Y0=10*log10(Y0/max(Y0)+)。 plot(Fn,fftshift(Y0),39。g39。,Fn,fftshift(Y),39。b39。,Fn,fftshift(S),39。r39。)。grid on。
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1