Notifications
Clear all

I give up. Stoopid maffs.

21 Posts
9 Users
0 Reactions
5 Views
(@micheal-daubyn-2-2-2)
Posts: 154
Estimable Member Registered
Topic starter
 

Hey guys.

Well, I'm stumped.

A friend of mine is currently studying at Uni and has come a question that needs answering. He turned to me for help. It's a graded question, so I am not playing with his actual numbers and won't pass on any numbers here.

My problem is this- I've never seen this problem before. I kept all my notes and study material but none of it covers this question.

It's a braced quad question, with two fixed points and two free points. It requires the variation of coordinates method to solve it. My problem is, the question itself gives angles + distances, whereas I learned how to solve it with directions + distances. Each angle has a different angular SD. The distance SDs are all equal.

His teacher isn't handing out the equation. What a world we live in.

Does anyone know what the equation is? I can't for the life of me figure out how to form the observation equation for the matrix. I got it to work with my own numbers in Starnet, then converted the angles to directions by de-propagating the angular SDs and got the same result in Starnet. But I can't reproduce it long-hand.

Would anyone happen to have the equation laying around? Has anyone seen this before?

Thanks,

Mick.

 
Posted : 13/05/2023 4:04 am
(@field-dog)
Posts: 1372
Noble Member Registered
 

Please post the question.

 
Posted : 13/05/2023 5:00 am
(@holy-cow)
Posts: 25292
 

Not grasping the true problem as this seems to be a fairly simple trigger-nom-uh-tree problem.

Using an American football expression, tell him to "PUNT".

 
Posted : 13/05/2023 8:37 am
(@bill93)
Posts: 9834
 

What are we to compute? Maybe best estimates of the positions of two points and their standard errors? How many measurements are given?

"A problem well stated is a problem half solved."~Charles Kettering

 
Posted : 13/05/2023 10:15 am
(@micheal-daubyn-2-2-2)
Posts: 154
Estimable Member Registered
Topic starter
 

Thanks for the replies guys. I'll paraphrase the question and get back to you. I'm trying to not cross any ethical lines here.

 
Posted : 13/05/2023 11:52 am
(@bill93)
Posts: 9834
 

The problem has to be linearized by using derivatives of angles and distances with respect to coordinates. Then you form the matrix equations.  Here are some notes that could help.

https://www.dropbox.com/s/esgmuxk2s0mqio4/LSQ_notes.doc?dl=0

 
Posted : 13/05/2023 12:00 pm
(@micheal-daubyn-2-2-2)
Posts: 154
Estimable Member Registered
Topic starter
 

The question itself:

"The following network was observed with both angles and distances. There are two fixed stations P and S.

Quad

Station

Easting

Northing

Notes

P

10000.000

10000.000

Fixed

S

11000.003

9499.995

Fixed

Q

10190.000

11000.000

Approx

R

11500.000

11410.000

Approx

 

Angle (At-From-To)

Obs Angle (DDD-MM-SS)

