|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
在做Matlab作业时，有一道题：设离散系统可由下列差分方程表示：y(n)-y(n-1)+0.9y(n-2)=f(n)，且f(n)=n^2。试分别绘制n=[0:30]时系统的冲激响应、阶跃响应和零状态响应。希望大神能够附上代码和注释。非常感谢！以下是本人写的代码，但是在运行的过程中出现了报错，我也不知道是什么原因：- n=0:30;%离散系统n的取值范围为0~30
- b=1;%右端系数
- a=[1,-1,0.9];%左端系数
- %构造离散系统方程：y(n)-y(n-1)+0.9y(n-2)=f(n)
- f=n^2;%激励函数
- h=impz(b,a,n);
- %离散系统的单位冲激响应h(k),impz参数中先写右端系数b，再写左端系数a，注意第3个量为序列的取值范围
- g=stepz(b,a,length(n));
- %离散系统的单位阶跃响应g(k),stepz参数中先写右端系数b，再写左端系数a，注意第3个量为序列的长度，调用length()函数
- y=filter(b,a,f);
- %离散系统的零状态响应,filter参数中先写右端系数b，再写左端系数a，第3个量为输入序列f
- subplot(1,2,1);%subplot函数没有画图功能,只是将窗口分割
- stem(n,h,'b');%离散系统的打印函数，自定义蓝色
- title('离散系统的单位冲激响应h(k)');
- xlabel('t'); %添加X轴的标签
- ylabel('y'); %添加Y轴的标签
- legend('h(k)'); %添加图例，即这条曲线是什么函数
- subplot(1,2,2);
- stem(n,g,'r');
- title('离散系统的单位阶跃响应g(k)');
- xlabel('t');
- ylabel('y');
- legend('g(k)');
- figure;
- stem(n,f);
- title('离散系统的零状态响应y_{zs}(t)');
- xlabel('t'); %添加X轴的标签
- ylabel('y'); %添加Y轴的标签
- legend('y_{zs}(t)'); %添加图例，即这条曲线是什么函数
( h% ]9 ]" Q/ R2 r+ i
- y( v) P5 d+ y8 w# p8 j1 U; I' B* }
; y9 c8 A* O- V; x以下是运行时出现的错误程序语句处：f=n^2;%激励函数
. G5 c8 E% M' b! z
命令行报错内容：
1 E! @ B5 u- \" M错误使用 ^8 P* d9 a+ L3 M( C: ^" j
用于对矩阵求幂的维度不正确。请检查并确保矩阵为方阵并且幂为标量。要执行按元素矩阵求幂，请使用 '.^'。$ v7 r5 f- z) u5 l8 x# C6 C; P
出错 Code2_1 (line 5)
2 ] e9 X: Y. H5 C' `) l+ q. l, Hf=n^2;%激励函数$ P+ P; z) x+ C: f
|
|
/1 
发表于 2020-4-16 14:07
收藏
淘帖
支持！
反对！