用matlab演示dsb,DSB-AM的MATLAB仿真

%DSB_AM调制

%题1求解

clear all;

clc;

echo on

t0 = 0.15;%定义信号持续时间

ts = 0.0001;%定义仿真时的信号采样率

fc = 250;%定义载波信号

snr = 20;%定义信噪比，单位dB

fs = 1/ts;%采样频率

df = 0.2;%定义频谱分辩力

t = [0:ts:t0];%时间矢量

snr_lin = 10^(snr/10);%线性的信噪比

%传输的信号

m = [ones(1,ceil(t0/(3*ts))),-2*ones(1,ceil(t0/(3*ts))),zeros(1,ceil(t0/(3*ts)+1))];

c = cos(2*pi*fc.*t);%载波

u = m.*c;%调制

[M,m,df1]=fftseq(m,ts,df);%传输的信号的频谱

M = M/fs;%频谱采样

[U,u,df1]=fftseq(u,ts,df);%调制信号频谱

U=U/fs;%频谱采样

[C,c,df]=fftseq(c,ts,df);%载波信号频谱

C = C/fs;%频率采样

%频谱矢量

f = [0:df1:df1*(length(m)-1)]-fs/2;

%信号的功率

signal_power=spower(u(1:length(t)));

%噪声的功率

noise_power=signal_power/snr_lin;

%噪声的标准差

noise_std=sqrt(noise_power);

%产生最后的噪声

n=noise_std*randn(1,length(u));

r=u+n;%将信号加到噪声上

[N,n,df1]=fftseq(n,ts,df);%对噪声进行fft

N = N/fs;%频率采样

[R,r,df1]=fftseq(r,ts,df);%对最后的信号进行fft

R = R/fs;%进行采样

%显示出没有噪声的情况

figure%显示第一张图片(第一题解)

subplot(3,2,1)%消息信号时域显示

plot(t,m(1:length(t)))

title('原始信号')

axis([0 0.15 -3 3])

subplot(3,2,2)%消息信号频域显示

plot(f,abs(fftshift(M)))

title('原始信号的频率')

axis([-500 500 0 0.15])

subplot(3,2,3)%调制载波时域显示

plot(t,c(1:length(t)))

title('载波信号')

axis([0 0.15 -3 3])

subplot(3,2,4)%调制载波频域显示

plot(f,abs(fftshift(C)))

title('载波的频率')

axis([-500 500 0 0.15])

subplot(3,2,5)%调制信号时域显示

plot(t,u(1:length(t)))

title('调制信号')

axis([0 0.15 -3 3])

subplot(3,2,6)%调制信号频域显示

plot(f,abs(fftshift(U)))

title('调制信号的频率')

axis([-500 500 0 0.15])

%噪声分析

figure%显示第二章图片(第二题解)

subplot(2,2,1)%噪声时域显示

plot(t,n(1:length(t)))

title('采样的噪声')

axis([0 0.15 -1 1])

subplot(2,2,2)%噪声频域显示

plot(f,abs(fftshift(N)))

axis([-500 500 0 0.05])

title('信号频谱')

subplot(2,2,3)%加噪调制信号时域显示

plot(t,r(1:length(t)))

title('信号与噪声')

axis([0 0.15 -3 3])

subplot(2,2,4)%加噪调制信号频域显示

plot(f,abs(fftshift(R)))

title('信号与噪声')

axis([-500 500 0 0.15])