Std Deviation (")

P-Q-R

35-39-50

4.5

P-R-S

69-35-21

6.5

S-P-Q

35-21-50

4.0

S-Q-R

42-48-52

3.5

Q-S-P

39-23-00

4.5

Q-R-P

118-24-21

7.0

R-S-P

32-13-41

3.5

R-P-Q

25-55-29

2.0

 

The following distances have been observed with a standard deviation of +/-0.015m +/- 4ppm

 

Distance (From-To)

Observed Distance (m)

P-Q

1019.91

P-R

2051.908

R-S

1964.765

S-Q

1700.081

R-Q

1360.161

 

Use the least squares method to determine the adjusted coordinates of Q and R. Note: you will need to iterate the solution using a convergence limit of 0.001."

So, that's the actual question itself, with the actual data he was given. As I said, I know how to do this when directions are given, but not angles.

To be clear, I'm trying to tread carefully here on someone else's behalf. So I'd rather not cross any lines. I'm only looking for the variation of coordinates formula that deals with angles and two unknown points.

Thank guys.

Mick.

 
Posted : 13/05/2023 12:03 pm
(@micheal-daubyn-2-2-2)
Posts: 154
Estimable Member Registered
Topic starter
 

@bill93 Thanks Bill. That doesn't address the problem at hand, but it's an interesting read. Anything that suggests we fudge figures to make it work is an interesting read!

Did that document come from Microsurvey?

 
Posted : 13/05/2023 12:15 pm
(@jon-payne)
Posts: 1595
Noble Member Registered
 

It's a braced quad question, with two fixed points and two free points. It requires the variation of coordinates method to solve it. My problem is, the question itself gives angles + distances, whereas I learned how to solve it with directions + distances.

My guess would be the professor is making the students work through several different ideas in one problem.  If there are two fixed points with the coordinates listed, as you show in your follow-up post, then there is also a fixed azimuth between them.  Would that not then mean there are azimuths from the fixed line PS (by adding observed angles) such that you would have directions (not just raw angles).

 

https://www.researchgate.net/publication/254250473_A_review_of_least_squares_theory_applied_to_traverse_adjustment

 
Posted : 13/05/2023 12:36 pm
(@bill93)
Posts: 9834
 

that document

I put that document together by reading stuff from various sources. I pretty much copied the procedures for getting derivatives for distances, changing to my own notation, and then developed the others.

That was my working paper as I wrote my own general 2-D LS program with some nice features like angles between lines that you don't know the vertex point (including parallels) and specified differences between lengths that are yet undetermined. Once you see how the basics work, such add-ons are easy. What wasn't easy was an attempt at graphics, and I've got a bug (probably memory leak) that crashes it sometimes.

 

 
Posted : 13/05/2023 12:44 pm
(@micheal-daubyn-2-2-2)
Posts: 154
Estimable Member Registered
Topic starter
 

@jon-payne I went down that road first. Finding a way to convert the angles to directions was my first choice, because I know the directions formula.

This would require starting with the azimuth on the P-S and swinging all the angular readings onto it. But the first swinging a second angle onto a first angle means swinging it onto an angle that has a SD that hasn't been solved yet. For example, if you take the azimuth P-S, then calculate the direction P-R by adding the negative angle to it, P-R ends up on the P-S azimuth. Of course. But the direction P-R still has the 6.5" SD attached to it that hasn't been dealt with yet.

If you then decide to swing line P-Q onto line P-R, you are adding a second angular SD (4.5") to a first angular SD (6.5") that hasn't been solved yet. That's fine too, because you can propagate those SDs.

That's all fine, and can (through the process of creating one giant mess) be figured out eventually. Rinse, lather, repeat. Then you end up using this equation:

Direction

But the professor doesn't want it done with that equation. He wants this:

Angles

Now, with the direction equation, or an angular equation that I convert to a direction equation, you work with a start point and an end point of each line. Then, depending on whether each start or end is fixed or free, you slot those numbers into your A matrix accordingly.

But when working with angles directly, you don't have a start point and end point. You have a left line and a right line. That's the part that is causing the problem, because the formula above is for resections, where your start point in unknown and you are relating your reading to two know points.

Similarly, that equation could be flipped to work for the sections of the quad where you dealing with intersections instead of resections. But what do you when, for example, you are at Q and reading to P, S and R? You are reading to known, known and unknown.

That's where I am stuck.

 
Posted : 13/05/2023 1:16 pm
(@jon-payne)
Posts: 1595
Noble Member Registered
 

But the direction P-R still has the 6.5" SD attached to it that hasn't been dealt with yet.

It has been too long since I have done one of these, but shouldn't the SDs be worked in as part of a weight matrix.

The following distances have been observed with a standard deviation of +/-0.015m +/- 4ppm

I think this would also similarly need to be taken into account within a weight matrix.

 
Posted : 13/05/2023 1:42 pm
(@micheal-daubyn-2-2-2)
Posts: 154
Estimable Member Registered
Topic starter
 

"but shouldn't the SDs be worked in as part of a weight matrix."

Yes they should, and do. But does the matrix know to propagate as the angles accumulate?

Regardless, it doesn't solve the issue of the requirement to work with angles rather than directions.

 
Posted : 13/05/2023 1:53 pm
(@dave-karoly)
Posts: 12001
 

I’m just a simple country land surveyor, not a mathematician, Jim!

 
Posted : 13/05/2023 2:25 pm
(@rover83)
Posts: 2346
Noble Member Registered
 

But the first swinging a second angle onto a first angle means swinging it onto an angle that has a SD that hasn't been solved yet. For example, if you take the azimuth P-S, then calculate the direction P-R by adding the negative angle to it, P-R ends up on the P-S azimuth. Of course. But the direction P-R still has the 6.5" SD attached to it that hasn't been dealt with yet.

 

Perhaps I am missing something, but if angles are converted to directions, the SD for each direction is independent of which direction everything is referenced to.

Angular SD is just the sum of two direction SDs, the square root of the sum of the squares of each direction SD.

If the foresight and backsight were observed the same manner (2DR, 4DR, etc for both FS and BS) then the a priori SD for each direction is the SD of the angle divided by the square root of 2.

 

So for the first angle P-Q-R (SD of 4.5") the equivalent direction SD for each direction (both P-Q and P-R) would be ~3.2", because adding two directions that each have that SD results in an angle between those two directions of 4.5" SD. Next angle P-R-S with SD of 6.5" translates to 4.6" SD directions, and so on.

Orienting each direction to a common starting azimuth doesn't change its individual a priori SD.

 

I really wish I hadn't left Ghilani's Adjustment Computations at work, I am sure the answer is in there somewhere. Maybe I need a second copy for my home office...

 
Posted : 14/05/2023 5:14 am
Page 1 / 2
Share: