확률의 안정 벡터

Question

안정 벡터가 발견되는 시간을 얻기 위해 scilab에서 프로그램을 만들려고합니다. (희망적으로는 matlab에도 동일하게 적용됩니다.) 말하자면, 제품 벡터와 매트릭스를 여러 번 만든 후에 그렇게되면 변화하지 않을 것입니다.

function [R]=vector_stable(v,m,i) 
    V=v*m; 
    if(V == v) then 
     R=i; 
     abort; 
    else 
     vector_stable(V,m,i+1); 
    end 
endfunction 

날, V 초기 벡터와 행렬의 제품이 조금 작은 것을 설명하자 경우 :

내가이 작업을 수행하는 가장 좋은 방법은 재귀 함수로 생각 때문에 나는 다음과 같은 코딩 결과는 벡터 매개 변수와 동일하므로 일어난 시간을 반환해야합니다. 그렇지 않으면 첫 번째 매개 변수와 동일한 함수를 결과와 함께 호출합니다. 그러나 나는 뒤에 오는 과실을 얻고있다

-->R=vector_stable(V,M,0) 
!--error 18 
: Too many names. 

나의 기능은 정확한가? 도와 줄수있으세요?

Answer 1

귀하의 기능이 올바르게 보이지 않습니다. 이것은 더 비슷할 수도 있습니다 :

function [R] = vector_stable(v, M, i) 

V = v*M; 
if(norm(V - v) < 0.001) 
    R = i; 
    return; 
else 
    R = vector_stable(V,m,i+1); 
end 

return 
end 

Answer 2

아마 "안정"은 "변화하지 않지만" "수렴"을 의미하지 않을까요? 그렇다면 V == v를 테스트하여 루프를 종료 할 수는 없습니다. 두 벡터 사이의 상대적인 차이를 살펴보고 작도가 작아지면 종료 할 수 있습니다. 1 % 또는 0.1 %.

V == v 대신 i == 10이면 (테스트 목적으로) 종료 될 때 오류가 발생합니까?