Matlab制作朱利表

朱利判据

在这里插入图片描述
在这里插入图片描述
其中
{ b n − k = a n − k − a n a 0 ∗ a k c n − k = b n − k − b n b 0 ∗ b k . . . q n − k = p n − k − p n p 0 ∗ p k \begin{cases} b_{n-k}=a_{n-k}-\frac{a_n}{a_0}*a_k\\ c_{n-k}=b_{n-k}-\frac{b_n}{b_0}*b_k\\ ...\\ q_{n-k}=p_{n-k}-\frac{p_n}{p_0}*p_k \end{cases} bnk=anka0anakcnk=bnkb0bnbk...qnk=pnkp0pnpk
系统稳定的充要条件就是 a 0 , b 0 , . . . , q 0 a_0,b_0,...,q_0 a0,b0,...,q0均为正。
可以用Matlab直接在工作区生成朱利表,需要先生成一个带交互式界面的M文件,它的代码主要依赖于Matlab本身的程序语言控制,M文件如下:

disp('利用朱利判据判断离散系统的稳定性')
disp('请输入特征方程的系数矩阵A')
A=input('A=');
n=size(A,2);%截取A的维数
for i=n:-1:1
    disp(A);%显示奇数行
    if i==1 
        if A(:)>0
            disp('系统稳定');
        else
            disp('系统不稳定');
        end
        break;
    end
    B=flip(A);
    disp(B);%显示偶数行
    if A(1,1)<=0
        disp('首元素非正!系统不稳定');
        break;
    end
    k=B(1,1)/A(1,1);
    A(1,:)=A(1,:)-k*B(1,:);
    A(:,i)=[];%A矩阵减维数
end

保存M文件命名为M1.m
执行的效果如下:
朱利判据效果图

  • 2
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

非线性光学元件

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值