tongxin.m

抑制载波双边带调幅和解调的实现, Matlab代码编写实现. 采用Matlab或者其它软件工具实现对信号进行抑制载波双边带调幅（DSB-SC）和解调

%将已知信号f(t)绘制出来
t=-2:0.001:2;%用0.001提高图象分辨率.
f=sinc(200*t);
subplot(231);
plot(t,f);
xlim([-0.1,0.1]);%将横坐标限制在[-0.1,0.1]显示图象.
xlabel('时间t/s');%修饰图象.
ylabel('幅度');
title('已知信号f(t)→时域关系');
grid on;%启动图象网格显示.


%将已知信号f(t)频谱绘制出来
fs=1000;%采样频率为1000Hz.
t1=-2:1/fs:2;
f=sinc(200*t1);
yk=fft(f,65536);%对信号f(t)做65536点的快速傅立叶变换.
yw=2*pi/40000*abs(fftshift(yk));%求频谱图的幅度值.
fw=[-32767:32768]/65536*fs;%求频谱图的横坐标，根据这个公式fw=1/dt*(-N/2+1:N/2)/N求解得到.
subplot(232);
plot(fw,yw);%画出信号f(t)的频谱图.
title('已知信号f(t)的频谱→频域关系');
xlabel('频率/Hz');
ylabel('幅度');
grid on;
xlim([-200,200]);%将横坐标限制在[-0.1,0.1]显示图象.


%将已知载波信号绘制出来
yc=cos(2*pi*200*t);
subplot(2,3,3);
plot(t,yc);
title('载波信号');
xlabel('时间：s');
ylabel('幅度');
grid on;
xlim([-0.1,0.1]);


%绘制已调信号
y4=sinc(t*200).*cos(2*pi*200*t);
subplot(2,3,4);
plot(t,y4);
title('已调信号');
xlabel('时间：s');
ylabel('幅度');
grid on;
xlim([-0.05,0.05]);


%绘制已调信号的频谱图
fs=1000;%已调信号频谱
yk=fft(y4,65536);%对已调信号y4做65536点的快速傅立叶变换.
yw=2*pi/40000*abs(fftshift(yk));%求频谱图的幅度值.
fw=[-32767:32768]/65536*fs;%求频谱图的横坐标，根据这个公式fw=1/dt*(-N/2+1:N/2)/N求解得到.
subplot(2,3,5);
plot(fw,yw);
title('已调信号的频谱');
xlabel('频率：hz');
ylabel('幅度');
grid on;
xlim([-400,400]);


[c,lags]=xcorr(y4,200); %%DSB信号自相关函数
figure             %200表示自相关函数时间т
subplot(2,1,1);   
plot(lags/fs,c);
title('DSB信号自相关函数');
xlabel('t');
ylabel('Rxx(t)');
grid on;


SDSBp=fft(c,65536); %%DSB功率谱
fs1=1000;  %已调信号频谱
yw=2*pi/4000*abs(fftshift(SDSBp));%频谱搬移
fw=[-32767:32768]/65536*fs1;
subplot(2,1,2);
plot(fw,yw);
title('DSB信号功率谱');
xlabel('w');
ylabel('Rxx(t)');
grid on;



y7=y4.*yc;  %%解调信号
figure;
subplot(2,1,1);
plot(t,y7);
title('解调信号');
xlabel('时间：s');
ylabel('幅度');
grid on;
xlim([-0.1,0.1]);


fs = 1000;%采样频率为1000Hz.
t=-2:1/fs:2;
fc = 200;
f=sinc(200*t);
y4=amod(f,fc,fs,'amdsb-sc');
figure;
subplot(211);
plot(t,y4);
 xlim([-0.1, 0.1]);
z=ademod(y4,fc,fs,'amdsb-sc');
subplot(212);
plot(t,z);
xlim([-0.1, 0.1]);