matlab 测试 数字二次混频

% test2
clear;
clf;
close all
Fs=800000;%采样频率800k
fz=80000;%载波频率80k

fz1=3000;%载波频率3k

fj=79000;%基波频率79k
i=5000;%采样点数
t=(1:i)/Fs;%时间
f=((1:i)-1)*Fs/i; %换算成实际的频率值
B=1;%基波的幅度
A=1;%直流分量
%%%%%%%%%%%%%%%%%%%
%产生基波
%%%%%%%%%%%%%%%%%%%
Mod_Sig=B*sin(2*pi*fj/Fs*(1:i));%基波

figure;plot(t,Mod_Sig);grid on;title('基波');xlabel('时间(S)');ylabel('幅度(V)');
sigI=cos(2*pi*fz/Fs*(1:i));%载波1/3*pi为初相
sigQ=sin(2*pi*fz/Fs*(1:i));%载波1/3*pi为初相
figure;plot(t(1:200),sigI(1:200));grid on;title('载波');xlabel('时间(S)');ylabel('幅度(V)');

fft1re=20*log10(abs(fft(sigI+sigQ*1i)));
figure;plot(f,fft1re);

sigSUM=sigI+sigQ;

fft2re=20*log10(abs(fft(sigSUM)));
figure;plot(f,fft2re);

sigI1=cos(2*pi*fz1/Fs*(1:i)).*sigSUM;%载波1/3*pi为初相
sigQ1=sin(2*pi*fz1/Fs*(1:i)).*sigSUM;%载波1/3*pi为初相

fft3re=20*log10(abs(fft(sigI1+sigQ1*1i)));
figure;plot(f,fft3re);