信号与系统Matlab实验作业

合集下载

信号与系统matlab实验报告

信号与系统matlab实验报告

信号与系统MATLAB实验报告实验目的本实验旨在通过MATLAB软件进行信号与系统的相关实验,探究信号与系统的特性与应用。

实验步骤1. 准备工作在正式进行实验之前,我们需要做一些准备工作。

首先,确保已经安装好MATLAB软件,并且熟悉基本的操作方法。

其次,准备好实验所需的信号与系统数据,可以是已知的标准信号,也可以是自己采集的实际信号。

2. 信号的生成与显示使用MATLAB编写代码,生成不同类型的信号。

例如,可以生成正弦信号、方波信号、三角波信号等。

通过绘制信号波形图,观察不同信号的特点和变化。

t = 0:0.1:10; % 时间范围f = 1; % 信号频率s = sin(2*pi*f*t); % 正弦信号plot(t, s); % 绘制信号波形图3. 系统的建模与分析根据实验需求,建立相应的系统模型。

可以是线性时不变系统,也可以是非线性时变系统。

通过MATLAB进行模型的建立和分析,包括系统的时域特性、频域特性、稳定性等。

sys = tf([1, 2], [1, 3, 2]); % 系统传递函数模型step(sys); % 绘制系统的阶跃响应图4. 信号与系统的运算对于给定的信号和系统,进行信号与系统的运算。

例如,进行信号的卷积运算、系统的响应计算等。

通过MATLAB实现运算,并分析结果的意义与应用。

x = [1, 2, 3]; % 输入信号h = [4, 5, 6]; % 系统响应y = conv(x, h); % 信号的卷积运算plot(y); % 绘制卷积结果的波形图5. 实验结果分析根据实验数据和分析结果,对实验进行结果总结与分析。

可以从信号的特性、系统的特性、运算结果等方面进行综合性的讨论和分析。

实验总结通过本次实验,我们学习了如何在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 称为抽样间隔,代表相邻样值间的时间间隔,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实验2

信号与系统matlab实验2

实验二 线性系统时域分析二、连续时间信号卷积请同学们利用MA TLAB 实现下述两个信号的卷积积分:解:1、f (t )=f1(t)*f2(t )= (t2/4+t/4+1/16)u (t+1/2)-(t2/4-t/2+1/4)u (t-1)-(t2/4+t/4-15/16)u (t-3/2)-( -t2/4+t/2+3/4)u (t-3)用MATLAB 绘制波形程序如下:t=-0.5:0.1:3;f=(t.^2/4+t./4+1/16).*u(t+1/2)-(t.^2/4-t./2+1/4).*u(t-1)-(t.^2/4+t./4-15/16).*u(t-3/2)-(-t.^2/4+t./2+3/4).*u(t-3);plot(t,f)title('f(t)=f1(t)*f2(t)')xlabel('t')ylabel('f(t)')图形如下:2、用MA TLAB 进行求解验证:程序如下:p=0.01;k1=-0.5:p:1;f1=0.*k1+1;k2=0:p:2;f2=0.5.*k2;[f,k]=sconv(f1,f2,k1,k2,p);图形如下:三、线性系统时域分析1、连续系统的冲激响应、阶跃响应及MATLAB实现解:(a)该系统的冲激响应h(t)=(2/√65)*sin(√65t/4)e-t/4,阶跃响应g(t)=(1/8)*[1-e-t/4cos(√65t/4)-(1/√65)*e t/4sin(√65t/4)]绘制波形的程序如下:syms tf1=sym('2/65^0.5*exp(-1/4*t)*sin(65^0.5/4*t)');f2=sym('1/8*(1-exp(-1/4*t)*cos(65^0.5/4*t)-1/65^0.5*exp(-1/4*t)* sin(65^0.5/4*t))');subplot(1,2,1)ezplot(f1,[0 8]);subplot(1,2,2)ezplot(f2,[0 8]);图形如下:(b)用MA TLAB进行验证:冲激响应程序如下:a=[2 1 8];b=[0 1];impulse(b,a);图形如下:阶跃响应程序如下:a=[2 1 8];b=[0 1];step(b,a);图形如下:2、LTI连续系统的求解解:(a)系统的零状态响应r(t)=(2e-t-2e-2t-te-2t)u(t)绘制波形的程序如下:f=sym('2*exp(-t)-2*exp(-2*t)-t*exp(-2*t)'); ezplot(f1,[0 5])图形如下:(b)用MA TLAB进行验证:程序如下:a=[1 4 4];b=[1 3];p=0.01;t=0:p:5;x=exp(-1*t);lsim(b,a,x,t);图形如下:。

东南大学信号与系统MATLAB实践第一次作业

东南大学信号与系统MATLAB实践第一次作业

