Short way to write this function?
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
José Antonio Torres
le 19 Mar 2019
Commenté : José Antonio Torres
le 19 Mar 2019
Hello,
Does anyone know of a solution to this problem?
function[y]=rayleigh_product(x,x2,x3,x4) % For up to 55 storm steps (threshold = 4 in the original case)
global SDSRp_ Tz_Tz_
y=(1-exp(-0.5*(x/SDSRp_(x2,x3,x4,1))^2))^(10800/Tz_Tz_(x4,2))...
*(1-exp(-0.5*(x/SDSRp_(x2,x3,x4,2))^2))^(10800/Tz_Tz_(x4,3))...
*(1-exp(-0.5*(x/SDSRp_(x2,x3,x4,3))^2))^(10800/Tz_Tz_(x4,4))...
*(1-exp(-0.5*(x/SDSRp_(x2,x3,x4,4))^2))^(10800/Tz_Tz_(x4,5))...
*(1-exp(-0.5*(x/SDSRp_(x2,x3,x4,5))^2))^(10800/Tz_Tz_(x4,6))...
*(1-exp(-0.5*(x/SDSRp_(x2,x3,x4,6))^2))^(10800/Tz_Tz_(x4,7))...
And it goes on till 55 parts. The code works now, but it would be useful to have the option of just doing the product of any number of functions. I.e, if I want to do the same analysys with a different data set.
0 commentaires
Réponse acceptée
Walter Roberson
le 19 Mar 2019
Modifié(e) : Walter Roberson
le 19 Mar 2019
y = prod(arrayfun(@(IDX) 1-exp(-0.5*(x./SDSRp_(x2,x3,x4,IDX).^(10800./Tz_Tz_(x4,IDX+1)))), 1:55));
This assumes that SDSRp_ and Tz_Tz_ are functions. If they were arrays then there is more potential for vectorization.
Is it practical to rewrite the two functions to accept a vector for the last parameter?
3 commentaires
Walter Roberson
le 19 Mar 2019
prod((1-exp(-0.5 *(x./SDSRp_(x2, x3, x4, 1:55)))).^(10800./Tz_Tz_(x4,1+(1:55))))
I think I had some brackets in the wrong place before.
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Matrix Indexing dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!