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])

本文来源：https://www.bwwdw.com/article/h9kd.html