信号与系统课后matlab作业
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
stem(n,abs(X),'file');
title('幅度谱')
xlabel('nw');
subplot(2,1,2);
stem(n,angle(X),'file');
title('相位谱')
M4-6(4)
x=[1,2,3,0,0];
X=fft(x,5);
subplot(2,1,1);
m=0:4;
sys=tf([2 1],[1 3 2]);
t=ts:dt:te;
x=exp(-3*t);
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];
title('x[-k]')
M3-1
（1）
ts=0;te=5;dt=0.01;
sys=tf([2 1],[1 3 2]);
t=ts:dt:te;
x=exp(-3*t);
y=lsim(sys,x,t);
plot(t,y)
xlabel('Time(sec)')
ylabel('y(t)')
（2）
ts=0;te=5;dt=1;
subplot(2,1,1)
stem(k,x,'file')
title('3倍内插')
t=-1:2;
y=[-2,-2,2,1];
subplot(2,1,2)
stem(t,y,'file')
Baidu Nhomakorabeatitle('3倍抽取')
axis([-3,4,-4,4])
（3）
k=-4:7;
x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1];
图象为：
(5)
A=2;w0=10*pi;phi=pi/6;
t=0:0.001:0.5;
xt=abs(A*sin(w0*t+phi));
plot(t,xt)
图象为：
(6)
A=1;w0=1;B=1;w1=2*pi;
t=0:0.001:20;
xt=A*cos(w0*t)+B*sin(w1*t);
plot(t,xt)
图象为：
(3)
A=10;a=-1;B=5;b=-2;
t=0:0.001:10;
xt=A*exp(a*t)-B*exp(b*t);
plot(t,xt)
图象为：
(4)
t=sym('t');
y=t*Heaviside(t);
ezplot(y,[-1,3]);
grid on
axis([-1 3 -0.1 3.1])
图象为：
M2-3
（1）functionyt=x2_3(t)
yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5);
（2）function yt=x2_3(t)
yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5);
t=0:0.001:6;
title('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')
（2）
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];
title('幅度谱')
xlabel('nw');
subplot(2,1,2);
stem(n,angle(X),'file');
title('相位谱')
（2）对于三角波信号的频谱是：Cn=- + +
n=-15:15;
X=sinc(n)-0.5*((sinc(n/2)).^2);
subplot(2,1,1);
h=impz(b,a,k);
stem(k,h,'file')
M4-1
（1）对于周期矩形信号的傅里叶级数cn =-1/2j*sin(n/2*pi)*sinc(n/2)
n=-15:15;
X=-j*1/2*sin(n/2*pi).*sinc(n/2);
subplot(2,1,1);
stem(n,abs(X),'file');
stem(m,real(X),'file');
title('X[m]实部')
subplot(2,1,2);
stem(m,imag(X),'file');
图象为：
(7)
A=4;a=-0.5;w0=2*pi;
t=0:0.001:10;
xt=A*exp(a*t).*cos(w0*t);
plot(t,xt)
图象为：
(8)
w0=30;
t=-15:0.001:15;
xt=cos(w0*t).*sinc(t/pi);
plot(t,xt)
axis([-15,15,-1.1,1.1])
subplot(2,1,1)
stem(k+2,x,'file')
title('x[k+2]')
subplot(2,1,2)
stem(k-4,x,'file')
title('x[k-4]')
（4）
k=-4:7;
x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1];
stem(-fliplr(k),fliplr(x),'file')
M2-1
(1)
t=-2:0.001:4;
T=2;
xt=rectpuls(t-1,T);
plot(t,xt)
axis([-2,4,-0.5,1.5])
图象为：
(2)
t=sym('t');
y=Heaviside(t);
ezplot(y,[-1,1]);
grid on
axis([-1 1 -0.1 1.1])
kx=-2:3;
h=[0.68,0.37,0.83,0.52,0.71];
kh=-1:3;
y=conv(x,h);
k=kx(1)+kh(1):kx(end)+kh(end);
stem(k,y,'file');
M3-8
k=0:30;
a=[1 0.7 -0.45 -0.6];
b=[0.8 -0.44 0.36 0.02];
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))
title('x(0.5t)')
axis([0,11,-2,3])
subplot(3,1,3)
plot(t,x2_3(2-0.5*t))
title('幅度谱')
xlabel('nw');
subplot(2,1,2);
stem(n,angle(X),'file');
title('相位谱')
M4-6(4)
x=[1,2,3,0,0];
X=fft(x,5);
subplot(2,1,1);
m=0:4;
sys=tf([2 1],[1 3 2]);
t=ts:dt:te;
x=exp(-3*t);
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];
title('x[-k]')
M3-1
（1）
ts=0;te=5;dt=0.01;
sys=tf([2 1],[1 3 2]);
t=ts:dt:te;
x=exp(-3*t);
y=lsim(sys,x,t);
plot(t,y)
xlabel('Time(sec)')
ylabel('y(t)')
（2）
ts=0;te=5;dt=1;
subplot(2,1,1)
stem(k,x,'file')
title('3倍内插')
t=-1:2;
y=[-2,-2,2,1];
subplot(2,1,2)
stem(t,y,'file')
Baidu Nhomakorabeatitle('3倍抽取')
axis([-3,4,-4,4])
（3）
k=-4:7;
x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1];
图象为：
(5)
A=2;w0=10*pi;phi=pi/6;
t=0:0.001:0.5;
xt=abs(A*sin(w0*t+phi));
plot(t,xt)
图象为：
(6)
A=1;w0=1;B=1;w1=2*pi;
t=0:0.001:20;
xt=A*cos(w0*t)+B*sin(w1*t);
plot(t,xt)
图象为：
(3)
A=10;a=-1;B=5;b=-2;
t=0:0.001:10;
xt=A*exp(a*t)-B*exp(b*t);
plot(t,xt)
图象为：
(4)
t=sym('t');
y=t*Heaviside(t);
ezplot(y,[-1,3]);
grid on
axis([-1 3 -0.1 3.1])
图象为：
M2-3
（1）functionyt=x2_3(t)
yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5);
（2）function yt=x2_3(t)
yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5);
t=0:0.001:6;
title('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')
（2）
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];
title('幅度谱')
xlabel('nw');
subplot(2,1,2);
stem(n,angle(X),'file');
title('相位谱')
（2）对于三角波信号的频谱是：Cn=- + +
n=-15:15;
X=sinc(n)-0.5*((sinc(n/2)).^2);
subplot(2,1,1);
h=impz(b,a,k);
stem(k,h,'file')
M4-1
（1）对于周期矩形信号的傅里叶级数cn =-1/2j*sin(n/2*pi)*sinc(n/2)
n=-15:15;
X=-j*1/2*sin(n/2*pi).*sinc(n/2);
subplot(2,1,1);
stem(n,abs(X),'file');
stem(m,real(X),'file');
title('X[m]实部')
subplot(2,1,2);
stem(m,imag(X),'file');
图象为：
(7)
A=4;a=-0.5;w0=2*pi;
t=0:0.001:10;
xt=A*exp(a*t).*cos(w0*t);
plot(t,xt)
图象为：
(8)
w0=30;
t=-15:0.001:15;
xt=cos(w0*t).*sinc(t/pi);
plot(t,xt)
axis([-15,15,-1.1,1.1])
subplot(2,1,1)
stem(k+2,x,'file')
title('x[k+2]')
subplot(2,1,2)
stem(k-4,x,'file')
title('x[k-4]')
（4）
k=-4:7;
x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1];
stem(-fliplr(k),fliplr(x),'file')
M2-1
(1)
t=-2:0.001:4;
T=2;
xt=rectpuls(t-1,T);
plot(t,xt)
axis([-2,4,-0.5,1.5])
图象为：
(2)
t=sym('t');
y=Heaviside(t);
ezplot(y,[-1,1]);
grid on
axis([-1 1 -0.1 1.1])
kx=-2:3;
h=[0.68,0.37,0.83,0.52,0.71];
kh=-1:3;
y=conv(x,h);
k=kx(1)+kh(1):kx(end)+kh(end);
stem(k,y,'file');
M3-8
k=0:30;
a=[1 0.7 -0.45 -0.6];
b=[0.8 -0.44 0.36 0.02];
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))
title('x(0.5t)')
axis([0,11,-2,3])
subplot(3,1,3)
plot(t,x2_3(2-0.5*t))