calculate the function in vector.

2 次查看（过去 30 天）
JaeSung Choi
JaeSung Choi 2017-6-12
评论： JaeSung Choi 2017-6-13
I want to calculate my 'poly' function for domain of linspace(0,1,100) so I tried ---------------------------------
%make poly function
function [output] = poly(input)
output= ([input^0 input^1 input^2 input^3 input^4 input^5]*transpose([1.0000 1.0001 0.4991 0.1703 0.0349 0.0139]) )
end
----------------------------------
x = linsapce(0,1,100)
poly(x)
----------------------------------
but it doesn't work. I found that for sin(x) it does. I want to know what's different between to func. and how to solve the problem.
  1 个评论
KSSV
KSSV 2017-6-12
What is the input you used? You have to take care of element by element operations.

请先登录，再进行评论。

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

采纳的回答

Andrei Bobrov
Andrei Bobrov 2017-6-13
编辑：Andrei Bobrov 2017-6-13
function [output] = AsPolyvalForJaeSung(input)
output = bsxfun(@power,input(:),0:5)*[1.0000;1.0001;0.4991;0.1703;0.0349;0.0139];
end
  1 个评论
JaeSung Choi
JaeSung Choi 2017-6-13
Oh my god, you exactly catched what i wanted. Thank you very much!!

请先登录，再进行评论。

更多回答（2 个）

KSSV
KSSV 2017-6-12
For
input = linspace(0,1,100) ;
In the line
output= ([input.^0 input.^1 input.^2 input.^3 input.^4 input.^5]*transpose([1.0000 1.0001 0.4991 0.1703 0.0349 0.0139]) )
The size of term in square braces would be 1X600 where as the term transpose i.e second term got only 6X1 terms. How you expect them to multiply? You need to rethink on your code.
  1 个评论
JaeSung Choi
JaeSung Choi 2017-6-12
That's what I'm in problem. I want to derive y = [poly(0) poly(0.01) ...... poly(1)] (i.e. calculate for each domain) As for 'sin' function If we take x = linspace(0,1,100) y = sin(x) then y = [sin(0) sin(0.01) sin(0.02)..... sin(1)] I want to do same for my own function

请先登录，再进行评论。

Torsten
Torsten 2017-6-12
output= ([(input.').^0 (input.').^1 (input.').^2 (input.').^3 (input.').^4 (input.').^5]*([1.0000 1.0001 0.4991 0.1703 0.0349 0.0139]).'
Best wishes
Torsten.
  2 个评论
JaeSung Choi
JaeSung Choi 2017-6-12
I've already tested for the same code. Thanks for your answer but that's not what I needed.
Torsten
Torsten 2017-6-13
??
According to your question, I think this is exactly what you needed.
Best wishes
Torsten.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 MATLAB 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!