设y（t）=(A+M*cos2pi*fm*t)*cos2pi*fc*t 其中M=2,fm=1000hz,A=2,fc=10000hz,取样率为100000，取样间隔为0.00001,补全以下程序 参考程序： dt=1e-5; T=3*1e-3; t=0:dt:T; fs=1/dt; %抽样频率 df=fs/length(t); %频率分辨率 f=-fs/2:df:fs/2-df; %频率绘制区间 M=2; fm=1e3; A=2 fc=1e4 input=M*cos(2*pi*fm*t); INPUT=fft(input); %绘制调制信号频谱 %fft 函数实现从时域到频域的转换（fft—快速傅里叶变换） carrier=cos( ); envelop_AM=input+（ ）; output=( ).*carrier; CARRIER=fft( ); %绘制载波频谱 OUTPUT=fft( ); %绘制已调信号频谱 figure; subplot(3,1,1);plot(t,input); xlabel('时间/s');ylabel('调制信号'); ……. subplot(3,1,1); plot(f,abs(fftshift(INPUT))); % fft 变换的结果需要使用 fftshift 修正偏移 title('基带信号频谱'); % abs 函数为取修正结果的绝对值 …….

dt=1e-5;
T=3*1e-3;
t=0:dt:T;
fs=1/dt; %抽样频率
df=fs/length(t); %频率分辨率
f=-fs/2:df:fs/2-df; %频率绘制区间
M=2;
fm=1e3;
A=2;
fc=1e4;
input=M*cos(2*pi*fm*t);
INPUT=fft(input); %绘制调制信号频谱
%fft 函数实现从时域到频域的转换（fft—快速傅里叶变换）
carrier=cos(2*pi*fc*t);
envelop_AM=input+(A*M/2)*cos(2*pi*fm*t); %计算包络信号
output=envelop_AM.*carrier;
CARRIER=fft(carrier); %绘制载波频谱
OUTPUT=fft(output); %绘制已调信号频谱

figure;
subplot(3,1,1);plot(t,input);
xlabel('时间/s');ylabel('调制信号');
title('时域基带信号');

subplot(3,1,2);plot(t,envelop_AM);
xlabel('时间/s');ylabel('包络信号');
title('时域包络信号');

subplot(3,1,3);plot(t,output);
xlabel('时间/s');ylabel('调制后信号');
title('时域已调信号');

figure;
subplot(3,1,1); plot(f,abs(fftshift(INPUT))); % fft 变换的结果需要使用 fftshift 修正偏移
xlabel('频率/Hz');ylabel('调制信号振幅');
title('基带信号频谱'); % abs 函数为取修正结果的绝对值

subplot(3,1,2); plot(f,abs(fftshift(CARRIER)));
xlabel('频率/Hz');ylabel('载波振幅');
title('载波频谱');

subplot(3,1,3); plot(f,abs(fftshift(OUTPUT)));
xlabel('频率/Hz');ylabel('已调信号振幅');
title('已调信号频谱');