基于MATLAB的信号与系统试验指导编程练习

合集下载

信号与系统实验指导书信号与系统matlab实验

信号与系统实验指导书信号与系统matlab实验

信号与系统实验指导书信号与系统matlab实验信号与系统实验指导书一、实验目的1、掌握用Matlab绘制波形图的方法,学会常见波形图的绘制。

2、掌握用Matlab编写函数的方法3、通过对周期信号和非周期信号的观察,加深对周期信号的理解。

二、实验内容1、实验原理与计算实例1.1 绘制波图的基本函数 Matlab是一种基于矩阵和数组的编程语言,它将所有的变量都看成矩阵。

它不仅有强大的计算功能,还有各种各样的画图功能。

这里主要介绍信号与系统分析中常见的几个Matlab函数,包括Matlab提供的内部函数和自定义函数。

我们可以在命令窗口中每次执行一条Matlab语句;或者生成一个程序,存为M文,供以后执行;或是生成一个函数,在命令窗口中执行。

下面介绍几个基本函数。

(1)单位阶跃函数 M文名:u.m%单位阶跃函数(连续或离散)%调用格式 y=u(t)产生单位阶跃函数 function y=u(t) y=(t>=0)(2)门函数 M文名:rectplus.m,是Matlab的内部函数。

调用格式 y=rectplus(t)产生高度为1,宽度为1的门函数调用格式y=rectplus(t,W) 产生高度为1,宽度为W的门函数(3)三角脉冲函数 M文名:tripuls.m,是Matlab的内部函数。

调用格式 y=tripuls(t) 产生高度为1,宽度为1的三角脉冲函数调用格式 y=tripuls(t,w) 产生高度为1,宽度为w的三角脉冲函数调用格式 y=tripuls(t,w,s)产生高度为1,宽度为w的三角脉冲函数,-1<s<1。

当s=0时,为对称三角形;当S=-1时,为三角形顶点左边。

(4)抽样函数 M文名:Sa.m %抽样函数(连续或者离散)% 高度为1 % 调用格式 y=Sa(t),产生高度为1,第一个过零点为π function f=Sa(t)f=sinc(t./pi) %sinc(t)=sin(πt)/(πt)是MATLAB函数(5)符号函数 M文名:sign.m是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 称为抽样间隔,代表相邻样值间的时间间隔,1s 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 ;% ƵÂÊ·Ö±ðΪ5Hz,10Hz,20Hz,30Hz f1=5;xa = cos(2*pi*f1*t) ; subplot(1, 2, 1) ;plot(t, xa) ;axis([0, max(t), min(xa), max(xa)]) ;xlabel('t(s)') ;ylabel('Xa(t)') ;line([0, max(t)],[0,0]) ; subplot(1, 2, 2) ;stem(t, xa, '.') ;line([0, max(t)], [0, 0]) ;axis([0, max(t), min(xa), max(xa)]) ;xlabel('n') ;ylabel('X(n)') ;频率越高,图像更加密集。

matlab软件仿真实验(信号与系统)(1)

matlab软件仿真实验(信号与系统)(1)

matlab软件仿真实验(信号与系统)(1)《信号与系统实验报告》学院:信息科学与⼯程学院专业:物联⽹⼯程姓名:学号:⽬录实验⼀、MATLAB 基本应⽤实验⼆信号的时域表⽰实验三、连续信号卷积实验四、典型周期信号的频谱表⽰实验五、傅⽴叶变换性质研究实验六、抽样定理与信号恢复实验⼀MATLAB 基本应⽤⼀、实验⽬的:学习MATLAB的基本⽤法,了解 MATLAB 的⽬录结构和基本功能以及MATLAB在信号与系统中的应⽤。

⼆、实验内容:例⼀已知x的取值范围,画出y=sin(x)的图型。

x=0:0.05:4*pi;y=sin(x);plot(y)例⼆计算y=sin(π/5)+4cos(π/4)例三已知z 取值范围,x=sin(z);y=cos(z);画三维图形。

z=0:pi/50:10*pi;x=sin(z);y=cos(z);plot3(x,y,z)xlabel('x')ylabel('y')zlabel('z')例四已知x的取值范围,⽤subplot函数绘图。

参考程序:x=0:0.05:7;y1=sin(x);y2=1.5*cos(x);y3=sin(2*x);y4=5*cos(2*x);subplot(2,2,1),plot(x,y1),title('sin(x)')subplot(2,2,2),plot(x,y2),title('1.5*cos(x)')subplot(2,2,3),plot(x,y3),title('sin(2*x)')subplot(2,2,4),plot(x,y4),title('5*cos(2*x)')连续信号的MATLAB表⽰1、指数信号:指数信号Ae at在MATLAB中可⽤exp函数表⽰,其调⽤形式为:y=A*exp(a*t) (例取 A=1,a=-0.4)参考程序:A=1;a=-0.4;t=0:0.01:10;ft=A*exp(a*t);plot(t,ft);grid on;2、正弦信号:正弦信号Acos(w0t+?)和Asin(w0t+?)分别由函数cos和sin表⽰,其调⽤形式为:A*cos(w0t+phi) ;A*sin(w0t+phi) (例取A=1,w0=2π,?=π/6) 参考程序:A=1;w0=2*pi; phi=pi/6; t=0:0.001:8;ft=A*sin(w0*t+phi);plot(t,ft);grid on ;3、抽样函数:抽样函数Sa(t)在MATLAB中⽤sinc函数表⽰,其定义为:sinc(t)=sin(πt)/( πt)其调⽤形式为:y=sinc(t)参考程序:t=-3*pi:pi/100:3*pi;ft=sinc(t/pi);plot(t,ft);grid on;4、矩形脉冲信号:在MATLAB中⽤rectpuls函数来表⽰,其调⽤形式为:y=rectpuls(t,width),⽤以产⽣⼀个幅值为1,宽度为width,相对于t=0点左右对称的矩形波信号,该函数的横坐标范围由向量t决定,是以t=0为中⼼向左右各展开width/2的范围,width的默认值为1。

信号与系统及matlab实验5到8

信号与系统及matlab实验5到8

