數(shù)字下變頻仿真

1、數(shù)字下變頻仿真原理：信號采樣的頻譜調(diào)頻信號：因為f0=30MHz，整體向左平移30MHz。帶寬為5MHz通過仿真得到的時域波形和頻域波形，如下圖所示。clc;clear all;close all;f0=30e6; 中心頻率B=5e6; 帶寬T=30e-6; 脈沖寬度fs=40e6; 采樣頻率N=T*fs; 采樣點數(shù)K=B/T; 頻率變換率ts=1/fs; 采樣周期t=-T/2:ts:T/2-ts;x=cos(2*pi*(f0*t+K*t.2/2);figure(1);title('時域波形');xlabel('point ');figure(2);plot(a

2、bs(fft(x);title('頻域波形'); xlabel('point');I路信號和Q路信號：ddc_i = x.*cos(2*pi*f0*(1:N)/fs); I路信號ddc_q = -x.*sin(2*pi*f0*(1:N)/fs); Q路信號figure(3);subplot(211);plot(t,ddc_i);grid;title('I路波形');subplot(212);plot(t,ddc_q);grid;title('Q路波形');FIR濾波器的設(shè)置在command Windows里輸入fdatool進入濾

3、波器參數(shù)設(shè)置的界面Hd=fir;h,f=freqz(Hd,512);figure(4);plot(f,20*log10(abs(h);gridtitle('FIR低通濾波器的幅頻特性響應(yīng)')hd=Hd.Numerator;fir_i=conv(ddc_i,hd);fir_q=conv(ddc_q,hd);figure(5);subplot(211);plot(fir_i);gridtitle('FIR低通濾波器后的I路信號');subplot(212);plot(fir_q);gridtitle('FIR低通濾波器后的Q路信號');s=fir_i

4、+j*fir_q;figure(6);subplot(211);sf=fftshift(abs(fft(s);plot(sf);grid;title('信號濾波后的頻譜')subplot(212);plot(20*log10(sf/max(sf);grid title('歸一化頻譜')clc;clear all;close all;f0=30e6;B=5e6;T=30e-6;fs=40e6;N=T*fs;K=B/T;ts=1/fs; t=-T/2:ts:T/2-ts;x=cos(2*pi*(f0*t+K*t.2/2);figure(1);plot(x);titl

5、e('時域波形');xlabel('point');figure(2);plot(abs(fft(x);title('頻域波形'); xlabel('point'); ddc_i = x.*cos(2*pi*f0*(1:N)/fs); ddc_q = -x.*sin(2*pi*f0*(1:N)/fs); figure(3);subplot(211);plot(t,ddc_i);grid;title(' I路波形');subplot(212);plot(t,ddc_q);grid;title(' Q路波形&#

6、39;); Hd=fir;h,f=freqz(Hd,512);figure(4);plot(f,20*log10(abs(h);gridtitle(' FIR低通濾波器的幅頻特性響應(yīng)')xlabel('w/pi');hd=Hd.Numerator;fir_i=conv(ddc_i,hd);fir_q=conv(ddc_q,hd);figure(5);subplot(211);plot(fir_i);gridtitle(' FIR低通濾波器后的I路信號');subplot(212);plot(fir_q);gridtitle(' FIR低通濾波器后的Q路信號'); s=fir_i+j*fir_q;figure(6);subplot(211);sf=fftshift(abs(fft(s);plot(