<信号与系统MATLAB实践> 练习一实验一二. 熟悉简单的矩阵输入1.实验代码>>A=[1,2,3;4,5,6;7,8,9]实验结果A =1 2 34 5 67 8 93.实验代码>>B=[9,8,7;6,5,4;3,2,1]C=[4,5,6;7,8,9;1,2,3]实验结果:B =9 8 76 5 43 2 1C =4 5 67 8 91 2 3 4.>> AA =1 2 34 5 67 8 9>> BB =9 8 76 5 43 2 1C =4 5 67 8 91 2 3三. 根本序列运算1.>>A=[1,2,3],B=[4,5,6]A =1 2 3B =4 5 6 >> C=A+BC =5 7 9 >> D=A-BD =-3 -3 -3 >> E=A.*BE =4 10 18 >> F=A./BF =>> G=A.^B1 32 729 >> stem(A)>>stem(B)>> stem(D)>> stem(F)再举例:>> a=[-1,-2,-3] a =-1 -2 -3 >> b=[-4,-5,-6]b =-4 -5 -6 >> c=a+bc =-5 -7 -9 >> d=a-bd =3 3 3 >> e=a.*be =4 10 18 >> f=a./bf =>> g=a.^bg =>> stem(a)>> stem(c)>> stem(e)>> stem(g)2. >>t=0:0.001:10f=5*exp(-t)+3*exp(-2*t);plot(t,f)ylabel('f(t)');xlabel('t');title('(1)');>> t=0:0.001:3;f=(sin(3*t))./(3*t);plot(t,f)ylabel('f(t)');xlabel('t');title('(2)');>> k=0:1:4;f=exp(k); 1 1.52 2.53 3.54 4.550102030405060四. 利用MATLAB求解线性方程组2.>>A=[1,1,1;1,-2,1;1,2,3]b=[2;-1;-1]x=inv(A)*bA =1 1 11 -2 11 2 3b =2-1-1x =4.>> A=[2,3,-1;3,-2,1;1,2,1]b=[18;8;24]x=inv(A)*bA =2 3 -13 -2 11 2 1b =18824x =468实验二二.1.>> k=0:50x=sin(k);stem(x)xlabel('k');ylabel('sinX');title('sin(k)ε(k)');2.>> k=-25:1:25x=sin(k)+sin(pi*k); stem(k,x)xlabel('k');ylabel('f(k)');title('sink+sinπk');3.>> k=3:50x=k.*sin(k);stem(k,x)xlabel('k');ylabel('f(k)');title('ksinkε(k-3)');4.%函数function y=f1(k)if k<0y=(-1)^k;else y=(-1)^k+(0.5)^k; end%运行代码for k=-10:1:10;y4(k+11)=f1(k);endk=-10:1:10;stem(k,y4);xlabel('k');ylabel('f(k)');title('4');七.2.>> f1=[1 1 1 1];f2=[3 2 1];conv(f1,f2)ans =3 5 6 6 3 1 3.函数定义:function [r]= pulse( k )if k<0r=0;elser=1;endend运行代码for k=1:10f1(k)=pulse(k);f2(k)=(0.5^k)*pulse(k);endconv(f1,f2)结果ans =Columns 1 through 10 Columns 11 through 20 Columns 21 through 30 Columns 31 through 394for i=1:10f1(i)=pulse(i);f2(i)=((-0.5)^i)*pulse(i); endconv(f1,f2)结果ans =Columns 1 through 10 Columns 11 through 20 Columns 21 through 30 Columns 31 through 39实验三2.clear;x=[1,2,3,4,5,6,6,5,4,3,2,1];N=0:11;w=-pi:0.01:pi;m=length(x);n=length(w);for i=1:nF(i)=0;for k=1:mF(i)=F(i)+x(k)*exp(-1j*w(i)*k);endendF=F/10;subplot(2,1,1);plot(w,abs(F),'b-');xlabel('w');ylabel('F');title('幅度频谱');grid subplot(2,1,2);plot(w,angle(F),'b-');xlabel('w');X=fftshift(fft(x))/10;subplot(2,1,1);hold on;plot(N*2*pi/12-pi,abs(X),'r.');legend('DIFT算法','DFT算法');subplot(2,1,2);hold on;plot(N*2*pi/12-pi,angle(X),'r.');xlabel('w');ylabel('相位');title('相位频谱');grid三.1.function y=fun1(x)if((-pi<x) && (x<0))y=pi+x;elseif ((0<x) && (x<pi))y=pi-x;elsey=0endclear allclcfor i=1:1000g(i)=fun1(2/1000*i-1);w(i)=(i-1)*0.2*pi;endfor i=1001:10000g(i)=0;w(i)=(i-1)*0.2*pi;endG=fft(g)/1000;subplot(1,2,1);plot(w(1:50),abs(G(1:50)));xlabel('w');ylabel('G');title('DFT幅度频谱'); subplot(1,2,2);plot(w(1:50),angle(G(1:50)))xlabel('w');ylabel('Fi');title('DFT相位频谱');0102030400.511.522.53wGDFT 幅度频谱010203040-3.5-3-2.5-2-1.5-1-0.5wF iDFT 相位频谱2.function y=fun2(x) if x<1 && x>-1 y=cos(pi*x/2); elsey=0; endfor i=1:1000g(i)=fun2(2/1000*i-1); w(i)=(i-1)*0.2*pi; endfor i=1001:10000 g(i)=0;w(i)=(i-1)*0.2*pi; endG=fft(g)/1000; subplot(1,2,1);plot(w(1:50),abs(G(1:50)));xlabel('w');ylabel('G');title('幅度频谱');subplot(1,2,2);plot(w(1:50),angle(G(1:50)))xlabel('w');ylabel('Fi');title('相位频谱');0102030400.10.20.30.40.50.60.7wGDFT 幅度频谱010203040-4-3-2-1123wF iDFT 相位频谱3.function y=fun3(x) if x<0 && x>-1 y=1;elseif x>0 && x<1 y=-1; elsey=0 endfor i=1:1000g(i)=fun3(2/1000*i-1); w(i)=(i-1)*0.2*pi; endfor i=1001:10000 g(i)=0;w(i)=(i-1)*0.2*pi;G=fft(g)/1000; subplot(1,2,1);plot(w(1:50),abs(G(1:50)));xlabel('w');ylabel('G');title('DFT 幅度频谱'); subplot(1,2,2);plot(w(1:50),angle(G(1:50)))xlabel('w');ylabel('Fi');title('DFT 相位频谱');0102030400.10.20.30.40.50.60.70.8wGDFT 幅度频谱010203040-4-3-2-1123wF iDFT 相位频谱练习二实验六一.用MA TLAB 语言描述如下系统,并求出极零点、 1.>> Ns=[1]; Ds=[1,1];sys1=tf(Ns,Ds) 实验结果: sys1 =-----s + 1>> [z,p,k]=tf2zp([1],[1,1])z =Empty matrix: 0-by-1p =-1k =12.>>Ns=[10]Ds=[1,-5,0]sys2=tf(Ns,Ds)实验结果:Ns =10Ds =1 -5 0sys2 =10---------s^2 - 5 s>>[z,p,k]=tf2zp([10],[1,-5,0]) z =Empty matrix: 0-by-1p =5k =10二.系统的系统函数如下,用MATLAB描述如下系统。

信号与系统 MATLAB综合实验

信号与系统    MATLAB综合实验

