matlab离散正弦调频变换函数,MATLAB实验二(修改).doc

MATLAB实验二(修改)

实验二 信号的表示及其基本运算

一、实验目的

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=