FM信号的调制与解调
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
探究FM信号的调制与解调
方波的调制：
clear all
ts=0.001; %信号抽样时间间隔
t=0:ts:10-ts; %时间向量
fs=1/ts; %抽样频率
df=fs/length(t); %fft的频率分辨率
msg=square(4*t);
msg2=reshape(msg.',1,length(t));
Pm=fft(msg2)/fs; %求消息信号的频谱
f=-fs/2:df:fs/2-df;
subplot(2,1,1)
plot(t,fftshift(abs(Pm)))
title('消息信号频谱')
int_msg(1)=0; %消息信号积分
for ii=1:length(t)-1
int_msg(ii+1)=int_msg(ii)+msg2(ii)*ts;
end
kf=50;
fc=250; %载波频率
Sfm=cos(2*pi*fc*t+2*pi*kf*int_msg); %调频信号
Pfm=fft(Sfm)/fs; % FM信号频谱
subplot(2,1,2)
plot(f,fftshift(abs(Pfm))) % 画出已调信号频谱
title('FM信号频谱')
Pc=sum(abs(Sfm).^2)/length(Sfm) %已调信号功率
Ps=sum(abs(msg2).^2)/length(msg2) %消息信号功率
fm=50;
betaf=kf*max(msg)/fm % 调制指数
W=2*(betaf+1)*fm % 调制信号带宽
012345678910
02
4
6
消息信号频谱
-500-400-300-200-100010020030040050001
2
3
FM 信号频谱
已调信号的功率：Pc =0.5000
消息信号的功率：Ps =1
调制指数：betaf =1
调制信号的带宽：W =200
正弦信号的FM 调制：
clear all
ts=0.001; %信号抽样时间间隔
t=0:ts:10-ts; %时间向量
fs=1/ts; %抽样频率
df=fs/length(t); %fft 的频率分辨率
msg=sin(2*pi*t);
msg2=reshape(msg.',1,length(t));
Pm=fft(msg2)/fs; %求消息信号的频谱
f=-fs/2:df:fs/2-df;
subplot(2,1,1)
plot(t,fftshift(abs(Pm)))
title('消息信号频谱')
int_msg(1)=0; %消息信号积分
for ii=1:length(t)-1
int_msg(ii+1)=int_msg(ii)+msg2(ii)*ts;
end
kf=50;
fc=250; %载波频率
Sfm=cos(2*pi*fc*t+2*pi*kf*int_msg); %调频信号
Pfm=fft(Sfm)/fs; % FM 信号频谱
subplot(2,1,2)
plot(f,fftshift(abs(Pfm))) % 画出已调信号频谱
title('FM 信号频谱')
Pc=sum(abs(Sfm).^2)/length(Sfm) %已调信号功率
Ps=sum(abs(msg2).^2)/length(msg2) %消息信号功率
fm=50;
betaf=kf*max(msg)/fm % 调制指数
W=2*(betaf+1)*fm % 调制信号带宽
012345678910
01
2
3
4
5
消息信号频谱
-500-400-300-200-1000100200300400500
00.5
1
FM 信号频谱
Pc = 0.5000 Ps = 0.5000
betaf =1 W = 200
已调信号的功率：Pc =0.5007
消息信号的功率：Ps =0.4975
调制指数：betaf =1
调制信号的带宽：W =200
锯齿波FM 调制：
clear all
ts=0.001; %信号抽样时间间隔
t=0:ts:10-ts; %时间向量
fs=1/ts; %抽样频率
df=fs/length(t); %fft的频率分辨率msg=sawtooth(8*pi*t);
msg2=reshape(msg.',1,length(t));
Pm=fft(msg2)/fs; %求消息信号的频谱f=-fs/2:df:fs/2-df;
subplot(2,1,1)
plot(t,fftshift(abs(Pm)))
title('消息信号频谱')
int_msg(1)=0; %消息信号积分
for ii=1:length(t)-1
int_msg(ii+1)=int_msg(ii)+msg2(ii)*ts;
end
kf=50;
fc=250; %载波频率
Sfm=cos(2*pi*fc*t+2*pi*kf*int_msg); %调频信号
Pfm=fft(Sfm)/fs; % FM信号频谱subplot(2,1,2)
plot(f,fftshift(abs(Pfm))) % 画出已调信号频谱title('FM信号频谱')
Pc=sum(abs(Sfm).^2)/length(Sfm) %已调信号功率
Ps=sum(abs(msg2).^2)/length(msg2) %消息信号功率
fm=50;
betaf=kf*max(msg)/fm % 调制指数
W=2*(betaf+1)*fm