信号与系统MATLAB综合实验一、实验目的:1、学习MATLAB语言的编程方法及熟悉MATLAB指令。

2、掌握连续时间信号的卷积运算方式,分析建立信号波形间的联系。

3、通过使用MATLAB函数研究线性时不变离散时间系统的时域特性,以加深对线性时不变离散时间系统的时不变性的理解。

二、实验仪器1、计算机2、MATLAB 软件三、实验原理一个离散时间系统是将输入序列变换成输出序列的一种运算。

若以T[•]表示这种运算,则一个离散时间系统可由图1-1来表示,即→∙→(1-1)x n T y n()[]()图1-1 离散时间系统离散时间系统中最重要的、最常用的是“线性时不变系统”。

时不变系统系统的运算关系T[•]在整个运算过程中不随时间(也不随序列的先后)而变化,这种系统称为时不变系统(或称移不变系统)。

这个性质可用以下关系表示:若输入)(ny,则将输入序列移动任意位后,其输出序列除了跟着x的输出为)(n移位外,数值应保持不变,即若)ynm[mT--(m为任意整数)=(xn(()]()][nT=,则)yxn满足以上关系的系统就称为时不变系统。

四、实验内容及结论1、连续时间系统的时域分析已知微分方程: )(2)(3)(2)(3)(t f t f t y t y t y +'=+'+'',1)0(-='-y , 2)0(=-y 若激励信号为)()(t u t f =,利用阶跃响应函数step(sys,t) 求解画波形;利用零状态响应函数lsim 求解画波形;利用卷积函数求解画波形;比较结果。

程序如下:dt=0.001;t1=0:dt:10;f1=-1*exp(-t1)+4*exp(-2*t1);t2=t1;f2=u(t2);f=conv(f1,f2);f=f*dt;t3=0:dt:20;subplot(311)plot(t3,f);xlabel('时间(t)');ylabel('y(t)');title('零状态响应(卷积法)');b=[3 2];a=[1 3 2];sys=tf(b,a);t=0:0.01:10;x=stepfun(t,0);y=lsim(sys,x,t);subplot(312)plot(t,y);xlabel('时间(t)');ylabel('y(t)');title('零状态响应(阶跃函数求法)');sys=tf(b,a);t=0:0.1:10;y=step(sys,t);subplot(313)plot(t,y);xlabel('时间t)');ylabel('y(t)');title('阶跃响应');结论:上述三种方法求得的都是输入为阶跃函数时候的零状态响应,也为阶跃响应,通过图形我们可以看出,利用卷积法求出的零状态和另外两种方法求出的零状态响应图形有一点差别,三者在0到10区间上响应都一致,而利用卷积法求的响应却在下面的区间内发生了变化,我试图修改程序,无论怎么改,发现只要调用了卷积函数,求得的图形就像上述的卷积法求的图形一样,不得解。

《信号与系统》matlab仿真实验

《信号与系统》matlab仿真实验

《信号与系统》matlab仿真实验综合实验一《信号与系统》的MATLAB仿真实验一.实验目的1.熟悉MA TLAB软件平台、工具箱、高效的数值计算及符号计算功能。

2.熟悉MATLAB软件的信号处理编程方法和结果的可视化3.了解数字信号处理的计算机仿真方法4.进一步加深对信号与系统的基本原理、方法及应用的理解。

二.实验软件MATLAB 6.5 界面三.实验内容1.基本信号的表示及可视化2.连续信号的时域运算与时域变换3.线性系统的时域分析及Matlab实现4.连续时间信号的频域分析及Matlab实现四.实验原理方法及相关MATLAB函数1.基本信号的表示及可视化1.1 连续时间信号(1)表示出连续信号f(t)=Sa(t)=sin(t)/tMatlab命令如下:t=-10:1.5:10;%向量t时间范围t=t1:p:t2,p为时间间隔f=sin(t)./t;plot(t,f); %显示该信号的时域波形title(‘f(t)=Sa(t)’);xlabel(‘t’)axis([-10,10,-0.4,1.1])注:改变p可使信号波形更平滑,同学们可以试一试。

(2)单位阶跃信号定义阶跃函数function f=Heaviside(t)f=(t>0)调用阶跃函数t=-1:0.01:3;f=Heaviside(t)plot(t,f);axis([-1,3,-0.2,1.2]);(2)单位冲击信号 (t)定义冲击函数functionchongji(t1,t2,t0)dt=0.01;t=t1:dt:t2;n=length(t);x=zeros(1,n);x(1,(-t0-t1)/dt+1)=1/dt;stairs(t,x);axis([t1,t2,0,1.2/dt])title('单位冲击信号δ(t)')调用的是chongji(-1,5,0);可以试着给别的t1,t2,t0.1.2离散时间信号(1)单位阶跃序列ε(k)定义阶跃序列function jyxulie(k1,k2,k0)k=k1:-k0-1;kk=-k0:k2;n=length(k);nn=length(kk);u=zeros(1,n); %k0前信号赋值为零uu=ones(1,nn);%k0后信号赋值为一stem(kk,uu,’filled’)hold onstem(k,u,’filled’)holdofftitle(‘单位阶跃序列’)axis([k1 k20 1.5])调用阶跃序列函数jyxulie(-2,6,0)(3)单位序列δ(k)定义单位序列函数functiondwxulie(k1,k2,k0)k=k1:k2;n=length(k);f=zeros(1,n);f(1,-k0-k1+1)=1;stem(k,f,’filled’)axis([k1,k2,0,1.5])title(‘单位序列δ(k)’)调用单位序列函数dwxulie(-3,5,0)2.连续信号的时域运算与时域变换运算、变换的符号运算方法:相加、相乘、移位、反折、尺度变换、倒相已知信号)]2()2([)21()(--+⨯+=ttttfεε,用matlab求f(t+2),f(t-2),f(-t),f(2t),-f(t),并绘出时域波形。

信号与系统及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中信号的运算

信号与系统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

信号与系统部分01.分别用MATLAB 表示并绘出下列连续时间信号的波形:2()(2)()t f t e u t -=-02.分别用MATLAB 表示并绘出下列连续时间信号的波形:[]()cos()()(4)2tf t u t u t π=--03.分别用MATLAB 表示并绘出下列离散时间信号的波形:()12()()kf k u k =-04.分别用MATLAB 表示并绘出下列离散时间信号的波形:[]()()(8)f t k u k u k =--05.已知信号f (t)的波形如下图所示,试用MATLAB 绘出满足下列要求的信号波形。

()f t -06.已知信号f (t)的波形如下图所示,试用MATLAB 绘出满足下列要求的信号波形。

()f ata =0.507.已知信号f (t)的波形如下图所示,试用MATLAB 绘出满足下列要求的信号波形。

()f at ,a =208.已知信号f (t)的波形如下图所示,试用MATLAB 绘出满足下列要求的信号波形。

(0.51)f t +09.已知两信号1()(1)()f t u t u t =+-,2()()(1)f t u t u t =--,求卷积积分12()()()g t f t f t =*。

10.已知两信号1()()f t tu t =,20()()0t tt te u t f t t e-≥⎧=⎨<⎩,求卷积积分12()()()g t f t f t =*。

11.已知{}{}12()1,1,1,2,()1,2,3,4,5f k f k ==,求两序列的卷积和。

12.已知描述系统的微分方程如下,试用理论分析并计算系统的单位冲激响应h(t),并用MATLAB 绘出系统单位冲激响应的波形,验证结果是否相同。

''()4'()4()'()3()y t y t y t f t f t ++=+13.已知描述系统的微分方程如下,试用解析方法求系统的单位冲激响应h(t),并用MATLAB 绘出系统单位冲激响应的波形,验证结果是否相同。

信号与系统matlab实验习题3 绘制典型信号及其频谱图

信号与系统matlab实验习题3  绘制典型信号及其频谱图

绘制典型信号及其频谱图答案在下面四个常用信号及其傅里叶变换式如表1所示。

(1)绘制单边指数信号及其频谱图的MATLAB程序如下: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)|');figure;max_logF=max(abs(F));plot(w,20*log10(abs(F)/max_logF));xlabel('\omega');ylabel('|F(\omega)| indB');figure;plot(w,angle(F));xlabel('\omega');ylabel('\phi(\omega)');请更改参数,调试此程序,绘制单边指数信号的波形图和频谱图。

观察参数a 对信号波形及其频谱的影响。

注:题目中阴影部分是幅频特性的对数表示形式,单位是(dB),请查阅相关资料,了解这种表示方法的意义及其典型数值对应的线性增益大小。

(2)绘制矩形脉冲信号、升余弦脉冲信号和三角脉冲信号的波形图和频谱图,观察并对比各信号的频带宽度和旁瓣的大小。

(3)更改参数,调试程序,绘制单边指数信号的波形图和频谱图。

观察参数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 ,对信号波形及其频谱的影响。

(完整版)信号与系统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进行信号与系统分析实验

【实验1】利用matlab 进行信号的时域分析 (1)指数信号 >>A=1; >> a=-0.4;>> t=0:0.01:10;>> ft=A*exp(a*t); >> plot(t,ft); >> grid;>> axis([0 10 -0.1 1.1]; >> xlabel('t') >> ylabel('ft')(2)正弦信号 >> A=1; >> w0=2*pi; >> phi=pi/6; >> t=0:0.01:3; >> ft=A*sin(w0*t+phi); >> plot(t,ft); >> grid;>> axis([0 3 -1.1 1.1]); >> xlabel('t') >> ylabel('ft')()t f t Ae α=()sin()f t A t ωϕ=+>>x=linspace(-20,20); >> y=sinc(x/pi); >> plot(x,y);>> grid; >> axis([-21 21 -0.5 1.1]); >> xlabel('x') >> ylabel('y')(4)矩形脉冲信号 >> t=0:0.001:4; >> T=1;>> ft=rectpuls(t-2*T,2*T); >> plot(t,ft); >> grid;>> axis([-1 5 -0.1 1.1]); >> xlabel('t') >>ylabel('ft')t t t Sa t f )sin()()(==)]()([)()(10τετετ+-+==t t A t G t ffunction ft=heaviside(t) ft=(t>0); >> t=-1:0.001:3; >> ft=heaviside(t); >> plot(t,ft); >> grid;>> axis([-1 3 -0.1 1.1]); >> xlabel('t') >>ylabel('ft')(6)复指数信号的时域波形 >> t=0:0.1:60;>> f=exp(-0.1*t).*sin(2/3*t); >> plot(t,f); >> grid;>> axis([0 60 -1 1]); >> xlabel('Time(sec)') >>ylabel('f(t)')⎩⎨⎧<>=)0(0)0(1)(t t t ε)32sin()(1.0t e t f t -=(7)加入随机噪声的正弦波>> t=0:0.001:50;>> y=sin(2*pi*50*t);>> s=y+randn(size(t)); >> subplot(2,1,1);>> plot(t(1:100),y(1:100)); >> grid;>> subplot(2,1,2);>> plot(t(1:100),s(1:100)); >>grid;(8)周期矩形波>> A=1;>> t=0:0.0001:5;>> y=A*square(2*pi*t,20); >> plot(t,y);>> grid;>> axis([0 5 -1.5 1.5]);(9)信号的基本运算>> syms t;>>f=sym('(t/2+1)*(heaviside(t+2)-heaviside(t-2))');>>subplot(3,2,1),ezplot(f,[-3,3]);>>grid;>> y1=subs(f,t,t+2);>> subplot(3,2,2),ezplot(y1,[-5,1]);>> title('f(t+2)');>> grid;>> y2=subs(f,t,t-2);>> subplot(3,2,3),ezplot(y2,[-1,5]);>> title('f(t-2)');>> grid;>> y3=subs(f,t,-t);>> subplot(3,2,4),ezplot(y3,[-3,3]);>> title('f(-t)');>> grid;>> y4=subs(f,t,2*t);>> subplot(3,2,5),ezplot(y4,[-2,2]);>> title('f(2t)');>> grid;例1求系统y ”(t )+2y ’(t )+100y (t )=10f (t )的零状态响应,已知f (t )=(sin2πt ) ε(t )。

利用matlab设计仿真信号与系统作业

利用matlab设计仿真信号与系统作业

Q1-2:以Q1_2为文件名存盘,产生实门信号)(2t g 和信号t t g t f π10cos )()(2=。

要求在图形中加上网格线,并使用函数axis()控制图形的时间范围在-2~2秒之间。

然后执行该程序,保存所的图形。

syms t;gt=sym(heaviside(t+1)-heaviside(t-1));ft=gt*cos(10*pi*t);ezplot(ft,[-2,2]);grid on;axis([-2,2,-2,2]);title('ft=gt*cos(10*pi*t),gt=u(t+2)-u(t-2)');xlabel('t');Q1-4:根据示例程序的编程方法,编写一个MATLAB 程序,以Q1_5为文件名存盘,由给定信号x(t) = e -0.5t u(t)求信号y(t) = x(1.5t+3),并绘制出x(t) 和y(t)的图形。

编写的程序Q1_5如下:syms txt=exp((-0.5*t))*heaviside(t);yt=xt*(1.5*t+3);subplot(1,2,1),ezplot(xt,[-2,6]),title(' xt=exp((-0.5*t))*u(t)')subplot(1,2,2),ezplot(yt,[-2,6]),title(' yt=xt*(1.5*t+3)')信号x(t)的波形图 和 信号y(t) = x(1.5t+3) 的波形Q1-6编写程序Q1_8,用Matlab 的方法计算并绘制由如下微分方程表示系统的冲激响应和阶跃响应,并分为上下两个子图绘制在一个图中。

)()(5.0)(2)(3)('22t f t f t y dt t dy dtt y d +=++num=[0.5,1];den=[1,3,2];t=0:0.01:8;subplot(2,1,1) ,impulse(num,den,8) subplot(2,1,2) ,step(num,den,8)Q1-7:做如下总结:1、信号与系统分析,就是基于信号的分解,在时域中,信号主要分解成:连续时间信号和离散时间信号2、写出卷积的运算步骤,并谈谈你对卷积的一些基本性质的理解。

北京交通大学《信号与系统》 课后matlab作业

北京交通大学《信号与系统》 课后matlab作业

Matlab课后作业1.M2-1(1)Matlab程序:t=-5:0.01:5;x=(t>0)-(t>2);plot(t,x);axis([-5,5,-2,2]);仿真结果:(8)Matlab程序:t=-10:0.01:10;pi=3.14;x=sin(pi*t)./(pi*t).*cos(30*t);plot(t,x);仿真结果:M2-2Matlab程序:t=-2:0.001:2;x=(t>-1)-(t>0)+2*tripuls(t-0.5,1,0); plot(t,x);axis([-2,2,-2,2]);仿真结果:M3-3(1)function yt=f(t)yt=t.*(t>0)-t.*(t>=2)+2*(t>=2)-3*(t>3)+(t>5); (2)Matlab程序:t=-10:0.01:11;subplot(3,1,1);plot(t,f(t));title('x(t)');axis([-1,6,-2,3]);subplot(3,1,2);plot(t,f(0.5*t));axis([-1,11,-2,3]);title('x(0.5t)');subplot(3,1,3);plot(t,f(2-0.5*t));title('x(2-0.5t)');axis([-9,5,-2,3]);仿真结果:M2-9(1)Matlab程序:k=-4:7;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; stem(k,x);仿真结果:(2)Matlab程序:k=-12:21;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; N=length(x);y=zeros(1,3*N-2);y(1:3:end)=x;stem(k,y);仿真结果:Matlab程序:k=-1:3;x=[0,0,-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; x1=x(1:3:end);stem(k-1,x1);仿真结果:(3)Matlab程序:k=-6:5;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; stem(k,x);仿真结果:程序>> k=-2:9;>> x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; >> stem(k,x);结果程序>> k=-4:7;>> x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; >> xk=fliplr(x);>> k1=-fliplr(k);>> stem(k1,xk);结果M3-1(1)程序>> ts=0;te=5;dt=0.01; >> sys=tf([2 1],[1 3 2]); >> t=ts:dt:te;>> x=exp(-3*t).*(t>=0); >> y=lsim(sys,x,t);>> plot(t,y);>> xlabel('Time(sec)') >> ylabel('y(t)')结果(2)程序>> ts=0;te=5;dt=0.0001; >>sys=tf([2 1],[1 3 2]); >>t=sys:dt:te;>>x=exp(-3*t).*(t>=0); >>y=lsim(sys,x,t);>>plot(t,y);>>xlabel('Time(sec)') >>ylabel('y(t)')结果M3-4>> x=[0.85,0.53,0.21,0.67,0.84,0.12]; >> k1=-2:3;>> h=[0.68,0.37,0.83,0.52,0.71];>> k2=-1:3;>> y=conv(x,h);>> k=(k1(1)+k2(1)):(k1(end)+k2(end)); >> stem(k,y)结果M6-1(1)>> num=[16 0 0];>> den=[1 5.6569 816 2262.7 160000]; >> [r,p,k]=residue(num,den)得r =0.0992 - 1.5147i0.0992 + 1.5147i-0.0992 + 1.3137i-0.0992 - 1.3137ip =-1.5145 +21.4145i-1.5145 -21.4145i-1.3140 +18.5860i-1.3140 -18.5860ik =[]所以可得 X(s)=j s j j s j j s j 5860.183140.13137.10992.05860.183140.13137.10992.04145.215145.15147.10992.021.4145j -1.5145s j 5147.1-0992.0++--+-++-++++++x(t)=3.0108e-1.5145tcos(21.4145t-1.5054)u(t)+2.635e-1.314tcos(18.586t+1.6462)u(t ) (2)X(s)=)2552^)(5(2^+++s s s s解:>> num=[1 0 0 0];den=conv([1 5],[1 5 25]);[r,p,k]=residue(num,den)[angle,mag]=cart2pol(real(r),imag(r))得r =-5.0000 + 0.0000i-2.5000 - 1.4434i-2.5000 + 1.4434ip =-5.0000 + 0.0000i-2.5000 + 4.3301i-2.5000 - 4.3301ik =1angle =3.1416-2.61802.6180mag =5.00002.88682.8868所以X(s)=3301.45.24434.15.23301.45.24434.15.25s 5.0-1j s j j s j +++-+-+--+++x(t)=δ(t)+5e-5tu(t)+5.7736e-2.5tcos(4.3301t-2.618)u(t)M6-2程序>> t=0:0.1:10;>> y1=(2.5*exp(-t)-1.5*exp(-3*t)).*(t>=0);>> y2=((1/3)+2*exp(-t)-(5/6)*exp(-3*t)).*(t>=0);>> y=((1/3)+(9/2)*exp(-t)-(7/3)*exp(-3*t)).*(t>=0);>> plot(t,y1,'r-',t,y2,'g--',t,y,'b-')>> xlabel('Time');>> legend('零输入响应','零状态响应','完全响应')结果M6-5>> num=[1 2];>> den=[1 2 2 1];>> sys=tf(num,den);>> pzmap(sys)>> num=[1 2];den=[1 2 2 1];[r,p,k]=residue(num,den) [angle,mag]=cart2pol(real(r),imag(r))1.0000 + 0.0000i-0.5000 - 0.8660i-0.5000 + 0.8660ip =-1.0000 + 0.0000i-0.5000 + 0.8660i-0.5000 - 0.8660ik =[]angle =-2.09442.0944mag =1.00001.00001.0000所以H(s)=866.05.0866.05.0866.05.0866.05.01s 1j s j j s j +++-+-+--++系统冲激响应h(t)=e-tu(t)+2e-0.5tcos(0.866t-2.0944)u(t)>> num=[1 2];>> den=conv([1 0],[1 2 2 1]);>> [r,p,k]=residue(num,den)r =-1.0000 + 0.0000i-0.5000 + 0.8660i-0.5000 - 0.8660i2.0000 + 0.0000ip =-1.0000 + 0.0000i-0.5000 + 0.8660i-0.5000 - 0.8660i0.0000 + 0.0000ik =[][angle,mag]=cart2pol(real(r),imag(r))angle =3.14162.0944-2.0944mag =1.00001.00001.00002.0000所以Y(s)=s j s j j s j 2866.05.0866.05.0866.05.0866.05.0-1s 1-+++--+-++++ 系统阶跃响应y(t)=e-tu(t)+2e-0.5tcos(0.866t+2.0944)u(t)因为系统的冲激响应h(t)=e-tu(t)-1.00001e-0.5tcos(0.866t)u(t)+1.73205e-0.5tsin(0.866t)u(t) 所以系统的频率响应H(j ω)=5.0)866.0(866.05.0)866.0(866.05.0)866.0(5.05.0)886.0(5.01j 1j j j j ++--+--++-+++++ωωωωω。

信号与系统matlab仿真实验报告

信号与系统matlab仿真实验报告

题号
1、对于一般的正弦信号、复指数信号、指数信号能够画出其波形图,分析其有
无周期,有周期的给出周期值;
2、能够画出任意输入信号的时域波形图;
注:请区别CT和DT信号。

3、能够画出信号的频谱图,包括幅度谱和相位谱;
4、能够实现信号卷积,并画出卷积后信号的波形图;
1.我把前四个要求集成在了一个主面板里,通过面板platform调用四个不同的功能。

2.第一个分界面general_signal用于实现题目1,由用户选择信号类型并输入信号相关参数。

3.第二个分界面any_signal用于实现题目1、2,同时显示出信号的时域波形以及频谱图、相位图。

4.第三、四个界面signal_ct_conv、signal_dt_conv分别用于实现连续和离散的卷积,对应题目4。

信号与系统matlab实验与答案

信号与系统matlab实验与答案

产生离散衰减正弦序列xn 0.8nsin πn, 0 n 10,并画出其波形图。

4n=0:10;x=sin(pi/4*n).*0.8.^n;stem(n,x);xlabel( 'n' );ylabel( 'x(n)' );用MATLAB生成信号sincat t0,a和t0都是实数,4 t 10,画波形图。

观察并分析a和t0的变化对波形的影响。

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的正弦波可表示为x a t cos2π,对其进行等间隔抽样,得到ft的离散样值序列可表示为xnx a t ,其中T称为抽样间隔,代表相邻tnT样值间的时间间隔,f s1表示抽样频率,即单位时间内抽取样值的个数。

T抽样频率取f s40Hz,信号频率f分别取5Hz,10Hz,20Hz和30Hz。

请在同一张图中同时画出连续信号x a t t和序列xn nT的波形图,并观察和对比分析样值序列的变化。

可能用到的函数为plot,stem,holdon。

fs=40;t=0:1/fs:1;%?μ?ê·?±e?a5Hz,10Hz,20Hz,30Hzf1=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实验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 西电版)实验2 常用离散时间

信号与系统实验(MATLAB 西电版)实验2  常用离散时间

实验2 常用离散时间信号的实现 图 2.3 数值法生成的单位阶跃序列
实验2 常用离散时间信号的实现
4) MATLAB
clf; c=-(1/12)+(pi/6)*i; K=2; n=0:40; x=K*exp(c*n); subplot(2,1,1);
实验2 常用离散时间信号的实现
stem(n,real(x)); ylabel(′幅值f(k)′); title(′实部′); subplot(2,1,2); stem(n,imag(x)); xlabel(′时间(k)′); ylabel(′幅值f(k)′); title(′虚部′); 用数值法生成的复指数序列如图2.4
实验2 常用离散时间信号的实现 图 2.6 数值法生成的正弦序列
实验2 常用离散时间信号的实现
7) 单位斜坡序列 MATLAB
clf; k1=-1; k2=20; k0=0; n=[k1:k2]; if k0>=k2 x=zeros(1,length(n)); elseif (k0<k2)&(k0>k1)
实验2 常用离散时间信号的实现 图 2.10 数值法生成的幅值调制序列
实验2 常用离散时间信号的实现
11) MATLAB clf; R=51; d=0.8*(rand(1,R)-0.5); % m=0:R-1; s=2*m.*(0.9.^m); % x=s+d; %
实验2 常用离散时间信号的实现
实验2 常用离散时间信号的实现 图 2.8 数值法生成的随机序列
实验2 常用离散时间信号的实现
9) 扫频正弦序列 MATLAB
n=0:100; a=pi/2/100; b=0; arg=a*n.*n + b*n; x=cos(arg); clf; stem(n,x); axis([0,100,-1.5,1.5]); grid; axis; title(′扫频正弦序列′) xlabel(′k′); ylabel(′f(k)′); 用数值法生成的扫频正弦序列如图2.9

