Solve a cubic equation using MATLAB code
显示 更早的评论
I have a cubic equation whose coefficients are varying according to a parameter say w in the following manner:
a=2/w;
b=(3/w+3);
c=(4/(w-9))^3;
d=(5/(w+6))^2;
a*(x^3)+b*(x^2)+c*x+d=0
I want to solve the above equation using a m-file not in the command window. Is it possible ? I am fairly new to MATLAB, so help would be appreciated.
采纳的回答
Matt Fig
2011-2-28
From your comments, it looks like you want this instead:
function [] = solve_cubic()
w=10:2:30;
a=4*((9*w)-7);
b=(-2)*((18*w)-16).*(w+1);
c=3*((w+1).^2).*((3*w)-1);
d=(-3)*((1+w).^3);
f = zeros(1,length(w));
for ii = 1:numel(w)
g = roots([a(ii),b(ii),c(ii),d(ii)]);
f(ii) = g(g<1);
end
ptr=2./(1+w);
u=(1+w)./2;
l=(u.*(1-ptr.*f).^3)./3;
plot(w,l);
5 个评论
Bhagat
2011-2-28
Bhagat
2011-2-28
编辑：Walter Roberson
2020-9-13
Walter Roberson
2011-2-28
Your a, b, c, and d are all row vectors. [a,b,c,d] is then a row vector of size equal to the sum of the sizes of a, b, c, d. Your w vector has 11 elements, so your a, b, c, and d are all 11 long, so the constructed vector is 11+11+11+11 = 44 elements long. And then roots() will calmly tell you the 43 roots of that 44'th degree polynomial.
I told you above that roots() must be applied to each [a(K),b(K),c(K),d(K)] combination individually.
Read the documentation for roots(): it clearly says that the input must be a row vector. No possibility of a matrix of rows to be solved independently is mentioned.
If you choose to believe that roots() might accept a matrix, then go ahead and try
f = roots([a(:),b(:),c(:),d(:)])
as that will at least have one row per value of w.
Also: You are throwing away some of the solutions for f. How many solutions will you have left? If for each w, you do not happen to have exactly one real-valued root with value less than 1, then the size of f will not match the size of the ptr matrix that you are trying to multiply by.
In order to use your formula for "l", for each value of f, you must use the value of w that was used to generate that f. I warned you above about having 0, 1, or 3 roots per w...
Matt Fig
2011-2-28
You keep posting the same erroneous code, but did you try to copy what I wrote and see if it gives you the values you expect? I ran it and it got no errors. If you don't want it as a function, simply delete the first line; it isn't hard to do!
Ubong
2014-9-11
Seen the final response by Matt, that was just so nice . Still have questions for Matt though
1. what is the consequence of not adding the f = zeros(1,length(w)); line? I found that I can still get my row vector, f, without adding it
2. What do I do, if g had returned both real and complex numbers for every ii, and I need to get only real numbers from g
更多回答（4 个）
Bruno Luong
2011-2-26
编辑：Image Analyst
2018-11-16
Please download this FEX,
w = linspace(-3,3,10);
a = 2./w;
b = (3./w+3);
c = (4./(w-9)).^3;
d = (5./(w+6)).^2;
X = CardanRoots([a(:) b(:) c(:) d(:)])
Each row of X contains three solutions of the respective cubic equation. Filter out the complex solutions if you don't need them.
3 个评论
Ann Ytterberg
2018-11-16
Where is this CardonRoots.m?
Image Analyst
2018-11-16
It's the first hit from Google: CardanRoots
Bruno Luong
2018-11-16
Click on the link it should work
Walter Roberson
2011-2-26
EDIT: corrected missing syms, added more detail
function sols = solve_cubic(a, b, c, d)
syms x
sols = solve(a*x^3 + b*x^2 + c*x + d);
end
The inputs to this can include symbolic expressions.
The outputs of this will be symbolic numeric radicals if the inputs are all numeric, but might include symbolic RootOf expressions if any of the arguments are symbolic. To convert symbolic numbers to double precision floating point numbers, use double() on the output.
5 个评论
Bhagat
2011-2-26
编辑：Walter Roberson
2020-9-13
Bhagat
2011-2-26
Walter Roberson
2011-2-26
Sorry I missed a statement before.
Walter Roberson
2011-2-26
Don't use the apostrophes, use the syms command.
Bhagat
2011-2-26
Walter Roberson
2011-2-26
Provided that the input values are scalar numeric:
function sols = solve_cubic(a, b, c, d)
sols = roots([a,b,c,d]);
end
EDIT: provide more detail
The solutions to this will be a column vector of double precision floating point numbers (possibly complex.)
13 个评论
Bhagat
2011-2-26
Walter Roberson
2011-2-26
编辑：Walter Roberson
2020-9-13
Wvals = 1:exp(1):pi^3; %chose appropriate vals
nW = length(Wvals);
out1 = nan(1,nW);
out2 = nan(1,nW);
out3 = nan(1,nW);
for K = 1 : nW
w = Wvals(K);
a=2/w; b=(3/w+3); c=(4/(w-9))^3; d=(5/(w+6))^2;
these_sols = solve_cubic(a,b,c,d);
these_sols = these_sols(imag(these_sols) == 0);
L = length(these_sols);
if L > 0
out1(K) = these_sols(1);
end
if L > 1
out2(K) = these_sols(2);
end
if L > 2
out3(K) = these_sols(3);
end
end
plot(Wvals, out1, '.', Wvals, out2, '.', Wvals, out3, '.');
This will plot up to 3 points per w value, as there can be up to 3 real solutions in general. (I would have to look further to see whether 3 real solutions are actually possible.)
Walter Roberson
2011-2-26
There is an algebraic theorem that any cubic in real coefficients has either one or three real roots, never 0 or 2. You can use that theorem to simplify the above code slightly.
I tried some values myself, and found that indeed for most values of w, there is only one real root. There is a region, though, that is somewhere within -10 to +10, in which there are 3 real roots for each value of w.
Bhagat
2011-2-26
Walter Roberson
2011-2-26
find(Matrix < ParticularNumber)
Bhagat
2011-2-28
Bhagat
2011-2-28
Walter Roberson
2011-2-28
Matrix(Matrix < ParticularNumber)
Bhagat
2011-2-28
Bhagat
2011-2-28
编辑：Walter Roberson
2020-9-13
Bhagat
2011-2-28
Walter Roberson
2011-2-28
编辑：Walter Roberson
2020-9-13
You _are_ putting in w as a constant value: it's value is the entire list 10:2:30 . Using a list like that is not varying w.
Although it is possible to slightly change your a, b, c, and d so that values for each of the w would be calculated, the roots() call will not accept a vector of values, and must be called for each individual a, b, c, d combination.
w0 = 10:2:30;
for K = 1 : length(w0)
w = w0(K);
%calculate a, b, c, d here, and calculate the roots and select those in the proper range and continue on through to the calculation of "l". Do _not_ do the plot() at this time though.
L{K} = l;
end
Now you can plot of w0(K) against L{K}, but keep in mind that you might have no real roots or one real root or three real roots for any individual w, so you will have to build appropriate plotting code.
Bhagat
2011-3-1
Awa Kologo
2020-9-13
0 个投票
- Solving for the roots of a third order polynomial (i.e. finding a value that will make the expression equal zero) may require tedious algebra to do by hand but can be solved easily by a computer using an iterative approach. For example, finding the roots of the expression: , (ie. a value of x so that the equation is satisfied) is time consuming to do by hand. However, plugging in a guess for and then modifying that guess until a tolerance is met gives . Write a MATLAB script that solves exponential equations of the form where are constants that the user inputs in a 1 x 4 row vector: [a b c d]. The program will first check whether or not the input is a 1 x 4 matrix of numerical entries throw an error if it is not. The program will solve for one root of the polynomial iteratively to 6 decimal places and print out the value. You can take this to mean that your program should keep running for as long as Your program will need to make an initial guess for in order to solve iteratively. In most cases, the will give a reasonable first guess. You can use the roots() function to check your answer.
1 个评论
Walter Roberson
2020-9-13
This does not appear to be an Answer to the Question that was asked.
类别
在 帮助中心 和 File Exchange 中查找有关 Common Operations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