实验五 matlab 运算基础一、 实验目的1、熟悉启动和退出matlab 的方法2、熟悉matlab 命令窗口的组成3、掌握建立矩阵的方法4、掌握matlab 各种表达式的书写规则以及常用函数的使用5、掌握matlab 关系运算和逻辑运算二、 实验内容1 求下列表达式的植,然后显示matlab 工作空间的使用情况并保存全部变量 1)21)85sin(21e z o += 2))1ln(2122x x z ++=,其中⎥⎦⎤⎢⎣⎡-+=545.0212i x 3))22arctan(3DBC E A z ππ+=,其中A=2.1,B=-4.5,C=6,D=3.5,E=-5, 2 已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=6821945753412A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=038147196B ;求下列表达式的值 1)A+6*B和A+B-22)A*B和B*A3)A/B和A\B4)[A,B]和[A([1,3],:);B^2].5)A^3和A.^33设矩阵A和B,⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=25242322212019181716151413121110987654321A ,⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--=11134079423096171603B 1)求他们的乘积2)将矩阵右下角23⨯子矩阵赋给D3)查看matlab 工作空间的使用情况.4 设矩阵⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=25242322212019181716151413121110987654321A ,取出A的前两列构成矩阵B,取出矩阵A的前两行构成矩阵C,转置B构成矩阵D,计算A*B,C<D,C&D,C|D,~C|~D 5求下列表达式的值,然后显示MA TLAB 工作空间的使用情况并保存全部变量1))3.0sin(213.03.0+-=-a e e z aa ,其中a=0.3,9.2,8.2,,8.2,9.2,0.3 --- 提示:a 可利用冒号表达式生成向量,求各点函数值时用点乘运算。

信号与系统实验指导书(Matlab版)宁波大学.

信号与系统实验指导书(Matlab版)宁波大学.
e.单位冲击信号可看作是宽度为 ,幅度为 的矩形脉冲,即t=t1处的冲击信号为
画出 ,t1=1的单位冲击信号。
f.画出复指数信号 当 (0<t<10)的实部和虚部的波形图。
2、信号的基本运算(相加、相乘、反折、移位、尺度变换)
a.画出教材(上册)中P13图1-16、1-17(取 ).
b.利用符号函数subs画出教材(上册)中P11图1-13(a)(b)(c)(d),并与P38习题1-4进行对比。
b.取样速率Fs=1000Hz,给出X2(jw)曲线
3、对上面产生的x1(n)、x2(n)序列,采用内插函数sinc重构xa(t)。
实验六
一、实验目的
掌握求系统频率响应的方法;掌握信号经系统的频域分析方法。
二、实验内容
1、设 ,利用freqs函数画出系统幅频特性曲线和相频特性曲线。
2、设 ,利用freqs函数画出系统幅频特性曲线和相频特性曲线。设激励信号为 ,画出稳态响应 的波形。
掌握常用离散信号、线性时不变系统的单位样值响应和零状态响应,理解卷积概念。
二、实验内容
1、常用离散信号的表示(单位序列、单位阶跃序列、正弦序列、实指数序列、复指数序列)
编写程序来产生下列基本脉冲序列:
a.单位脉冲序列,起点n0,终点nf,在ns处有一单位脉冲。
b.单位阶跃序列,起点n0,终点nf,在ns前为0,在ns后为1。
3、信号的奇偶分解
利用符号函数subs画出教材(上册)中P40习题1-18(c)(d)的波形并画出信号的奇分量和偶分量的波形。
实验二
一、实验目的
掌握利用Matlab工具箱求解线性时不变系统的冲激响应、阶跃响应和零状态响应,理解卷积概念。
二、实验内容

基于MATLAB的信号与系统实验教程

基于MATLAB的信号与系统实验教程

基于MATLAB的信号与系统实验教程第一部分 MATLAB基础第1章 MATLAB环境1.1 MATLAB界面图1.1 MATLAB主界面图1.2 Workspace图1.3 MATLAB.m文件编辑窗口界面1.2 文件类型图1.4 设置路径图1.5 例1-1运行结果1.3 系统和程序控制指令1.4 练习第2章 数据类型与数学运算2.1 数值、变量和表达式2.1.1 数值的记述2.1.2 变量命名规则2.1.3 运算符和表达式2.2 数组、矩阵及其运算2.2.1 复数和复数矩阵2.2.2 数组和矩阵的运算2.2.3 特殊矩阵(Specialized matrices)2.3 关系和逻辑运算2.4 练习第3章 数值计算与符号计算3.1 线性代数与矩阵分析3.1.1 线性代数3.1.2 特征值分解3.1.3 奇异值分解3.1.4 矩阵函数3.2 线性方程组求解3.2.1 确定性线性方程组求解3.2.2 线性最小二乘问题的方程求解3.3 数据分析函数图3.1 例3-4运行结果3.4 符号计算图3.2 数值型与符号型数据转换关系3.5 练习第4章 绘图4.1 基本绘图指令4.1.1 plot的基本调用格式图4.1 例4-1运行结果4.1.2 stem: 离散数据绘制(火柴杆图)图4.2 例4-2运行结果4.1.3 polar: 极坐标图图4.3 例4-3运行结果4.2 各种图形标记、控制指令图4.4 例4-4运行结果4.2.1 图的创建与控制4.2.2 轴的产生与控制4.2.3 分格线(grid)、坐标框(box)、图保持(hold)4.2.4 图形标志4.3 其他常用绘图指令4.3.1 其他类型图的绘制图4.5 例4-5运行结果图4.6 例4-6运行结果简易绘图指令图4.7 例4-7运行结果4.4 练习第5章 SIMULINK5.1 SIMULINK的基本使用方法图5.1 Simulink Library Browser窗口图5.2 Pulse Generator模块的参数设置5.2 SIMULINK模型概念及基本模块介绍图5.4 SIMULINK模型的一般结构5.2.1 常用的sources——信号源模块5.2.2 常用的sinks——信号显示与输出模块图5.5 示波器纵坐标设置对话框图5.6 示波器属性对话框5.2.3 math operations——数学运算单元模块5.2.4 continuous——连续系统模块5.2.5 discrete——离散系统模块5.3 SIMULINK模型的仿真5.3.1 仿真参数设置图5.7 仿真设置对话框5.3.2 建立子系统图5.8 例5-2的SIMULINK模型图5.9 例5-2的子系统模型图5.10 例5-2仿真输出波形5.4 练习第6章 M函数和工具箱6.1 M函数6.2 工具箱图6.1 演示程序中的工具箱(Toolbox)使用帮助6.3 练习第7章 MATLAB实用技术遴选7.1 图形用户界面设计7.1.1 设计原则与设计步骤7.1.2 界面与控件介绍图7.1 标准菜单样式7.1.3 GUI实例分析。

信号与系统MATLAB实验-实验二 Matlab中信号的运算

信号与系统MATLAB实验-实验二 Matlab中信号的运算

1、运用funtool对f(x)=sin(x)/x分别进行信号的尺度变换f(2x)、f(0.5x)和信号的移位运算f(x+1)、f(x-1)操作以及f(0.5x+1),分别记录相应波形。

f(x)=sin(x)/x f(x+1)f(2x) f(x-1)f(0.5x) f(0.5x+1)2、已知两连续时间信号如下图所示,1)写出信号的函数表达式,并计算f(t)=f1(t)* f2(t)的解析表达式; 2)用MATLAB 求f(t)=f1(t)* f2(t),并绘出f(t)的时域波形图。

(设定取样时间间隔为dt )【实验思考】:通过不断改变dt 的取值并对比所得到的实验效果,观察当取样时间dt 为多大时,函数conv_cs()的计算结果就是连续时间卷积f(t)=f1(t)* f2(t)的较好近似结果?3、已知两连续时间信号如下图所示,1)写出信号的函数表达式,并计算f(t)=f1(t)* f2(t)的解析表达式;2)用MATLAB 求f(t)=f1(t)* f2(t),并绘出f(t)的时域波形图。

(设定取样时间间隔为dt)【实验思考】:不断改变dt的取值并对比实验效果,当取样时间dt为多大时,函数conv_cs()的计算结果就是连续时间卷积f(t)=f1(t)* f2(t)的较好近似结果?clear alldt = 0.01;t1 = -3:dt:3;f1 = 2*(u(t1+1) - u(t1-1));figure;stairs(t1,f1);hold allgrid ont2 = -3:dt:3;f2 = u(t2+2)-u(t2-2);stairs(t2,f2)[fn, tn] = conv_cs(f1, t1, f2, t2, dt);plot(tn, fn)grid onlegend('f1', 'f2', 'f1*f2')。

信号与系统实验--matlab-软件部分

信号与系统实验--matlab-软件部分

