您现在的位置：毕业论文 >> 论文 >> 正文

通信原理仿真实验系统设计与实现第9页

更新时间：2009-6-5: 来源：毕业论文

通信原理仿真实验系统设计与实现第9页
subplot(7,1,7)
plot(f,abs(fftshift(R))) %作出频域总信号曲线
title('Signal and noise spectrum'); xlabel('Frequency')

%.抑制双边代调制
clear
echo on
t0=2; %信号持续时间
ts=0.001; %抽样时间间隔
fc=100;    %载波频率
fs=1/ts;
df=0.3;    %频率分辨力
t=[-t0/2:ts:t0/2]; %定义时间序列
%以下三句为定义信号序列
x=sin(200*t);
m=x./(200*t);
m(1001)=1;    %避免产生无穷大的值
c=cos(2*pi*fc.*t); %载波
u=m.*c;         %抑制载波调制
[M,m,df1]=fftseq(m,ts,df); %付里叶变换
M=M/fs;
[U,u,df1]=fftseq(m,ts,df);
U=U/fs;    %频率压缩
f=[0:df1:df1*(length(m)-1)]-fs/2;
%pause
clf
subplot(4,2,1)
plot(t,m(1:length(t)))    %作出为调信号的波形
axis([-0.4,0.4,-0.5,1.1])
xlabel('时间');   title('未调信号')
subplot(4,2,3)
plot(t,c(1:length(t)))
axis([-0.1,0.1,-1.5,1.5]);xlabel('时间')
title('载波');
subplot(4,2,2)
plot(t,u(1:length(t)))
axis([-0.2,0.2,-1,1.2]);
xlabel('时间')
title('已调信号');
subplot(4,1,3)
plot(f,abs(fftshift(M)))
xlabel('频率');    title('未调信号的频谱')
subplot(4,1,4)
plot(f,abs(fftshift(U)))
title('已调信号的频谱'); xlabel('频率')

%单边带调制 ssb.m
clear
echo on
t0=2;%信号持续时间
ts=0.001;%抽样时间
fc=100;%载波频率
fs=1/ts;
df=0.3; %频率分辩力
t=[-t0/2:ts:t0/2];
%以下四句定义信号序列
x=sin(200*t);   m=x./(200*t);
m(1001)=1;      m=m.*m;
c=cos(2*pi*fc.*t);%定义载波同向分量
b=sin(2*pi*fc.*t);%定义载波正交分量
v=m.*c+imag(hilbert(m)).*b;%计算出下边带调幅分量
u=m.*c-imag(hilbert(m)).*b;%计算出上边带调幅分量
[M,m,df1]=fftseq(m,ts,df);%付里叶变换
M=M/fs;
[U,u,df1]=fftseq(u,ts,df);
U=U/fs;
[V,v,df1]=fftseq(v,ts,df);
V=V/fs;
f=[0:df1:df1*(length(m)-1)]-fs/2;
   clf
subplot(5,2,1)
plot(t,m(1:length(t))); axis([-0.1,0.1,-0.5,1.1])
xlabel('时间'); title('未调信号')
subplot(5,2,2)
plot(t,c(1:length(t))); axis([-0.1,0.1,-1.5,1.5])
xlabel('时间'); title('载波')
subplot(5,2,3)
plot(t,u(1:length(t)));
axis([-0.2,0.2,-1,1.2]); xlabel('时间');
title('上边带已调信号')
subplot(5,2,4)
plot(t,v(1:length(t)))
axis([-0.2,0.2,-1,1.2]); xlabel('时间')
title('下边带已调信号');
subplot(5,1,3); plot(f,abs(fftshift(M)))
xlabel('频率'); title('未调信号的频谱')
subplot(5,1,4)
plot(f,abs(fftshift(U)))
title('上边带已调信号的频谱'); xlabel('频率');
subplot(5,1,5); plot(f,abs(fftshift(V)))
title('下边带已调信号的频谱'); xlabel('频率')

%频率调制
echo on
t0=0.15;   %信号持续时间
tz=0.0005; %抽样频率
fc=200;    %载波频率
kf=50;     %调制系数
fz=1/tz;
t=[0:tz:t0]; %定义时间序列
df=0.25;
%定义信号序列
m=[ones(1,t0/(3*tz)),-2*ones(1,t0/(3*tz)),zeros(1,t0/(3*tz)+1)];
int_m(1)=0;        %对M积分，以便后面的调频
for i=1:length(t)-1
    int_m(i+1)=int_m(i)+m(i)*tz;
end
[M,m,df1]=fftseq(m,tz,df);
M=M/fz;
f=[0:df1:df1*(length(m)-1)]-fz/2;
u=cos(2*pi*fc*t+2*pi*kf*int_m); %调制信号调制在载波上
[U,u,df1]=fftseq(u,tz,df);
U=U/fz; %pause;
subplot(4,1,1)
plot(t,m(1:length(t)))   %作出信号波形
axis=([0 0.15 -2.1 2.1]); xlabel('时间')
title('未调信号');   subplot(4,1,2)
plot(t,u(1:length(t)))   %作出调频信号波形
axis=([0 0.15 -2.1 2.1]); xlabel('时间')
title('调频信号');    % pause
subplot(4,1,3)
plot(f,abs(fftshift(M))) %作出原信号频谱
xlabel('频率'); title('信号的频谱')
subplot(4,1,4)
plot(f,abs(fftshift(U))) %作出调频信号频谱