信号与系统——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求解线性方程组。

信号与系统matlab实验6

信号与系统matlab实验6

实验六 拉普拉斯变换及其逆变换一、实验目的(1)掌握连续系统及信号拉普拉斯变换概念(2)掌握利用MATLAB 绘制系统零极点图的方法 (3)掌握利用MATLAB 求解拉普拉斯逆变换的方法二、实验内容1、求解下述信号的拉普拉斯变换,并利用MATLAB 绘制拉普拉斯变换的曲面图:(1))()2cos()(t u t t f = (2))()sin()(2t u t e t f t -=2、已知连续时间信号)()(3t u e t f t -=,试求出该信号的拉普拉斯变换)(s F 和傅立叶变换)(ωj F ,用MATLAB 绘出拉普拉斯变换曲面图)(s F 及幅频曲线)(ωj F ,观察曲面图在虚轴上的剖面图,并将其与幅频曲线相比较,分析频域与复频域的对应关系。

3、已知信号的拉普拉斯变换如下所示,试用MATLAB 绘制曲面图,观察拉普拉斯变换零极点分布对曲面图的影响。

(1))3)(2()4)(1()(++++=s s s s s s F(2)44)(22+-=s s s F4、试用MATLAB 求下列信号的拉普拉斯逆变换(1)ss s s s s F 6545)(232++++=(2)1221)(23+++=s s s s F三、实验程序及结果1、求解下述信号的拉普拉斯变换,并利用MATLAB绘制拉普拉斯变换的曲面图:(1))tuf(tt2()cos()程序如下:clf;a=-0.5:0.08:0.5;b=-3.5:0.08:3.5;[a,b]=meshgrid(a,b);d=ones(size(a));c=a+i*b; %确定绘制曲面图的复平面区域c=c.*c;c=c+4.*d;c=1./c;c=abs(c); %计算拉普拉斯变换的样值surf(a,b,c); %绘制曲面图title('cos(2t)u(t)的拉氏变换曲面图');colormap(hsv);结果如图:(2))tef t-=tu((2t)sin()程序如下:clf;a=-3:0.06:-1;b=-2:0.06:2;[a,b]=meshgrid(a,b);d=ones(size(a));c=a+i*b; %确定绘制曲面图的复平面区域c=c+2;c=c.*c;c=c+d;c=1./c;c=abs(c); %计算拉普拉斯变换的样值surf(a,b,c); %绘制曲面图title('exp(-2t)sin(t)u(t)的拉氏变换曲面图');colormap(hsv);结果如图:2、已知连续时间信号)()(3t u e t f t -=,试求出该信号的拉普拉斯变换)(s F 和傅立叶变换)(ωj F ,用MATLAB 绘出拉普拉斯变换曲面图)(s F 及幅频曲线)(ωj F ,观察曲面图在虚轴上的剖面图,并将其与幅频曲线相比较,分析频域与复频域的对应关系。

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

实验一 典型连续时间信号和离散时间信号一、实验目的掌握利用Matlab 画图函数和符号函数显示典型连续时间信号波形、典型时间离散信号、连续时间信号在时域中的自变量变换。

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

2)画出复指数信号()()j t f t e σω+=当0.4, 8σω==(0<t<10)时的实部和虚部的波形图。

5)单位冲击信号可看作是宽度为∆,幅度为1/∆的矩形脉冲,即t=t 1处的冲击信号为2、典型离散信号的表示(单位样值序列、单位阶跃序列、实指数序列、正弦序列、复指数序列)编写函数产生下列序列:1)单位脉冲序列,起点n 0,终点n f ,在n s 处有一单位脉冲。

