|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
本帖最后由 eecsseudl 于 2013-4-29 10:01 编辑
教材上用m文件编程:
求LTI系统的零状态响应,已知二阶连续系统的常微分方程表示为:y''+2y'+8y=u
其中y'表示y对时间t求一次导数,u=3t+cos(0.1t);
求该系统的零状态响应。课本上的编程是这样的:
a=input('多项式分母系数向量a=');
b=input('多项式分子系数向量b=');
t=input('输入时间序列 t=');
u=input('输入序列 u=');
tf=t(end);
dt=tf/(length(t)-1);
%用极点留数法求冲击响应
[r,p,k]=residue(b,a);
h=r(1)*exp(p()*t)+r(2)*exp(p(2)*t);
subplot(2,1,1),plot(t,h);grid;
%求u和h的卷积,得输出y(t)
y=conv(u,h)*dt;
subplot(2,1,2),plot(t,y(1:length(t)));grid
上述编程中,小弟我搞不明白为什么y=conv(u,h)*dt,而不是y=conv(u,h),y=conv(u,h)*dt又代表着什么意义呢?请各位大侠指教!
|
|