How to find only positive root of a polynomial

16 ビュー (過去 30 日間)
Atom
Atom 2013 年 4 月 16 日
How to find only positive root of a polynomial equation x^4+7*x^2-A=0 where A is varying from 1:.1:3. If rr is the positive real root, then find rr/(rr+1) for each case.
for A=1:.1:3
poly = [1 0 7 0 A];
R = roots(poly);
if isa(R,'complex') && (R<=0)
continue;
else
R/(R+1)
end
Please correct the code.
  2 件のコメント
Matt J
Matt J 2013 年 4 月 16 日
What do you mean by "without solving a polynomial equation"? A root is, by definition, the solution to such a problem.
Atom
Atom 2013 年 4 月 16 日
Yes. you are right. Please ignore the words "without solving".

サインインしてコメントする。

採用された回答

Matt J
Matt J 2013 年 4 月 16 日
編集済み: Matt J 2013 年 4 月 16 日
By the quadratic formula, the largest solution for x^2 is
x^2 = (-7+ sqrt(49+4*A))/2
For positive A, this will always be positive. You can then get a positive root for x by doing
x = sqrt( (-7+ sqrt(49+4*A))/2 )
I don't know if this satisfies your requirement "without solving a polynomial equation". It doesn't seem possible that you meant this literally (see my comment above).
  2 件のコメント
Atom
Atom 2013 年 4 月 16 日
Please ignore the words "without solving". Please modify the above code so that I can use it for general one.
Matt J
Matt J 2013 年 4 月 16 日
編集済み: Matt J 2013 年 4 月 16 日
for A=1:.1:3
b=poly(3);
a=poly(1);
R = sqrt( (-b+ sqrt(b^2+4*A*a))/2/a );
R/(R+1)
end

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangePolynomials についてさらに検索

タグ

Community Treasure Hunt

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

Start Hunting!

Translated by