Calling multiple .m files into separate function

2019 3 19
0 个回答
更新时间：2019 3 19
3 次查看（30 天）

0 个投票

function [B, H] = IsStable(polynomial)
if polynomial ~= AllNonZero(polynomial(false))
B=0;
H=[];
elseif polynomial ~= AllSameSign(polynomial(false))
B=0;
H=[];
else
H=HurwitzMatrix(polynomial);
pm = length(H);
for i=1:pm
minor(i)=det(polynomial(1:i,1:i));
end
if minor(i)>0
B=1;
else
B=0;
end
end
I am trying to calculate the stability of a Hurwitz Matrix by checking if the principle minors are more than zero. Any sugestions to improve the above code as it is not working when I call three separate .m files e.g. AllNonZero, AllSameSign and HurwitzMatrix? Thank you.

7 个评论
显示 5更早的评论 隐藏 5更早的评论

Walter Roberson
Walter Roberson 2019-3-19
What are you passing in as the first parameter? If you are not passing in an object or function handle, then polynomial(false) is a logical indexing request to return none of the elements of polynomial so polynomial(false) is likely to be empty.
Mughees Asif
Mughees Asif 2019-3-19
The aim is to enter a row vector into the function, which, will cross check with the other three .m files before displaying 1 if the principle minors are more than zero and 0 if the principle minors are less than 0. When the row vector is entered into the function, it should check if all the signs are same, if any number is 0 and display the associated Hurwitz matrix. But for this file, I need it to check and display all three results from the files and then check for the principle minors.
Mughees Asif
Mughees Asif 2019-3-19
function IsStable(polynomial)
if AllNonZero(polynomial) == false
H=[]
elseif AllSameSign(polynomial) == false
H=[]
else
HurwitzMatrix(polynomial);
pm = length(polynomial);
for i=1:pm
minor(i)=det(polynomial(1:i,1:i));
end
if minor(i)>0
B=1
else
B=0
end
end
I tried changing the parameter but i get the error,
Index in position 1 exceeds array bounds (must not exceed 1).
Error in IsStable (line 10)
minor(i)=det(polynomial(1:i,1:i));
Adam
Adam 2019-3-19
You should presize your array before the for loop:
minor = zeros(1,pm);
Mughees Asif
Mughees Asif 2019-3-19
function IsStable(polynomial)
if AllNonZero(polynomial) == false
H = [];
elseif AllSameSign(polynomial) == false
H = [];
else
HurwitzMatrix(polynomial);
pm = length(polynomial);
minor = zeros(1,pm);
for i = 1:pm
minor(i) = det(polynomial(1:i,1:i));
if minor(i) > 0
B(i) = 1
else
B(i) = 0
end
end
end
I am still getting the same error,
Index in position 1 exceeds array bounds (must not exceed 1).
Error in IsStable (line 11)
minor(i) = det(polynomial(1:i,1:i));
Adam
Adam 2019-3-19
So, as Walter asked, what are you passing in as the parameter?
Walter Roberson
Walter Roberson 2019-3-19
If AllNonZero is false then isn't the matrix all 0 or at least a row of 0?

请先登录，再进行评论。

请先登录，再回答此问题。

回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Matrices and Arrays 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by