Discussion :

[cyberkamikaz]

Bonjour,
Je ne vois pas comment résoudre ce problème. Voici le code source:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
function x = householder(A, b)
 
	[m,n] = size(A);
	Q = eye(m);
	for k=1:n
		s=((norm(A(k:m,k))).^2);
		if A(k,k) < 0
			s=-s;
		endif
		v(1:k-1) = 0;
		v(k) = A(k,k)+s;
		A(k,k) = -s;
		v(k+1:m)=A(k+1:m,k);
		A=zeros(k+1:m,k);		
		p=s*v(k);
 
		for j = k+1:n
			t=(v(k:m)'*A(k:m,j))/p;	
			A(k:m,j)=A(k:m,j)-t*v(k:m);
		endfor
		t=((v(k:m)'*b(k:m))/p);
		b(k:m)=b(k:m)-(t*v(k:m));
 
	for j=k+1:n
			t=v(k:m)'*Q(k:m,j)/p;
			Q(k:m,j)=Q(k:m,j)-t*v(k:m);
		endfor
	endfor
	for i=n:-1:1
		s=A(i,i+1:n)*x(i+1:n)
		x(i)=(b(i)-s)/A(i,i)
	endfor
	R=triu(A);
	Q=Q';
endfunction

Et voici l'erreur :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
octave:3> x=householder(A,b)
warning: implicit conversion from range to real scalar
error: A(I): Index exceeds matrix dimension.
error: called from:
error:   /home/Desktop/householder.m at line 18, column 5

Merci d'avance.

[Jean Dumoncel]

Bonjour,

Pour comprendre le message d'erreur :
Index exceeds matrix dimensions.

je ne sais pas ce que tu cherches à faire, mais apparemment, tu travailles par rapport aux dimensions d'un tableau A, et tu modifies les dimensions de ce tableau en faisant

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

A=zeros(k+1:m,k);

[cyberkamikaz]

Merci pour ton aide

ci dessous le m-file correct fonctionel si quelqu'un aurras besoin d'une fonction qui resoud les Ax=b par la methode Householder.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
function x = householder(A, b)
 
	[m,n] = size(A);
	Q = eye(m);
	for k=1:n
		s=(norm(A(k:m,k)));
		if A(k,k) < 0
			s=-s;
		endif
		v(1:k-1) = 0;
		v(k) = A(k,k)+s;
		A(k,k) = -s;
		v(k+1:m)=A(k+1:m,k);
		A(k+1:m,k)=0;		
		p=s*v(k);
 
		for j = k+1:n
			t=v(k:m)*A(k:m,j)/p;	
			A(k:m,j)=A(k:m,j)'-t*v(k:m);
		endfor
		t=((v(k:m)*b(k:m))/p);
		b(k:m)=b(k:m)'-(t*v(k:m));
 
	for j=k+1:n
			t=v(k:m)*Q(k:m,j)/p;
			Q(k:m,j)=Q(k:m,j)'-t*v(k:m);
		endfor
	endfor
	x = zeros(n,1);
	for i=n:-1:1
    	suma=0;
   		for j=(i+1):n
        		suma+=A(i,j)*x(j);
    		endfor
    	x(i)=(b(i)-suma)/A(i,i);
    	endfor
endfunction

Ps: Desole pour les fautes d'ortographe je suis pas francais.