信号与系统课后matlab作业Word格式.docx
《信号与系统课后matlab作业Word格式.docx》由会员分享，可在线阅读，更多相关《信号与系统课后matlab作业Word格式.docx（24页珍藏版）》请在冰豆网上搜索。
ezplot（y,[-1,3]）;
axis（[-13-0.13.1]）
（5）
A=2;
w0=10*pi;
phi=pi/6;
0.5;
xt=abs（A*sin（w0*t+phi））;
（6）
A=1;
w0=1;
B=1;
w1=2*pi;
20;
xt=A*cos（w0*t）+B*sin（w1*t）;
（7）
A=4;
a=-0.5;
w0=2*pi;
xt=A*exp（a*t）.*cos（w0*t）;
（8）
w0=30;
t=-15:
15;
xt=cos（w0*t）.*sinc（t/pi）;
axis（[-15,15,-1.1,1.1]）
M2-3
（1）functionyt=x2_3（t）
yt=（t）.*（t>
=0&
t<
=2）+2*（t>
=2&
=3）-1*（t>
=3&
=5）;
（2）functionyt=x2_3（t）
yt=（t）.*（t>
6;
subplot（3,1,1）
plot（t,x2_3（t））
title（'
x（t）'
）
axis（[0,6,-2,3]）
subplot（3,1,2）
plot（t,x2_3（0.5*t））
x（0.5t）'
axis（[0,11,-2,3]）
subplot（3,1,3）
plot（t,x2_3（2-0.5*t））
x（2-0.5t）'
axis（[-6,5,-2,3]）
图像为：
M2-9
（1）
k=-4:
7;
xk=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1];
stem（k,xk,'
file'
k=-12:
21;
x=[-3,0,0,-2,0,0,3,0,0,1,0,0,-2,0,0,-3,0,0,-4,0,0,2,0,0,-1,0,0,4,0,0,1,0,0,-1];
subplot（2,1,1）
stem（k,x,'
3倍内插'
t=-1:
2;
y=[-2,-2,2,1];
subplot（2,1,2）
stem（t,y,'
3倍抽取'
axis（[-3,4,-4,4]）
x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1];
stem（k+2,x,'
x[k+2]'
stem（k-4,x,'
x[k-4]'
stem（-fliplr（k）,fliplr（x）,'
x[-k]'
M3-1
ts=0;
te=5;
dt=0.01;
sys=tf（[21],[132]）;
t=ts:
dt:
te;
x=exp（-3*t）;
y=lsim（sys,x,t）;
plot（t,y）
xlabel（'
Time（sec）'
ylabel（'
y（t）'
dt=1;
y=lsim（sys,x,t）
y=
0
0.6649
-0.0239
-0.0630
-0.0314
-0.0127
从
（1）
（2）对比我们当抽样间隔越小时数值精度越高。
M3-4
x=[0.85,0.53,0.21,0.67,0.84,0.12];
kx=-2:
3;
h=[0.68,0.37,0.83,0.52,0.71];
kh=-1:
y=conv（x,h）;
k=kx
（1）+kh
（1）:
kx（end）+kh（end）;
stem（k,y,'
M3-8
k=0:
30;
a=[10.7-0.45-0.6];
b=[0.8-0.440.360.02];
h=impz（b,a,k）;
stem（k,h,'
M4-1
（1）对于周期矩形信号的傅里叶级数cn=-1/2j*sin（n/2*pi）*sinc（n/2）
n=-15:
X=-j*1/2*sin（n/2*pi）.*sinc（n/2）;
subplot（2,1,1）;
stem（n,abs（X）,'
幅度谱'
nw'
subplot（2,1,2）;
stem（n,angle（X）,'
相位谱'
（2）对于三角波信号的频谱是：
Cn=-
+
X=sinc（n）-0.5*（（sinc（n/2））.^2）;
M4-6（4）
x=[1,2,3,0,0];
X=fft（x,5）;
m=0:
stem（m,real（X）,'
X[m]实部'
stem（m,imag（X）,'
X[m]虚部'
M4-7（3）
x=0.5.^k;
subplot（3,1,1）;
x[k]'
X=fft（x,10）;
subplot（3,1,2）;
9;
subplot（3,1,3）;
M5-2
0.05:
2.5;
T=1;
xt1=rectpuls（t-0.5,T）;
subplot（2,2,1）
plot（t,xt1）
x（t1）'
axis（[0,2.5,0,2]）
xt2=tripuls（t-1,2）;
subplot（2,2,2）
plot（t,xt2）
x（t2）'
xt=xt1+xt2.*cos（50*t）;
subplot（2,2,[3,4]）
figure;
b=[10000];
a=[1,26.131,341.42,2613.1,10000];
w=linspace（0,2*pi,200）;
[H,w]=freqs（b,a,w）;
plot（w,abs（H））;
幅度曲线'
plot（w,angle（H））;
相位曲线'
sys=tf（[10000],[126.131341.422613.110000]）;
yt1=lsim（sys,xt,t）;
plot（t,yt1）;
y（t1）'
yt2=lsim（sys,xt.*cos（50*t）,t）;
plot（t,yt2）;
y（t2）'
M6-1
num=[41.6667];
den=[13.744425.760441.6667];
[r,p,k]=residue（num,den）
r=-0.9361-0.1895i-0.9361+0.1895i1.8722
p=-0.9361+4.6237i-0.9361-4.6237i-1.8722
k=[]
num=[1600];
den=[15.65698162262.7160000];
r=0.0992-1.5147i0.0992+1.5147i-0.0992+1.3137i-0.0992-1.3137i
p=-1.5145+21.4145i-1.5145-21.4145i-1.3140+18.5860i-1.3140-18.5860i
num=[1000];
den=conv（[15],[1525]）;
r=-2.5000-1.4434i-2.5000+1.4434i-5.0000
p=-2.5000+4.3301i-2.5000-4.3301i-5.0000
k=1
num=[833.3025];
den=conv（[14.112328.867],[19.927928.867]）;
r=-2.4819+1.0281i-2.4819-1.0281i2.4819-5.9928i2.4819+5.9928i
p=-2.0562+4.9638i-2.0562-4.9638i-4.9640+2.0558i-4.9640-2.0558i
M6-2
sys=tf（b,a）;
x=1*（t>
0）+0*（t<
=0）;
y1=lsim（sys,x,t）;
plot（t,y1）;
零状态响应'
[ABCD]=tf2ss（b,a）;
sys=ss（A,B,C,D）;
zi=[12];
y=lsim（sys,x,t,zi）;
plot（t,y）;
完全响应'
y2=y-y1;
plot（t,y2）;
零输入响应'
M6-5
a=[1221];
b=[12];
pzmap（sys）
冲击响应
impulse（sys）
阶跃响应
step（sys）
频率响应
w=linspace（-10,10,20000）;
H=freqs（b,a,w）;
plot（w,abs（H））
M7-1
num=[216445632];
den=[33-1518-12];
[r,p,k]=residuez（num,den）
r=-0.01779.4914-3.0702+2.3398i-3.0702-2.3398i
p=-3.23611.23610.5000+0.8660i0.5000-0.8660i
k=-2.6667
num=[4-8.86-17.9826.74-8.04];
den=[1-210665];
r=1.0849+1.3745i1.0849-1.3745i0.9769-1.2503i0.9769+1.2503i
p=2.0000+3.0000i2.0000-3.0000i-1.0000+2.0000i-1.0000-2.0000i
k=-0.1237
M7-2
a=[2-13];
b=[2-1];
y1=filter（b,a,x）;
stem（k,y1）
y2=conv（x,h）;
N=length（y2）;
stem（0:
N-1,y2）
完全响应
y=y1+y2
M7-3
zplane（num,den）
由图可知极点不是都在单位圆内，所以系统不是稳定系统。
num=[4-8.68-17.9826.74-8.04];
升级会员 