다항식 선에 대한 방정식을 계산하는 방법

Question

나는 다항식 y = 0.230 + -0.046*x + -0.208*x^2을가집니다. (X, Y)에서 다른 선을 자르는이 선에 대한 수직선을 계산하고 싶습니다.

Answer 1

%Example data 
    x=0:0.1:10; 
    y = 0.230 + -0.046*x + -0.208*x.^2 ;  
    plot(x,y); 

%Find the tangent and normals at all points 
    dy = [0 diff(y)./diff(x)]; 
    py = -1./dy; 

%Choose a point 
    n = 60; 
    X = x(n); 
    Y = y(n); 
    hold on 
    plot(X, Y, 'or') 

%Find the equation of the straight line normal to that point. You can do this in one step (yn = py(n)*(x - X) + Y) but I've done it in two to illustrate where this comes from 
    c = Y - py(n)*X; 
    yn = py(n)*x + c; 
    plot(x, yn, 'g') 

Answer 2

다른 방법으로는 그리 어려운 것은 아니지만 분석 결과를 계산하는 것입니다. 는 (당신은을위한 상징적 인 도구 상자를 사용할 수 있지만, 당신의 머리에 앉아 NN 할 것) :

%Example data 
x=0:0.1:10; 
y = 0.230 + -0.046*x + -0.208*x.^2 ; 
plot(x,y); 

%Find the tangent and normals at all points (*edited*) 
slope = (-0.046 + -2*0.208*x); 
py = -1./slope;   % <-- modified from Dan's expression 
          %  to use analytical derivative 


%Choose a point 
n = 60; 
X = x(n); 
Y = y(n); 
hold on 
plot(X, Y, 'or') 

% Copying @Dan: Find the equation of the straight line normal to that point. You can do this in one step (yn = py(n)*(x - X) + Y) but I've done it in two to illustrate where this comes from 
c = Y - py(n)*X; 
yn = py(n)*x + c; 
plot(x, yn, 'g') 
axis tight equal 

당신이 정말로 직교 곡선을 가지고 있는지 axis equal 또한이 예에서는 좋은 생각이다 사용. 의