Register to reply

Need help with Fortran

Share this thread:
sue132
#1
Mar14-12, 06:02 AM
P: 14
Hi,

I'm trying to generate a set of random vectors with binary values. I need to generate them such that each time, one new vector gets added to the existing matrix. I have been using call random_number(ranval) for the generation, but I haven't been able to get linearly independent vectors.

(I need to orthogonalize them, so the vectors have to be linearly independent)

Any help would be much appreciated. Any info on random number generators would also be useful.

Thanks,
S.

(P.S. I am currently working with Fortran77)
Phys.Org News Partner Science news on Phys.org
Hoverbike drone project for air transport takes off
Earlier Stone Age artifacts found in Northern Cape of South Africa
Study reveals new characteristics of complex oxide surfaces
Mark44
#2
Mar14-12, 08:27 PM
Mentor
P: 21,216
Quote Quote by sue132 View Post
Hi,

I'm trying to generate a set of random vectors with binary values. I need to generate them such that each time, one new vector gets added to the existing matrix. I have been using call random_number(ranval) for the generation, but I haven't been able to get linearly independent vectors.
You don't provide enough context for me to get very specific, but I'll take a stab at it.

How big are the vectors? And when you say "binary values" I'm assuming you mean the vectors consist of 0s and 1s - just checking.

It might be that you can't do what you're trying to do. Let's say for the sake of simplicity that your vectors have 4 elements, such as <0, 1, 1, 1>. If you generate a vector with four values, that vector and the one you just generated might not form a linearly independent set of vectors. In any case, you'll never get more than 4 vectors for a linearly independent set (with n = 4).

The same limitation exists for higher dimension vectors. If n = 25, you won't be able to get more than 25 vectors in your matrix.

Maybe you know all this, but since I don't know your mathematical background, I felt I should mention it.

Quote Quote by sue132 View Post

(I need to orthogonalize them, so the vectors have to be linearly independent)

Any help would be much appreciated. Any info on random number generators would also be useful.

Thanks,
S.

(P.S. I am currently working with Fortran77)
AlephZero
#3
Mar14-12, 09:40 PM
Engineering
Sci Advisor
HW Helper
Thanks
P: 6,929
Your question seems to be about generating Hadamard matrices, or a generalization of them. http://en.wikipedia.org/wiki/Hadamard_matrix

sue132
#4
Mar14-12, 11:14 PM
P: 14
Need help with Fortran

Quote Quote by AlephZero View Post
Your question seems to be about generating Hadamard matrices, or a generalization of them. http://en.wikipedia.org/wiki/Hadamard_matrix
I think that is what my end result would be. I have to generate a matrix with random vectors, and then when i do the orthogonalization, I think I'll get a Hadamard matrix as the result.
Thanks for the info
sue132
#5
Mar14-12, 11:23 PM
P: 14
Quote Quote by Mark44 View Post
You don't provide enough context for me to get very specific, but I'll take a stab at it.

How big are the vectors? And when you say "binary values" I'm assuming you mean the vectors consist of 0s and 1s - just checking.

It might be that you can't do what you're trying to do. Let's say for the sake of simplicity that your vectors have 4 elements, such as <0, 1, 1, 1>. If you generate a vector with four values, that vector and the one you just generated might not form a linearly independent set of vectors. In any case, you'll never get more than 4 vectors for a linearly independent set (with n = 4).

The same limitation exists for higher dimension vectors. If n = 25, you won't be able to get more than 25 vectors in your matrix.

Maybe you know all this, but since I don't know your mathematical background, I felt I should mention it.
So, if I have 10 as my vector length, I should be able to get at least 9 linearly independent vectors, right? I'm currently using +1's and -1's as my binary values.

What I'm doing at present is adding one vector at a time to my existing matrix using a random number generator (my random numbers are between 0 and 1, so I made anything lesser than 0.5 a -1, and the rest are +1's.) Should I change my random number generator, or would it be a better idea to ignore a dependent vector and generate a new vector in its place?

Thanks for the reply
Mark44
#6
Mar15-12, 01:23 AM
Mentor
P: 21,216
Why do you have -1 entries in your vectors? There's nothing wrong with it, but when most people think binary, they think 0 and 1. Having a 0 instead of a -1 might make the calculations easier.

If your vectors are length 10, you should be able to come up with 10 lin. independent vectors.

What I would do is generate a vector, and then figure out if that vector and the ones I already had form a linearly independent set.

I don't know any good algorithms for that, but at each step (with N vectors, where 2 <= N <= 10), you need to solve N equations for N unknowns. If the only solution is c1 = c2 = ... = cN = 0, then the vectors are linearly independent.

If the vector I just came up with turned out to be a linear combination of the others (hence the set of vectors is now linearly dependent) I would discard it, and try again.
sue132
#7
Mar15-12, 02:21 AM
P: 14
Quote Quote by Mark44 View Post
Why do you have -1 entries in your vectors? There's nothing wrong with it, but when most people think binary, they think 0 and 1. Having a 0 instead of a -1 might make the calculations easier.
I'm using -1's as they are more convenient than zeros for my further calculations.

What I would do is generate a vector, and then figure out if that vector and the ones I already had form a linearly independent set.
I have been trying to do that. As I generate each vector, I'm checking its projection with the previous vectors. If it's in the same direction as any of those, then I get a projection value of 0, and hence it will be linearly dependent. Hope I'm right in doing it this way.

If the vector I just came up with turned out to be a linear combination of the others (hence the set of vectors is now linearly dependent) I would discard it, and try again.
This is what I'm trying to do now.

Thanks again for your patience.

(P.S. I'm still fairly new to Fortran, so am having some trouble finding my way around)
sue132
#8
Mar16-12, 05:17 AM
P: 14
I figured it out. Thanks


Register to reply

Related Discussions
Fortran 77 help making an empty array (or blank list if they exist in fortran) Programming & Computer Science 5
FORTRAN 95: New to fortran, want to learn how to input a function Programming & Computer Science 1
MPI in Fortran Programming & Computer Science 2
Accessing Fortran Modules within a Fortran library from Fortran Programming & Computer Science 0
Fortran v.s. Visual Fortran Programming & Computer Science 0