Defining an anonymous function with 3*(N+1) variables

2 visualizaciones (últimos 30 días)
Anthony Gurunian
Anthony Gurunian el 8 de Abr. de 2021
Comentada: Anthony Gurunian el 10 de Abr. de 2021
I need to define a function
f = @(X)...
where X has 3*(N+1) variables in the form: X(:,1), X(:,2), ... X(:,3*(N+1)) . The problem is that the expression for the function is very long. However, there is a pattern which I want to take advantage of.
The function is
Notice that the function can be expressed in terms of summations and products. I want to take advantage of this fact so that I don't have to manually type out the products and sums. Is there any way of doing this using for loops?
  4 comentarios
Cris LaPierre
Cris LaPierre el 8 de Abr. de 2021
Perhaps I'm getting hung up on symantics here, but it appears you have a single variable, X. Is it correct to say that your array X has 3*(N+1) columns?
Anthony Gurunian
Anthony Gurunian el 8 de Abr. de 2021
Yes that is correct. I'm trying to use the syntax required here: https://www.mathworks.com/matlabcentral/fileexchange/53477-monte-carlo-integration
" f is a vectorized function a la f=@(X)sin(X(:,1))+cos(X(:,2))"

Iniciar sesión para comentar.

Respuesta aceptada

David Hill
David Hill el 8 de Abr. de 2021
f=prod(exp(-beta*(diff(x(1:n+1)).^2+diff(x(n+2:2*n+2)).^2+diff(x(2*n+3:3*n+3)).^2)))*sum((x(1:n+1)-x(1:n+1)').^2+(x(n+2:2*n+2)-x(n+2:2*n+2)').^2+(x(2*n+3:3*n+3)-x(2*n+3:3*n+3)).^2,'all');
  7 comentarios
David Hill
David Hill el 9 de Abr. de 2021
Editada: David Hill el 9 de Abr. de 2021
Is there a different f for each row? If not, I don't understand your equation. If so, then just do a loop with my equation. How can you not know the number of rows? X must be defined before this equation is processed.
for k=1:size(x,1)
f(k)=prod(exp(-beta*(diff(x(k,1:n+1)).^2+diff(x(k,n+2:2*n+2)).^2+diff(k,x(2*n+3:3*n+3)).^2)))*...
sum((x(k,1:n+1)-x(k,1:n+1)').^2+(x(k,n+2:2*n+2)-x(k,n+2:2*n+2)').^2+...
(x(k,2*n+3:3*n+3)-x(k,2*n+3:3*n+3)).^2,'all');
end
Or,
for k=1:size(x,1)
X=reshape(x(k,:),[],3);
f(k)=prod(exp(-beta*(sum(diff(X,1).^2),2)))*sum((X(:,1)-X(:,1)').^2+...
(X(:,2)-X(:,2)').^2+(X(:,3)-X(:,3)').^2,'all');
end
Anthony Gurunian
Anthony Gurunian el 10 de Abr. de 2021
I ended up writing a custom function f = Rg(X) and called it using integralN_mc(@(X) Rg(X), ...)

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Argument Definitions en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by