图5 SSB频谱
4 所有函数的使用
三角波函数sawtooth：
调用格式为：x = sawtooth(t, width)
功能：产生一个周期为2π、幅度在-1到+1之间的周期性三角波信号。其中width表示最大幅度出现的位置：即在一个周期内，信号从t=0到width×2π时函数值从-1到+1线性增加，而从width×2π到2π又是从+1到-1线性下降。width取值在0 ~ 1之间。 若x = sawtooth(Ωt, width)，则对应的周期为2π/Ω。
abs(x)： 纯量的绝对值或向量的长度
Abs函数返回一个数的绝对值。用法为：result = Abs(number) 其中result是number参数的绝对值。
filter ： 一文数字滤波
filter(fb,fa,mo)，这里fb,fa分别为滤波器的上下限截止频率，而mo为滤波器的输入信号。
其他简单函数：
cla 清除当前坐标轴
clc 清除命令窗口显示
clf 清除当前图形窗口
demo 运行MATLAB演示程序
function MATLAB函数表达式的引导符
grid 给图形加网格线
gtext 在鼠标指定的位置加文字说明
help 启动联机帮助文件显示
hold 当前图形保护模式
imag 求取虚部函数
length 查询向量的文数
linspace 构造线性分布的向量
logspace 构造等对数分布的向量
pi 圆周率π
plot 线性坐标图形绘制
quit 退出MATLAB环境
real 求取实部函数
size 查询矩阵的文数
sqrt 平方根函数
stem 函数序列柄状图形绘制
subplot 将图形窗口分成若干个区域
title 给图形加标题
xlabel 给图形加x坐标说明
ylabel 给图形加y坐标说明
5 MATLAB仿真源程序
5.1调制信号为三角波
Fs=100000; %抽样频率
t=[0:1/Fs:0.01];
m=1;
y=sawtooth(300*2*pi*t); %调制信号
Y=fft(y);
Y=abs(Y(1:length(Y)/2+1)); %调制信号频谱
frqY=[0:length(Y)-1]*Fs/length(Y)/2;
Fc=30000;
y1=amod(m*y,Fc,Fs,'amssb');%ssb解调
n=awgn(y1,30); %已调信号加噪声
Y1=fft(y1);
Y1=abs(Y1(1:length(Y1)/2+1));
frqY1=[0:length(Y1)-1]*Fs/length(Y1)/2;
y2=ademod(n,Fc,Fs,'amssb'); %ssb信号解调
Y2=fft(y2);
Y2=abs(Y2(1:length(Y2)/2+1));
frqY2=[0:length(Y2)-1]*Fs/length(Y2)/2;%解调信号频谱
subplot(3,2,1);
plot(t,y); title('调制信号')
subplot(3,2,2);
plot(frqY,Y);title('调制信号频谱')
axis([0 3000 0 max(Y)]);
subplot(3,2,3);
plot(t,y1); title('已调信号')
subplot(3,2,4);
plot(frqY1,Y1); title('已调信号频谱')
subplot(3,2,5);
plot(t,y2);title('解调信号')
subplot(3,2,6);
plot(frqY2,Y2); title('解调信号频谱')
axis([0 3000 0 max(Y2)]);