matlab部分程序
1. 用MATLAB命令产生如下信号,并绘出波形图 (1) 2e?1.5tu(t)
t=-2*pi:pi/50:2*pi; a=t>=0;
f(a)=2*exp(-1.5*t(a)); f(~a)=0; plot(t,f)
(2)2sin(2?t??4)
t=-2*pi:pi/50:2*pi; y=2*sin(2*pi*t+pi/4); plot(t,y)
(3) f(t)???2?00?t?1 t?0,t?1t=-5:0.01:5;
y=2*rectpuls(t-0.5,1); plot(t,y) axis equal
2. 已知信号f(t)的波形如图所示,试用MATLAB命令画出f(t?2)、 f(3t)、f(?t)、f(?3t?2)的波形图。 f(t)1?201t
function y=p(t) y1=(t<0&t>=-2); y2=(-t+1).*(t<=1&t>=0); y=y1+y2
t=-5:0.01:5; y=p(t-2); plot(t,y) axis equal
t=-5:0.01:5; y=p(3*t); plot(t,y) axis equal
t=-5:0.01:5; y=p(-t); plot(t,y) axis equal
t=-5:0.01:5; y=p(-3*t-2); plot(t,y) axis equal
3.已知描述系统的微分方程和激励信号如下 y??(t)?4y?(t)?4y(t)?f?(t)?3f(t),f(t)?e?tu(t)
要求:(1)从理论上求解系统的冲激响应和零状态响应,并根据求解结果用MATLAB绘制其时域波形;
(2)分别用MATLAB的impulse()函数和lsim()函数绘制系统的冲激响应和零状态响应,验证(1)中的结果。 clf;
a=[1 4 4]; b=[0 1 3]; subplot(2,1,1); impulse(b,a) t=0:0.01:3.5;
h=(1+t).*exp(-2*t); subplot(2,1,2); plot(t,h) axis([0,3.5,0,1])
figure(2) t=0:0.01:3.5; x=exp(-t); subplot(2,1,1); lsim(b,a,x,t) y=2*exp(-t)-(2+t).*exp(-2*t); subplot(2,1,2); plot(t,y) axis([0,3.5,0,1])
感谢您的阅读,祝您生活愉快。