%plot(t,f7); %plot(t,f9); plot(t,f11); xlabel('N=11') subplot(3,2,1); plot(t,f9); xlabel('N=9') subplot(3,2,2); plot(t,f7); xlabel('N=7') subplot(3,2,3); plot(t,f5); xlabel('N=5') subplot(3,2,4); plot(t,f3) xlabel('N=3') subplot(3,2,5); plot(t,f1) xlabel('基波') subplot(3,2,6) 运行结果如下:
plot(t,x2); axis([0 20 -3 3]) subplot(1,3,3); plot(t,x3); axis([0 20 -3 3]) 3、一个程序分别对横轴和纵轴给出锯齿波和三角波的标注。 clear all; w=0.5; t1=0:0.0001:1; t=10*pi*t1; y=sawtooth(t); y1=sawtooth(t,w); subplot(2,1,1) plot(t,y) xlabel('锯齿波'); subplot(2,1,2) plot(t,y1) xlabel('三角波')
1 1 1 [sin t sin 3t sin 5t sin nt 3 5 n n 1,3,5
用 Matlab 编程,依次取一次谐波;一三次谐波的合成;一三五次谐波的合成; 一三五七次谐波的合成…。观察谐波叠加后的波形。 基本程序如下: t0=0,tn=10,dt=0.001;t=t0:dt:tn; w=0.5*pi; f1=(4/pi).*sin(w.*t); f3=(4/pi).*(sin(w.*t)+(1/3)*sin(3.*w.*t)); f5=(4/pi).*(sin(w.*t)+(1/3)*sin(3.*w.*t)+(1/5)*sin(5.*w.*t)); f7=(4/pi).*(sin(w.*t)+(1/3)*sin(3.*w.*t)+(1/5)*sin(5.*w.*t)+(1/7)*sin (7.*w.*t)); f9=(4/pi).*(sin(w.*t)+(1/3)*sin(3.*w.*t)+(1/5)*sin(5.*w.*t)+(1/7)*sin (7.*w.*t)+(1/9)*sin(9.*w.*t)); f11=(4/pi).*(sin(w.*t)+(1/3)*sin(3.*w.*t)+(1/5)*sin(5.*w.*t)+(1/7)*si n(7.*w.*t)+(1/9)*sin(9.*w.*t)+(1/11)*sin(11.*w.*t)) %plot(t,f1); %plot(t,f3); %plot(t,f5);

基于MATLAB的信号与系统试验指导编程练习

基于MATLAB的信号与系统试验指导编程练习

2连续时间信号在MATLAB中的表示2-1.利用MATLAB命令画出下列连续信号的波形图(1)>> t=0:0.01:3;>> ft=2*cos(3*t+pi/4);>> plot(t,ft),grid on;>> axis([0 3 -2.2 2.2]);>> title('2cos(3t+pi/4)')(2)>> t=0:0.01:3;>> ft=2-exp(-t);>> plot(t,ft),grid on;>> title('(2-exp(-t))u(t)')(3)>> t=-1:0.01:1;>> ft=t.*(uCT(t)-uCT(t-1));>> plot(t,ft),grid on>> axis([-1 1 -0.2 1.2]);>> title('t[u(t)-u(t-1)]')(4)>> t=-1:0.01:3;>> ft=(1+cos(pi*t)).*(uCT(t)-uCT(t-2));>> plot(t,ft),grid on>> axis([-1 3 -0.2 2.2]);>> title('[1+cos(pi*t)][u(t)-u(t-2)]')2-2.利用MATLAB命令画出下列复信号的实部、虚部、模和辐角(1)>> t=0:0.01:3;>> ft=2+exp(i*(pi/4)*t)+exp(i*(pi/2)*t);>> subplot(2,2,1);plot(t,real(ft));title('实部');axis([0 3 0 4]);grid on;>> subplot(2,2,2);plot(t,imag(ft));title('虚部');axis([0 3 0 2]);grid on;>> subplot(2,2,3);plot(t,abs(ft));title('模');axis([0 3 0 4]);grid on;>> subplot(2,2,4);plot(t,angle(ft));title('相角');axis([0 3 0 2]);grid on;(2)t=0:0.01:3;>> ft=2*exp(i*(t+pi/4));>> subplot(2,2,1);plot(t,real(ft));title('实部');axis([0 3 0 2]);grid on;');axis([0 3 0 2]);grid on; 虚部>> subplot(2,2,2);plot(t,imag(ft));title('>> subplot(2,2,3);plot(t,abs(ft));title('模');axis([0 3 0 4]);grid on;>> subplot(2,2,4);plot(t,angle(ft));title('相角');axis([0 3 0 4]);grid on;2-3.利用MATLAB命令产生幅度为1、周期为1、占空比为0.5的一个周期矩形脉冲信号>> t=-0.5:0.01:3;>> ft=square(2*pi*t,50);>> plot(t,ft);grid on;axis([-0.5 3 -1.2 1.2]);>> title('幅度为1、周期为1、占空比0.5的周期举行脉冲信号')3连续时间信号在MATLAB中的运算3-1.试用MATLAB命令绘出以下信号的波形图(1)>> syms x t;>> t=-1:0.01:1;>> x=exp(-t).*sin(10*pi*t)+exp(-0.5*t).*sin(9*pi*t);>> plot(t,x)(2)>> syms x t;>> t=-1:0.01:1;>> x=sinc(t).*cos(10*pi*t);>> plot(t,x)3-2.已知连续时间信号f(t)的波形如图3-6所示,试用MATLAB 命令画出下列信号的波形图先画出图3-6:>> t=-2:0.01:2;>>f=(-t-1).*(-uCT(t+2)+uCT(t+1))+uCT(t+1)+uCT(t)-uCT(t-1)-(t-1).*(uCT(t-1)-uCT(t-2))-uCT(t-2);>> plot(t,f)>> axis([-4 4 -1 2])>> title('图3-6')>> t=-2:0.01:2;>> f1=funct2(t-1);>> f2=funct2(2-t);>> f3=funct2(2*t+1);>> f4=funct2(4-t/2);>> f5=(funct2(t)+funct2(-t)).*uCT(t);>> subplot(231);plot(t,f1);grid on;title('f(t-1)');axis([-3 3 -1 2]);>> subplot(232);plot(t,f2);grid on;title('f(2-t)');axis([-3 3 -1 2]);>> subplot(233);plot(t,f3);grid on;title('f(2t-1)');axis([-3 3 -1 2]);>> subplot(234);plot(t,f4);grid on;title('f(4-t/2)');axis([-3 3 -1 2]);>> subplot(235);plot(t,f5);grid on;title('(f(t)+f(-t))u(t)');axis([-3 3 -1 2]);3-3.试用MATLAB命令绘出如图3-7所示信号的偶分量和奇分量>> t=0:0.01:2;>> f=(uCT(t)-uCT(t-2)).*(-t+1);>> plot(t,f);title('图3-7')>> f1=fliplr(f);>> fe=(f+f1)/2;fo=(f-f1)/2;>> subplot(211),plot(t,fe);grid on>> title('fe')>> subplot(212),plot(t,fo);grid on;title('fo')4连续时间信号的卷积计算命令绘出下列信号的卷积积分的时域4-1用MATLAB 波形图>>dt=0.001;t1=-0.5:dt:3.5;>> f1=uCT(t1)-uCT(t1-2);>> t2=t1;>> f2=uCT(t2)+uCT(t2-1)-uCT(t2-2)-uCT(t2-3);>> [t,f]=ctsconv(f1,f2,t1,t2,dt);6周期信号的傅里叶级数及频谱分析6-1已知周期三角信号如图6-5所示,试求出该信号的傅里叶级数,利用MATLAB编程实现其各次谐波的叠加,并验证其收敛性。

(完整版)信号与系统Matlab实验作业

(完整版)信号与系统Matlab实验作业

(完整版)信号与系统Matlab实验作业实验一典型连续时间信号和离散时间信号一、实验目的掌握利用Matlab 画图函数和符号函数显示典型连续时间信号波形、典型时间离散信号、连续时间信号在时域中的自变量变换。

二、实验内容1、典型连续信号的波形表示(单边指数信号、复指数信号、抽样信号、单位阶跃信号、单位冲击信号)1)画出教材P28习题1-1(3) ()[(63)(63)]t f t e u t u t =----的波形图。

