Save outputs of nested for loops in MATLAB -

-

i have following codes wish have output matrix rpp of (10201,3). run code (which takes bit long) check matrix size of rpp , see (1,3), tried many things couldn't find proper way. logic of codes take 6 values (contain 4 constant values , 2 variable values (chosen 101 values)) , make calculation 3 different i1 , store every output vector of 3 in matrix (101*101 (pairs of 2 variable values)) rows , 3 (for each i1) columns.

i appreciate help

    vp1=linspace(3000,3500,101);     vp2=3850;     rho1=2390;     rho2=2510;     vs1=linspace(1250,1750,101);     vs2=2000;     i1=[10 25 40];       rpp = zeros(length(vp1)*length(vs1),length (i1));      n=1:length(vp1)*length(vs1)     m=1:length (i1)     l=1:length(vp1)         k=1:length(vs1)              p=sin(i1)/vp1(l);            i2=asin(p*vp2);            j1=asin(p*vs1(k));            j2=asin(p*vs2);            a=rho2*(1-2*vs2^2*p.^2)-rho1*(1-2*vs1(k).^2*p.^2);            b=rho2*(1-2*vs2^2*p.^2)+2*rho1*vs1(k)^2*p.^2;            c=rho1*(1-2*vs1(k)^2*p.^2)+2*rho2*vs2^2*p.^2;            d=2*(rho2*vs2^2-rho1*vs1(k)^2);            e=b.*cos(i1)./vp1(l)+c.*cos(i2)/vp2;            f=b.*cos(j1)./vs1(k)+c.*cos(j2)/vs2;            g=a-d*(cos(i1)/vp1(l)).*(cos(j2)/vs2);            h=a-d*(cos(i2)/vp2).*(cos(j1)/vs1(k));            d=e.*f+g.*h.*p.^2;            rpp=((b.*(cos(i1)/vp1(l))-c.*cos((i2)/vp2)).*f-(a+d*((cos(i1)/vp1(l))).*(cos(j2)/vs2)).*h.*p.^2)./d         end     end     end     end

the way use loop wrong. you're running calculation length(vp1)*length(vs1) * length (i1) * length(vp1) * length(vs1) times. here's correct way. changed l lll won't confuse number 1. in each iteration of first loop, you're running length(vs1) times, , need assign result (a 1x3 array) rpp using row number specified k+(lll-1)*length(vp1).

for lll=1:length(vp1)     k=1:length(vs1)        p=sin(i1)/vp1(lll);        i2=asin(p*vp2);        j1=asin(p*vs1(k));        j2=asin(p*vs2);        a=rho2*(1-2*vs2^2*p.^2)-rho1*(1-2*vs1(k).^2*p.^2);        b=rho2*(1-2*vs2^2*p.^2)+2*rho1*vs1(k)^2*p.^2;        c=rho1*(1-2*vs1(k)^2*p.^2)+2*rho2*vs2^2*p.^2;        d=2*(rho2*vs2^2-rho1*vs1(k)^2);        e=b.*cos(i1)./vp1(lll)+c.*cos(i2)/vp2;        f=b.*cos(j1)./vs1(k)+c.*cos(j2)/vs2;        g=a-d*(cos(i1)/vp1(lll)).*(cos(j2)/vs2);        h=a-d*(cos(i2)/vp2).*(cos(j1)/vs1(k));        d=e.*f+g.*h.*p.^2;        rpp(k+(lll-1)*length(vp1),:)=((b.*(cos(i1)/vp1(lll))-c.*cos((i2)/vp2)).*f-(a+d*((cos(i1)/vp1(lll))).*(cos(j2)/vs2)).*h.*p.^2)./d;     end end

Comments

Post a Comment

Popular posts from this blog

dns - How To Use Custom Nameserver On Free Cloudflare? -

-
i have problem. i own small hosting company. trying add cloudflare site free cdn. problem when sign have change nameservers nameservers. have nameserver dns set through godaddy, how add them cloudflare? main site uses company's nameservers whenever add them site goes offline. when go cloudflare's dns records click on dropdown, click ns, , there's 2 boxes 'name' , 'nameserver'. mean input 'ns1' 'name' , 'ns1.domain.com' 'nameserver'? input server's ip addresses? is there way of accomplishing free version of cloudflare? that impossible. option have custom name servers available on business , enterprise. see https://www.cloudflare.com/plans . also, might mistaken. must change name servers in domain register settings, not on cloudflare.
Read more

python - Pygame screen.blit not working -

-
in program creating in python/pygame, screen.blit not working. in other programs have made, screen.blit works. error: traceback (most recent call last): file "c:...\main.py", line 46, in <module> screen.blit(sight.image, sight.rect) attributeerror: 'int' object has no attribute 'blit' the error in question caused screen variable being allocated integer instead of pygame.surface() object.perhaps accidentally rewrote variable somewhere? from can make of doing, seems stating position of blit full rect instead of coordinate. when do: screen.blit(sight.image,sight.rect) what program seeing is: screen.blit(sight.image,(sight.rect.x,sight.rect.y,sight.rect.width,sight.rect.height)) i assuming don't intend this. remember rect object consists of position, width , height. to correct code, put (sight.rect.x,sight.rect.y) location goes instead this: screen.blit(sight.image,(sight.rect.x,sight.rect.y)) the .x , .y corres
Read more

c# - Web API response xml language -

-
i having simple post url clients send post request. http://mydomain.com/requests/request post data: <request> <orderid>1e008921</orderid> <ordername>mc1</ordername> </request> to respond sending them xml <srequest> <requestedby>client 1 </requestedby> <requestname>test name</requestname> <requeststatus>success</requeststatus> <srequest> now client wants above xml in different language (say french ) how handle ? how send response language preference in web api ? specially post scenario ? (append language string ? or best practices ?) i suggest either include in url: http://mydomain.com/fr/requests/request or use http header accept-language . read more in answer to: good way changing language resources dynamically request
Read more