信号与系统matlab实验报告
一、实验目的
1.初步了解matlab6.5应用环境;
2.练习利用matlab6.5中的命令或程序描述基本信号;
3.练习利用matlab6.5中的命令或程序实现信号的运算;
4.练习matlab6.5中的绘图方法;
5.初步练习M-File的建立、保存、调用。
二、实验原理与方法
各种信号、运算形式的数学描述及相应的命令或程序。
3.画出以下各序列在给定区间的图形。
1)x(n)=2δ(n+2)-δ(n-4)+δ(n),-5≤n≤5
2)x(n)=n[u(n)-u(n-10)]+10e-0.3(n-10)[u(n-10)-u(n-20)],0≤n≤20
3)x(n)=2-cos(0.4πn)-cos(0.8πn),0≤n≤50
4.令x(n)={1,2,3,4,5,6,7,6,5,4,3,2,1},画出下列序列的图形
1)x1(n)=2x(n-5)-3x(n-4)
2)x2(n)=x(3-n)-x(n)x(n-2)
3)xe(n)和xo(n)---通过编写求xe(n)、xo(n)的M-file来实现(需要给出测试函数)xe(n)和xo(n)分别为x(n)的偶部和奇部
解:
3
1)x(n)=2δ(n+2)-δ(n-4)+δ(n),-5≤n≤5
程序:
clear
clc
n1=-5;n2=5;
n=n1:n2;
x0=[n==0];
x2=[(n+2)==0];
x4=[(n-4)==0];
x2=2x2;
x=x2-x4+x0;
stem(n,x);xlabel(''n'');ylabel(''x'');
holdon
图像:
2)x(n)=n[u(n)-u(n-10)]+10e-0.3(n-10)[u(n-10)-u(n-20)],0≤n≤20
程序:
N=20;
n=0:N;
x1=[n>=0];
x2=[(n>=10)];
x3=[(n-20)>=0];
x4=exp(-0.3).^(n-10);
x=x1-x2+10.x4.(x2-x3);
stem(n,x);xlabel(''n'');ylabel(''x'');
holdon
图像:
3)x(n)=2-cos(0.4πn)-cos(0.8πn),0≤n≤50
程序:
N=50;
n=0:N;
x=2-cos(0.4pin)-cos(0.8pin);
stem(n,x);xlabel(''n'');ylabel(''x'');
holdon
图像:
4:
x1(n)=2x(n-5)-3x(n-4)
程序:
x=[1234567654321];
n=0:12;
[y1,n1]=sigshift(x,n,5);
[y2,n2]=sigshift(x,n,4);
x1=2y1-3y2;
stem(n,x);xlabel(''n'');ylabel(''x'');
holdon
图像:
2)x2(n)=x(3-n)-x(n)x(n-2)
程序:
x=[1234567654321];
n=0:12;
[y1,n1]=sigshift(x,n,-3);
[y2,n2]=sigfold(y1,n1);
[y3,n3]=sigshift(x,n,2);
[y4,n4]=sigmult(x,n,y3,n3);
[x2,n5]=sigadd(y2,n2,-y4,n4);
stem(n5,x2);xlabel(''n'');ylabel(''x'');
holdon
图像:
3)xe(n)和xo(n)---通过编写求xe(n)、xo(n)的M-file来实现(需要给出测试函数)xe(n)和xo(n)分别为x(n)的偶部和奇部
程序:
奇部函数:
function[y1,m1]=jishu(x,n)
fori=1:2:n
ifi<=n
m1((i+1)/2)=(i+1)/2;
y1((i+1)/2)=x(i);
end
end
偶部函数:
function[y1,m1]=oushu(x,n)
fori=2:2:n
ifi<=n
m1((i)/2)=(i)/2;
y1((i)/2)=x(i);
end
end
测试程序:
x=[1,2,3,4,5,6,7,6,5,4,3,2,1];n0=1:13;n=13;
[y1,m1]=jishu(x,n);
[y2,m2]=oushu(x,n);
subplot(311);stem(n0,x);xlabel(''n'');ylabel(''x(n)'');title(''x(n)'');
subplot(312);stem(m1,y1);xlabel(''n'');ylabel(''x(n)'');title(''x(n)'');
subplot(313);stem(m2,y2);xlabel(''n'');ylabel(''x(n)'');title(''x(n)'');
图像
一、实验目的
1.进一步了解matlab6.5应用环境;
2.练习利用matlab6.5中的命令或程序对离散信号和离散时间系统进行分析---单位样值响应、频率响应、零极点增益分析等。
3.进一步练习matlab6.5中的绘图方法;
4.进一步练习matlab6.5中信号的描述及运算方法;
5.进一步练习M-File的建立、保存、调用。
二、实验原理与方法
1.对离散信号和离散时间系统进行分析的有关公式。
7.已知一因果的LTI系统的系统函数为
H(z)=(8-4z-1+11z-2-2z-3)/(1-1.25z-1+0.75z-2-0.125z-3)。
1)绘制它的幅频相频特性;
2)绘制系统的脉冲响应序列;
3)绘制输入x(n)={1,1,1,1}时系统的输出y(n)。
解:
程序:
b=[8-411-2];a=[1-1.250.75-0.125];
x=[zeros(1,10),1111,zeros(1,15)];
n=-10:18;
figure(1);
freqz(b,a,128);title(''上图为幅度响应下图为频率响应'');
figure(2);
impz(b,a,50);title(''系统的脉冲响应序列'');
figure(3);
y=filter(b,a,x);
subplot(121);stem(n,x);title(''输入序列'');
subplot(122);stem(n,y);title(''输出序列'');
图像:
|
|