dsb-sc.m

DSB-SC信号的处理

%DSB-SC信号的处理
close all
clear all
%采样点数的设置
k=input('采样点数N=2^k，请输入k的值，默认为15:');
if isempty(k),
    k=15;
end;
N=2^k;%采样点数
dt=0.01;%时域取点间隔
df=1.0/(N*dt);%频域取点间隔
T=N*dt;%时域截断区间
Bs=N*df/2;%频域截断区间
t=linspace(-T/2,T/2,N);%生成时域采样点
f=linspace(-Bs,Bs,N);%生成频域采样点
fm=1;%基带信号频率
fc=6*fm;%载波频率
m=sin(2*pi*fm*t);%基带信号的产生
c=sin(2*pi*fc*t);%载波信号的产生
s=m.*c;%DSB-SC信号的产生
S=t2f(s,dt);%傅立叶变换
figure(1)%开启一个编号为1的绘图窗口
%设置窗口1左上角的位置在距屏幕左侧100像素、下侧300像素的地方，长为400像素，宽为300像素
set(1,'position',[100,300,400,300])
plot(t,s)%画时域图
%设置窗口的各种参数
title('时域图')
xlabel('t(ms)')
ylabel('s(t)')
axis([-4,4,min(s),max(s)])
grid on
figure(2)%开启一个编号为2的绘图窗口
%设置窗口2左上角的位置在距屏幕左侧500像素、下侧300像素的地方，长为400像素，宽为300像素
set(2,'position',[500,300,400,300])
plot(f,abs(S))%画频域图
title('频域图')
xlabel('f(kHz)')
ylabel('S(f)')
axis([-10,10,min(abs(S)),max(abs(S))])
grid on