function y=u(t) y=t>=0; t=-3:0.01:3;f='exp(t)*(u(6-3*t)-u(-6-3*t))'; ezplot(f,t); grid on;2)画出复指数信号()()j t f t e σω+=当0.4, 8σω==(0<t<10)时的实部和虚部的< p="">波形图。

t=0:0.01:10;f1='exp(0.4*t)*cos(8*t)'; f2='exp(0.4*t)*sin(8*t)'; figure(1) ezplot(f1,t); grid on; figure(2) ezplot(f2,t); grid on;t=-10:0.01:10; f='sin(t)/t'; ezplot(f,t); grid on;t=0:0.01:10;f='(sign(t-3)+1)/2'; ezplot(f,t);grid on;5)单位冲击信号可看作是宽度为?,幅度为1/?的矩形脉冲,即t=t 1处的冲击信号为11111()()0 t t t x t t t otherδ??<<+?=-=画出0.2?=, t 1=1的单位冲击信号。

t=0:0.01:2;f='5*(u(t-1)-u(t-1.2))'; ezplot(f,t); grid on;axis([0 2 -1 6]);2、典型离散信号的表示(单位样值序列、单位阶跃序列、实指数序列、正弦序列、复指数序列)编写函数产生下列序列:1)单位脉冲序列,起点n0,终点n f,在n s处有一单位脉冲。

《基于MATLAB的信号与系统实验指导》编程练习资料

《基于MATLAB的信号与系统实验指导》编程练习资料

《基于MATLAB的信号与系统实验指导》编程练习资料MATLAB是一种强大的计算机计算和编程语言,它可以用于分析、模拟和可视化大量复杂的信号与系统过程,为电子工程师和数字系统工程师提供了全面而又直观的分析和设计工具。

本文旨在为我们介绍一些基于MATLAB的信号与系统实验练习,帮助工程师们更好地理解MATLAB在信号与系统实验中的应用。

首先,我们介绍一种基于 MATLAB 的信号 and 系统实验练习:连续时间信号的分析与处理。

我们以常见的音频信号为例,利用MATLAB的图形用户界面(GUI)和应用程序工具(APPS),可以很方便地构建信号读取、分析及显示工具。

同时,MATLAB的内置函数可以用来进行连续时间信号的傅立叶变换,以及在频域中进行简单的滤波运算,使得信号进行滤波和处理更加方便。

再者,MATLAB也可以用于离散时间信号分析。

因为MATLAB运行结果快速稳定,效率高,所以,它可以很好地用于离散时间信号的微分,积分,以及自回归运算等操作。

此外,MATLAB的典型离散信号工具箱也支持基于序列的采样处理,可以自动识别可能的故障状态,并对采样信号进行FFT、PCA及时域,频域及谱图模型分析,这样可以方便快捷地得到信号采样结果。

最后,MATLAB可以用于实时信号处理实验练习,具有很高的灵活性和可扩展性。

MATLAB提供了一个丰富而完善的实时数据获取、图形显示以及实时算法处理框架。

用户可以轻松设计实时信号处理的仿真实验,在计算机与实体之间的互动中,以信号作为中间载体,模拟实时时延等等影响信号传输的参数,为实时信号处理提供直观的可视化指导。

基于MATLAB的信号与系统实验指导可以帮助工程师更好地理解信号与系统实验的处理过程以及MATLAB的运用,从而更好地实现系统的设计和优化。

基于Matlab的信号与系统实验指导2

基于Matlab的信号与系统实验指导2

基于Matlab的信号与系统实验指导2基于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 π+=的实部、虚部、模和辐⾓。

长江大学信号与系统matlab实验答案

长江大学信号与系统matlab实验答案

实验1 信号变换与系统非时变性质的波形绘制●用MA TLAB画出习题1-8的波形。

●用MA TLAB画出习题1-10的波形。

Eg 1.8代码如下:function [y]=zdyt(t) %定义函数zdyty=-2/3*(t-3).*(heaviside(-t+3)-heaviside(-t));endt0=-10;t1=4;dt=0.02;t=t0:dt:t1;f=zdyt(t);y=zdyt(t+3);x=zdyt(2*t-2);g=zdyt(2-2*t);h=zdyt(-0.5*t-1);fe=0.5*(zdyt(t)+zdyt(-t));fo=0.5*(zdyt(t)-zdyt(-t));subplot(7,1,1),plot(t,f);title('信号波形的变化')ylabel('f(t)')grid;line([t0 t1],[0 0]);subplot(7,1,2),plot(t,y);ylabel('y(t)')grid;line([t0 t1],[0 0]);subplot(7,1,3),plot(t,x);ylabel('x(t)')grid;line([t0 t1],[0 0]);subplot(7,1,4),plot(t,g);ylabel('g(t)')grid;line([t0 t1],[0 0]);subplot(7,1,5),plot(t,h);ylabel('h(t)')grid;line([t0 t1],[0 0]);subplot(7,1,6),plot(t,fe);ylabel('fe(t)')grid;line([t0 t1],[0 0]);subplot(7,1,7),plot(t,fo);ylabel('fo(t)')grid;line([t0 t1],[0 0]);xlabel('Time(sec)')结果:Eg1.10代码如下:function [u]=f(t) %定义函数f(t) u= heaviside(t)-heaviside(t-2); endfunction [u] =y(t) %定义函数y(t)u=2*(t.*heaviside(t)-2*(t-1).*heaviside(t-1)+(t-2).*heaviside(t-2)); endt0=-2;t1=5;dt=0.01; t=t0:dt:t1; f1=f(t); y1=y(t); f2=f(t)-f(t-2); y2=y(t)-y(t-2); f3=f(t)-f(t+1); y3=y(t)-y(t+1);subplot(3,2,1),plot(t,f1); title('激励——响应波形图') ylabel('f1(t)')grid;line([t0 t1],[0 0]);-10-8-6-4-2024012信号波形的变化f (t)-10-8-6-4-2024012y (t)-10-8-6-4-2024012x (t)-10-8-6-4-2024012g (t)-10-8-6-4-2024012h (t)-10-8-6-4-202400.51f e (t)-10-8-6-4-2024-101f o (t)Time(sec)subplot(3,2,2),plot(t,y1); ylabel('y1(t)')grid;line([t0 t1],[0 0]); subplot(3,2,3),plot(t,f2); ylabel('f2(t)')grid;line([t0 t1],[0 0]); subplot(3,2,4),plot(t,y2); ylabel('y2(t)')grid;line([t0 t1],[0 0]); subplot(3,2,5),plot(t,f3); ylabel('f3(t)')grid;line([t0 t1],[0 0]); subplot(3,2,6),plot(t,y3); ylabel('y3(t)')grid;line([t0 t1],[0 0]); xlabel('Time(sec)')结果:实验2 微分方程的符号计算和波形绘制上机内容用MA TLAB 计算习题2-1,并画出系统响应的波形。

信号与系统实验指导书(matlab软件仿真)

