dsb.m

对于dsb信号波形的仿真。包括时域波形和频域波形图

close all	%关闭所有的窗口
clear all	%清除所有的变量
k=input('取样点数=2^k,k=[10]');
if isempty(k),
   k=10;
end;
N=2^k;
dt=0.001;%仿真系统对信号波形的分辨率
df=1.0/(N*dt);
T=N*dt;
Bs=N*df/2;
t=linspace(-T,T,N);
f=linspace(-Bs,Bs,N);
fm=1;%fm=1kHz
x1=2*pi*fm*t;
m=sin(x1);
fc=6*fm;%fc=6fm
x2=2*pi*fc*t;
c=5*sin(x2);
s=m.*c;%dsb-sc信号产生
S=fftshift(fft(s)).*dt;%傅立叶变换
subplot(2,1,1)%开窗口
plot(t,s);
grid on		%显示网格
%设置坐标范围
axis([-2^10*dt/2,2^10*dt/2,min(s),max(s)])	
xlabel('t(ms)')	%给X轴加标签
ylabel('s(t)')	%给Y轴加标签
title('dsb-sc信号时域图')
%画频谱曲线
subplot(2,1,2)
plot(f,abs(S))
grid on
axis([-4*(fm+fc),4*(fm+fc),min(abs(S)),max(abs(S))])
xlabel('f(kHz)')
ylabel('S(f)')
title('dsb-sc信号频域图')