一. 生成一个基频为20Hz ，幅度为±1的周期方波信号；
（1）试确定该信号的周期？
（2）画出该信号5个周期内的波形；画出该信号的幅度频谱图；
（3）该信号谱线间隔和频带宽度是多少，对它进行时域采样，采样频率Fs 至少为多少？
（4）设计一个IIR/FIR 低通滤波器，滤去该周期信号中80Hz 以后的频率成分，要求通带截止频率80Hz ，阻带截止频率90Hz ，通带纹波小于3dB ，阻带衰减大于50dB ,用matlab 命令画出滤波器的频率响应，以及滤波后的信号波形和频谱。
答：1.周期为0.05秒；
3.谱线间隔为20Hz, 带宽为40Hz ，采样频率至少为80Hz 。
Fs=1000;
N=256; % 时域/频率采样点
t=0:1/Fs:(N-1)/Fs;
x=square(2*pi*20*t); %产生幅度为[-1,1]，频率为20Hz 的方波
subplot(2,1,1) %方波绘图
plot(t,x); title('5个周期的矩形波')
axis([0,0.1,-1.5,1.5])
T=0.05; %求方波的频谱
t1=0:1/Fs:(T-1/Fs);
xn=square(2*pi*20*t1);
Xk=fft(xn,N);
Xk=fftshift(abs(Xk));
Xkl=Xk/N;
subplot(2,1,2)
plot(Fs/N*[0:1:N-1],Xkl)
title('方波信号的频谱')
[n,Wn]=buttord(80/500,90/500,3,50); %确定butterworth 滤波器的阶次 butterworth 低通滤波
%器：通带截止频率80Hz ，阻带截止频率90Hz ，通带纹波小于3dB ，阻带衰减大于50dB ,
%以上参数可以变化
[b,a]=butter(n,Wn); %确定滤波器传递函数
y=filter(b,a,x); %滤波
figure %滤波器频率特性
freqs(b,a); title('低通滤波器频率特性') ;
figure
subplot(2,1,1)
plot(t,y) ;title('矩形波通过低通滤波器后的波形') ;axis([0,0.1,-1.5,1.5])
subplot(2,1,2)
N=512; %可适当增减
fy=fftshift(fft(y,N)); %矩形波通过低通滤波器后的频谱
f=(-N/2:(N/2-1))*2; f=f*Fs/N;
plot(f,abs(fy)/N) ;title('矩形波通过低通滤波器后的频谱')
二.设有连续信号()cos(240)f t t π=??，
（1）利用DFT 计算该连续时间信号的频谱，采样频率至少应取多少？80HZ
（2）以采样频率Hz F s 1000=对信号)(t f 采样得)(t f a ，请画出两个周期内的离散 信号)(t f a 的波形；[要求横
坐标限定到0-0.3s ，纵坐标在-1.4—1.5之间]
（3）若以频率为100Hz 的余弦信号对信号)(t f 进行幅度调制，采样频率仍选取Hz F s 1000=，试画出调制信号()()*cos(2200)y t f t t π=??的波形图和幅度频谱图。（频谱图只显示-1000Hz 到1000HZ 范围，频谱图纵坐标要求是dB 为单位，单独画一幅图）（30分）
（4）对y(t)进行解调：1()()*cos(2200)*cos(2200)y t f t t t ππ=????,设计一个Butter-Worth 低通滤波器，要求从y 1(t)中恢复出原连续时间信号f(t)，确定滤波器的参数，采样频率为1000Hz ；画出滤波器的幅频特性曲线及y 1(t)经过滤波器后的输出信号的时域波形。（要求横坐标是Hz 为单位的频率）。
clear all
clc
Fs=1000;
Fc=200;
N=512;
n=0:1:N-1;
t=n/Fs;
fa=cos(2*pi*20*t);
y=cos(2*pi*200*t).*fa;%调制
stem(n/Fs,fa); %画调制信号
title('采样信号 fa(t)')
axis([0,0.3,-1.4,1.5])
figure %画已调制信号波形和频谱
subplot(2,1,1)
plot(y);title('已调制信号')
subplot(2,1,2)
yk=fftshift(fft(y,N));
f=(-N/2:(N/2-1))*2;
f=f*Fs/N;
plot(f,abs(yk))
title('调制信号的频谱')
[n,Wn]=buttord(60/500,100/500,3,50);%确定butterworth 滤波器阶次
[b,a]=butter(n,Wn);%确定滤波器传递函数
figure
freqs(b,a);%滤波器频率特性
title('低通滤波器频率特性')
y1=y.*cos(2*pi*200*t) %相干解调
y=filter(b,a,y1); %滤波
figure
plot(t,y);
title('调制信号y （t ）通过滤波器后的波形')
三.已知信号)1002sin()152sin()102sin()(t t t t x ?+?+?=πππ，
（1）若对信号进行数字处理，采样频率理论上至少为多少Hz ？200Hz
（2）该信号的频率分辨率是多少Hz ？在500Hz 采样频率下，若想区分三个正弦分量的频谱，采样点数至少为多少？128
（3）在500Hz 采样频率绘制信号的时域波形，绘出信号的幅度频谱图。
（4）设计一个IIR 或FIR 滤波器，使)(t x 经过滤波器后，滤除)200
sin(t π。画出滤波器的频率响应特性图，画出滤波器的输出信号时域波形图。
clear all
clc
fs=500;
n=[0:127];
x=cos(2*pi*10*n/fs)+cos(2*pi*15*n/fs)+cos(2*pi*100*n/fs);
x1=sin(200*pi*n/fs);
subplot(2,1,1);plot(n/fs,x);
title('signal x(t)')
Xk=fft(x);
Xk=fftshift(abs(Xk));
subplot(2,1,2);
plot(n/fs,Xk);
title('信号幅度频谱图');
[n,Wn]=buttord(10/500,100/500,3,50);%确定butterworth 滤波器阶次
[b,a]=butter(n,Wn);%确定滤波器传递函数
figure
freqs(b,a);%滤波器频率特性
title('低通滤波器频率特性')
y=filter(b,a,x1);%滤波
figure
n=[0:127];
fs=500;
plot(n/fs,y)
title('x1通过低通滤波器后的波形')
四,设有连续信号
(t)cos(2 5.51000)cos(261000)cos(291000)a x t t t πππ=??+??+??
(1) 若利用DFT 计算该连续时间信号频谱,采样频率应该取多少?18000HZ
(2) 为区分上述的三个频率分量,至少应该采集多少个样本?128
(3) 一采样频率Fs=32kHz 对信号Xa(t)采样,画出Xa(t)的波形
(4) 一采样频率Fs=32kHz 对信号Xa(t)采样,采样数据长度N=128点,做N=128点的DFT,并画出幅频特性.
(5) 设计一个滤波器,滤除信号Xa(t)中频率为6kHz 和9kHz 的两个分量,画出滤波器的幅频特性曲线
(6) 画出滤波器的输出信号的时域波形图.
fs=32000;
n=[0:127];
x=cos(2*pi*5.5*1000*n/fs)+cos(2*pi*6*1000*n/fs)+cos(2*pi*9*1000*n/fs);
subplot(2,1,1);plot(n/fs,x);title('signal x(t)')
Xk=fft(x,128);
magXk=abs(Xk);
k1=0:127;w1=fs/256*k1;
subplot(2,1,2);plot(w1,magXk);title('16点DTFT的幅值');
xlabel('w(单位：Hz)')
[n,Wn]=buttord(5500/16000,6000/16000,3,50);%确定butterworth滤波器阶次[b,a]=butter(n,Wn);%确定滤波器传递函数
figure
freqs(b,a);%滤波器频率特性
y=filter(b,a,x);%滤波
figure
n=[0:127];
fs=32000;
plot(n/fs,y)
title('xa(t)通过低通滤波器后的波形')
五,产生一个周期为0.01秒,幅值为-1到1的方波信号,画出信号的时域波形图;
(1)该信号的谱线间隔是多少Hz?/频带宽度是多少Hz?
(2)若利用DFT计算该连续时间信号的频谱,采样频率至少应取多少?
(3)用Matlab画出该信号的幅度频谱图;
(4)设计一个滤波器,滤除该信号的基波分量;画出滤波器的频响特性;
(5)画出滤波器的输出信号的时域波形;
答(1)谱线间隔是100Hz;频带宽度是200Hz’
(2)采样频率至少应取400Hz
clear all
clc
fs=1000;
t=0:1/fs:0.1;
x=square(2*pi*100*t);%产生幅度为[-1，+1]的方波
plot(t,x)
title('周期矩形波')
axis([0,0.1,-1.5,1.5])
%求周期矩形波的幅度频谱
T=0.01;
t1=0:1/fs:(T-1/fs);
x1=square(2*pi*100*t1);
Xk=fft(x1);
Xk=(1/10)*(fftshift(abs(Xk)));
figure
stem((-fs/10)*(-5:4),Xk)
title('周期矩形波的幅度频谱图')
[n,Wn]=buttord(150/500,200/500,3,50);%确定butterworth滤波器阶次
[b,a]=butter(n,Wn);%确定滤波器传递函数
freqs(b,a);%滤波器频率特性
title('低通滤波器频率特性')
y=filter(b,a,x);%滤波
figure
plot(t,y)
title('矩形波通过低通滤波器后的波形')
axis([0,0.1,-1.5,1,5])
六
1.生成一个基频为10Hz，幅度为0到1，占空比为50%的周期矩形方波信号，画出该信号3个周期内的波形（要
求横坐标限定在0到0.4秒之间，纵坐标限定在-0.3到+1.3之间）；
2.选取合适的采样频率，画出该信号的幅度频谱图，只显示-500Hz到+500Hz范围；
3.设计一个带阻滤波器，要求该滤波器滤除掉频率为30Hz的三次谐波成分，用MATLAB命令只画出滤波器的幅度
频谱响应特性曲线，要求纵坐标单位为dB。
解：
clc
clear all
fs=1000;
N=256;
t=0:1/fs:0.3;
x=(square(2*pi*10*t)+1)/2; %产生方波信号
subplot(2,1,1)
plot(t,x);title('周期方波信号')
axis([0,0.4,-0.3,1.3])
T=0.1; %求方波的频谱
t1=0:1/fs:(T-1/fs);
Xk=fft(x,N);
Xk=fftshift(abs(Xk));
Xk1=Xk/N;
f=(-N/2:(N/2-1));
f=f*fs/N;
subplot(2,1,2)
plot(f,Xk1)
title('方波信号的频谱')
%设计带阻滤波器
wp=[30,35]*2/fs;
ws=[20,40]*2/fs;
Rp=3;
Rs=50;
[N,wn]=buttord(wp,ws,Rp,Rs);
[b,a]=butter(N,wn,'stop');
figure
freqs(b,a);
title('低通滤波器频率特性')
[H,W]=freqz(b,a);
n=0:511;
%plot((fs/2)/512*n,abs(H));grid on
plot((fs/2)/512*n,10*log10(abs(H)));grid on
title('滤波器幅度频率响应曲线')
七设有连续信号()cos(220)f t t π=??,
1. 以Fs=800Hz 对信号()f t 采样得离散信号()a f t ，画出()a f t 4个周期内的离散波形（要求横坐标限定在0到0.3
秒之间，纵坐标限定在-1.4到1.5之间）；
2. 若以频率为100Hz 的余弦信号对信号()f t 进行幅度调制，选取合适的采样频率Fs ，试画出调制信号
1()()cos(2100)y t f t t π=???的连续时域波形图和幅度频谱图（频谱图只显示-800Hz 到+800Hz 范围，频谱图纵坐标要求是dB 为单位，单独画一幅图）。
3. 对1()y t 进行解调：2()()cos(2100)cos(2100)y t f t t t ππ=??????，设计一个Butter-Worth 带通滤波器，
要求从2()y t 中恢复出原连续时间信号()f t ，确定滤波器的参数，采样频率1600Hz ；只画出滤波器的幅度频谱图（要求横坐标是Hz 为单位的频率）。
解：
%第一题
clc
clear
fs=800;
t=0:1/fs:0.2;
f=cos(2*pi*20*t);
figure
stem(t,f);
axis([0,0.3,-1.4,1.5])
title('离散信号波形')
xlabel('t');ylabel('f(t)')
%第二题
fs=1600;
fc=100;
t=0:1/fs:0.2;
f=cos(2*pi*20*t);
y1=modulate(f,fc,fs,'Amdsb-sc');
figure
plot(t,y1)
title('调制信号波形')
xlabel('t/s')
ylabel('y(t)/V')
Y1=fftshift(abs(fft(y1,512)));
w=-256:255;
w=w*fs/512;
figure
plot(w,Y1);
plot(w,10*log10(Y1)); %以dB 为单位
title('调制信号频谱图')
xlabel('w/Hz');ylabel('Y(jw)/dB')
%设计带通滤波器
fs=1600;
wp=[20,25]*2/fs;
ws=[10,30]*2/fs;
Rp=3;Rs=45;
[N,wn]=buttord(wp,ws,Rp,Rs);
[b,a]=butter(N,wn,'bandpass');
figure
freqz(b,a);
title('带通滤波器的频谱响应')
y1=f.*cos(2*pi*200*t) %相干解调
f=filter(b,a,y1); %滤波
figure
plot(t,f);
title('调制信号y（t）通过滤波器后的波形')
八
的周期方波信号；画出该信号五个周期的时域波形；
1.生成一个周期为0.02s，幅度为1
2.画出该信号的幅度频谱图；
3.设计一个butterworth低通滤波器，滤波器的参数自行确定，滤去该周期信号中50Hz以后的频率成分，用matlab
命令画出滤波器的频率响应特性曲线。
4.画出该信号经过滤波后的输出信号的波形和频谱图。
解：
%第一题
clc
clear
fs=1000;
t=0:1/fs:0.1;
x=square(2*pi*50*t); %产生方波信号
figure
plot(t,x)
title('时域波形')
axis([0,0.1,-1.2,1.2])
%第二题
N=512;
n=0:N-1;
X=fftshift(fft(x,512));
figure
plot(-fs/2+fs/512*n,abs(X/N))
title('幅度频谱图')
%低通滤波器
fs=200;
wp=50*2/fs;
ws=45*2/fs;
Rp=3;Rs=45;
[N,wn]=buttord(wp,ws,Rp,Rs);
[b,a]=butter(N,wn,'low');
figure
[H,W]=freqz(b,a);
k=0:511;
plot((fs/2)/512*k,abs(H));grid
title('低通滤波器的频率响应')
%滤波
N=4*(fs/10);
y=filter(b,a,x);
figure
plot(t,y)
title('滤波后的波形')
Y=fftshift(fft(y,512));
yk=Y/N;
figure
plot(-fs/2+fs/512*k,abs(yk))
title('滤波后的频谱')
产生离散衰减正弦序列()π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) ;
实验一 基本信号的产生与运算 一、 实验目的 学习使用MATLAB 产生基本信号、绘制信号波形、实现信号的基本运算。 二、 实验原理 MATLAB 提供了许多函数用于产生常用的基本信号：如阶跃信号、脉冲信号、指数信号、正弦信号和周期方波等等。这些信号是信号处理的基础。 1、 利用MATLAB 产生下列连续信号并作图。 （1）51),1(2)(<<---=t t u t x （2）300),3 2sin()(3.0<<=-t t e t x t （3）1.01.0,3000cos 100cos )(<<-+=t t t t x （4）2000),8.0cos()1.0cos()(<<=t t t t x ππ 答：（1）、 >> t=-1:0.02:5; >> x=(t>1); >> plot(t,-2*x); >> axis([-1,5,-3,1]); >> title('杨婕婕 朱艺星'); >> xlabel('x(t)=-2u(t-1)'); （2）、 >> t=0:0.02:30; >> x=exp(-0.3*t).*sin(2/3*t); >> plot(t,x);
>> title('杨婕婕朱艺星'); >> xlabel('x(t)=exp(-0.3*t).*sin(2/3*t)'); 因为原函数在t=15后x(t)取值接近于零，所以将横坐标改成0到15，看得更清晰 axis([0,15,-0.2,0.6]); （3）>> t=-0.1:0.01:0.1;x=cos(100*t)+cos(3000*t);plot(t,x); >> title('杨婕婕朱艺星'); >>xlabel('x=cos(100*t)+cos(3000*t)'); 因为t的间隔取太大，以至于函数不够准确，缩小t的间隔： t=-0.1:0.002:0.2;x=cos(100*t)+cos(3000*t); plot(t,x);title('杨婕婕') >> t=-0.1:0.0001:0.1; x=cos(100*t)+cos(3000*t); >> plot(t,x);title('杨婕婕朱艺星'); >> xlabel('x=cos(100*t)+cos(3000*t)');
《信号与系统》实验指导书 张建奇骆崇编写 浙江工业大学之江学院信息工程分院 2012年2月
目录 实验一MATLAB的基本使用 (1) 一、实验目的 (1) 二、实验原理 (1) 三、实验内容与要求 (8) 四、实验报告 (9) 实验二时域波形的MATLAB实现 (10) 一、实验目的 (10) 二、预习要求 (10) 三、实验原理 (10) 四、实验内容与要求 (18) 五、实验报告 (19) 实验三用MATLAB对系统时域分析 (20) 一、实验目的 (20) 二、预习要求 (20) 三、实验原理 (20) 四、实验内容与要求 (29)
实验一MATLAB的基本使用 一、实验目的 1、了解和掌握MATLAB的基本操作 2、了解MATLAB的库函数 3、会用MATLAB进行简单的操作。 二、实验原理 1、界面操作 MATLAB是“MATrix LABoratory”的缩写（矩阵实验室），它是由美国Mathworks公司于1984年正式推出的一种科学计算软件，由于其强大的功能，在欧美的一些大学里MATLAB已经成为许多诸如数字信号处理、自动控制理论等高级教程的主要工具软件，同时也成为理工科学生，必须掌握的一项基本技能。 当需要运行程序时，只需选择桌面上（或开始）中的MATLAB6.5应用程序图标即可 通常情况下，MATLAB的工作环境主要由一下几个窗口组成： 命令窗口（Command Window）
工作区间浏览器（Workspace） 历史命令窗口（Command History） 图形窗口(Figure) 文本编辑窗口(Editor) 当前路径窗口（Current Directory） MATLAB的命令窗与命令操作 当用户使用命令窗口进行工作时，在命令窗口中可以直接输入相应的命令，系统将自动显示信息。 例如在命令输入提示符“>>”后输入指令： >>t=[1,2,3;4,5,6;7,8,9]; 按回车键(Enter)后，系统即可完成对变量t的赋值。 MATALB提供了非常方便的在线帮助命令(help)，它可提供各个函数的用法指南，包括格式、参数说明、注意事项及相关函数等内容。 2、图形窗 MATLAB图形窗(Figure)主要用于显示用户所绘制的图形。 通常，只要执行了任意一种绘图命令，图形窗就会自动产生。
一、设计目的 通过运用MATLAB对函数进行Z域分析和单边带信号的调制与解调，使我们进一步加深对MATLAB的认识和运用，以实现以下目的： 1.本次试验进一步熟悉了MATLAB软件的使用方法及相关的操作。 2.对Z变换及其反变换函数在MATLAB中的调用有了掌握。 3.理论与实际的仿真相结合，更直观的看到结果。 4.观察了单边带信号调制与解调后的图像，加深认识。 二、设计原理 MATLAB是The MathWorks公司在1984年推出的一种商品化软件，它提供了大量丰富的应用函数，并且具有扩充的开放性结构。目前，该软件包涵盖了控制系统应用、数字信号处理、数字图像处理、通讯、神经网络、小波理论分析、优化与统计、偏微分方程、动态系统实时仿真等多学科专业领域。 其中单边带调制信号是将双边带信号中的一个边带滤掉而形成的。根据方法的不同，产生单边带调制信号的方法有：滤波和相移法。 由于滤波法在技术上比较难实现所以在此我们将用相移法对单边带调制与解调系统进行讨论与设计。 三、设计内容和MATLAB图像
1、数字系统的响应 源代码如下： b=[0 1 2 1 0]; a=[1 -0.5 0 0.3 -0.005]; subplot(421);zplane(b,a); title('系统的零极点图'); subplot(422);impz(b,a,21); title('单位脉冲响应'); subplot(423);stepz(b,a,21); title('单位阶跃响应');
N=21;n=0:N-1; x=exp(-n); x0=zeros(1,N); y0=[1,-1]; xi=filtic(b,a,y0); y1=filter(b,a,x0,xi); xi0=filtic(b,a,0); y2=filter(b,a,x,xi0); y3=filter(b,a,x,xi); [h w]=freqz(b,a,21); subplot(424);stem(n,y1); title('零输入响应');grid on; subplot(425);stem(n,y2); title('零状态响应');grid on; subplot(426);stem(n,y3); title('系统的全响应');grid on; subplot(427);plot(w,abs(h)); title('幅频特性曲线');grid on; subplot(428);plot(w,angle(h)); title('相频特性曲线');grid on;
《信号与系统》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 （） plot(t,f) %用plot 函数绘制连续函数 title('f=exp(t)') %用title 函数设置图形的名称 axis([-5,5,-1,100]) %用axis 函数规定横纵坐标的围 a=2时： t=-5:0.01:5 f=exp(2*t) % 调用指数函数exp （） plot(t,f) title('f=exp(2*t)') axis([-5,5,-1,100]) a=-2时： t=-5:0.01:5 f=exp(-2*t) plot(t,f) title('f=exp(-2*t)') axis([-5,5,-1,100]) （4）)()(t R t f = t=-5:0.01:5 f=rectpuls(t,2) % 用rectpuls(t,a)表示门函数，默认以零点为中心，宽度为a plot(t,f) title('f=R(t)') axis([-5 5 -0.5 1.5]) （5）)()(t Sa t f ω= ω=1时： t=-20:0.01:20 f=sin(t)./t % 调用正弦函数sin （）,并用sin （t ）./t 实现抽 样函数 plot(t,f)
实验1 用MATLAB 进行信号频谱分析 一、实验目的 ㈠ 初步掌握MATLAB 产生常用离散时间信号的编程方法。 ㈡ 学习编写简单的FFT 算法程序，对离散信号进行幅频谱分析。 ㈢ 观察离散时间信号频谱的特点。 二、实验原理 ㈠ 常用的离散时间信号 在 MATLAB 语言主要是研究离散信号的。常用的离散信号有： 1．单位取样序列 ???≠==000 1)(n n n δ 2．单位阶跃序列 ?? ?<≥=0 01 )(n n n u 3．实指数序列 R a n a n x n ∈?=;)( 4．复指数序列 n e n x n j ?=+)(0)(ωσ 5．正(余)弦序列 )c o s ()(0θω+=n n x n ? 6．周期序列 n N n x n x ?+=)()( ㈡ 离散信号的产生 离散信号的图形显示使用stem 指令。 在 MATLAB 中的信号处理工具箱中,主要提供的信号是离散信号。
由于MATLAB 对下标的约定为从1开始递增，例如x=[5,4,3,2,1,0],表示x(1)=5,x(2)=4,X(3)=3… 因此要表示一个下标不由1开始的数组x(n)，一般应采用两个矢量，如 n=[-3，-2，-1，0，l ，2，3，4，5]； x=[1，-l ，3，2，0，4，5，2，1]； 这表示了一个含9个采样点的矢量：X(n)={x(-3)，x(-2)，x(-1)，x(0)，x(1)，x(2)，x(3)，x(4)，x(5)}。 1．单位取样序列 ?? ?≠==δ0 001)(n n n 这一函数实现的方法有二： 方法一：可利用MATLAB 的zeros 函数。 x=zeros(1，N)； %建立一个一行N 列的全零数组 x(1)=1； %对X （1）赋1 方法二：可借助于关系操作符实现 n=1:N; x=[n==1]; %n 等于1时逻辑关系式结果为真，x=1；n 不等于1时为假，x=0 如要产生 ?????≤<<=≤≤=-δ2 0210 100)(10)(n n n n n n n n n n n n 则可采用MATLAB 实现： n=n1:n2； x=[(n-n0)==0]；%n=n0时逻辑关系式结果为真，x=1；n ≠n0时为
第一章课程设计内容及要求 1. 直流电动机的机械特性仿真； 2. 直流电动机的直接起动仿真； 3. 直流电动机电枢串联电阻启动仿真； 4. 直流电动机能耗制动仿真； 5．直流电动机反接制动仿真； 6. 直流电动机改变电枢电压调速仿真； 7. 直流电动机改变励磁电流调速仿真。 要求：编写M文件，在Simulink环境画仿真模型原理图，用二维画图命令画仿真结果图或用示波器观察仿真结果，并加以分析
第二章直流电动机的电力拖动仿真绘制 1）直流电动机的机械特性仿真 clear; U_N=220;P_N=22;I_N=115; n_N=1500;R_a=;R_f=628; Ia_N=I_N-U_N/R_f; C_EPhi_N=(U_N-R_a*Ia_N)/n_N; C_TPhi_N=*C_EPhi_N; Ia=0;Ia_N; n=U_N/C_EPhi_N-R_a/(C_EPhi_N)*Ia; Te=C_TPhi_N*Ia; P1=U_N*Ia+U_N*U_N/R_f; T2_N=9550*P_N/n_N; figure(1); plot(Te,n,'.-'); xlabel('电磁转矩Te/'); ylabel('转矩n/rpm'); ylim([0,1800]); figure(2); plot(Te,n,'rs'); xlabel('电磁转矩Te/'); ylabel('转矩n/rpm');
hold on; R_c=0; for coef=1:;; U=U_N*coef; n=U/C_EPhi_N-(R_a+R_c)/(C_EPhi_N*C_TPhi_N)*Te; plot(Te,n,'k-'); str=strcat('U=',num2str(U),'V'); s_y=1650*coef; text(50,s_y,str); end figure(3); n=U_N/C_EPhi_N-(R_a+R_c)/(C_EPhi_N*C_TPhi_N)*Te; plot(Te,n,'rs'); xlabel('电磁转矩Te/'); ylabel('转矩n/rpm'); hold on; U=U_N;R_c=; for R_c=0::; n=U/C_EPhi_N-(R_a+R_c)/(C_EPhi_N*C_TPhi_N)*Te; plot(Te,n,'k-'); str=strcat('R=',num2str(R_c+R_a),'\Omega'); s_y=400*(4-R_c*; text(120,s_y,str);
习题三 绘制典型信号及其频谱图 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 )');
实验一 基本信号的产生与运算 一、 实验目的 学习使用MATLAB 产生基本信号、绘制信号波形、实现信号的基本运算。 二、 实验原理 MATLAB 提供了许多函数用于产生常用的基本信号：如阶跃信号、脉冲信号、指数信号、正弦信号和周期方波等等。这些信号是信号处理的基础。 1、 利用MATLAB 产生下列连续信号并作图。 （1）51),1(2)(<<---=t t u t x （2）300),3 2 sin()(3.0<<=-t t e t x t （3）1.01.0,3000cos 100cos )(<<-+=t t t t x （4）2000),8.0cos()1.0cos()(<<=t t t t x ππ 答：（1）、 >> t=-1:0.02:5; >> x=(t>1); >> plot(t,-2*x); >> axis([-1,5,-3,1]); >> title('杨婕婕 朱艺星'); >> xlabel('x(t)=-2u(t-1)');
（2）、 >> t=0:0.02:30; >> x=exp(-0.3*t).*sin(2/3*t); >> plot(t,x); >> title('杨婕婕朱艺星'); >> xlabel('x(t)=exp(-0.3*t).*sin(2/3*t)');
因为原函数在t=15后x(t)取值接近于零，所以将横坐标改成0到15，看得更清晰 axis([0,15,-0.2,0.6]);
（3）>> t=-0.1:0.01:0.1;x=cos(100*t)+cos(3000*t);plot(t,x); >> title('杨婕婕朱艺星'); >>xlabel('x=cos(100*t)+cos(3000*t)'); 因为t的间隔取太大，以至于函数不够准确，缩小t的间隔：t=-0.1:0.002:0.2;x=cos(100*t)+cos(3000*t); plot(t,x);title('杨婕婕')
MATLAB基础与应用简明教程 张明等编著 北京航空航天大学出版社(2001.01) MATLAB软件环境是美国New Mexico大学的Cleve Moler博士首创的，全名为MATrix LABoratory（矩阵实验室）。它建立在20世纪七八十年代流行的LINPACK（线性代数计算）和ESPACK（特征值计算）软件包的基础上。LINPACK和ESPACK软件包是从Fortran语言开始编写的，后来改写为C语言，改造过程中较为复杂，使用不便。MA TLAB是随着Windows环境的发展而迅速发展起来的。它充分利用了Windows环境下的交互性、多任务功能语言，使得矩阵计算、数值运算变得极为简单。MA TLAB语言是一种更为抽象的高级计算机语言，既有与C语言等同的一面，又更为接近人的抽象思维，便于学习和编程。同时，它具有很好的开放性，用户可以根据自己的需求，利用MA TLAB提供的基本工具，灵活地编制和开发自己的程序，开创新的应用。 本书重点介绍了MA TLAB的矩阵运算、符号运算、图形功能、控制系统分析与设计、SimuLink仿真等方面的内容。 Chap1 MATLAB入门与基本运算 本章介绍MATLAB的基本概念，包括工作空间；目录、路径和文件的管理方式；帮助和例题演示功能等。重点介绍矩阵、数组和函数的运算规则、命令形式，并列举了可能得到的结果。由于MA TLAB的符号工具箱是一个重要分支，其强大的运算功能在科技领域有特殊的帮助作用。 1.1 MATLAB环境与文件管理 1.2 工作空间与变量管理 1.2.1 建立数据 x1=[0.2 1.11 3]; y1=[1 2 3;4 5 6]建立一维数组x1和二维矩阵y1。分号“；”表示不显示定义的数据。 MATLAB还提供了一些简洁方式，能有规律地产生数组： xx=1:10 %xx从1到10，间隔为1 xx=-2:0.5:1 %xx从-2到1，间隔为0.5 linespace命令等距离产生数组，logspace在对数空间中等距离产生数组。对于这一类命令，只要给出数组的两端数据和维数就可以了。 xx=linespace(d1,d2,n) %表示xx从d1到d2等距离取n个点 xx=logspace(d1,d2,n) %表明xx从10d1到10d2等距离取n个点 1.2.2 who和whos命令 who: 查看工作空间中有哪些变量名 whos: 了解这些变量的具体细节 1.2.3 exist命令 查询当前的工作空间内是否存在一个变量，可以调用exist()函数来完成。 调用格式：i=exist(…A?); 式中，A为要查询的变量名。返回的值i表示A存在的形式： i=1 表示当前工作空间内存在一个变量名为A的矩阵； i=2 表示存在一个名为A.m的文件； i=3 表示MATLAB的工作路径下存在一个名为A.mex的文件；
实验一 名称：连续时间信号分析 姓名：王嘉琦 学号：201300800636 班级：通信二班 一、实验目的 （一）掌握使用Matlab 表示连续时间信号 1、学会运用Matlab 表示常用连续时间信号的方法 2、观察并熟悉常用信号的波形和特性 （二）掌握使用Matlab 进行连续时间信号的相关运算 1、学会运用Matlab 进行连续时间信号的时移、反褶和尺度变换 2、学会运用Matlab 进行连续时间信号微分、积分运算 3、学会运用Matlab 进行连续时间信号相加、相乘运算 4、学会运用Matlab 进行连续时间信号卷积运算 二、实验条件 Matlab 三、实验内容 1、利用Matlab 命令画出下列连续信号的波形图。 （1）)4/3t (2cos π+ 代码： k=2;w=3;phi=pi/4; t=0:0.01:3; ft=k*cos(w*t+phi); plot(t,ft),grid on; axis([0,3,-2.2,2.2]) title('余弦信号')
（2） )t (u )e 2(t -- 代码： k=-1;a=-1; t=0:0.01:3; ft=2-k*exp(a*t); plot(t,ft),grid on axis([0,3,2,3]) title('指数信号')
（3）)]2()(u )][t (cos 1[--+t u t π 代码： k=1;w=pi;phi=0; t=0:0.01:2; ft=1+k*cos(w*t+phi); plot(t,ft),grid on; axis([0,3,0,2]) title('余弦信号')
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光滑得多。
S i m u l i n k仿真环境基础学习Simulink是面向框图的仿真软件。 7.1演示一个Simulink的简单程序 【例7.1】创建一个正弦信号的仿真模型。 步骤如下： (1)在MATLAB的命令窗口运行simulink命令，或单击工具栏中的图标，就可以打开Simulink模块库浏览器(SimulinkLibraryBrowser)窗口，如图7.1所示。
图7.1Simulink界面 (2)单击工具栏上的图标或选择菜单“File”——“New”——“Model”，新建一个名为“untitled”的空白模型窗口。 (3)在上图的右侧子模块窗口中，单击“Source”子模块库前的“+”(或双击Source)，或者直接在左侧模块和工具箱栏单击Simulink下的Source子模块库，便可看到各种输入源模块。 (4)用鼠标单击所需要的输入信号源模块“SineWave”(正弦信号)，将其拖放到的空白模型窗口“untitled”，则“SineWave”模块就被添加到untitled窗口；也可以用鼠标选中“SineWave”模块，单击鼠标右键，在快捷菜单中选择“addto'untitled'”命令，就可以将“SineWave”模块添加到untitled窗口，如图7.2所示。
(5) Scope ”模块(示波器)拖放到“untitled ”窗口中。 (6)在“untitled ”窗口中，用鼠标指向“SineWave ”右侧的输出端，当光标变为十字符时，按住鼠标拖向“Scope ”模块的输入端，松开鼠标按键，就完成了两个模块间的信号线连接，一个简单模型已经建成。如图7.3所示。 (7)开始仿真，单击“untitled ”模型窗口中“开始仿真”图标 ，或者选择菜单“Simulink ”——“Start ”，则仿真开始。双击“Scope ” 模块出现示波器显示屏，可以看到黄色的正弦波形。如图7.4所示。 图7.2Simulink 界面
MATLAB在电力电子技术中的应用 目录 MATLAB在电力电子技术中的应用 (1) MATLAB in power electronics application (2) 目录 (4) 1绪论 (6) 1.1关于MATLAB软件 (6) 1.1.1MATLAB软件是什么 (6) 1.1.2MATLAB软件的特点和基本操作窗口 (7) 1.1.3MATLAB软件的基本操作方法 (10) 1.2电力电子技术 (12) 1.3MATLAB和电力电子技术 (13) 1.4本文完成的主要内容 (14) 2MATLAB软件在电路中的应用 (15) 2.1基本电气元件 (15) 2.1.1基本电气元件简介 (15) 2.1.2如何调用基本电器元件功能模块 (17) 2.2如何简化电路的仿真模型 (19) 2.3基本电路设计方法 (19) 2.3.1电源功能模块 (19) 2.3.2典型电路设计方法 (20) 2.4常用电路设计法 (21) 2.4.1ELEMENTS模块库 (21) 2.4.2POWER ELECTRONICS模块库 (22) 2.5MATLAB中电路的数学描述法 (22) 3电力电子变流的仿真 (25) 3.1实验的意义 (25) 3.2交流-直流变流器 (25)
3.2.1单相桥式全控整流电路仿真 (26) 3.2.2三相桥式全控整流电路仿真 (38) 3.3三相交流调压器 (53) 3.3.1无中线星形联结三相交流调压器 (53) 3.3.2支路控制三角形联结三相交流调压器 (59) 3.4交流-交流变频电路仿真 (64) 3.5矩阵式整流器的仿真 (67)
实验一、MATLAB 编程基础及典型实例 一、实验目的 （1）熟悉MATLAB 软件平台的使用； （2）熟悉MATLAB 编程方法及常用语句； （3）掌握MATLAB 的可视化绘图技术； （4）结合《信号与系统》的特点，编程实现常用信号及其运算。 示例一：在两个信号进行加、减、相乘运算时，参于运算的两个向量要有相同的维数，并且它们的时间变量范围要相同，即要对齐。编制一个函数型m 文件，实现这个功能。function [f1_new,f2_new,n]=duiqi(f1,n1,f2,n2) a=min(min(n1),min(n2)); b=max(max(n1),max(n2)); n=a:b; f1_new=zeros(1,length(n)); f2_new=zeros(1,length(n)); tem1=find((n>=min(n1))&(n<=max(n1))==1); f1_new(tem1)=f1; tem2=find((n>=min(n2))&(n<=max(n2))==1); f2_new(tem2)=f2; 四、实验内容与步骤 （2）绘制信号x(t)=)3 2sin(2t e t ?的曲线，t 的范围在0~30s ，取样时间间隔为0.1s 。t=0:0.1:30; y=exp(-sqrt(2)*t).*sin(2*t/3); plot(t,y);
（3）在n=[-10:10]范围产生离散序列：?? ?≤≤?=Other n n n x ,033,2)(,并绘图。n=-10:1:10; z1=((n+3)>=0); z2=((n-3)>=0); x=2*n.*(z1-z2); stem(n,x);（4）编程实现如下图所示的波形。 t=-2:0.001:3; f1=((t>=-1)&(t<=1)); f2=((t>=-1)&(t<=2)); f=f1+f2; plot(t,f); axis([-2,3,0,3]);
基带信号眼图实验——matlａb 仿真
————————————————————————————————作者：————————————————————————————————日期： ?
数字基带信号的眼图实验——matla ｂ仿真 一、实验目的 1、掌握无码间干扰传输的基本条件和原理，掌握基带升余弦滚降系统的实现方法； 2、通过观察眼图来分析码间干扰对系统性能的影响，并观察在输入相同码率的ＮRZ 基带信号下，不同滤波器带宽对输出信号码间干扰大小的影响程度； 3、熟悉MATL ＡB 语言编程。 二、实验预习要求 1、复习《数字通信原理》第七章７.1节——奈奎斯特第一准则内容； 2、复习《数字通信原理》第七章7．2节——数字基带信号码型内容； 3、认真阅读本实验内容,熟悉实验步骤。 三、实验原理和电路说明 １、基带传输特性 基带系统的分析模型如图3-1所示,要获得良好的基带传输系统，就应该 () n s n a t nT δ-∑() H ω() n s n a h t nT -∑基带传输抽样判决 图3-１?基带系统的分析模型 抑制码间干扰。设输入的基带信号为()n s n a t nT δ-∑，s T 为基带信号的码元周期,则经过基 带传输系统后的输出码元为 ()n s n a h t nT -∑。其中 1 ()()2j t h t H e d ωωωπ +∞ -∞ = ? ?（3-1） 理论上要达到无码间干扰,依照奈奎斯特第一准则,基带传输系统在时域应满足： 10()0,s k h kT k =?=? ? ， 为其他整数 ?? ?(３－２) 频域应满足：
MATLAB的实验仿真 目录 实验一MATLAB在控制系统模型建立与仿真中的应用 (1) 实验二典型系统的时域响应分析 (13) 实验三线性控制系统的根轨迹与频域分析 (17) 实验四线性系统的校正 (22) 附录一 MATLAB6.5 控制系统工具箱函数和结构化的控制语句 (30) 附录二 SIMULINK 基本模块介绍 (34)
实验一 MATLAB 在控制系统模型建立与仿真中的应用 一、 MATLAB 基本操作与使用 1. 实验目的 1) 掌握MATLAB 仿真软件的安装及启动，熟悉 MATLAB 工作环境平台。 2) MATLAB 命令窗口，包括工具条以及菜单选项的使用；MATLAB 语言的基本规定，包括数值的表示、变量命名规定、基本运算符、预定义变量以及表达式等。 3) MATLAB 图形绘制功能、M 文件程序设计和线性控制系统传递函数模型的建立等。 2. 实验仪器 PC 计算机一台，MATLAB 软件1套 3. 实验内容 1) MATLAB 的启动 这里介绍MATLAB 装入硬盘后，如何创建MATLAB 的工作环境。 方法一 MATLAB 的工作环境由matlab.exe 创建，该程序驻留在文件夹matlab\bin\ 中。它的图标是 matlab 。只要从<我的电脑>或<资源管理器>中去找这个程序，然后双击此图标，就会自动创建如图1所示的MATLAB6.5 版的工作平台 。 图1 在英文Windows 平台上的MATLAB6.5 MATLAB 工作平台 方法二 假如经常使用MATLAB ，则可以在Windows 桌面上创建一个MATLAB 快捷方式图标。具体办法为： 把<我的电脑>中的 matlab 图标用鼠标点亮，然后直接把此图标拖到Windows 桌面上即可。此后，直接双击Windows 桌面上的matlab 图标，就可建立图1所示的 MATLAB 工作平台。 2) MATLAB 工作环境平台 桌面平台是各桌面组件的展示平台，默认设置情况下的桌面平台包括 6 个窗口，具体如下： ① MATLAB 窗口 Command Window
《信号与系统》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 （） plot(t,f) %用plot 函数绘制连续函数 title('f=exp(t)') %用title 函数设置图形的名称 axis([-5,5,-1,100]) %用axis 函数规定横纵坐标的范围 a=2时： t=-5:0.01:5 f=exp(2*t) %调用指数函数exp （） plot(t,f) title('f=exp(2*t)') axis([-5,5,-1,100]) a=-2时： t=-5:0.01:5 f=exp(-2*t) plot(t,f) title('f=exp(-2*t)') axis([-5,5,-1,100]) （4）)()(t R t f = t=-5:0.01:5 f=rectpuls(t,2) %用rectpuls(t,a)表示门函数，默认以零点为中心，宽度为a plot(t,f) title('f=R(t)') axis([-5 5 -0.5 1.5]) （5）)()(t Sa t f ω= ω=1时：
基于Matlab 的信号与系统实验指导 实验一 连续时间信号在Matlab 中的表示 一、实验目的 1、学会运用Matlab 表示常用连续时间信号的方法 2、观察并熟悉这些信号的波形和特性 二、实验原理及实例分析 1、信号的定义与分类 2、如何表示连续信号？ 连续信号的表示方法有两种；符号推理法和数值法。 从严格意义上讲，Matlab 数值计算的方法不能处理连续时间信号。然而，可利用连续信号在等时间间隔点的取样值来近似表示连续信号，即当取样时间间隔足够小时，这些离散样值能被Matlab 处理，并且能较好地近似表示连续信号。 3、Matlab 提供了大量生成基本信号的函数。如： （1）指数信号：K*exp(a*t) （2）正弦信号：K*sin(w*t+phi)和K*cos(w*t+phi) （3）复指数信号：K*exp((a+i*b)*t) （4）抽样信号：sin(t*pi) 注意：在Matlab 中用与Sa(t)类似的sinc(t)函数表示，定义为：)t /()t (sin )t (sinc ππ= （5）矩形脉冲信号：rectpuls(t,width) （6）周期矩形脉冲信号：square(t,DUTY)，其中DUTY 参数表示信号的占空比
DUTY%，即在一个周期脉冲宽度（正值部分）与脉冲周期的比值。占空比默认为0.5。 （7）三角波脉冲信号：tripuls(t, width, skew)，其中skew 取值范围在-1~+1之间。 （8）周期三角波信号：sawtooth(t, width) （9）单位阶跃信号：y=(t>=0) 三、实验内容 1、验证实验内容 直流及上述9个信号 2、程序设计实验内容 （1）利用Matlab 命令画出下列连续信号的波形图。 （a ）)4/3t (2cos π+ （b ） )t (u )e 2(t -- （c ）)]2()(u )][t (cos 1[--+t u t π （2）利用Matlab 命令画出复信号) 4/t (j 2e )t (f π+=的实部、虚部、模和辐角。 四、实验报告要求 1、格式：实验名称、实验目的、实验原理、实验环境、实验内容、实验思考等 2、实验内容：程序设计实验部分源代码及运行结果图示。