Login
Password
Sources on this Page

> Headlines by Category

 Home / Science / Mathematics

You are using the plain HTML view, switch to advanced view for a more complete experience.

Reading gromacs gro file in MATLAB

This function will read all the coordinates (x,y,z) from a *.gro file into an array in matlab.
The resulting N x 4 array will have the atomic number in first column followed by x,y,z in others. Make sure you have the gro file in the same directory.
You can add a line of code to include the type of atom but then take care of variable type conflict !

twelvepulseinductionfurnace.mdl

this is model of induction furnace with twelve pulse rectifire which has lower total harmonica distrotion effect in network.

Automated phenotyping of mouse behavior

This package contains tools for automating the scoring of several common behavior tasks used by the neuroscience community. Automated behavior tasks include: barnes maze, open field, zero maze, y-maze, fear conditioning, social interaction, spatial object recognition and morris water maze. Please visit our website: www.seas.upenn.edu/~molneuro/autotyping.html for detailed user guide, instructional videos, sample experiments and accompanying manuscript in Frontiers in Behavioral Neuroscience describing the use of this software.
To get started, add mmread and Inpaint_nans to your matlab path. Type "autotyping" at the command window and select the appropriate behavior experiment you'd like to analyze. Batch processing is recommended.

Step #11 in the Human Conquest of Space
In my previous post I talked about the magical quality of an orbit: Each time a spacecraft settles into a permanent path around a new object, humanity has taken one more step in venturing off this little blue world of ours and becoming colonizers of the universe. When the Rosetta reached Comet Churyumov-Gerasimenko on August 6 it marked the 11th celestial body we’ve orbited: Step #11 in the human conquest of space. Going from zero to 11 has taken just under 57 years; the anniversary of Sp
Re: Random Vectors with Fixed Sum / Naturals
>
> for i = 1:P
> s = sort(randperm(n+L1-1,n-1));
> v(i,:) = diff([0 s n+L1])-1;
> end
>
> which WORKS, but not sure if there's a way to do "P" all at once.

RANDPERM can't be vectorized in a straight manner. There are different ways, but the good choice depends on typical values of P, n, L1 you are dealing with.

Bruno
Re: Selective Randomized Matrix
On 10/01/2014 2:21 PM, Bruno Luong wrote:
> dpb <none@non.net> wrote in message <m0hhn4$mfn$1@dont-email.me>...
>
>>
>> x(randi(numel(x),n,1))=v;
>>
>
> OP should be warned the effect when randi returns duplicated values.
> RANDPERM is perhaps better option.

Good point, Bruno...won't happen often, but _is_ possible...

--
Re: fprintf standard error
On 10/01/2014 11:54 AM, Eric Simon wrote:
...

>> The facility is there to be able to separate the two but again, afaik,
>> within an interactive Matlab session there's no way to actually make
>> the redirection.
>
>
> Thanks!

Hmmm...on reflection--

I've not tried it; not sure what (if anything) might happen if you
started a Matlab session with the redirection

   matlab >&> filename

OK, tried w/ a new session--as expected, didn't redirect standard error
from command window.

That might be worth a query to TMW official support on whether it's
doable, and if so, how, 'cuz afaict it's not documented if so. Only
thing I can find is that the default device/handle is open.

--
ybus.m

This is to represent the Ybus power system

symbolic computation of large equation
Hi

Let's say I have a equation like this:
(((L/r)-r(4*L))*(Qr)+L-((r+m)^2)/4*L)/(((Qr)+r+m))=0

for
Qr we have to put
Qr=((l*r)/(2*L))-(2*l*L/r)-a+sqrt(((16*L^2*l^2/r^2)-8*l^2+l^2*r^2/L^2 -16*l*L+16*l*yo+16*a*l*L/r -4*a*r*l/L)/2)


so we have ...
(((L/r)-r(4*L))*(((l*r)/(2*L))-(2*l*L/r)-a+sqrt(((16*L^2*l^2/r^2)-8*l^2+l^2*r^2/L^2 -16*l*L+16*l*yo+16*a*l*L/r -4*a*r*l/L)/2))+L-((r+m)^2)/4*L)/(((((l*r)/(2*L))-(2*l*L/r)-a+sqrt(((16*L^2*l^2/r^2)-8*l^2+l^2*r^2/L^2 -16*l*L+16*l*yo+16*a*l*L/r -4*a*r*l/L)/2))+r+m))=0

how to using symbolic computation determine variable "a" from equation above???

My all problem look like that:

K1=(L/r)-(r/(4*L));

Kc=-(Xc+a)/(2*l);

K3=(((L/r)-r(4*L))*(Qr)+L-((r+m)^2)/4*L)/(((Qr)+r+m));

Qr=((l*r)/(2*L))-(2*l*L/r)-a+sqrt(((16*L^2*l^2/r^2)-8*l^2+l^2*r^2/L^2 -16*l*L+16*l*yo+16*a*l*L/r -4*a*r*l/L)/2);

abs(((k1)-(kc))/(1+(k1)*(kc)))==abs(((k3)-(kc))/(1+(k3)*(kc)))


using solve (equation, a)
Matlab give error
Re: combination of arrays (matrix of strings)
Bruno Luong" <b.luong@fogale.findmycountry> wrote in message <m0gni7$nbt$1@newscl01ah.mathworks.com>...
> If you insist on using KRON (why?) :
> [u, ~, J] = unique(label);
> u(kron(reshape(J,size(label)),ones(1,3)))
>
> You have many other methods that works on indexing:
>
> http://www.mathworks.com/matlabcentral/newsreader/view_thread/260143
>
> % Bruno

Thank you Bruno!
Well I don't insist but that was just a way I already knew that. By the way thank you for the code and introducing the other methods.
I'd like to draw your attention to the output matrix which should have been as:

'N' 'N' 'N' 'P' 'P' 'P' 'F' 'F' 'F'
'N' 'P' 'F' 'N' 'P' 'F' 'N' 'P' 'F'

but in your code is:

'N' 'N' 'N' 'P' 'P' 'P' 'F' 'F' 'F'
'N' 'N' 'N' 'P' 'P' 'P' 'F' 'F' 'F'

Am I doin' anything wrong?
Post Selected Items to:

Showing 10 items of about 1000

home  •   advertising  •   terms of service  •   privacy  •   about us  •   contact us  •   press release design by Popshop •   Official PR partner B2BLogger.com •   © 1999-2014 NewsKnowledge