%%?功能簡介：?獲取用戶輸入的特征方程的系數，
%%????????????構造勞斯判據，進而判斷該系統是否穩定。
%%????????????適用于任意階系統。
%%?作????者：?朱其斯
%%?版????本：?matlab?R2011b

%%?獲取相關數據
%?獲取用戶輸入的特征方程系數
num?=?input（‘Please?input?the?particular?function（example:[1?2?3?4]）:?‘）;
%?獲取特征方程的階數（即length-1）
length?=?max（size（num））;
%?len表示除去勞斯表前兩行后，所余行數，用于確定嵌套for循環中for的次數
len?=?length?-?2;
%?cow表示除去全零列后的列數（為了去掉全零列）
if?mod（length?2）???%?若length為奇數，則cow取第二行列數
????cow?=?（（length?-?1）?/?2）?+?1;
else?????????%?若length為偶數，則cow取第一或第二行列數
????cow?=?length?/?2;
end
%?構造length*length的全零方陣，用于存放勞斯表數據
lsb?=?zeros（length）;

%%?構造勞斯表的第一行，從第一個數開始，每隔一個取一個
k?=?0;?%?用于記錄勞斯表行數，0表示第一行
j?=?1;?%?用于記錄lsb方陣的序號，從1開始，豎著數
for?i?=?1:2:length???????
????lsb（j）?=?num（i）;???%?第一行數據
????j?=?j?+?length;????%?第一行數據所在的序號
end

%%?構造勞斯表的第二行，從第二個數開始，每隔一個取一個
k?=?k?+?1;?%?勞斯表第一行記錄完畢，行數加一，開始記錄第二行
j?=?k?+?1;?%?j=2表示勞斯表第二行第一個數
for?i?=?2:2:length
????lsb（j）?=?num（i）;???%?第二行數據
????j?=?j?+?length;????%?第二行數據所在的序號
end

%%?構造勞斯表的第三行以及以后各行（通過勞斯判據的計算法則獲取數據）
for?m?=?1:l