Matlab仿真AM已调信号相干解调过程

clc;
clear all;

SimTime = 0.1;    % 总的仿真时间（单位：秒）
DispTime = 0.001; % 显示的时间范围（单位：秒）
fs = 450000;      % 系统总的采样频率（单位：Hz）
Ts = 1/fs;        % 采样时间间隔
t = 0:Ts:SimTime; % 时间矢量
N = length（t）;    % 总采样点数
A=1;              %调制信号的幅度因子
m=0.5;            %调制指数

df = fs/N;        % 计算频率分辨率（事实上df=1/SimTime）
f = -fs/2:df:fs/2-df; % 定义频率矢量（频谱图的横坐标）

%被调制信号m（t）=sin6000πt
mt = sin（6000*pi*t）;  
figure（1）;
plot（tmt）;grid on;
axis（[0 DispTime -1.1 1.1]）;
xlabel（‘时间（s）‘）;ylabel（‘电压值（V）‘）;
title（‘被调制信号m（t）=sin6000πt‘）;

figure（2）;
fmt = fft（mt）;
fmt = fftshift（fmt）;
fmt = abs（fmt）;
plot（ffmt）;grid on;
minF = min（abs（fmt））;
maxF = max（abs（fmt））; 
axis（[-4000 4000 minF maxF]）;
xlabel（‘频率（Hz）‘）;ylabel（‘频谱幅度值‘）;
title（‘和信号幅频特性‘）;

%载波信号c（t）=cos90000πt
ct = cos（90000*pi*t）;  
figure（3）;
plot（tct）;grid on;
axis（[0 DispTime -1.1 1.1]）;
xlabel（‘时间（s）‘）;ylabel（‘电压值（V）‘）;
title（‘载波信号m（t）=cos90000πt‘）;

%AM已调制信号Sam（t）= A*[1+m*m（t）]*c（t）
samt= A*（1+m*mt）.*ct;  
figure（4）