马赫曾德儿调制器的matlab程序

初学成果

function [ output_args ] = mz_double（ input_args ）
%双臂调制器双边带调制等等（输入两个射频信号）
%   Detailed explanation goes here
tau=1e-13; 
delta_t=tau;
F=1/delta_t; 
N=1e5;
T=delta_t*（N-1）;
t=linspace（-T/2T/2N）; 
f=linspace（-F/2F/2N）;
omega=2*pi*f;

f_RF=1e10;omega_RF=2*pi*f_RF;

V_pi=6;
%V_DC1=10;V_DC2=7;V_RF=1;Phi_1=pi;Phi_2=pi/3;%双边带调制
V_DC1=10;V_DC2=4;V_RF=6;Phi_1=pi;Phi_2=pi/4;%载波抑制调制
%V_DC1=10;V_DC2=7;V_RF=1;Phi_1=pi;Phi_2=pi/2;%单边带调制，抑制上边带
%V_DC1=10;V_DC2=7;V_RF=0.5;Phi_1=pi/2;Phi_2=pi;%单边带调制，抑制下边带
%V_DC1=10;V_DC2=10;V_RF=3;Phi_1=pi;Phi_2=0;%所有奇数边带被抑制
%V_DC1=10;V_DC2=4;V_RF=3;Phi_1=pi;Phi_2=0;%所有偶数边带被抑制
%V_DC1=10;V_DC2=4;V_RF=3;Phi_1=0;Phi_2=pi;%所有偶数边带被抑制（相同）

E_out1_t=0.5*exp（1i*pi*（V_DC1+V_RF*cos（omega_RF*t+Phi_1））/V_pi）;
E_out1_omega=fftshift（fft（E_out1_t）/N*2）;
E_out2_t=0.5*exp（1i*pi*（V_DC2+V_RF*cos（omega_RF*t+Phi_2））/V_pi）;
E_out2_omega=fftshift（fft（E_out2_t）/N*2）;
E_out_t=E_out1_t+E_out2_t;
E_out_omega=fftshift（fft（E_out_t）/N*2）;
i_t=abs（E_out_t）.^2;
i_omega=fftshift（fft（i_t）/N*2）;


subplot（611）;
plot（t（E_out1_t））;axis（[-1e-10 1e-10 -1 1]）;
subplot（612）;
plot（fabs（E_out1_omega））;axis（[-1e11 1e11 -1 1]）;
subplot（613）;
plot（t（E_out2_t））;axis（[-1e-10 1e-10 -1 1]）;
subplot（614）;
plot（fabs（E_out2_omega））;axis（[-1e11 1e11 -1 1]）;
subplot（615）;
%plot（t（E_out_t））;axis（[-1e-10 1e-10 -1 1]）;
plot（fabs（E_out_omega））;axis（[-1e11 1e11 -1 1]）;
subplot（616）;
%plot（t（i_t））;axis（[-1e-10 1e-10 -1 1]）;
plot（fabs（i_omega））;axis（[-1e11 1e11 0 1]）;


return;
 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-12-02 13:43  tiaozhi\
     文件        1619  2017-11-15 20:31  tiaozhi\mz_double.m
     文件        1413  2017-12-02 10:19  tiaozhi\mz_double_s.m
     文件        1195  2017-11-17 13:59  tiaozhi\mz_kongzai.m
     文件        1797  2017-12-02 10:07  tiaozhi\mz_single.m