create equation from polyfit

Question

0 个投票

Hi,

Just got matlab installed. You can only get that far in Excel ;)

The problem I have is that I am not completely into the syntax in matlab yet so I would really appreciate some help.

I have a data set in an csv file which represent a curve. I read them from matlab into two arrays. x's and y's. Then I use polyfit and find a good fit.

The issue I have now is that how do I use the results from polyfit as an equation. First I assumed that I had define a variable x using 'syms x' and then simply use eq1=polyval(p,x); (where p is the array with the constants from the polyfit). This doesnt work. Now I am thinking simply using P*[x^n x^n-1 .. 1]' but this is a bit of a workaround since then I have to create the variable x vector. Is there a simpler way?

Regards

Jakob

2 个评论
显示无隐藏无

Stephen23 2016-2-16

编辑：Stephen23 2016-2-16

What do you want to do with the equation? If you just need to plot the data (or do some other operations) then you do not need the equation, but can use the output of polyfit: the polyfit documentation has examples of how to do this.

Why exactly do you need this equation?

Jakob 2016-2-16

I would like to solve it for different values. In this case x is time and the curve is rising from y=20 to y=60. An example would be that I would like to know how long time it takes to reach y=30.

For this I would use solve(eqn==30,x) and then use this result later.

请先登录，再进行评论。

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

Follow Question

Answer 1

Stephen23 2016-2-16

编辑：Stephen23 2016-2-16

在 MATLAB Online 中打开

0 个投票

A pure numeric solution is very easy using polyval and fzero (no syms needed):

>> X = 1:9;
>> Y = X.^2 + 8*rand(1,numel(X))-4;
>> P = polyfit(X,Y,2);
>> T = fzero(@(x)polyval(P,x)-30,X([1,end]))
T =
       5.3825

And it is easy to plot too:

>> plot(X,Y,'or', X,polyval(P,X),'xb')
>> hold on
>> plot(T([1,1]),get(gca,'Ylim'),'g')
>> plot(X([1,end]),[30,30],'g')
>> legend({'original','polyfit','solved'})

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Jakob 2016-2-16

编辑：Jakob 2016-2-16

ok - thats look nice.

I have though now accepted that maybe I should just use the data set as they are and then dont use polyfit.

Instead I simply use:

[c index]=min(abs(y1-30)) ;

x1(index)

This give the time value that I am interested in.

请先登录，再进行评论。

create equation from polyfit

2 个评论
显示无隐藏无

回答（1 个）

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

类别

标签

Community Treasure Hunt

create equation from polyfit

2 个评论 显示 无 隐藏 无

回答（1 个）

1 个评论 显示 -1更早的评论 隐藏 -1更早的评论

类别

标签

另请参阅

Community Treasure Hunt

2 个评论
显示无隐藏无

1 个评论
显示 -1更早的评论隐藏 -1更早的评论