信号与系统实验报告.doc

文档介绍：信号与系统MATLAB实验报告
实验二离散时间信号与系统
实验目的
了解信号处理的基本操作
熟悉一些常用的序列及其应用
]
我们所接触的信号大多为连续信号,为使之便于处理,往往要对其进行采样,对信号抽样并保证其能完全恢复,对抽样频率有一定的限制。
基本的离散序列的定义如下:
单位采样序列
单位阶跃序列
实指数序列
复数指数序列
正余弦序列
周期序列

用MATLAB实现函数impseq(n0,n1,n2),使函数实现,。
程序:
function [x,n]=impseq(n0,n1,n2)
if n1>n2
error('n1 should smaller than or equal to n2.');
end
n=n1:n2;
x=zeros(1,n2-n1+1);
if n0>=n1 & n0<=n2
x(n0-n1+1)=1;
end
运行结果:
>>[x n]=impseq(2,-1,4)
x = 0 0 0 1 0 0
n = -1 0 1 2 3 4
用MATLAB实现函数stepseq(n0,n1,n2),使函数实现u(n-n0),。
程序:
function [x,n]=stepseq(n0,n1,n2)
if n1>n2
error('n1 should smaller than or equal to n2.');
end
n=n1:n2;
x=zeros(1,n2-n1+1);
if n0<=n1
x=x+1;
elseif n0<=n2
for i=(n0-n1+1):n2-n1+1
x(i)=1;
end
end
运行结果:
>>[x n]=stepseq(2,-1,4)
x = 0 0 0 1 1 1
n = -1 0 1 2 3 4

a.               
程序:
function [x,n]=seq3a
x=0:10;
for n=0:10
x(n+1)=power(,n);
end
n=0:10;
运行结果:
>>seq3a
ans =
Columns 1 through 9

Columns 10 through 11

b.           
程序
function [x,n]=seq3b
x=0:10;
for n=0:10;
x(n+1)=exp(n*(2+3j));
end
n=0:10;
运行结果:
>> seq3b
ans =
+008 *
Columns 1 through 4
- + - - +
Columns 5 through 8
- - + - - +
Columns 9 through 11
- - + -
c.
程序
function [x,n]=seq3c
x=0:10;
for n=0:10;
x(n+1)=3*cos(pi*n*+1/3*pi)+2*sin(*pi*n);
end
n=0:10;
运行结果:
>> seq3c
ans =
Columns 1 through 8
- - - - - -
Columns 9 through 11
- - -
“+”实现信号相加,但和的长度必须相等。如果序列长度不等,或者长度虽然相等但采样的位置不同,就不能运用
“+”了。试用MATLAB写出任意序列相加的函数sigadd,其格式如下:
程序:
function [y,n]=sigadd(x1,n1,x2,n2)
if n1(1)<=n2(1);
nsmall=n1;
nlarge=n2;
xsmall=x1;
xlarge=x2;
else
nsmall=n2;
n