Double intégrale de la fonction d'objectif

Question

Je souhaite évaluer l'intégrale double de ma fonction d'objectif (nommée myfunction, voir ci-dessous) en utilisant la fonction de construction integral2.

J'ai essayé d'exécuter ce script;

f = @(r,theta) myfunction(r,theta); 
integral2(f,0,2,0,2*pi); 

où myfunction est la fonction suivante:

function fkt=myfunction(r,theta) 
x=r.*cos(theta); 
y=r.*sin(theta); 
los(:,1)=x; 
los(:,2)=y; 
norm = (sum(sqrt(los(:,1).^2 + los(:,2).^2))); 
fkt=norm*r; 
end 

Je fais l'intégrale en coordonnées polaires, que la raison pour laquelle FKT = norme * r.

Matlab me donne le message d'erreur suivant:

>> untitled2 
Subscripted assignment dimension mismatch. 

Error in myfunction (line 8) 
los(:,1)=x; 

Je ne peux pas comprendre, quel est le problème.

Answer 1

Il y a deux choses qui peuvent être améliorées:

1. los est non défini
2. los(:,1) est une colonne en x est une ligne, de sorte que la cession a à l'échec.

Vous pouvez corriger cela en définissant los et modifier votre affectation. Par exemple:

los = NaN(numel(r), 2); 
los(:,1) = x'; 
los(:,2) = y'; 

Mais, pourquoi avez-vous besoin de la variable los? Il suffit de l'enlever et l'erreur sera disparu:

function fkt=myfunction(r,theta) 
x=r.*cos(theta); 
y=r.*sin(theta); 
norm = (sum(sqrt(x.^2 + y.^2))); 
fkt=norm*r; 
end 

Best,

Answer 2

x est une matrice et vous essayez de l'affecter à un vecteur de colonne.