Short way to write this function?
2 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
José Antonio Torres
il 19 Mar 2019
Commentato: José Antonio Torres
il 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 Commenti
Risposta accettata
Walter Roberson
il 19 Mar 2019
Modificato: Walter Roberson
il 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 Commenti
Walter Roberson
il 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.
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Matrix Indexing in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!