信号与系统实验指导书(matlab软件仿真)

信号与系统实验指导书(MATLAB仿真)目录实验一MATLAB 基本应用 (2)实验二信号的时域表示 (7)实验三连续信号卷积 (11)实验四典型周期信号的频谱表示 (18)实验五傅立叶变换性质研究 (23)实验六离散信号分析 (26)实验七离散系统的Z域分析 (29)Matlab相关符号及函数说明 (37)实验一MATLAB 基本应用一、实验目的:学习MATLAB的基本用法,了解 MATLAB 的目录结构和基本功能以及MATLAB在信号与系统中的应用。

二、实验内容:例一已知x的取值范围,画出y=sin(x)的图型。

参考程序:x=0:0.05:4*pi;y=sin(x);plot(y)例二计算y=sin(π/5)+4cos(π/4)例三已知z 取值范围,x=sin(z);y=cos(z);画三维图形。

z=0:pi/50:10*pi;x=sin(z);y=cos(z);plot3(x,y,z)xlabel('x')ylabel('y')zlabel('z')例四已知x的取值范围,用subplot函数绘图。

参考程序:x=0:0.05:7;y1=sin(x);y2=1.5*cos(x);y3=sin(2*x);y4=5*cos(2*x);subplot(2,2,1),plot(x,y1),title('sin(x)')subplot(2,2,2),plot(x,y2),title('1.5*cos(x)')subplot(2,2,3),plot(x,y3),title('sin(2*x)')subplot(2,2,4),plot(x,y4),title('5*cos(2*x)')连续信号的MATLAB表示1、指数信号:指数信号Ae at在MATLAB中可用exp函数表示,其调用形式为:y=A*exp(a*t) (例取A=1,a=-0.4)参考程序:A=1;a=-0.4;t=0:0.01:10;ft=A*exp(a*t);plot(t,ft);grid on;注:grid on是一个函数,表示在画图的时候添加网格线。

信号与系统——Matlab部分习题

信号与系统——Matlab部分习题

课程设计(论文)课程名称:信号与系统课程论文题目:信号与系统——MATLAB课程设计姓名:系:专业:物理学(电子信息工程方向)年级:大学二年级学号:指导教师:职称:2017年 5 月22日目录摘要---------------------------------------------------------------------- 1关键词-------------------------------------------------------------------- 1一、M2-7 --------------------------------------------------------------- 21、问题重述--------------------------------------------------------------- 22、问题分析--------------------------------------------------------------- 23、仿真程序与仿真结果----------------------------------------------------- 2(1)仿真程序---------------------------------------------------------------------------------------------------------- 2频率为262Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 2频率为294Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 3频率为330Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 3频率为349Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 4频率为392Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 4频率为440Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 5频率为494Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 5频率为524Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 6结果分析:--------------------------------------------------------------------------------------------------------------- 6(2)仿真程序:------------------------------------------------------------------------------------------------------- 6仿真结果:--------------------------------------------------------------------------------------------------------------- 7结果分析:--------------------------------------------------------------------------------------------------------------- 7(3)仿真程序:------------------------------------------------------------------------------------------------------- 7仿真结果:--------------------------------------------------------------------------------------------------------------- 7结果分析:--------------------------------------------------------------------------------------------------------------- 8二、M3-3 ------------------------------------------------------------------ 81、问题重述--------------------------------------------------------------- 82、问题分析--------------------------------------------------------------- 83、仿真程序与仿真结果----------------------------------------------------- 8(1)仿真程序:------------------------------------------------------------------------------------------------------- 8仿真结果:--------------------------------------------------------------------------------------------------------------- 8结果分析:--------------------------------------------------------------------------------------------------------------- 9(2)仿真程序:------------------------------------------------------------------------------------------------------- 9仿真结果:--------------------------------------------------------------------------------------------------------------- 9结果分析:--------------------------------------------------------------------------------------------------------------- 9三、M4-3 ----------------------------------------------------------------- 101、问题重述:------------------------------------------------------------ 102、问题分析-------------------------------------------------------------- 103、仿真程序与仿真结果---------------------------------------------------- 10(1)仿真程序:----------------------------------------------------------------------------------------------------- 10仿真结果:------------------------------------------------------------------------------------------------------------- 10结果分析:------------------------------------------------------------------------------------------------------------- 10(2)仿真程序:----------------------------------------------------------------------------------------------------- 11仿真结果:------------------------------------------------------------------------------------------------------------- 11结果分析:------------------------------------------------------------------------------------------------------------- 12(3)女声变男声:-------------------------------------------------------------------------------------------------- 12仿真程序:------------------------------------------------------------------------------------------------------------- 12仿真结果:------------------------------------------------------------------------------------------------------------- 12结果分析:------------------------------------------------------------------------------------------------------------- 124、自主学习内容---------------------------------------------------------- 125、阅读文献-------------------------------------------------------------- 126、发现问题-------------------------------------------------------------- 127、问题探究-------------------------------------------------------------- 13信号与系统——MATLAB课程设计摘要:1.(1)掌握信号的建模,了解基本信号及其特点;(2)掌握基本信号的叠加及播放;(3)掌握MA TLAB对语音信号的读取与播放及其时域波形分析。

matlab与信号系统实验题库(1)

matlab与信号系统实验题库(1)

matlab与信号系统实验题库(1)1.设,试⽤MATLAB画出该系统的幅频特性和相频特性。

w=0:0.025:15;b=[1];a=[0.08,0.4,1];H=freqs(b,a,w);subplot(2,1,1);plot(w,abs(H));grid;title('H(jw)的幅频特性');subplot(2,1,2);plot(w,angle (H));grid;title('H(jw)的相频特性');2. 绘制出y(t)=ε(t-1)- ε(t-4)的幅频特性曲线。

syms t wGt=sym('heaviside(t-1)-heaviside(t-4)');Fw=fourier(Gt,t,w);FFw=maple('convert',Fw,'piecewise');FFP=abs(FFw);ezplot(FFP,[-10*pi 10*pi]);grid;axis([-10*pi 10*pi 0 3.5])3.若某连续系统的输⼊为e(t),输出为r(t),系统的微分⽅程为:,y’(0_)=5,y(0_)=0;,求全响应。

a=[1 5 6];b=[0 3 2];%定义a,b[A B C D]=tf2ss(b,a);%定义A,B,C,,D的返回值t=0:0.01:5;%定义时间范围f=exp(-2*t); %定义函数X0=[5,0];%系统状态变量X=[x1,x2,…..xn]'在t=0时刻的初值。

xlabel('t');%定义横标题title('全响应y(t)');%定义总横标题grid on;%⽹格4.求的拉⽒变换式,并⽤MATLAB绘制拉⽒变换在s平⾯的三维曲⾯图。

syms t s %定义符号变量ft=sym('2*exp(-t)*Heaviside(t)+5*exp(-3*t)*Heaviside(t)'); %定义时间函数f(t)的表达式Fs=laplace(ft) %求f(t)的拉⽒变换式F(s)FFss=abs(Fs); %求出F(s)的模ezmesh(FFss); %画出拉⽒变换的⽹格曲⾯图ezsurf(FFss); %画出带阴影效果的三维曲⾯图colormap(hsv); %设置图形中多条曲线的颜⾊顺序5.求的拉⽒变换式,并⽤MATLAB绘制拉⽒变换在s平⾯的三维曲⾯图。

(完整版)信号与系统matlab实验3

(完整版)信号与系统matlab实验3

