课程设计任务书
学生姓名：毛丽娟专业班级：通信0906
指导教师：尹勇工作单位：信息工程学院
1.题目: SSB信号的仿真分析
初始条件
①MATLAB软件
②数字信号处理与图像处理基础知识
要求完成的主要任务:
调制信号：分别为300Hz正弦信号和三角波信号；载波频率：30kHz；解调方式：同步解调；
要求：画出以下三种情况下调制信号、已调信号、解调信号的波形、频谱以及解调器输入输出信噪比的关系曲线；
1）调制信号幅度=0.8×载波幅度；2）调制信号幅度=载波幅度；
3）调制信号幅度=1.5×载波幅度。
参考书目
时间安排
第1周，安排任务（鉴主15楼实验室）
第1-17周，仿真设计（鉴主13楼计算机实验室）
第18周，完成（答辩，提交报告，演示）
指导教师签名： 2011年月日
系主任（或责任教师）签名：2011年月日
目录
摘要 (1)
Abstract (2)
1 SSB调制原理 (3)
2 函数的使用 (4)
3 SSB调制的实现 (5)
3.1 调制信号为正弦信号 (5)
3.1.1 调制信号幅度=0.8×载波幅度 (5)
3.1.2 调制信号幅度=载波幅度 (8)
3.1.3 调制信号幅度=1.5×载波幅度 (12)
3.2 调制信号为三角波信号 (16)
3.2.1 调制信号幅度=0.8×载波幅度 (16)
3.2.2 调制信号幅度=载波幅度 (19)
3.2.3 调制信号幅度=1.5×载波幅度 (23)
4 心得体会 (27)
5 参考文献 (28)
摘要
新版的MATLAB增强了图形处理功能，并在WINDOWS环境下运行。现今，MATLAB的发展已大大超出了“矩阵实验室”的范围，它的配备了涉及到自动控制、信息处理、计算机仿真等种类繁多的工具箱（Tool Box），这些工具箱有数理统计、信号处理、系统辨识、最优化、稳健等等。本次课程设计主要利用MATLAB 软件对通信原理中的模拟信号进行进行抑制双边带调制的仿真分析，即SSB信号调试与解调的仿真分析。设计中主要是对SSB已调信号进行时域和频域分析。
Abstract
The new version of MATLAB enhanced graphics processing functions, and in WINDOWS environment. Today, MATLAB development has gone far beyond the "matrix laboratory" the scope, it is equipped with a related to automatic control, information processing, computer simulation, such as a wide variety of toolbox (Tool Box), a few of these toolbox of mathematical statistics, signal processing, system identification, optimization, sound and so on. The curriculum design using MATLAB software, the main principle of communication of analog signals to inhibit bilateral band modulation simulation analysis, namely, SSB signal demodulation debugging and simulation analysis. Design mainly SSB signals have had time and frequency domain analysis.
1 SSB调制原理
与标准幅度调制相比，单边带调制（SSB）对于频谱和输出功率的利用率更高。尽管很少用于数据传送，SSB仍广泛地用于HF和VHF低端的语音通讯。双边带调制信号包含有两个完全相同的基带信号，即上、下边带。由于两个边带含的信息相同，因而从信息传输角度考虑，传送一个边带同样可以达到信息传输的目的,本设计只考虑上边带信号。单边带调制，就是通过某种办法，只传送一个边带的调制方法。
设调制信号为单频信号f(t)=A
m cosω
m
t，载波为c(t)=cosω
c
t，则调制后的
双边带时域波形为：
S DSB (t)＝A
m
cosω
m
t cos t＝[A
m
cos(ω
c
+ω
m
)t+ A
m
cos(ω
c
-ω
m
)t] /2
保留上边带，波形为：
S USB (t)＝[A
m
cos(ω
c
+ω
m
)t]/2＝A
m
(cosω
c
tcosω
m
t-sinω
c
tsinω
m
t) /2
保留下边带，波形为：
S LSB (t)＝[A
m
cos(ω
c
-ω
m
)t]/2＝A
m
(cosω
c
tcosω
m
t+sinω
c
tsinω
m
t) /2
上两式中的第一项与调制信号和载波信号的乘积成正比，称为同相分量；而第二项的乘积则是调制信号与载波信号分别移相90°后相乘的结果，称为正交分量。
SSB调制框图如下
图1 SSB调制框图
2 函数的使用
三角波函数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演示程序
grid 给图形加网格线
gtext 在鼠标指定的位置加文字说明
hold 当前图形保护模式
length 查询向量的维数
linspace 构造线性分布的向量
logspace 构造等对数分布的向量
pi 圆周率π
plot 线性坐标图形绘制
subplot 将图形窗口分成若干个区域
title 给图形加标题
xlabel 给图形加x坐标说明
ylabel 给图形加y坐标说明
3 SSB调制的实现
3.1 调制信号为正弦信号
3.1.1 调制信号幅度=0.8×载波幅度
程序为：
Fs=100000; %抽样频率
t=[0:1/Fs:0.01]; %抽样间隔
Fc=30000; %载波频率
m=cos(300*2*pi*t); %调制信号
X=fft(m);
X=abs(X(1:length(X)/2+1)); %调制信号频谱
frqX=(0:length(X)-1)*Fs/length(X)/2
sm = modulate(m,Fc,Fs,'amssb'); %对信号进行调制
Y=fft(sm);
Y=abs(Y(1:length(Y)/2+1)); %已调信号频谱
frqY=(0:length(Y)-1)*Fs/length(Y)/2
k=sm+awgn(sm,1,0); %对已调信号加信噪比为1db的噪声
Y1=fft(k); Y1=abs(Y1(1:length(Y1)/2+1)); %加入噪声后的频谱
frqY1=(0:length(Y1)-1)*Fs/length(Y1)/2;
sn=demod(k,Fc,Fs,'amssb'); %ssb信号解调
Z=fft(sn);
Z=abs(Z(1:length(Z)/2+1));
frqZ=[0:length(Z)-1]*Fs/length(Z)/2; %解调信号频谱
figure(1);
plot(t,m); title('调制信号')
figure
plot(frqX,X);title('调制信号频谱') axis([0 3000 0 max(X)]); figure
plot(t,sm);title('已调信号') figure
plot(frqY,Y); title('已调信号频谱') figure
plot(frqY1,Y1);title('已调信号加噪声的频谱') figure
plot(t,sn);title('解调信号') figure
plot(frqZ,Z); title('解调信号频谱') axis([0 3000 0 max(Z)]);
00.0010.0020.0030.0040.0050.0060.0070.0080.0090.01
-0.8
-0.6-0.4-0.200.20.40.60.8调制信号
图3-1（a) 调制信号时域图
050010001500200025003000
50100150200250300350400调制信号频谱
图3-1(b) 调制信号频谱图
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
0.01
-1-0.8-0.6-0.4-0.200.20.40.60.81已调信号
图3-1(c) 已调信号时域图
00.51 1.52 2.53 3.54 4.5
5
x 10
4
50
100
150
200
250
300
350
已调信号频谱
图3-1(d) 已调信号频谱图
00.51 1.52 2.53 3.54 4.5
5
x 10
4
100
200300400500600700800已调信号加噪声的频谱
图3-(e) 已调信号加噪声的频谱图
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
0.01
-2-1.5-1-0.500.511.522.53解调信号
图3-1(f) 解调信号时域图
050010001500200025003000
50100150200250300350400解调信号频谱
图3-1(g) 解调信号频谱图
3.1.2 调制信号幅度=载波幅度
程序为：
Fs=100000; %抽样频率 t=[0:1/Fs:0.01]; %抽样间隔 Fc=30000; %载波频率 m=cos(300*2*pi*t); %调制信号 X=fft(m);
X=abs(X(1:length(X)/2+1)); %调制信号频谱
frqX=(0:length(X)-1)*Fs/length(X)/2
sm = modulate(m,Fc,Fs,'amssb'); %对信号进行调制
Y=fft(sm);
Y=abs(Y(1:length(Y)/2+1)); %已调信号频谱
frqY=(0:length(Y)-1)*Fs/length(Y)/2
k=sm+awgn(sm,1,0); %对已调信号加信噪比为1db的噪声
Y1=fft(k); Y1=abs(Y1(1:length(Y1)/2+1)); %加入噪声后的频谱
frqY1=(0:length(Y1)-1)*Fs/length(Y1)/2;
sn=demod(k,Fc,Fs,'amssb'); %ssb信号解调
Z=fft(sn);
Z=abs(Z(1:length(Z)/2+1));
frqZ=[0:length(Z)-1]*Fs/length(Z)/2; %解调信号频谱
figure(1);
plot(t,m); title('调制信号')
figure
plot(frqX,X);title('调制信号频谱')
axis([0 3000 0 max(X)]);
figure
plot(t,sm);title('已调信号')
figure
plot(frqY,Y); title('已调信号频谱')
figure
plot(frqY1,Y1);title('已调信号加噪声的频谱')
figure
plot(t,sn);title('解调信号')
figure
plot(frqZ,Z); title('解调信号频谱')
axis([0 3000 0 max(Z)]);
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
0.01
-1-0.8-0.6-0.4-0.200.20.40.60.8
1调制信号
图3-2(a) 调制信号时域图
500
1000
1500
2000
2500
3000
050100150200250300350400450500调制信号频谱
图3-2(b) 调制信号频谱图
00.001
0.0020.0030.0040.0050.0060.0070.0080.0090.01
-1.5
-1
-0.5
0.5
1
1.5
已调信号
图3-2(c) 已调信号时域图
00.511.522.533.544.55
x 10
4
50100150200250300350400450已调信号频谱
图3-2(d) 已调信号频谱图
00.511.522.533.544.55
x 10
4
100200300400500600700800900已调信号加噪声的频谱
图3-2(e) 已调信号加噪声的频谱图
0.0010.00
20.0030.0040.0050.0060.0070.0080.0090.01
-3-2
-1
1
2
3
4
解调信号
图3-2(f) 解调信号时域图
500
1000
1500
2000
2500
3000
050100150200250300350400450500解调信号频谱
图3-2(g) 解调信号频谱图
3.1.3 调制信号幅度=1.5×载波幅度
程序为：
Fs=100000; %抽样频率 t=[0:1/Fs:0.01]; %抽样间隔 Fc=30000; %载波频率 m=1.5*cos(300*2*pi*t); %调制信号 X=fft(m);
X=abs(X(1:length(X)/2+1)); %调制信号频谱 frqX=(0:length(X)-1)*Fs/length(X)/2
sm = modulate(m,Fc,Fs,'amssb'); %对信号进行调制 Y=fft(sm);
Y=abs(Y(1:length(Y)/2+1)); %已调信号频谱 frqY=(0:length(Y)-1)*Fs/length(Y)/2
k=sm+awgn(sm,1,0); %对已调信号加信噪比为1db 的噪声
Y1=fft(k); Y1=abs(Y1(1:length(Y1)/2+1)); %加入噪声后的频谱 frqY1=(0:length(Y1)-1)*Fs/length(Y1)/2; sn=demod(k,Fc,Fs,'amssb'); %ssb 信号解调
Z=fft(sn);
Z=abs(Z(1:length(Z)/2+1));
frqZ=[0:length(Z)-1]*Fs/length(Z)/2; %解调信号频谱 figure(1);
plot(t,m); title('调制信号') figure
plot(frqX,X);title('调制信号频谱') axis([0 3000 0 max(X)]); figure
plot(t,sm);title('已调信号') figure
plot(frqY,Y); title('已调信号频谱') figure
plot(frqY1,Y1);title('已调信号加噪声的频谱') figure
plot(t,sn);title('解调信号') figure
plot(frqZ,Z); title('解调信号频谱') axis([0 3000 0 max(Z)]);
00.0010.0020.0030.0040.0050.0060.0070.0080.0090.01
-1.5
-1
-0.5
0.5
1
1.5
调制信号
图3-3(a)调制信号时域图
050010001500200025003000
100
200300400500600700调制信号频谱
图3-3(b) 调制信号频谱图
00.0010.0020.0030.0040.0050.0060.0070.0080.0090.01
-2
-1.5-1-0.500.511.52已调信号
图3-3(c) 已调信号时域图
00.511.522.533.544.55
x 10
4
100
200
300
400
500
600
700
已调信号频谱
图3-3(d) 已调信号频谱图
00.51 1.52 2.53 3.54 4.5
5
x 10
4
200
400
600
800
1000
1200
1400
已调信号加噪声的频谱
图3-3(e) 已调信号加噪声的频谱图
0.0010.0020.0030.0040.0050.0060.0070.0080.0090.01
-4-3
-2
-1
1
2
3
解调
信号
图3-3(f) 解调信号时域图
050010001500200025003000
100
200300400500600700解调信号频谱
图3-3(g) 解调信号频谱图
3.2 调制信号为三角波信号
3.2.1 调制信号幅度=0.8×载波幅度
程序为：
Fs=100000; %抽样频率
t=[0:1/Fs:0.01]; %抽样间隔
Fc=30000; %载波频率
m=0.8*sawtooth(300*2*pi*t); %调制信号
X=fft(m);
X=abs(X(1:length(X)/2+1)); %调制信号频谱
frqX=(0:length(X)-1)*Fs/length(X)/2
sm = modulate(m,Fc,Fs,'amssb'); %对信号进行调制
Y=fft(sm);
Y=abs(Y(1:length(Y)/2+1)); %已调信号频谱
frqY=(0:length(Y)-1)*Fs/length(Y)/2
k=sm+awgn(sm,1,0); %对已调信号加信噪比为1db的噪声
Y1=fft(k); Y1=abs(Y1(1:length(Y1)/2+1)); %加入噪声后的频谱
frqY1=(0:length(Y1)-1)*Fs/length(Y1)/2;
sn=demod(k,Fc,Fs,'amssb'); %ssb信号解调
Z=fft(sn);
Z=abs(Z(1:length(Z)/2+1));
frqZ=[0:length(Z)-1]*Fs/length(Z)/2; %解调信号频谱
figure(1);
plot(t,m); title('调制信号')
figure
plot(frqX,X);title('调制信号频谱')
axis([0 3000 0 max(X)]);
figure
plot(t,sm);title('已调信号') figure
plot(frqY,Y); title('已调信号频谱') figure
plot(frqY1,Y1);title('已调信号加噪声的频谱') figure
plot(t,sn);title('解调信号') figure
plot(frqZ,Z); title('解调信号频谱') axis([0 3000 0 max(Z)]);
00.0010.0020.0030.0040.0050.0060.0070.0080.0090.01
-0.8
-0.6-0.4-0.200.20.40.6
0.8调制信号
图3-4(a) 调制信号时域图
050010001500200025003000
50
100
150
200
250
调制信号频谱
图3-4(b) 调制信号频谱图
00.0010.0020.0030.0040.0050.0060.0070.0080.0090.01
-3
-2
-1
1
2
3
已调信号
图3-4(c) 已调信号时域图
0.5
1
1.5
2
2.5
3
3.5
4
4.5
5
x 10
4
020406080100120140160180200已调信号频谱
图3-4(d) 已调信号频谱图
00.51 1.52 2.53 3.54 4.5
5
x 10
4
50100150200250300350400已调信号加噪声的频谱
图3-4(e) 已调信号加噪声的频谱图
课题一： 连续时间信号和系统时域分析及MATLAB实现 课题要求： 深入研究连续时间信号和系统时域分析的理论知识。利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能，实现连续时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB实现常用连续时间信号的时域波形（通过改变参数，分析其时域特性）。 1、单位阶跃信号， 2、单位冲激信号， 3、正弦信号， 4、实指数信号， 5、虚指数信号， 6、复指数信号。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘， 4、微分， 5、积分 三、用MATLAB实现信号的时域变换（参数变化，分析波形变化） 1、反转， 2、使移（超时，延时）， 3、展缩， 4、倒相， 5、综合变化 四、用MATLAB实现信号简单的时域分解 1、信号的交直流分解， 2、信号的奇偶分解 五、用MATLAB实现连续时间系统的卷积积分的仿真波形 给出几个典型例子，对每个例子，要求画出对应波形。 六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。 给出几个典型例子，四种调用格式。 七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题二： 离散时间信号和系统时域分析及MATLAB实现。 课题要求： 深入研究离散时间信号和系统时域分析的理论知识。利用MATLAB强大的图
形处理功能、符号运算功能以及数值计算功能，实现离散时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB绘制常用信号的时域波形（通过改变参数分析其时域特性） 1、单位序列， 2、单位阶跃序列， 3、正弦序列， 4、离散时间实指数序列， 5、离散时间虚指数序列， 6、离散时间复指数序列。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘。 三、用MATLAB实现信号的时域变换（参数变化，分析波形的变化） 1、反转， 2、时移（超时，延时）， 3、展缩， 4、倒相。 四、用MATLAB实现离散时间系统卷积和仿真波形 给出几个典型例子，对每个例子要求画出e(k)，h(k)，e(i)，h(i)，h(-i)，Rzs(k)波形。 五、用MATLAB实现离散时间系统的单位响应，阶跃响应的仿真波形 给出几个典型例子，四中调用格式。 六、用MATLAB实现离散时间系统对实指数序列信号的零状态响应的仿真波形 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题三： 连续时间信号傅里叶级数分析及MATLAB实现。 课题要求： 深入研究连续时间信号傅里叶级数分析的理论知识，利用MATLAB强大的图形处理功能，符号运算功能以及数值计算功能，实现连续时间周期信号频域分析的仿真波形。 课题内容： 一、用MATLAB实现周期信号的傅里叶级数分解与综合 以周期矩形波信号为例，绘出包含不同谐波次数的合成波形，观察合成波形与原矩形 波形之间的关系及吉布斯现象。
产生离散衰减正弦序列()π0.8sin 4n x n n ?? = ??? , 010n ≤≤，并画出其波形图。 n=0:10; x=sin(pi/4*n).*0.8.^n; stem(n,x);xlabel( 'n' );ylabel( 'x(n)' ); 用MATLAB 生成信号()0sinc at t -, a 和0t 都是实数,410t -<<，画波形图。观察并分析a 和0t 的变化对波形的影响。 t=linspace(-4,7); a=1;
t0=2; y=sinc(a*t-t0); plot(t,y); t=linspace(-4,7); a=2; t0=2; y=sinc(a*t-t0); plot(t,y);
t=linspace(-4,7); a=1; t0=2; y=sinc(a*t-t0); plot(t,y);
三组对比可得a 越大最大值越小，t0越大图像对称轴越往右移 某频率为f 的正弦波可表示为()()cos 2πa x t ft =，对其进行等间隔抽样，得到的离散样值序列可表示为()()a t nT x n x t ==，其中T 称为抽样间隔，代表相邻样值间的时间间隔，1 s f T = 表示抽样频率，即单位时间内抽取样值的个数。抽样频率取40 Hz s f =，信号频率f 分别取5Hz, 10Hz, 20Hz 和30Hz 。请在同一张图中同时画出连续信号()a x t t 和序列()x n nT 的波形图，并观察和对比分析样值序列的变化。可能用到的函数为plot, stem, hold on 。 fs = 40; t = 0 : 1/fs : 1 ; % ?μ?ê·?±e?a5Hz,10Hz,20Hz,30Hz f1=5; xa = cos(2*pi*f1*t) ; subplot(1, 2, 1) ;
设计出一套完整的系统，对信号进行频谱分析和滤波处理； 1．产生一个连续信号，包含低频，中频，高频分量，对其进行采样，进行频谱分析，分别设计三种高通，低通，带通滤波器对信号进行滤波处理，观察滤波后信号的频谱。 2.采集一段含有噪音的语音信号（可以录制含有噪音的信号，或者录制语音后再加进噪音信号），对其进行采样和频谱分析，根据分析结果设计出一合适的滤波器滤除噪音信号。 %写上标题 %设计低通滤波器： [N,Wc]=buttord() %估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc [a,b]=butter(N,Wc); %设计Butterworth低通滤波器 [h,f]=freqz(); %求数字低通滤波器的频率响应 figure(2); % 打开窗口2 subplot(221); %图形显示分割窗口 plot(f,abs(h)); %绘制Butterworth低通滤波器的幅频响应图 title(巴氏低通滤波器''); grid; %绘制带网格的图像 sf=filter(a,b,s); %叠加函数S经过低通滤波器以后的新函数 subplot(222); plot(t,sf); %绘制叠加函数S经过低通滤波器以后的时域图形 xlabel('时间(seconds)'); ylabel('时间按幅度'); SF=fft(sf,256); %对叠加函数S经过低通滤波器以后的新函数进行256点的基—2快速傅立叶变换 w= %新信号角频率 subplot(223); plot()); %绘制叠加函数S经过低通滤波器以后的频谱图 title('低通滤波后的频谱图'); %设计高通滤波器 [N,Wc]=buttord() %估算得到Butterworth高通滤波器的最小阶数N和3dB截止频率Wc [a,b]=butter(N,Wc,'high'); %设计Butterworth高通滤波器 [h,f]=freqz(); %求数字高通滤波器的频率响应 figure(3); subplot(221); plot()); %绘制Butterworth高通滤波器的幅频响应图 title('巴氏高通滤波器'); grid; %绘制带网格的图像 sf=filter(); %叠加函数S经过高通滤波器以后的新函数 subplot(222); plot(t,sf); ;%绘制叠加函数S经过高通滤波器以后的时域图形 xlabel('Time(seconds)'); ylabel('Time waveform'); w; %新信号角频率 subplot(223);
现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 (一)Ansoft公司的仿真工具 现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 Ansoft的信号完整性工具采用一个仿真可解决全部设计问题： SIwave是一种创新的工具，它尤其适于解决现在高速PCB和复杂IC封装中普遍存在的电源输送和信号完整性问题。 该工具采用基于混合、全波及有限元技术的新颖方法，它允许工程师们特性化同步开关噪声、电源散射和地散射、谐振、反射以及引线条和电源/地平面之间的耦合。该工具采用一个仿真方案解决整个设计问题，缩短了设计时间。 它可分析复杂的线路设计，该设计由多重、任意形状的电源和接地层，以及任何数量的过孔和信号引线条构成。仿真结果采用先进的3D图形方式显示，它还可产生等效电路模型，使商业用户能够长期采用全波技术，而不必一定使用专有仿真器。 (二)SPECCTRAQuest Cadence的工具采用Sun的电源层分析模块： Cadence Design Systems的SpecctraQuest PCB信号完整性套件中的电源完整性模块据称能让工程师在高速PCB设计中更好地控制电源层分析和共模EMI。 该产品是由一份与Sun Microsystems公司签署的开发协议而来的，Sun最初研制该项技术是为了解决母板上的电源问题。 有了这种新模块，用户就可根据系统要求来算出电源层的目标阻抗；然后基于板上的器件考虑去耦合要求，Shah表示，向导程序能帮助用户确定其设计所要求的去耦合电容的数目和类型；选择一组去耦合电容并放置在板上之后，用户就可运行一个仿真程序，通过分析结果来发现问题所在。 SPECCTRAQuest是CADENCE公司提供的高速系统板级设计工具，通过它可以控制与PCB layout相应的限制条件。在SPECCTRAQuest菜单下集成了一下工具： （1）SigXplorer可以进行走线拓扑结构的编辑。可在工具中定义和控制延时、特性阻抗、驱动和负载的类型和数量、拓扑结构以及终端负载的类型等等。可在PCB详细设计前使用此工具，对互连线的不同情况进行仿真，把仿真结果存为拓扑结构模板，在后期详细设计中应用这些模板进行设计。 （2）DF/Signoise工具是信号仿真分析工具，可提供复杂的信号延时和信号畸变分析、IBIS 模型库的设置开发功能。SigNoise是SPECCTRAQUEST SI Expert和SQ Signal Explorer Expert进行分析仿真的仿真引擎，利用SigNoise可以进行反射、串扰、SSN、EMI、源同步及系统级的仿真。 （3）DF/EMC工具——EMC分析控制工具。 （4）DF/Thermax——热分析控制工具。 SPECCTRAQuest中的理想高速PCB设计流程： 由上所示，通过模型的验证、预布局布线的space分析、通过floorplan制定拓朴规则、由规
课程设计任务书 学生姓名：董航专业班级：电信1006班 指导教师：阙大顺，李景松工作单位：信息工程学院 课程设计名称：Matlab应用课程设计 课程设计题目：Matlab运算与应用设计5 初始条件： 1.Matlab6.5以上版本软件； 2.课程设计辅导资料：“Matlab语言基础及使用入门”、“Matlab及在电子信息课程中的应 用”、线性代数及相关书籍等； 3.先修课程：高等数学、线性代数、电路、Matlab应用实践及信号处理类相关课程等。 要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 1.课程设计内容：根据指导老师给定的7套题目，按规定选择其中1套完成； 2.本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析， 针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结。具体设计要求包括： ①初步了解Matlab、熟悉Matlab界面、进行简单操作； ②MATLAB的数值计算：创建矩阵矩阵运算、多项式运算、线性方程组、数值统计； ③基本绘图函数：plot, plot3, mesh, surf等，要求掌握以上绘图函数的用法、简单图形 标注、简单颜色设定等； ④使用文本编辑器编辑m文件，函数调用； ⑤能进行简单的信号处理Matlab编程； ⑥按要求参加课程设计实验演示和答辩等。 3.课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括： ①目录； ②与设计题目相关的理论分析、归纳和总结； ③与设计内容相关的原理分析、建模、推导、可行性分析； ④程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结； ⑤课程设计的心得体会（至少500字）； ⑥参考文献（不少于5篇）； ⑦其它必要内容等。 时间安排：1.5周（分散进行） 参考文献： [1](美)穆尔,高会生,刘童娜,李聪聪．MA TLAB实用教程(第二版) . 电子工业出版社,2010. [2]王正林,刘明．精通MATLAB(升级版) .电子工业出版社,2011. [3]陈杰. MA TLAB宝典(第3版) . 电子工业出版社,2011. [4]刘保柱,苏彦华,张宏林. MATLAB 7.0从入门到精通(修订版) . 人民邮电出版社,2010. 指导教师签名：年月日 系主任（或责任教师）签名：年月日
实验报告 实验名称MATLAB仿真实验 课程名称信号分析与处理 院系部: 专业班级：学生姓名：学号:同组人:实验台号:指导教师：成绩：实验日期:2015-11-29
实验一信号的产生与运算 1.单位阶跃信号 （1）源程序 t=-0.5:0.01:1.5; u=stepfun(t,0); u1=stepfun(t,0.5); figure(1) plot(t,u);axis([-0.5 1.5 -0.2 1.2]);title('单位阶跃信号波形'); figure(2) plot(t,u1);axis([-0.5 1.5 -0.2 1.2]);title('延迟单位阶跃信号波形'); （2）实验结果
2.单位冲激信号 （1）源程序 clear;clc; t=-1:0.001:1; for i=1:3 dt=1/(i^4); x=(1/dt)*((t>=-(1/2*dt))-(t>=(1/2*dt))); subplot(1,3,i); stairs(t,x); end （2）实验结果
3.抽样信号 （1）源程序 clear;clc; t=-20:0.01:20; x=sinc(t/pi); plot(t,x); title('抽样信号'); （2）实验结果
4.单位样值序列（1）源程序 clear;clc; n1=input('n1='); n2=('n2='); n=n1:n2; k=length(n); x1=zeros(1,k); x1(1,-n1+1)=1 subplot(1,2,1); stem(n,x1,'filled') （2）实验结果
自控系统仿真软件课程设计报告 MATLAB 设计题目：牛顿摆球 姓名： 学号： 院系： 班级：1203 指导教师： 2014年12月20日
一．课程设计目的 1、熟悉课程设计的基本流程； 2、掌握MATLAB语法结构及调试方法； 3、熟悉MATLAB函数调用，熟练二维画图; 4、掌握MATLAB语言在控制方面的运用； 5、学会用MATLAB进行基本仿真； 6、掌握MATLAB编程技巧，提高编程水平。 二．系统分析 1.题目的描述： （1）牛顿摆球原理描述 五个质量相同的球体由吊绳固定，彼此紧密排列。当摆动最右侧的球并在回摆时碰撞紧密排列的另外四个球，最左边的球将被弹出，并仅有最左边的球被弹出。当然此过程也是可逆的，当摆动最左侧的球撞击其它球时，最右侧的球会被弹出。当最右侧的两个球同时摆动并撞击其他球时，最左侧的两个球会被弹出。同理相反方向同样可行，并适用于更多的球。 为了更接近现实，在这里我将考虑重力及空气阻力的影响，摆球将不会永无止境的运动下去，由于外界因素的影响，摆球运动一段时间后将回归静止状态。（2）通过MATLAB动画程序制作软件，实现下述过程 当运行程序时，把最右边的小球拉到一定的高度放下，让其碰撞其余四个小球，仅让最左边的小球被弹出，当最左边小球回摆碰撞其它球时，最右边小球又被弹出，如此循环。由于是非理想条件下，摆球的摆动幅度会随摆动次数的增加越来越小，直到静止。 时间停顿两秒，把右边两小球一起拉到一定高度放下，让其碰撞其余三个球，同样仅让左边两球被弹出，当球回摆再次碰撞时，最右边两球又被同时弹出，如此循环，因为外界因素的影响，最终五个球都会静止下来。 （3）整个实验看似简单，但要在MATLAB上完成这样一个动画过程，还是需要下点功夫，克服困难的。经过自己的努力，终于实现了整个过程，这也是一种不小的收获。 2.设计要求： （1）能够实现有阻尼摆动，即摆幅随摆动次数增加越来越小，直到静止。（2）能够让摆球弧线摆动。 三．系统设计 1.系统设计过程 （1）通过函数axis建立坐标系 （2）在坐标系范围内通过函数line画各个支架 （3）通过函数title添加标题“动量守恒实验”、函数text添加标注“牛顿摆球” （4）通过函数line画出五个球，并设定其初始位置，颜色,大小,线条的擦拭方式
《信号与系统》MATLAB实验报告 院系：专业： 年级：班号： 姓名：学号： 实验时间： 实验地点：
实验一 连续时间信号的表示及可视化 实验题目： )()(t t f δ=；)()(t t f ε=；at e t f =)(（分别取00<>a a 及）； )()(t R t f =；)()(t Sa t f ω=；)2()(ft Sin t f π=（分别画出不同周期个数 的波形）。 解题分析： 以上各类连续函数，先运用t = t1: p:t2的命令定义时间范围向量，然后调用对应的函数，建立f 与t 的关系，最后调用plot （）函数绘制图像，并用axis （）函数限制其坐标范围。 实验程序： （1） )()(t t f δ= t=-1:0.01:3 %设定时间变量t 的范围及步长 f=dirac(t) %调用冲激函数dirac （） plot(t,f) %用plot 函数绘制连续函数 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的范围 （2） )()(t t f ε= t=-1:0.01:3 %设定时间变量t 的范围及步长 f=heaviside(t) %调用阶跃函数heaviside （） plot(t,f) %用plot 函数绘制连续函数 title('f(t)=heaviside(t)') %用title 函数设置图形的名称 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的范围 （3） at e t f =)( a=1时： t=-5:0.01:5 %设定时间变量t 的范围及步长 f=exp(t) %调用指数函数exp （）
第23卷第3期湖南理工学院学报(自然科学版)Vol.23 No.3 2010年9月 Journal of Hunan Institute of Science and Technology (Natural Sciences) Sep. 2010信号的频谱分析及MATLAB实现 张登奇, 杨慧银 (湖南理工学院信息与通信工程学院, 湖南岳阳 414006) 摘 要: DFT是在时域和频域上都已离散的傅里叶变换, 适于数值计算且有快速算法, 是利用计算机实现信号频谱分析的常用数学工具. 文章介绍了利用DFT分析信号频谱的基本流程, 重点阐述了频谱分析过程中误差形成的原因及减小分析误差的主要措施, 实例列举了MATLAB环境下频谱分析的实现程序. 通过与理论分析的对比, 解释了利用DFT分析信号频谱时存在的频谱混叠、频谱泄漏及栅栏效应, 并提出了相应的改进方法. 关键词: MA TLAB; 频谱分析; 离散傅里叶变换; 频谱混叠; 频谱泄漏; 栅栏效应 中图分类号: TN911.6 文献标识码: A 文章编号: 1672-5298(2010)03-0029-05 Analysis of Signal Spectrum and Realization Based on MATLAB ZHANG Deng-qi, YANG Hui-yin (College of Information and Communication Engineering, Hunan Institute of Science and Technology, Yueyang 414006, China) Abstract:DFT is a Fourier Transform which is discrete both in time-domain and frequency-domain, it fits numerical calculation and has fast algorithm, so it is a common mathematical tool which can realize signal spectrum analysis with computer. This paper introduces the basic process of signal spectrum analysis with DFT, emphasizes the causes of error producing in spectrum analysis process and the main ways to decrease the analysis error, and lists the programs of spectrum analysis based on MATLAB. Through the comparison with the theory analysis, the problems of spectrum aliasing, spectrum leakage and picket fence effect are explained when using DFT to analyze signal spectrum, and the corresponding solution is presented. Key words:MATLAB; spectrum analysis; DFT; spectrum aliasing; spectrum leakage; picket fence effect 引言 信号的频谱分析就是利用傅里叶分析的方法, 求出与时域描述相对应的频域描述, 从中找出信号频谱的变化规律, 以达到特征提取的目的[1]. 不同信号的傅里叶分析理论与方法, 在有关专业书中都有介绍, 但实际的待分析信号一般没有解析式, 直接利用公式进行傅里叶分析非常困难. DFT是一种时域和频域均离散化的傅里叶变换, 适合数值计算且有快速算法, 是分析信号的有力工具. 本文以连续时间信号为例, 介绍利用DFT分析信号频谱的基本流程, 重点阐述频谱分析过程中可能存在的误差, 实例列出MATLAB 环境下频谱分析的实现程序. 1 分析流程 实际信号一般没有解析表达式, 不能直接利用傅里叶分析公式计算频谱, 虽然可以采用数值积分方法进行频谱分析, 但因数据量大、速度慢而无应用价值. DFT在时域和频域均实现了离散化, 适合数值计算且有快速算法, 是利用计算机分析信号频谱的首选工具. 由于DFT要求信号时域离散且数量有限, 如果是时域连续信号则必须先进行时域采样, 即使是离散信号, 如果序列很长或采样点数太多, 计算机存储和DFT计算都很困难, 通常采用加窗方法截取部分数据进行DFT运算. 对于有限长序列, 因其频谱是连续的, DFT只能描述其有限个频点数据, 故存在所谓栅栏效应. 总之, 用DFT分析实际信号的频谱, 其结果必然是近似的. 即使是对所有离散信号进行DFT变换, 也只能用有限个频谱数据近似表示连续频 收稿日期: 2010-06-09 作者简介: 张登奇(1968? ), 男, 湖南临湘人, 硕士, 湖南理工学院信息与通信工程学院副教授. 主要研究方向: 信号与信息处理
《MATLAB课程设计》报告题目：基于MATLAB的DSB调制与解调分析专业班级: 通信1104班 学生姓名： 指导教师：
MATLAB课程设计任务书 学生姓名：专业班级： 指导教师：工作单位： 题目: 基于MATLAB的DSB调制与解调分析 设计内容和要求 DSB信号的仿真分析 调制信号：分别为300Hz正弦信号和矩形信号；载波频率：30kHz； 解调：同步解调； 要求：画出以下三种情况下调制信号、已调信号、解调信号的波形、频谱以及解调器输入输出信噪比的关系曲线； 1）调制信号幅度=×载波幅度；2）调制信号幅度=载波幅度； 3）调制信号幅度=×载波幅度； 时间安排 2013年12月25日：复习DSB的原理，初步构想设计的流程。 2013年12月26日至28日：程序编写及调试。 2013年12月29日：写报告。 指导教师签名：年月日
目录
摘要 调制在通信系统中有十分重要的作用。通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于传播的已调信号，而且它对系统的传输有效性和传输的可靠性有着很大的影响，调制方式往往决定了一个通信系统的性能。MATLAB软件广泛用于数字信号分析，系统识别，时序分析与建模，神经网络、动态仿真等方面有着广泛的应用。本课题利用MATLAB软件对DSB 调制解调系统进行模拟仿真，分别利用300HZ正弦波和矩形波，对30KHZ正弦波进行调制，观察调制信号、已调信号和解调信号的波形和频谱分布，并在解调时引入高斯白噪声，对解调前后信号进行信噪比的对比分析，估计DSB调制解调系统的性能。 Abstract Modulation in communication systems have an important role. Through the modulation, not only can move the spectrum, the modulated signal spectrum move to the desired position, which will convert into a modulated signal suitable for transmission of modulated signals, and that its transmission system, the effectiveness and reliability of transmission has a great impact, the modulation method is often decided on a communication system performance. MATLAB software is widely used in digital signal analysis, system identification, time series analysis and modeling, neural networks, dynamic simulation have a wide range of applications. This topic using MATLAB software DSB modulation and demodulation system simulation, use, respectively, 300HZ sine wave and rectangular wave, sine wave modulation of the 30KHZ observed modulated signal modulated signal and demodulate the signal waveform and spectrum distribution, and in the solution white Gaussian noise introduced when adjusted for demodulating the signal-noise ratio before and after the comparative analysis, it is estimated DSB modulation and demodulation performance of the system.
第一类：单位转换 1.长度单位换算的设计与实现 2.面积单位换算的设计与实现 3.体积单位换算的设计与实现 4.容积单位换算的设计与实现 5.质量单位换算的设计与实现 6.时间单位换算的设计与实现 7.温度单位换算的设计与实现 7.压强单位换算的设计与实现 8.角度单位换算的设计与实现 8.功率单位换算的设计与实现 第二类：曲线绘制 1.直线的自动绘制和相关计算 2.椭圆的自动绘制和相关计算 3.双曲线的自动绘制和相关计算 4.抛物线的自动绘制和相关计算 5.心脏线的自动绘制和相关计算 6.渐开线的自动绘制和相关计算 7.滚圆线的自动绘制和相关计算 8.三叶玫瑰线的自动绘制和相关计算9.四叶玫瑰线的自动绘制和相关计 10.阿基米德螺线的自动绘制和相关计算第三类：曲面绘制 1.球面的自动绘制和相关计算 2.椭球面的自动绘制和相关计算 3.单叶双曲面的自动绘制和相关计算 4.双叶双曲面的自动绘制和相关计算 5.抛物面的自动绘制和相关计算 6.双曲抛物面的自动绘制和相关计算 7.双曲柱面的自动绘制和相关计算 8.椭圆柱面的自动绘制和相关计算 9.抛物柱面的自动绘制和相关计算 10.圆锥面的自动绘制和相关计算 第四类：线性回归 1.男士身高体重相关计算经验公式 2.女士身高体重相关计算经验公式 3.男士胖瘦等级的确定 4.女士胖瘦等级的确定 5.男士身高脚长相关计算经验公式 6.女士身高脚长相关计算经验公式 7.父子身高相关性研究 8.母子身高相关性研究 9.父女身高相关性研究 10.母女身高相关性研究 第五类：学习成绩 1.期末总评自动计算的设计与实现 2.成绩等级自动评定的设计与实现 3.成绩分段自动统计的设计与实现 4.成绩分布折线自动绘制的设计与实现 5.成绩自动统计分析的设计与实现 6.试卷分布自动分析的设计与实现 7.试卷难度自动分析的设计与实现 8.考试成绩名次自动生成的设计与实现
习题三 绘制典型信号及其频谱图 1.更改参数，调试程序，绘制单边指数信号的波形图和频谱图。观察参数a对信号波形 及其频谱的影响。 程序代码： close all; E=1;a=1; t=0:0.01:4; w=-30:0.01:30; f=E*exp(-a*t); F=1./(a+j*w); plot(t,f);xlabel('t');ylabel('f(t)'); figure; plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|';
E=1，a=1,波形图频谱图更改参数E=2，a=1;
更改参数a，对信号波形及其频谱的影响。（保持E=2）上图为a=1图像 a=2时
a=4时 随着a的增大，f(t)曲线变得越来越陡，更快的逼近0，而对于频谱图，随着a增大，图像渐渐向两边张开，峰值减小，陡度减小，图像整体变得更加平缓。 2.矩形脉冲信号 程序代码： close all; E=1;tao=1; t=-4:0.1:4; w=-30:0.1:30;
f=E*(t>-tao/2&tao/2)+0*(t<=-tao/2&t>=tao/2); F=(2*E./w).*sin(w*tao/2); plot(t,f);xlabel('t');ylabel('f(t)'); figure; plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|') ; figure; plot(w,20*log10(abs(F))); xlabel('\omega');ylabel('|F(\omega)| in dB'); figure; plot(w,angle(F));xlabel('\omega');ylabel('\phi(\omega )');
SI 五款信号完整性仿真工具介绍 （一）Ansoft公司的仿真工具 现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB 设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，An soft公司的仿真工具能够从三维场求解的角度出发，对PCB 设计的信号完整性问题进行动态仿真。 Ansoft 的信号完整性工具采用一个仿真可解决全部设计问题： Slwave是一种创新的工具，它尤其适于解决现在高速PCB和复杂IC封装中普遍存在的电源输送和信号完整性问题。 该工具采用基于混合、全波及有限元技术的新颖方法，它允许工程师们特性化同步开关噪声、电源散射和地散射、谐振、反射以及引线条和电源/地平面之间的耦合。该工具采用一个仿真方案解决整个设计问题，缩短了设计时间。 它可分析复杂的线路设计，该设计由多重、任意形状的电源和接地层，以及任何 数量的过孔和信号引线条构成。仿真结果采用先进的3D 图形方式显示，它还可产生等效电路模型，使商业用户能够长期采用全波技术，而不必一定使用专有仿 （二）SPECCTRAQuest Cade nee的工具采用Sun的电源层分析模块: Cade nee Design System 的SpeeetraQuest PCB信号完整性套件中的电源完整性模块据称能让工程师在高速PCB设计中更好地控制电源层分析和共模EMI 。 该产品是由一份与Sun Microsystems公司签署的开发协议而来的，Sun最初研制该项技术是为了解决母板上的电源问题。 有了这种新模块，用户就可根据系统要求来算出电源层的目标阻抗；然后基于板上的器件考虑去耦合要求，Shah表示，向导程序能帮助用户确定其设计所要求的去耦合电容的数目和类型；选择一组去耦合电容并放置在板上之后，用户就可运行一个仿真程序，通过分析结果来发现问题所在。 SPECCTRAQuest是CADENCE公司提供的高速系统板级设计工具，通过它可以控制与PCB layout相应的限制条件。在SPECCTRAQuest菜单下集成了一下工具： （1）SigXplorer 可以进行走线拓扑结构的编辑。可在工具中定义和控制延时、特性阻抗、驱动和负载的类型和数量、拓扑结构以及终端负载的类型等等。可在
Matalab课后作业 学院：电气信息工程及其自动化 班级： 学号： 姓名： 完成日期： 2012年12月23日
1、 matlab 软件主要功能是什么？电气工程及其自动化专业本科生主要用到哪 些工具箱，各有什么功能？ 答：（1）主要功能：工业研究与开发； 数学教学，特别是线性代数；数值分析和科学计算方面的教学与研究；电子学、控制理论和物理学等工程和科学学科方面的教学与研究； 经济学、化学和生物学等计算问题的所有其他领域中的教学与研究；符号计算功能；优化工具；数据分析和可视化功能；“活”笔记本功能；工具箱；非线性动态系统建模和仿真功能。 （2）常用工具箱： （a ） MATLAB 主工具箱：扩充matlab 的数值计算、符号运算功能、图形建模仿真功能、文字处理功能以及与硬件实时交互功能。 （b ）符号数学工具箱：符号表达式、符号矩阵的创建；符号可变精度求解；因式分解、展开和简化；符号代数方程求解；符号微积分；符号微分方程。 （c ） SIMULINK 仿真工具箱： Simulink 是用于动态系统和嵌入式系统的多领域仿真和基于模型的设计工具。对各种时变系统，包括通讯、控制、信号处理、视频处理和图像处理系统，Simulink 提供了交互式图形化环境和可定制模块库来对其进行设计、仿真、执行和测试。 （d ）信号处理工具箱：数字和模拟滤波器设计、应用及仿真；谱分析和估计；FFT 、DCT 等 变换；参数化模型。 （e ）控制系统工具箱：连续系统设计和离散系统设计；状态空间和传递函数以及模型转换；时域响应（脉冲响应、阶跃响应、斜坡响应）；频域响应（Bode 图、Nyquist 图）；根轨迹、极点配置。 2、设y=23e t 4-sin(43t+3 ),要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。 程序如下： clc clear x=0:0.01:1.5; y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3); y1=diff(y); subplot(2,1,1) plot(x,y) subplot(2,1,2) plot(x(1:150),y1) 曲线如下图所示：
2016-2017学年第一学期 信号与系统实验报告 班级： 姓名： 学号： 成绩： 指导教师：
实验一常见信号的MATLAB 表示及运算 一．实验目的 1．熟悉常见信号的意义、特性及波形 2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二．实验原理 信号一般是随时间而变化的某些物理量。按照自变量的取值是否连续，信号分为连续时间信号和离散时间信号，一般用()f t 和()f k 来表示。若对信号进行时域分析，就需要绘制其波形，如果信号比较复杂，则手工绘制波形就变得很困难，且难以精确。MATLAB 强大的图形处理功能及符号运算功能，为实现信号的可视化及其时域分析提供了强有力的工具。 根据MATLAB 的数值计算功能和符号运算功能，在MATLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MATLAB 语句表示出信号后，就可以利用MATLAB 中的绘图命令绘制出直观的信号波形了。下面分别介绍连续时间信号和离散时间信号的MATLAB 表示及其波形绘制方法。 1.连续时间信号 所谓连续时间信号，是指其自变量的取值是连续的，并且除了若干不连续的点外，对于一切自变量的取值，信号都有确定的值与之对应。从严格意义上讲，MATLAB 并不能处理连续信号。在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t 的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 说明：plot 是常用的绘制连续信号波形的函数。 严格说来，MATLAB 不能表示连续信号，所以，在用plot()命令绘制波形时，要对自变量t 进行取值，MATLAB 会分别计算对应点上的函数值，然后将各个数据点通过折线连接起来绘制图形，从而形成连续的曲线。因此，绘制的只是近似波形，而且，其精度取决于t 的取样间隔。t 的取样间隔越小，即点与点之间的距离越小，则近似程度越好，曲线越光滑。例如：图1-1是在取样间隔为p=0.5时绘制的波形，而图1-2是在取样间隔p=0.1时绘制的波形，两相对照，可以看出图1-2要比图1-1光滑得多。
信号与系统仿真实验报告
实验一 (1)()t δ Function-M 文件 function [x,t]=dirac(t1,t2,t0) %y=dirac(t-t0),t1
分析：在新版的Matlab 函数库中有自带的阶跃函数，调用方法为f=heaviside(t)，这里为了方便画位移后0()t t ε-的图像，故自定义了一个阶跃函数。 (3)指数 ①a=1； >> f=sym('exp(t)'); >> ezplot(f,[-3,3]) >> xlabel('时间t') >> ylabel('函数f （x ）') ②a=-1； f=sym('exp((-1)*t)'); >> ezplot(f,[-3,3]) >> xlabel('时间t') >> ylabel('函数f （x ）') 图a ）a=1时的指数信号图像 图b ）a=-1时的指数函数图像 分析：y=sym （‘f （x ）’）是用了符号运算法 （4）(),5N R t N = >> t=-1:0.001:10; >> y=heaviside(t,0)-heaviside(t,5); >> plot(t,y) >> axis([0,10,-0.2,1.2]) 分析：采用两个跳变点不等的阶跃函数相减得到一个矩形函数的方法生成的门函数。