Plotting a function (Horner Schema)
显示 更早的评论
I am sorry if this question has already been asked, but I am totally new to Matlab and Programming but I need it for my class at university... So the task is to programm the "Horner Schema". This is what I got:
function y = horner(a,x)
n = length(a)-1;
y = zeros(length(x));
for i = n:-1:1
y = y(i)*x+a(i);
end
end
It seems toi work and gives back no errors to me but I don't know, if it gives back the right array.
Anyway, my problem lies within the fact that we should also analyse 
in the interval [-1,1] and I have absolutely no idea how to do this.
in the interval [-1,1] and I have absolutely no idea how to do this.I am glad for any help that you can give me.
5 个评论
John D'Errico
2020-7-9
What does analyze mean to you? For example, I tried to get that polynomial to lay down on a couch and discuss its feelings, but I've gotten nothing from it.
Luisa Maria Korb
2020-7-9
Luisa Maria Korb
2020-7-9
madhan ravi
2020-7-9
Ask your teacher/professor what they want you to do, seems like you’re not sure what you have to do.
采纳的回答
John D'Errico
2020-7-9
编辑：John D'Errico
2020-7-9
While polynomic psychological analysis seems well beyond my talents, Horner's rule is pretty easy. :) And since you made a credible effort, I'll look at your code. The solution is easier than you think. You were pretty close, though with a couple of significant mistakes.
It looks like you want to evaluate a polynomial with coeffiicents provided in a, at a vector of points in x. So you preallocated y. That is good. And you looped over the coefficients in a. But first, what does Horner's rule mean?
Given a polynomial represented by coefficients in a, such that
a(n+1)*x^n + a(n)*x^n-1 + a(n-1)*x^n-2 + ... + a(1)
then we can write the polynomial as
a(1) + x*(a(2) + x*(a(3) + x*(a(4) + x*(..... ))))))
Horner's method starts in the very deepest part of the parens in that expression. That is, it starts with
y = a(n+1)
Then you multiply by x, and add the preceding coefficient. Repeat until you are done.
(I will point out that this is the opposite sequence the coefficients are stored for tools like polyfit and polyval., but that is not truly relevant here.)
I've carefully called the function myhorner, since there is also a function called horner in the symbolc toolbox. You can name it as you wish.
function y = myhorner(a,x)
% Horner's method to evaluate a polynomial
% a contains coefficient of the polynomial, stored in increasing order of the power of x.
% x may be a scalar, vector, or array of any size or shape.
% what degree polynomial is this?
n = length(a)-1;
% preallocate y to be the same shape and size as x, but
% initialized to contain copies of a(n+1). repmat serves
% this purpose this perfectly.
y = repmat(a(n+1),size(x));
for i = n:-1:1
% Note use of .* to multiply by x. Recall that y is potentially a vector
% or array, of the same shape and size as x. You wish to multiply every
% element of y by the corresponding element of x.
y = y.*x + a(i);
end
end
I've added comments in there to explain what I did, however the code itself is still very much yours in its flow and methodology.
Finally, some people might recommend not using i as a loop index. I sort of agee with them, as it then makes creating complex numbers more difficult at some point for you. At the same time, I come from an old school of fortran coders where i was always the loop index in all loops, so I am willing to forgive this use. :)
A quick test of the code is in order. First, does it do what we expected in terms of symbolic manipulations?
a = sym('a',[1,5])
a =
[ a1, a2, a3, a4, a5]
syms x
y = myhorner(a,x)
y =
a1 + x*(a2 + x*(a3 + x*(a4 + a5*x)))
expand(y)
ans =
a5*x^4 + a4*x^3 + a3*x^2 + a2*x + a1
Yes. That is indeed perfect.
A numerical test seems important too.
myhorner([1 2 3],[-1 0 1 5])
ans =
2 1 6 86
You can check those values. Or, we can flip the coefficients around, and use polyval as a check.
polyval([3 2 1],[-1 0 1 5])
ans =
2 1 6 86
And, finally, we can plot the function with those coefficients.
fplot(@(x) myhorner([1 2 3],x),[-5,5])
That looks quite reasonably quadratic to me.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Text Data Preparation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