一、实现下述周期信号的傅立叶级数分解与合成(a )首先,推导出求解0a ,n a ,nb 的公式,计算出前10次系数; (b )利用MATLAB 求解0a ,n a ,n b 的值,其中n a ,nb 求解前10次系数,并给出利用这些系数合成的信号波形。

解:(a)110220[sign(t) - sign(t - 1)]0.25Ta dt ==⎰ 112202[sign(t) - sign(t - 1)][cos()]Tn n t a dt Tπ=⋅⎰ 112202[sign(t) - sign(t - 1)][sin()]T n n t b dt T π=⋅⎰ 程序:function [A_sym,B_sym]=CTFShchsymsyms t n k xT=4;if nargin<4;Nf=10;endif nargin<5;Nn=32;endx=time_fun_x(t);A0=int(x,t,0,1)/T;As=int(2*x*cos(2*pi*n*t/T)/T,t,0,1);Bs=int(2*x*sin(2*pi*n*t/T)/T,t,0,1);A_sym(1)=double(vpa(A0,Nn));for k=1:NfA_sym(k+1)=double(vpa(subs(As,n,k),Nn));B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn));endif nargout==0c=A_sym;disp(c);d=B_sym;disp(d);t=-8:0.01:9;f1=c(1)+c(2).*cos(2*pi*1*t/T)+d(2).*sin(2*pi*1*t/T); f2=c(3).*cos(2*pi*2*t/T)+d(3).*sin(2*pi*2*t/T); f3=c(4).*cos(2*pi*3*t/T)+d(4).*sin(2*pi*3*t/T); f4=c(5).*cos(2*pi*4*t/T)+d(5).*sin(2*pi*4*t/T); f5=c(6).*cos(2*pi*5*t/T)+d(6).*sin(2*pi*5*t/T);f6=c(7).*cos(2*pi*6*t/T)+d(7).*sin(2*pi*6*t/T);f7=c(8).*cos(2*pi*7*t/T)+d(8).*sin(2*pi*7*t/T);f8=c(9).*cos(2*pi*8*t/T)+d(9).*sin(2*pi*8*t/T);f9=c(10).*cos(2*pi*9*t/T)+d(10).*sin(2*pi*9*t/T);f10=c(11).*cos(2*pi*10*t/T)+d(11).*sin(2*pi*10*t/T);ff1=f1+f2+f3+f4+f5+f6+f7+f8+f9+f10;ff2=f1+f2+f3+f4+f5+f6+f7;ff3=ff2+f8;ff4=ff3+f9;subplot(2,2,1)plot(t,ff1),hold ony=time_fun_e(t) %µ÷ÓÃÁ¬Ðøʱ¼äº¯Êý-ÖÜÆÚ¾ØÐÎÂö³å plot(t,y,'r:')title('ÖÜÆÚ¾ØÐ⨵ÄÐγɡª1+2+3+4+5+6+7+8+9+10´Îг²¨')axis([-4,4.5,-0.5,1.5])grid onsubplot(2,2,2)grid onplot(t,ff2),hold ony=time_fun_e(t)plot(t,y,'r:')title('ÖÜÆÚ¾ØÐ⨵ÄÐγɡª1+2+3+4+5+6+7´Îг²¨')axis([-4,4.5,-0.5,1.5])grid onsubplot(2,2,3)plot(t,ff3),hold ony=time_fun_e(t)plot(t,y,'r:')title('1+2+3+4+5+6+7+8´Îг²¨')axis([-4,4.5,-0.5,1.5])grid onsubplot(2,2,4)plot(t,ff4),hold ony=time_fun_e(t)plot(t,y,'r:')title('1+2+3+4+5+6+7+8+9´Îг²¨')axis([-4,4.5,-0.5,1.5])grid onendfunction x=time_fun_x(t)h=1;x1=sym('0.5+0.5*sign(t)')*h;x=x1-sym('(0.5+0.5*sign(t-1))')*h;%-------------------------------------------function y=time_fun_e(t)a=0.5;T=5;h=1;t=-8:0.01:9;e1=(1/2+1/2.*sign(t))-(1/2+1/2.*sign(t-1));e2=(1/2+1/2.*sign(t-4))-(1/2+1/2.*sign(t-5));e3=(1/2+1/2.*sign(t+4))-(1/2+1/2.*sign(t+3));y=e1+e2+e3;结果如下:A_sym =0.2500 0.3183 0.0000 -0.1061 -0.0000 0.0637 0.0000 -0.0455 -0.0000 0.0354 0.0000B_sym =0 0.3183 0.3183 0.1061 0.0000 0.0637 0.1061 0.0455 0.0000 0.0354 0.0637二、知周期为T=4的三角波,在第一周期(-2<t<2)内表示成:)(,试用MATLAB求该信号的傅立叶级数,并绘制它的频谱图。

matlab信号与系统实验程序

matlab信号与系统实验程序

%两个信号的相加t=0:0.1:4*pi;y1=sin(t);y2=cos(t);y=y1+y2;subplot(311)plot(t,y1);subplot(312)plot(t,y2);subplot(313)plot(t,y);%两个信号的相减t=0:0.1:4*pi;y1=sin(t);y2=cos(t);y=y2-y1;subplot(311)plot(t,y1);subplot(312)plot(t,y2);subplot(313)plot(t,y);%信号的卷积s=0.01;u1=0:s:6;u2=u1;x1=(u1>0);x2=(u1>2);x3=(u2>0);x4=(u2>3);f1=x1-x2;f2=x3-x4;f=conv(f1,f2);f=f*s;u0=u1(1)+u2(1);u3=length(u1)+length(u2)-2; u=u0:s:u3*s;subplot(311)plot(u1,f1);axis([-0.5,2.5,0,1.5]); subplot(312)plot(u2,f2);axis([-0.5,3.5,0,1.5]); subplot(313) plot(u,f);axis([0,6,0,3])%信号的采样s=0.01;t=0:s:8;TS=1/8;n=0:TS:8;y=sin(pi*t);yn=sin(pi*n);subplot(211)plot(t,y);subplot(212)plot(n,yn);stem(n,yn,'.')%信号的恢复:y(t)=sin2pi*t wm=pi*2;a=input('输入ws/wm=:');wc=wm;t0=2;t=-t0:0.01:t0;u1=(t>1);u2=(t>-1);u=u2-u1;x=(cos(pi*t)+1).*u;subplot(221);plot(t,x);axis([-2,2,-0.5,2.5]);title('x(t)');ws=a*wm;Ts=2*pi/ws;N=fix(t0/Ts);n=-N:N;nTs=n*Ts;un1=(nTs>1);un2=(nTs>-1);un=un2-un1;xs=(cos(pi*nTs)+1).*un; subplot(2,2,2);stem(n,xs,'.');title('x(n)');xr=zeros(1,length(t));L=length(-N:N);xa=xr;figure(2);stem(nTs,xs,'.');hold on;for i=1:L;m=(L-1)/2+1-i;xa=Ts*(wc)*xs(i)*sinc((wc)*(t+m*Ts)/pi)/pi;plot(t,xa,'b:');axis([-2,2,-0.5,2.5]);hold on;pausexr=xr+xa;endplot(t,xr,'r');axis([-2,2,-0.5,2.5]);hold on;figure(1);subplot(223);plot(t,xr,'r');axis([-2,2,-0.5,2.5]);title('xr(t)');error=abs(xr-x);subplot(2,2,4)plot(t,error);title('error');%冲激响应b=[1,3];a=[1,3,2];sys=tf(b,a);t=0:0.2:22;y=impulse(sys,t);plot(t,y);axis([0,22,-0.1,1.1])%阶跃响应b=[1,3];a=[1,3,2];sys=tf(b,a);t=0:0.2:22;y=step(sys,t);plot(t,y);axis([0,25,-0.1,1.6])%零状态响应b=[1,3];a=[1,3,2];sys=tf(b,a);t=0:0.01:9; x1=exp(-5*t);y=lsim(sys,x1,t);plot(t,y);axis([0,9,-0.1,0.2])%频率响应w=linspace(-5,5,200);b=[2];a=[1,1];H=freqs(b,a,w);subplot(2,1,1.5]);subplot(2,1,2);plot(w,angle(H));axis([-6,6,-2,2])%稳定性b=[1,-1];a=[1,3,2];zplane(b,a);legend('零点','极点')%离散信号相加a1=[1,3,5];n1=0:2;a2=[3,3,4,2];n2=-1:2;n=min([n1,n2]):max([n1,n2]);f1=zeros(1,length(n));f2=zeros(1,length(n));f1(find((n>=min(n1))&(n<=max(n1))==1))=a1; f2(find((n>=min(n2))&(n<=max(n2))==1))=a2; f=f1+f2;subplot(3,1,1)stem(n1,a1,'filled');axis([-2,3,0,5]);subplot(3,1,2)stem(n2,a2,'filled');axis([-2,3,0,5]);subplot(3,1,3)stem(n,f,'filled');axis([-2,3,0,9]);%离散信号相减a1=[1,3,5];n1=0:2;a2=[3,3,4,2];n2=-1:2;n=min([n1,n2]):max([n1,n2]);f1=zeros(1,length(n));f2=zeros(1,length(n));f1(find((n>=min(n1))&(n<=max(n1))==1))=a1; f2(find((n>=min(n2))&(n<=max(n2))==1))=a2; f=f1-f2;subplot(3,1,1)stem(n1,a1,'filled');axis([-2,3,0,5]);subplot(3,1,2)stem(n2,a2,'filled');axis([-2,3,0,5]);subplot(3,1,3)stem(n,f,'filled');axis([-2,3,-5,5]);%离散信号卷积f1=[1,3,5];k1=0:2;f2=[3,4,2];k2=0:2;f=conv(f1,f2);subplot(3,1,1)stem(k1,f1,'filled');axis([-1,3,0,6]);subplot(3,1,2)stem(k2,f2,'filled');axis([-1,3,0,5]);subplot(3,1,3)stem(0:length(f)-1,f,'filled');axis([-1,5,0,32]);%y(k)-3y(k-1)+3y(k-2)-y(k-3)=f(k)单位序列响应a=[1,-0.7,0.12];b=[2,-1];impz(b,a);axis([-1,10,-1,3])%y(k)-3y(k-1)+3y(k-2)-y(k-3)=f(k)单位阶跃响应a=[1,-0.7,0.12];b=[2,-1];stepz(b,a);axis([-1,10,0,3])%y(k)-3y(k-1)+3y(k-2)-y(k-3)=f(k)零状态响应a=[1,-0.7,0.12];b=[2,-1];n=0:20;f=(1/5).^n; y=filter(b,a,f);subplot(2,1,1)stem(n,f,'filled');axis([-1,21,-1,2]);subplot(2,1,2)stem(n,y,'filled');axis([-1,21,-1,3]);%H(Z)=e^3z+2e^2z+2e^z+1/e^3z频率响应响应w=-4*pi:8*pi/511:4*pi;b=[1 2 2 1];a=[1];h=freqz(b,a,w);subplot(2,1,1)plot(w/pi,abs(h));subplot(2,1,2)plot(w/pi,angle(h));。

