Discussion :

[junkie1986]

Bonjour,

j´ai écrit le code suivant pour résoudre une équation basique. Mais le fait est qu´il est extrêmement lent. Quelqu´un aurait-il une solution pour l´optimiser ?

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
function [ result ] = SolveEq( data )
%SOLVEEQ Summary of this function goes here
%   Detailed explanation goes here
syms x
[l,~] = size(data);
result = zeros(l,1);
 
for i = 1:1:l
result(i,1) = solve(sprintf('(x+50)/(x-50)=%.5d',data(i,1)));
end
 
end

D´autre part mon principal probleme est que le vecteur data contient un nombre complexe et si j´essaye ce code là :

Code :

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

result(i,1) = solve(sprintf('(x+50)/(x-50)=%8.3f + %8.3f*1i',real(data(i,1)),imag(data(i,1))));

j´obtiens l´erreur suivante:

Code :

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

1
2
3
4
5
??? Error using ==> solve>getEqns at 178
' (x+50)/(x-50)=   0.927 +    0.148*1i ' is not a valid expression or equation.
 
Error in ==> solve at 67
[eqns,vars] = getEqns(varargin{:});

La fonction solve ne peut elle résoudre des équations complexes ?


MATLAB Version 7.9.0.529 (R2009b)

[phryte]

Bonjour,
Exemple de calcul (4s) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
clear
data=10*ones(1,100)';
data(5,1)=(50+20i)';%*ones(1,100);
syms x
[l,~] = size(data);
result = zeros(l,1);
for i = 1:1:l
 result(i,1) = solve(sprintf('(x+50)/(x-50)=%.5d',data(i,1))); 
 % ou
%result(i,1) = solve(sprintf('(x+50)/(x-50)=%8.3f',real(data(i,1))+imag(data(i,1))*j));
end
 result

[junkie1986]

Bonjour phryte,

merci d´avoir pris le temps de me répondre. Le fait est que data contient environs 10k lignes et le calcul prend quelques minutes.
Je me demandais plus si on pouvait par exemple empêcher solve, paramétrer cette fonction pour lui demander de retourner moins de chiffres
significatifs après la virgule. Ou bien encore utiliser quelque chose de ce genre :

http://stackoverflow.com/questions/1...m-of-equations

Sinon je vais prendre mon mal en patience ....

[phryte]

Bonjour,
Si tu as MuPad avec Matlab tu peux essayer.