MATLAB实验二(修改)剖析.doc

实验二 信号的表示及其基本运算 一、实验目的 1、掌握连续信号及其MATLAB实现方法； 2、掌握离散信号及其MATLAB实现方法 3、掌握离散信号的基本运算方法，以及MATLAB实现 4 熟悉应用MATLAB实现求解系统响应的方法 例1.1：单位冲击信号的MATLAB实现程序如下： t1=-4; t2=4; t0=0; 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]); 任意函数 例1.2：用MATLAB画出如下表达式的脉冲序列 3 单位阶跃函数 例1.3：用MATLAB实现单位阶跃函数 clear all; t=-0.5:0.001:1; t0=0; u=stepfun(t,t0); plot(t,u) axis([-0.5 1 -0.2 1.2]) 斜坡函数 例1.4：用MATLAB实现g(t)=3(t-1) clear all; t=0:0.01:3; B=3; t0=1; u=stepfun(t,t0); n=length(t); for i=1:n u(i)=B*u(i)*(t(i)-t0); end plot(t,u) axis([-0.2 3.1 -0.2 6.2]) 5 抽样信号 抽样信号Sa(t)=sin(t)/t在MATLAB中用 sinc 函数表示。 定义为 t=-3*pi:pi/100:3*pi; ft=sinc(t/pi); plot(t,ft); grid on; axis([-10,10,-0.5,1.2]); %定义画图范围，横轴，纵轴 title(抽样信号) %定义图的标题名字 6 指数函数 例1.5：用MATLAB实现 7 正弦函数 例1.6：用MATLAB实现正弦函数f(t)=3cos(10πt+1) 8 虚指数信号 例 虚指数信号 调用格式是f=exp((j*w)*t) t=0:0.01:15; w=pi/4; X=exp(j*w*t); Xr=real(X); %取实部 Xi=imag(X); %取虚部 Xa=abs(X); %取模 Xn=angle(X); %取相位 subplot(2,2,1),plot(t,Xr),axis([0,15,-(max(Xa)+0.5),max(Xa)+0.5]), title(实部); subplot(2,2,3),plot(t,Xi),axis([0,15,-(max(Xa)+0.5),max(Xa)+0.5]), title(虚部); subplot(2,2,2), plot(t,Xa),axis([0,15,0,max(Xa)+1]),title(模); subplot(2,2,4),plot(t,Xn),axis([0,15,-(max(Xn)+1),max(Xn)+1]),title(相角); %subplot(m,n,i) 命令是建立m行n列画图窗口，并指定画图位置i 9 复指数信号 例 复指数信号 调用格式是f=exp((a+j*b)*t) t=0:0.01:3; a=-1;b=10; f=exp((a+j*b)*t); subplot(2,2,1),plot(t,real(f)),title(实部) subplot(2,2,3),plot(t,imag(f)),title(虚部) subplot(2,2,2),plot(t,abs(f)),title(模) subplot(2,2,4),plot(t,angle(f)),title(相角) （二）、离散信号及其MATLAB实现 单位冲激序列 例2.1：用MATLAB产生64点的单位冲激序列 clear all; N=64; x=zeros(1,N); x(1)=1; xn=0:N-1; stem(xn,x) axis([-1 65 0 1.1]) 任意序列 例2.2：用MATLAB画出如下表达式的脉冲序列 单位阶跃序列 例2.3：用MATLAB实现单位阶跃函数 斜坡序列 例2.4：用MATLAB实现g(n)=3(n-4)点数为32的斜坡序列 clear all; N=32; k=4 B=3; t0=1; x=[zeros(1,k) ones(1,N-k)]; for i=1:N x(i)=B*x(i)*(i-k); end xn=0:N-1; stem(xn,x)