信号与系统MATLAB实验全

信号与系统MATLAB实验全

实验篇 信号与系统实验指导实验一、MATLAB 编程基础及典型实例一、实验目的(1) 熟悉MATLAB 软件平台的使用; (2) 熟悉MATLAB 编程方法及常用语句; (3) 掌握MATLAB 的可视化绘图技术;(4) 结合《信号与系统》的特点,编程实现常用信号及其运算。

二、实验原理连续信号是指自变量的取值范围是连续的,且对于一切自变量的取值,除了有若干个不连续点以外,信号都有确定的值与之对应。

严格来说,MATLAB 并不能处理连续信号,而是用等时间间隔点的样值来近似表示连续信号。

当取样时间间隔足够小时,这些离散的样值就能较好地近似连续信号。

矩阵是MATLAB 进行数据处理的基本单元,矩阵运算是MATLAB 最重要的运算。

通常意义上的数量(也称为标量)在MATLAB 系统中是作为1×1的矩阵来处理的,而向量实际上是仅有一行或者一列的矩阵。

通常用向量表示信号的时间取值范围,如n = -5:5,但信号x(n)、向量n 本身的下标都是从1开始的,因此必须用一个与向量x 等长的定位时间变量n ,以及向量x ,才能完整地表示序列x(n)。

这一点详情可参考预备篇示例7的程序说明。

三、实验内容与步骤(1) 新建一个文件夹,以自己的汉语名字命名,以后就用该文件夹专门存放自己所编制的M 文件和产生的图形;将该文件夹设置成当前工作目录。

