fm.m

pcm编码设计和最佳量化器设计仿真程序。建议在matlab 6.1 或以上的版本

% 频率调制到MATLAB演示程序。消息信号
% 当 0 < t < t0/3为＋1, 当t0/3 < t < 2t0/3 时为－2，其它情况下为0。
echo off
clear all
close all
t0=.15; % 信号到持续时间
ts=0.0005; % 抽样间隔
fc=200; % 载波频率
kf=50; % 偏移常数
fs=1/ts; % 抽样频率
t=[0:ts:t0]; % 时间向量
df=0.25; % 所需频率分辨率
% message signal
x=ones(1,t0/(3*ts));
y=zeros(1,t0/(3*ts)+1);
m=[x,(-2)*x,y];
int_m(1)=0;
for i=1:length(t)-1  % m的积分
  int_m(i+1)=int_m(i)+m(i)*ts;
end
[M,m,df1]=fftseq(m,ts,df); % 傅立叶变换
M=M/fs; % 缩放
f=[0:df1:df1*(length(m)-1)]-fs/2; % 频率向量
u=cos(2*pi*fc*t+2*pi*kf*int_m); % 调制信号
[U,u,df1]=fftseq(u,ts,df); % 傅立叶变换 
U=U/fs; % 缩放
waitx=input('按任意键显示消息信号和调制信号');
figure(1)
subplot(2,1,1)
plot(t,m(1:length(t)))
axis([0 0.15 -2.1 2.1])
xlabel('Time')
title('The message signal')
subplot(2,1,2)
plot(t,u(1:length(t)))
axis([0 0.15 -2.1 2.1])
xlabel('Time')
title('The modulated signal')
waitx=input('按任意键显示消息信号和调制信号的功率谱');
figure(2)
subplot(2,1,1)
plot(f,abs(fftshift(M))) 
xlabel('Frequency')
title('Magnitude-spectrum of the message signal')
subplot(2,1,2)
plot(f,abs(fftshift(U))) 
title('Magnitude-spectrum of the modulated signal')
xlabel('Frequency')