2)单位阶跃序列,起点n 0,终点n f ,在n s 前序列值为0,在n s 后序列值为1。

对于1)、2)小题,最后以参数n 0= -10,n f =10,n s = -3为例,画出各自波形。

3)画出教材P21图1-26,即[][]n x n a u n =当a =1.2, 0.6, -1.5, -0.8的单边指数序列(-2≤4)画出教材P21图1-27,即00[]sin(), x n n =ΩΩ=的正弦序列(-7≤n ≤14)。

5)画出复指数序列/6[]j nx n eπ=和3[]j nx n e=的实部和虚部(-50≤n≤50)。

1)编写程序(函数),画出教材P10图1-13(a)即f(t)的波形(-6<t<6);2)利用1)中建立的函数,通过自变量替换方式依次画出图1-13(b)、(c)、(d)即实验二连续和离散时间LTI系统的响应及卷积一、实验目的掌握利用Matlab工具箱求解连续时间系统的冲激响应、阶跃响应,离散时间系统的单位样值响应,理解卷积概念。

二、实验内容1、连续时间系统的冲击响应、阶跃响应a. 利用impulse函数画出教材P44例2-15: LTI系统()3()2()dy ty t x tdt+=的冲击响应的波形。

b. 利用step函数画出教材P45例2-17: LTI系统1''()3'()2()'()2()2y t y t y t x t x t++=+的阶2、离散时间系统的单位样值响应利用impz函数画出教材P48例2-21: []3[1]3[2][3][]y n y n y n y n x n--+---=的单位3、连续时间信号卷积画出函数f1(t)=(1+t)[u(t)-u(t-1)]和f2(t)=u(t-1)-u(t-2)的图形,并利用附在后面的4、画出教材P60例2-28中h[n]、x[n]的图形(图2-14(a)(b)),并利用conv函数求出卷积x[n]*h[n]并画出图形(图2-14(f))。

