Forum de mathématiques - Bibm@th.net

\prod {\sum_{n=0}^infty} {\sum_{m=0}^infty} x(m)y(n-m))z^-n

j'ai essayé en Latex,mais j'arrive pas.
désolé

voila les deux codes l'un pour la méthode itérative et l'autre pour la corde ,est ce que c'est juste

function [x, error, total_iters, it_hist] = ...
                     pcgsol(x0, b, atv, params, pcv)
% Preconditioned Conjugate-Gradient solver
%
% C. T. Kelley, December 12, 1993
%
% This code comes with no guarantee or warranty of any kind.
%
% function [x, error, total_iters, it_hist]
%                    = pcgsol(x0, b, atv, params, pcv)
%
%
% Input:        x0=initial iterate
%               b=right hand side
%               atv, a matrix-vector product routine
%            atv must return Ax when x is input
%            the format for atv is
%                       function ax = atv(x)
%               params = two dimensional vector to control iteration
%                        params(1) = relative residual reduction factor
%                        params(2) = max number of iterations
%               pcv, a routine to apply the preconditioner
%                       if omitted, the identity is used.
%                       The format for pcv is
%                       function px = pcv(x).
%
% Output:       x=solution
%               error = vector of iteration residual norms
%               total_iters = number of iterations
%               it_hist (optional) = matrix of all iterations
%            useful for movies
%
%
%

%
% initialization
%
if nargout == 4 it_hist=[]; end
n=length(b); errtol = params(1); maxiters = params(2); error=[]; x=x0;
if nargout == 4; it_hist=[it_hist, x]; end
r=b - feval(atv, x);
if nargin == 4
    z=r;
else
    z = feval(pcv, r);
end
rho=z'*r;
tst=norm(r);
terminate=errtol*norm(b);
error=[error,tst];
it=1;
while((tst > terminate) & (it <= maxiters))
%
%
%
if(it==1)
    p = z;
else
    beta=rho/rhoold;
    p = z + beta*p;
%
% end if
%
end
w = feval(atv, p);
alpha=p'*w;
%
% Test here to make sure the linear transformation is positive definite.
% A non-positive value of alpha is a very bad sign.
%
if(alpha <= 0)
    [alpha, rho, it]
    error(' negative curvature ')
end
alpha=rho/alpha;
x=x+alpha*p;
if nargout == 4; it_hist=[it_hist, x]; end
r = r - alpha*w;
tst=norm(r);
rhoold=rho;
if nargin == 4
    z=r;
else
    z = feval(pcv, r);
end
rho=z'*r;
it=it+1;
error=[error,tst];
%
% end while
%
total_iters=it-1;
end

clc

t=0:0.01:2;

theta1=0.1*cos(sqrt(9.82/1)*t);

y=1-1*cos(theta1);

x=1*sin(theta1);

for i=1:length(t)

cla

hold on

plot([0 0],[-0.2 1],'b--')

fill([-0.1 -0.1 0.1 0.1],[1 1.1 1.1 1],[.8 .8 .8])

fill([0 x(i)],[1 y(i)],[0.8 0.8 0.8])

plot(x(i),y(i),'r.','MarkerSize',45)

axis([-0.5 0.5 -0.2 1.2]);

pause(0.01)

end

Bonsoir

j'attends toujours vos réponses,si quelqu'un possède ces trois programmes ;

merci

merci pour votre réponse

ma question est déja posée,les conditions pour lesquelles les méthodes convergent vers la solution,et comment conditionner la matrice pour que le système d'equations linéaires admet une solution, et finalement essayer des matrices non conditionnées dans le programme en c++ pour voir les difficultés ainsi que pour comparer la rapidité des quatres méthodes

merci encore une fois

bonsoir
ma question c est exactement cela:comment prouver ce resultat de racine p
merci

merci
vraiment j'ai oublié les propriétés des normes
pouvez vous me dire comment est disparu le terme de e^ia/2 pour obtenir 4 cos² (a/2)

merci
j'ai trouvé l'erreur,c 'est que j'ai pas renvoyé le resultat,vers une variable,

det=determinant(C,n);