信号与系统实验第八章.docx
- 文档编号:24812921
- 上传时间:2023-06-01
- 格式:DOCX
- 页数:14
- 大小:97.24KB
信号与系统实验第八章.docx
《信号与系统实验第八章.docx》由会员分享,可在线阅读,更多相关《信号与系统实验第八章.docx(14页珍藏版)》请在冰豆网上搜索。
信号与系统实验第八章
8.1
(2)
a=[142];
b=[103];
sys=tf(b,a);
p=0.01;
t=0:
p:
10;
f=exp(-t);
lsim(sys,f,t);
8.4
(2)
a=[11/2];
b=[12];
n=0:
20;
x=2*cos(n*pi/3);
y=filter(b,a,x)
subplot(2,1,1)
stem(n,x,'filled')
title('输入序列x(n)')
subplot(2,1,2)
stem(n,y,'filled')
title('输出序列x(n)')
y=
Columns1through11
2.00004.0000-1.0000-3.5000-3.25000.62503.68753.1563-0.5781-3.7109-3.1445
Columns12through21
0.57233.71393.1431-0.5715-3.7142-3.14290.57143.71433.1429-0.5714
8.6(3)
a=[111/4];
b=[1];
n=[0:
20]
[H,T]=impz(b,a,0:
20)
y=step(b,a,n)
subplot(2,1,1)
impz(b,a,0:
20)
title('单位序列响应')
subplot(2,1,2)
stem(n,y,'filled')
title('阶跃响应')
n=
Columns1through15
01234567891011121314
Columns16through21
151617181920
H=
1.0000
-1.0000
0.7500
-0.5000
0.3125
-0.1875
0.1094
-0.0625
0.0352
-0.0195
0.0107
-0.0059
0.0032
-0.0017
0.0009
-0.0005
0.0003
-0.0001
0.0001
-0.0000
0.0000
T=
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
y=
0
0.3608
1.0570
1.7687
2.3760
2.8508
3.2034
3.4564
3.6337
3.7556
3.8383
3.8937
3.9306
3.9549
3.9708
3.9812
3.9879
3.9923
3.9951
3.9969
3.9980
时间n<0时,系统单位序列响应h(n)=0,故系统为因果系统。
另外,系统单位序列响应h(n)是随时间增长而衰减的振荡序列,当n趋近于无穷时,h(n)趋近于4,即h(n)满足绝对可和条件,因此该系统是一个稳定系统。
8.7(a)
h=ones(1,4)
x=[1,2,1]
X=conv(h,x)
n1=0:
3
n2=-1:
1
n=(n1
(1)+n2
(1)):
(n1
(1)+n2
(1)+length(n1)+length(n2)-2)
subplot(3,1,1)
stem(n1,h,'filled')
title('h(n)')
xlabel('n')
axis([-1,4,0,1])
subplot(3,1,2)
stem(n2,x,'filled')
title('x(n)')
xlabel('n')
axis([-2,2,0,2])
subplot(3,1,3)
stem(n,X,'filled')
title('X(n)=h(n)*x(n)')
xlabel('n')
axis([-2,5,0,5])
h=
1111
x=
121
X=
134431
n1=
0123
n2=
-101
n=
-101234
8.8
function[x,n]=gghconv(x1,x2,n1,n2)
%计算序列x1与x2的卷积和实用函数
x=conv(x1,x2)%计算序列x1与x2的卷积和并赋值给向量x
ns=n1
(1)+n2
(1);%计算两序列卷积和序列x非零样值的起点位置
leg=length(x1)+length(x2)-2;%计算卷积和序列x非零样值的时域宽度
n=ns:
(ns+leg)%求出卷积和序列x非零样值的时间序号向量
subplot(2,2,1)
stem(n1,x1,'filled')
title('x1(n)')
xlabel('n')
subplot(2,2,2)
stem(n2,x2,'filled')
title('x2(n)')
xlabel('n')
subplot(2,2,3)
stem(n,x,'filled')
title('x(n)=x1(n)*x2(n)')
xlabel('n')
p=get(gca,'position');%获取当前坐标轴对象位置
p(3)=2.4*p(3);
set(gca,'position',p)%将子图3的横坐标范围扩为原来的2.4倍
(2)
n2=-2:
2;
x2=ones(1,5);
n3=0:
2;
x3=[321];
[x,n]=gghconv(x2,x3,n2,n3)
x=
3566631
n=
-2-101234
(3)
n1=0:
2;
x1=[321];
n2=0:
3;
x2=[1-11-1];
[x,n]=gghconv(x1,x2,n1,n2)
x=
3-12-2-1-1
n=
012345
8.9
function[f,t]=gggfconv(f1,f2,t1,t2)
%计算连续信号卷积积分实用函数
d=input('');
f=conv(f1,f2);%计算序列f1与f2的卷积和
f=f*d;%计算卷积积分信号f(t)离散样值
ts=t1
(1)+t2
(1);%计算序列f非零样值的起点位置
l=length(f1)+length(f2)-2;%计算卷积积分f的非零样值的宽度
t=ts:
d:
(ts+l*d);%确定卷积积分f非零样值的时间向量
subplot(2,2,1)
plot(t1,f1)
axis([min(t1),max(t1),min(f1)-min(f1)*0.2,max(f1)+max(f1)*0.2])
title('f1(t)')
xlabel('t')
subplot(2,2,2)
plot(t2,f2)
axis([min(t2),max(t2),min(f2)-abs(min(f2)*0.2),max(f2)+max(f2)*0.2])
title('f2(t)')
xlabel('t')
subplot(2,2,3)
plot(t,f)
axis([min(t),max(t),min(f)-min(f)*0.2,max(f)+max(f)*0.2])
p=get(gca,'position');
p(3)=2.4*p(3);
set(gca,'position',p)
title('f(t)=f1(t)*f2(t)')
xlabel('t')
(3)
t1=0:
0.01:
4;
f1=Heaviside(t1-1)-Heaviside(t1-3);
t2=-1:
0.01:
3;
f2=0.5*t2.*(Heaviside(t2)-Heaviside(t2-2));
[t,f]=gggfconv(f1,f2,t1,t2)
(5)
t1=-1:
0.01:
3;
f1=0.5*t1.*(Heaviside(t1)-Heaviside(t1-2));
t2=-3:
0.01:
3;
f2=0.5*(t2+2).*(Heaviside(t2+2)-Heaviside(t2))+0.5*(2-t2).*(Heaviside(t2)-Heaviside(t2-2));
%定义信号,向左平移2个单位后再翻转
[t,f]=gggfconv(f1,f2,t1,t2)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 信号 系统 实验 第八
![提示](https://static.bdocx.com/images/bang_tan.gif)