设 y(t ) =(A+ Mcos(2pifmt) ) cos2pifct ， 其中：M=2，，fm=1000hz，A=2，fc=10000hz。画出调制信号和已调信号频谱。为保证图形显 示效果，取样率 100000，相应的取样间隔 dt=1/100000 s。本例中新函数 fft 、fftshift。 参考程序： dt=1e-5; T=31e-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=Mcos(2pifmt); 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=31e-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); %绘制调制信号频谱

carrier=cos(2*pi*fc*t); %载波信号
envelop_AM=input+A; %包络线信号
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(f,abs(fftshift(INPUT))); % fft 变换的结果需要使用 fftshift 修正偏移
title('调制信号频谱');
subplot(3,1,3);
plot(t,output);
xlabel('时间/s');
ylabel('已调信号');
title('已调信号时域波形');
figure;
subplot(2,1,1);
plot(f,abs(fftshift(CARRIER)));
title('载波信号频谱');
subplot(2,1,2);
plot(f,abs(fftshift(OUTPUT)));
title('已调信号频谱');