How to force slope to be zero at a particular point using function PolyFit?

Question

Harshit Jain 2013-4-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/73085-how-to-force-slope-to-be-zero-at-a-particular-point-using-function-polyfit

Q1 = 1xn;

T = 1xn;

Finding a curve fitting (T,Q1) such that slope is zero at Q1(i), where i = 1,2,...,n

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Teja Muppirala 2013-4-22

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/73085-how-to-force-slope-to-be-zero-at-a-particular-point-using-function-polyfit#answer_83072

在 MATLAB Online 中打开

If you have LSQLIN in the Optimization Toolbox, this can be done with a little bit of effort, as described here http://www.mathworks.com/support/solutions/en/data/1-12BBUC/

% Making some random data
Q1 = 0:0.01:1;
T = cos(2.1*pi*Q1)+0.2*randn(size(Q1));
plot(Q1,T,'k.');
% Polyfit without constraints
order = 4;
C1 = polyfit(Q1,T,order);
hold on;
plot(Q1,polyval(C1,Q1))
V = bsxfun(@power,Q1(:),order:-1:0); % Make Vandermonde Matrix
% Make Constraints on the derivatives
Aleft = [(order:-1:1).*Q1(1).^(order-1:-1:0) 0]; 
Aright = [(order:-1:1).*Q1(end).^(order-1:-1:0) 0];
Aeq = [Aleft; Aright];
beq = [0;0]; %Value of the derivative is set to zero
% Call LSQLIN with options to prevent warnings
opts = optimset('lsqlin');
opts.LargeScale = 'off';
C2 = lsqlin(V,T,[],[],Aeq,beq,[],[],[],opts);
plot(Q1,polyval(C2,Q1),'r')
hold off;
legend({'Data' 'Polyfit' 'Constrained Polyfit'},'location','best');

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

Harshit Jain 2013-4-23

thanks for the answer. I have done it through spline. It works fine and is very easy.

请先登录，再进行评论。

Answer 2

Matt J 2013-4-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/73085-how-to-force-slope-to-be-zero-at-a-particular-point-using-function-polyfit#answer_83050

在 MATLAB Online 中打开

If the slope is zero at all i=1...n, it means you are fitting Q1 with a constant.

p=mean(Q1);

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

Harshit Jain 2013-4-22

Slope is not zero at all the points.

Eg - only at i = 1 and n.

请先登录，再进行评论。

Answer 3

Image Analyst 2013-4-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/73085-how-to-force-slope-to-be-zero-at-a-particular-point-using-function-polyfit#answer_83064

编辑：Image Analyst 2013-4-22

Your assignment says nothing about requiring the polyfit() function. There's no way to have the slope be 0 at "1" and "end" in general, unless you use Matt's solution. For specific functions, e.g. 4th order, you may luck out if your data happens to go in between the right elements, e.g. the two humps of a 4th order, but in general that won't happen. So, I'd probably use a spline. You might have to replicate the first and last value of the array to make sure that the slope is zero there.