mdsb.m

mdsb matlab

%显示模拟调调制信号的波形及解调方法DSB文件

%信源
clear all;
close all;
dt=0.001;       %时间采样间隔
fm=1;           %信源最高频率
fc=10;          %载波中心频率
T=5;            %信号时长
t=0:dt:T;
mt=sqrt(2)*cos(2*pi*fm*t);%信源
%N0=0.01                   %白噪边带功率谱密度

%DSP modulation
s_dsb=mt.*cos(2*pi*fc*t);
B=2*fm;
%noise=niose_nb(fc,B,N0,t);
%s_dsb=s_dsb+niose;
figure(1)
subplot(311);
plot(t,s_dsb);hold on;%画出DSB信号波形
plot(t,mt,'r--');
title('DSB调制信号')
xlabel('t');
%DSB demolution
rt=s_dsb.*cos(2*pi*fc*t);
rt=rt-mean(rt);
[f,rf]=T2F(t,rt);
[t,rt]=lpf(f,rf,2*fm);
subplot(312)
plot(t,rt);
hold on;
plot(t,mt/2,'r--');
title('相干解调后的信号波形与输入信号的比较');
xlabel('t');
subplot(313)
[f,sf]=T2F(t,s_dsb);%求调制信号的频谱
psf=(abs(sf).^2)/T;
plot(f,psf);
axis([-2*fc 2*fc 0 max(psf)]);
title('DSB信号功率谱');
xlabel('f');