实验三连续时间周期信号的傅里叶级数一、实验目的掌握连续时间周期信号的傅里叶级数的展开和合成,理解吉布斯现象,掌握周期矩形脉冲信号的频谱及脉冲宽度、周期对周期信号频谱的影响。

二、实验内容1、周期信号的傅里叶级数的展开和合成画出如下图对称方波(取E=1、T=1),并采用有限项傅里叶级数对原函数进行逼近,画出对称方波的1、3、5、7、9、11次谐波的傅里叶级数合成波形,观察吉布斯现象。

2、周期矩形脉冲信号的频谱a. 取E=1,?=1, 画出周期矩形脉冲(教材P83图3-6)的傅里叶级数的频谱(教材P83图3-7);b. 取E=1,?=1, 画出教材P85图3-8(a);c. 取E=1,?=1, 画出教材P85图3-8(c)。

实验四 非周期信号的频域分析一、实验目的理解非周期信号的频域分析方法,掌握典型信号的幅度谱和相位谱,理解信号的调制特性,掌握傅里叶变换的性质:尺度变换、时移、频移、卷积定理、对称性、微分特性。

二、实验内容1、利用符号函数fourier 和ifourier 求傅里叶变换和傅里叶逆变换。

a. 利用符号函数fourier 求教材P91双边指数信号||()a t f t e -=当a=3时的傅里叶变换表达式。

