clc
clear all
close all
df=0.2; %频率分辨力
t0=10; %定义t0信号的持续时间的值
tz=0.001; %定义抽样时间
fa=10; %定义载波频率
a=0.8; %定义调制系数
fz=1/tz;
t=[0:tz:t0]; %定义出抽样点数据
%定义信号m
m=sin(pi.*t);
c=cos(2*pi*fa.*t); %载波信号
m_n=m/max(abs(m));
[M,m,df1]=fftseq(m,tz,df); %傅立叶变换
M=M/fz; %功率缩放
f=[0:df1:df1*(length(m)-1)]-fz/2; %定义频率矢量
u=(1+a*m_n).*c; %将调制信号调制在载波上
[U,u,df1]=fftseq(u,tz,df); %对已调信号作傅立叶变换
U=U/fz; %功率缩放
[C,c,df1]=fftseq(c,tz,df); %对载波信号作傅立叶变换
C=C/fz; %功率缩放
%绘图
subplot(3,2,1)
plot(t,m(1:length(t))) %做出调制信号曲线
axis([0 5 -1.5 1.5]);
xlabel('时间')
title('原信号')
subplot(3,2,2)
plot(f,abs(fftshift(M))) %做出频域的调制信号
axis([-20 20 0 6])
xlabel('频率')
title('原信号频谱')
subplot(3,2,3)
plot(t,c(1:length(t))) %做出载波曲线
axis([0 5 -1.5 1.5]);
xlabel('时间')
title('载波信号')
subplot(3,2,4)
plot(f,abs(fftshift(C))) %做出频域的载波信号
axis([-20 20 0 6])
xlabel('频率')
title('载波信号频谱')
subplot(3,2,5)
plot(t,u(1:length(t))) %做出已调制信号曲线
axis([0 5 -2 2]);
xlabel('时间')
title('已调制信号')
subplot(3,2,6)
plot(f,abs(fftshift(U))) %做出频域的已调制信号
axis([-20 20 0 6])
xlabel('频率')
title('已调信号频谱')
web -browser http://www.ilovematlab.cn/thread-26843-1-1.html
