我正在尝试确定以下系统的稳定性: y(n-2)-4y(n-1)+(λ+6)y(n)-4y(n+1)+y(n+2)=λx(n).以下是我的代码:
function stable2(k)
k =5;
b = [0 0 k 0 0];
a = [1 -4 k+6 -4 1];
zplane(b,a);
p = roots(a);
pm = roundn(abs(p),-4);
if max(pm)>1 %double-sided
disp('System is stable');
else
disp('System is not stable');
end我的问题是:在我的代码中,我只能为λ输入一个值,比如1或1.5。如何将λ作为包含多个值的向量输入,从而可以同时证明具有不同λ的系统的稳定性?
发布于 2021-09-21 13:59:20
eigVec = [1, 2, 3];
for i = 1:length(eigVec)
chckStb(eigVec(i)) % Goes over all eigenvalues
end将你的稳定性检查器组织成这个部分下面的一个函数,它就完成了。
https://stackoverflow.com/questions/69263658
复制相似问题