信号与系统课后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];