武汉理工大学信号与系统实验答案Word文档格式.docx
- 文档编号:8253338
- 上传时间:2023-05-10
- 格式:DOCX
- 页数:12
- 大小:97.50KB
武汉理工大学信号与系统实验答案Word文档格式.docx
《武汉理工大学信号与系统实验答案Word文档格式.docx》由会员分享,可在线阅读,更多相关《武汉理工大学信号与系统实验答案Word文档格式.docx(12页珍藏版)》请在冰点文库上搜索。
y=exp(-0.5*(1.5*t+3)).*u(1.5*t+3);
subplot(211)
x(t)'
subplot(212)
plot(t,y)
y(t)=x(1,5*t+3)'
1-3
clear;
closeall;
t0=0;
t1=3.5;
dt=0.001;
t=t0:
t1;
x1=u(t)-u(t-2);
x2=u(t-1)-u(t-3);
y=dt*conv(x1,x2);
subplot(221)
plot(t,x1)
gridon
Signalx1(t)'
axis([t0,t1,0,1.5])
subplot(222)
plot(t,x2)
signalx2(t)'
2*t1;
Theconvolutionofx1(t)andx2(t)'
axis([2*t0,2*t1,0,2.5])
+
2-1
num=input('
Typeintherightcoeffcientvectorofdifferentialequation:
'
);
den=input('
Typeintheleftcoefficientvectorofdifferentialequation:
t=0:
0.01:
10;
x=input('
Typeintheexpressionoftheinputsignalx(t):
lsim(num,den,x,t)
运行过程:
8
[1,3,2]
(exp(-2*t)-exp(-3*t)).*u(t)
3-1
T=2;
dt=0.00001;
x1=(u(t+1)-u(t-1-dt)).*abs(t);
x=0;
form=-1:
1
x=x+u(t-m*T)-u(t-1-m*T-dt);
end
w0=2*pi/T;
N=10;
L=2*N+1;
fork=-N:
N;
ak(N+1+k)=(1/T)*x1*exp(-j*k*w0*t'
)*dt;
phi=angle(ak);
运行结果:
ak=
Columns1through6
0.0000-0.0000i-0.0025-0.0000i0.0000+0.0000i-0.0041+0.0000i
0.0000+0.0000i-0.0081-0.0000i
Columns7through12
0.0000-0.0000i-0.0225+0.0000i0.0000-0.0000i-0.2026+0.0000i
0.5000-0.2026-0.0000i
Columns13through18
0.0000+0.0000i-0.0225-0.0000i0.0000+0.0000i-0.0081+0.0000i
0.0000-0.0000i-0.0041-0.0000i
Columns19through21
0.0000-0.0000i-0.0025+0.0000i0.0000+0.0000i
3-2
b=input('
b='
a=input('
a='
[H,w]=freqs(b,a);
Hm=abs(H);
phai=angle(H);
hr=real(H);
Hi=imag(H);
plot(w,Hm)
Magnituderesponse'
Frequencyinrad/sec'
subplot(223)
plot(w,phai)
Phaseresponse'
plot(w,hr),
Realpartoffrequencyresponse'
subplot(224)
plot(w,Hi),
Imaginarypartoffrequencyresponse'
1.b=[1,0]
a=[1,1,25]
2.
b=[1,-1]
a=[1,1]
3.
b=[262]
a=[1,10,48,148,306,401,262]
4-1
a=-0:
0.1:
b=-20:
20;
[a,b]=meshgrid(a,b);
c=a+i*b;
c=(1-exp(-2*(c+eps)))./(c+eps);
c=abs(c);
mesh(a,b,c);
surf(a,b,c);
view(-60,20)
axis([0,5,-20,20,0,2]);
TheLaplacetransformoftherectangularpulse'
w=-20:
Fw=(2*sin(w+eps).*exp(i*(w+eps)))./(w+eps);
plot(w,abs(Fw))
TheFouriertransformoftherectangularpulse'
xlabel('
Frequencyw'
4-2
b=[0,1/16,0,1];
a=[1,2,2,1];
splane(b,a)
Thezero-polediagram'
附:
functionsplane(num,den)
p=roots(den);
q=roots(num);
p=p'
;
q=q'
x=max(abs([p,q]));
x=x+1;
y=x;
plot([-x,x],[0,0],'
:
holdon;
plot([0,0],[-y,y],'
plot(real(p),imag(p),'
x'
plot(real(q),imag(q),'
o'
zero-poleplot'
Realpart'
ylabel('
Imaginalpart'
axis([-x,x,-y,y]);
(a=0.6)(a=16)
(a=0.8)(a=4)
5-2
tmax=4;
tmax;
Ts=1/10;
ws=2*pi/Ts;
w0=20*pi;
dw=0.1;
w=-w0:
dw:
w0;
n=0:
1:
tmax/Ts;
x=exp(-4*t).*u(t);
xn=exp(-4*n*Ts);
plot(t,x);
Acontinuous-timesingnalx(t)'
Timet'
axis([0,tmax,0,1]);
stem(n,xn,'
.'
Thesampledversionx[n]ofx(t)'
Timeindexn'
axis([0,tmax/Ts,0,1]);
Xa=x*exp(-j*t'
*w)*dt;
X=0;
fork=-8.8;
X=X+x*exp(-j*t'
*(w-k*ws))*dt;
plot(w,abs(Xa))
Magnitudespctrumofx(t)'
gridon;
axis([-60,60,0,1.8*max(abs(Xa))])
plot(w,abs(X))
Magnitudespctrumofx(n)'
Frequencyinradians/s'
(Ts=1/2)(Ts=1/5)
(Ts=1/7)(Ts=1/10)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 武汉理工大学 信号 系统 实验 答案