(2) 绘制信号t)32sin(e x(t)t 2-=的曲线,t 的范围在0 ~ 30s ,取样时间间隔为0.1s.(3) 在n = [-10:10] 范围产生离散序列:⎩⎨⎧≤≤-=其余n0,3n 32n,x(n) ,并绘图。

四、实验报告要求整理并给出“实验内容与步骤”(2)、(3)的程序代码与产生的图形;并回答下面的问题。

(1) 在调用某一函数文件时,该文件中除了输入、输出变量外的其它变量在调用函数结束后是否还存在?这些变量是全局还是局部变量?(2) 设n = -10:0.2:20,你可以通过哪些方法查看向量n 的维数?经过关系运算y = (n >= 3)以后,y 的维数是多少?y 又等于什么?(3) 通过MATLAB 的帮助系统,学习fliplr 函数的功能和使用方法。

信号与系统——MATLAB基本实验

信号与系统——MATLAB基本实验

信号与系统——MATLAB基本实验《信号与系统MATLAB实践》第一次上机作业实验一、熟悉MATLAB基本操作三、基本序列运算1.数组的加减乘除和乘方运算A=[1 2 3];B=[4 5 6];C=A+B;D=A-B;E=A.*B;F=A./B;G=A.^B;subplot(2,4,1);stem(A)subplot(2,4,2);stem(B)subplot(2,4,3);stem(C)subplot(2,4,4);stem(D)subplot(2,4,5);stem(E)subplot(2,4,6);stem(F)subplot(2,4,7);stem(G)(2)t=0:0.001:10x=5*exp(-t)+3*exp(-2*t);plot(t,x)ylabel('f(t)');xlabel('t');title('(2)');(3)t=0:0.001:3x=exp(-t).*sin(2*pi*t); plot(t,x)ylabel('f(t)');xlabel('t');title('(3)');(4)t=0:0.001:3 x=sin(3*t)./(3*t); plot(t,x)ylabel('f(t)'); xlabel('t');title('(4)');(5)k=1:1:6 x=(-2).^(-k); stem(k) xlabel('k'); ylabel('f(k)'); title('(5)');(6)k=0:1:4 x=exp(k); stem(k) xlabel('k'); ylabel('f(k)'); title('(6)');(7)k=1:1:99 x=k;stem(k) xlabel('k'); ylabel('f(k)'); title('(7)');四、利用MATLAB求解线性方程组。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

2连续时间信号在MATLAB中的表示2-1.利用MATLAB命令画出下列连续信号的波形图(1)>> t=0:0.01:3;>> ft=2*cos(3*t+pi/4);>> plot(t,ft),grid on;>> axis([0 3 -2.2 2.2]);>> title('2cos(3t+pi/4)')(2)>> t=0:0.01:3;>> ft=2-exp(-t);>> plot(t,ft),grid on;>> title('(2-exp(-t))u(t)')(3)>> t=-1:0.01:1;>> ft=t.*(uCT(t)-uCT(t-1));>> plot(t,ft),grid on>> axis([-1 1 -0.2 1.2]);>> title('t[u(t)-u(t-1)]')(4)>> t=-1:0.01:3;>> ft=(1+cos(pi*t)).*(uCT(t)-uCT(t-2));>> plot(t,ft),grid on>> axis([-1 3 -0.2 2.2]);>> title('[1+cos(pi*t)][u(t)-u(t-2)]')2-2.利用MATLAB命令画出下列复信号的实部、虚部、模和辐角(1)>> t=0:0.01:3;>> ft=2+exp(i*(pi/4)*t)+exp(i*(pi/2)*t);>> subplot(2,2,1);plot(t,real(ft));title('实部');axis([0 3 0 4]);grid on;>> subplot(2,2,2);plot(t,imag(ft));title('虚部');axis([0 3 0 2]);grid on;>> subplot(2,2,3);plot(t,abs(ft));title('模');axis([0 3 0 4]);grid on;>> subplot(2,2,4);plot(t,angle(ft));title('相角');axis([0 3 0 2]);grid on;(2)t=0:0.01:3;>> ft=2*exp(i*(t+pi/4));>> subplot(2,2,1);plot(t,real(ft));title('实部');axis([0 3 0 2]);grid on;');axis([0 3 0 2]);grid on; 虚部>> subplot(2,2,2);plot(t,imag(ft));title('>> subplot(2,2,3);plot(t,abs(ft));title('模');axis([0 3 0 4]);grid on;>> subplot(2,2,4);plot(t,angle(ft));title('相角');axis([0 3 0 4]);grid on;2-3.利用MATLAB命令产生幅度为1、周期为1、占空比为0.5的一个周期矩形脉冲信号>> t=-0.5:0.01:3;>> ft=square(2*pi*t,50);>> plot(t,ft);grid on;axis([-0.5 3 -1.2 1.2]);>> title('幅度为1、周期为1、占空比0.5的周期举行脉冲信号')3连续时间信号在MATLAB中的运算3-1.试用MATLAB命令绘出以下信号的波形图(1)>> syms x t;>> t=-1:0.01:1;>> x=exp(-t).*sin(10*pi*t)+exp(-0.5*t).*sin(9*pi*t);>> plot(t,x)(2)>> syms x t;>> t=-1:0.01:1;>> x=sinc(t).*cos(10*pi*t);>> plot(t,x)3-2.已知连续时间信号f(t)的波形如图3-6所示,试用MATLAB 命令画出下列信号的波形图先画出图3-6:>> t=-2:0.01:2;>>f=(-t-1).*(-uCT(t+2)+uCT(t+1))+uCT(t+1)+uCT(t)-uCT(t-1)-(t-1).*(uCT(t-1)-uCT(t-2))-uCT(t-2);>> plot(t,f)>> axis([-4 4 -1 2])>> title('图3-6')>> t=-2:0.01:2;>> f1=funct2(t-1);>> f2=funct2(2-t);>> f3=funct2(2*t+1);>> f4=funct2(4-t/2);>> f5=(funct2(t)+funct2(-t)).*uCT(t);>> subplot(231);plot(t,f1);grid on;title('f(t-1)');axis([-3 3 -1 2]);>> subplot(232);plot(t,f2);grid on;title('f(2-t)');axis([-3 3 -1 2]);>> subplot(233);plot(t,f3);grid on;title('f(2t-1)');axis([-3 3 -1 2]);>> subplot(234);plot(t,f4);grid on;title('f(4-t/2)');axis([-3 3 -1 2]);>> subplot(235);plot(t,f5);grid on;title('(f(t)+f(-t))u(t)');axis([-3 3 -1 2]);3-3.试用MATLAB命令绘出如图3-7所示信号的偶分量和奇分量>> t=0:0.01:2;>> f=(uCT(t)-uCT(t-2)).*(-t+1);>> plot(t,f);title('图3-7')>> f1=fliplr(f);>> fe=(f+f1)/2;fo=(f-f1)/2;>> subplot(211),plot(t,fe);grid on>> title('fe')>> subplot(212),plot(t,fo);grid on;title('fo')4连续时间信号的卷积计算命令绘出下列信号的卷积积分的时域4-1用MATLAB 波形图>>dt=0.001;t1=-0.5:dt:3.5;>> f1=uCT(t1)-uCT(t1-2);>> t2=t1;>> f2=uCT(t2)+uCT(t2-1)-uCT(t2-2)-uCT(t2-3);>> [t,f]=ctsconv(f1,f2,t1,t2,dt);6周期信号的傅里叶级数及频谱分析6-1已知周期三角信号如图6-5所示,试求出该信号的傅里叶级数,利用MATLAB编程实现其各次谐波的叠加,并验证其收敛性。

6-2 试用MATLAB分析图6-5中周期三角信号的频谱。

当周期三角信号的周期和三角信号的宽度变化时,试观察分析其频谱的变化。

7 傅里叶变换及其性质7-1试用MATLAB命令求下列信号的傅里叶变换,并绘出其幅度谱和相位谱。

)2 ((1)解:(1)ft1=sym('sin(2*pi*(t-1))/(pi*(t-1))');>> Fw1=simplify(fourier(ft1));>> subplot(211)>> ezplot(abs(Fw1)),grid on>> title('幅度谱')>> phase=atan(imag(Fw1)/real(Fw1));>> subplot(212)>> ezplot(phase);grid on>> title('相位谱'))2(7-2.试用MATLAB命令求下列信号的傅里叶反变换,并绘出其时域信号图。

)2 (1()解:(1)>> syms t>> Fw=sym('10/(3+w*i)-4/(5+w*i)');>> ft=ifourier(Fw,t);>> ezplot(ft),grid on(2)>> syms t>> Fw2=sym('exp(-4*w^2)');>> ft2=ifourier(Fw2,t)ft2 =exp(-t^2/16)/(4*pi^(1/2))3.试用MATLAB数值计算方法求图7-8所示信号的傅里叶变换,并画出其频谱图。

解:4.已知两个门信号的卷积为三角波信号,试用MATLAB命令验证傅里叶变换的时域卷积定理。

解:将门函数先进行时域卷积运算,再将卷积后的结果做傅里叶变换,程序和结果如下:dt = 0.01; t = -2:dt:2.5;f1 = uCT(t+0.5)- uCT(t-0.5);f = conv(f1,f1)*dt;ft=sym('f');Fw = fourier(ft)Fw =2*i*pi*dirac(1,w)将一个门函数先进行傅里叶变换,再将结果与自身相乘,程序和结果如下:dt = 0.01; t = -2:dt:2.5;f1 = uCT(t+0.5)- uCT(t-0.5);ft=sym('f1');Fw = fourier(ft);Fw=Fw*FwFw =-4*pi^2*dirac(1,w)^2由此来验证傅里叶变换的时域卷积定理第8章连续时间LTI系统的频率特性及频域分析8.1试用MATLAB命令求图8-8所示电路系统的幅频特性和相频特性。

相关文档
最新文档