b. 利用符号函数ifourier 求教材P92第一个公式222()a F j a ωω=+当a=1时的傅里叶逆变换表达式。

c. 利用符号函数fourier 和ezplot 画出21()()2t f t e u t -=及其幅频谱。

>> F=sym('2/(1+w*w)');>> x=ifourier(F)x = exp(-x)*Heaviside(x)+exp(x)*Heaviside(-x)已知线性调制信号表示式如下:a. 0cos()cos()t t ωΩ;b. 0[1.5sin()]cos()t t ω+Ω式中09ω=Ω,试分别画出它们的波形图和频谱图。

3、傅里叶变换的性质(尺度变换、时移、频移、卷积定理、对称性) a. 设2()(1)(1)()f t u t u t G t =+--=,求1()(21)(21)()y t u t u t G t =+--=的频谱()Y j ω,并与()f t 的频谱()F j ω进行比较。

b. 画出21()()2t f t e u t -=、 2(0.4)11()(0.4)2t f t e u t --=-和2(0.4)21()(0.4)2t f t e u t -+=+的幅度谱和相位谱,观察信号时移对信号频谱的影响。

c. 画出()(1)(1)f t u t u t =+--、201()()j t f t f t e -=和202()()j t f t f t e =的频谱,进行相互比较。

d. 画出()(1)(1)f t u t u t =+--、()()*()y t f t f t =及其()F j ω、()()F j F j ωω∙和()Y j ω的图形,验证时域卷积定理。

e. 设()()f t Sa t =,已知信号()f t 的傅里叶变换为2()()[(1)(1)]F j G u u ωπωπωω==+--,求12()()f t G t π=的傅里叶变换1()F j ω,画出各自的图形,并验证对称性。

实验五连续信号的抽样和恢复一、实验目的理解模拟信号的抽样与重构过程,理解信号时域抽样对频域的影响,理解抽样定理。

二、实验内容设信号f(t)=Sa(t)=sin(t)/t,在抽样间隔分别为(1)T s=0.7?(令?m=1,?c=1.1?m)(2)T s=1.5?(令?m=1,?c=1.1?m)的两种情况下,对信号f(t)进行采样,试编写MATLAB程序代码,并绘制出抽样信号波形、由抽样信号得到的恢复信号波形。

实验六拉普拉斯变换一、实验目的掌握系统零极点求法, 理解其含义; 并能利用零极点分析系统的时域和频域特性; 掌握系统的复频域和频域之间的关系;掌握求系统频率响应的方法。

二、实验内容2、利用meshgrid 、mesh 、surf 函数画出信号f(t)= u(t)-u(t-2)的拉普拉斯变换的曲面图,观察曲面图在虚轴剖面上的曲线,并将其与信号傅里叶变换()F j ω绘制的振幅频谱进行比较。

3、画出22(3)(3)()(5)(10)s s H s s s -+=-+的曲面图,观察拉普拉斯变换的零极点。

4、利用roots函数求根,画出24324()2321sF ss s s s-=+-++和2325(45)()51630s s sF ss s s++=+++的零5、已知拉普拉斯变换324()4s F s s s +=+,利用residue 函数求其拉普拉斯逆变换。

>>a=[2 4];>>b=[1 0 4 0];>> [r p k]=residue(a,b)r =-0.5000 - 0.5000i-0.5000 + 0.5000i1.0000p =0 + 2.0000i0 - 2.0000ik =[]6、已知系统函数为324()32s H s s s s+=++,利用residue 函数求该系统的冲击响应h(t),并利用impulse 函数画出其时域波形,判断系统的稳定性。

>>b=[1,4];>>a=[1,3,2,0];>> [r p k]=residue(b,a)r =1-32p =-2-1k =[]>>impulse(b,a)7、设21()0.08()0.41H j j j ωωω=++,利用freqs 函数画出系统幅频特性曲线和相频特性实验七 离散系统的z 域分析 一、实验目的理解并掌握系统函数的概念; 掌握利用系统函数零极点分析系统的稳定性和频率特性,掌握序列的z 变换及其性质;掌握z 域系统表示和差分方程求解。

二、实验内容1、 利用residuez 函数对)3.0)(6.0()9.05.2()(---=z z z z z F (即)3.0)(6.0(9.05.2)(---=z z z z z F )进行部分分式展开。

2/(1-3/5/z)+1/2/(1-3/10/z)2、设某离散系统的系统函数为:541()341z H z z z +=-+,利用roots 函数求出系统的零极点,并画出系统的零极点图,判断系统是否稳定。

3、利用freqz函数画出离散系统0.5()zH zz-=的系统的幅频特性和相频特性曲线。

相关文档
最新文档