没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
二大地七参数,vb
大家可以互相交换学习
我的 QQ :35643195
下面是一个超人师兄编的:
Private Function 二大地七参数(a, b短, a1(), a2(), da, df) '椭球的长轴和曲率变化率已知,那么就剩余其它7个参数
e = Sqr(a * a - b短 * b短) / a
f = (a - b短) / a
n1 = UBound(a1()) / 5
n2 = UBound(a2()) / 3
Dim b1(), l1(), p()
ReDim b1(1 To 2 * n1, 1 To 7), l1(1 To 2 * n1), p(1 To 2 * n1)
For i = 1 To n1
l = hu(a1(5 * i - 3))
b = hu(a1(5 * i - 2))
n = bn(a, e, b)
m = bm(a, e, b)
b1(2 * i - 1, 1) = -Sin(l) / (n * Cos(b)) 'L的系数
b1(2 * i - 1, 2) = Cos(l) / (n * Cos(b))
b1(2 * i - 1, 3) = 0
b1(2 * i - 1, 4) = Tan(b) * Cos(l)
b1(2 * i - 1, 5) = Tan(b) * Sin(l)
b1(2 * i - 1, 6) = -1
b1(2 * i - 1, 7) = 0
l1(2 * i - 1) = hu(a1(5 * i - 3)) - hu(a1(5 * i - 1))
p(2 * i - 1) = 1
b1(2 * i, 1) = -Sin(b) * Cos(l) / m 'B的系数
b1(2 * i, 2) = -Sin(b) * Sin(l) / m
大家可以互相交换学习
我的 QQ :35643195
下面是一个超人师兄编的:
Private Function 二大地七参数(a, b短, a1(), a2(), da, df) '椭球的长轴和曲率变化率已知,那么就剩余其它7个参数
e = Sqr(a * a - b短 * b短) / a
f = (a - b短) / a
n1 = UBound(a1()) / 5
n2 = UBound(a2()) / 3
Dim b1(), l1(), p()
ReDim b1(1 To 2 * n1, 1 To 7), l1(1 To 2 * n1), p(1 To 2 * n1)
For i = 1 To n1
l = hu(a1(5 * i - 3))
b = hu(a1(5 * i - 2))
n = bn(a, e, b)
m = bm(a, e, b)
b1(2 * i - 1, 1) = -Sin(l) / (n * Cos(b)) 'L的系数
b1(2 * i - 1, 2) = Cos(l) / (n * Cos(b))
b1(2 * i - 1, 3) = 0
b1(2 * i - 1, 4) = Tan(b) * Cos(l)
b1(2 * i - 1, 5) = Tan(b) * Sin(l)
b1(2 * i - 1, 6) = -1
b1(2 * i - 1, 7) = 0
l1(2 * i - 1) = hu(a1(5 * i - 3)) - hu(a1(5 * i - 1))
p(2 * i - 1) = 1
b1(2 * i, 1) = -Sin(b) * Cos(l) / m 'B的系数
b1(2 * i, 2) = -Sin(b) * Sin(l) / m
b1(2 * i, 3) = Cos(b) / m
b1(2 * i, 4) = -Sin(l)
b1(2 * i, 5) = Cos(l)
b1(2 * i, 6) = 0
b1(2 * i, 7) = -n * e * e * Sin(b) * Cos(b) / m
l1(2 * i) = hu(a1(5 * i - 2)) - hu(a1(5 * i)) + (n * e * e * Sin(b) * Sin(b) * da / a + m * (2 - e * e * Sin(b) * Sin(b)) * Sin(b) * Cos(b) * df / (1 - f)) / m
p(2 * i) = 1
Next
Dim u(), nni(), xx()
ReDim u(1 To 7), nni(1 To 28), xx(1 To 7)
nni() = 法系数(b1(), p())
u() = 常数项(b1(), p(), l1())
nni() = ni(nni(), 2 * t)
xx() = bx(nni(), u())
Dim a3()
ReDim a3(1 To 3 * n2)
For i = 1 To n2
l = hu(a2(3 * i - 1))
b = hu(a2(3 * i))
n = bn(a, e, f)
m = bm(a, e, b)
a3(2 * i - 1) = (-Sin(l) * xx(1) / Cos(b) + Cos(l) * xx(2) / Cos(b)) / n + Tan(b) * Cos(l) * xx(4) + Tan(b) * Sin(l) * xx(5) - xx(6)
a3(2 * i) = (-Sin(b) * Cos(l) * xx(1) - Sin(b) * Sin(l) * xx(2) + Cos(b) * xx(3)) / m - Sin(l) * xx(4) + Cos(l) * xx(5) - n * e * e * Sin(b) * Cos(b) * xx(7) / m + (n * e * e * Sin(b) * Cos(b) * da / a + m * (2 - e * e * Sin(b) * Sin(b)) * Sin(b) * Cos(b) * df / (1 - f)) / (m + n)
a3(2 * i - 1) = l + a3(2 * i - 1)
a3(2 * i) = b + a3(2 * i)
a3(2 * i - 1) = du(a3(2 * i - 1))
a3(2 * i - 1) = Format(a3(2 * i - 1), "0.000000")
a3(2 * i) = du(a3(2 * i))
a3(2 * i) = Format(a3(2 * i), "0.000000")
Next
b1(2 * i, 4) = -Sin(l)
b1(2 * i, 5) = Cos(l)
b1(2 * i, 6) = 0
b1(2 * i, 7) = -n * e * e * Sin(b) * Cos(b) / m
l1(2 * i) = hu(a1(5 * i - 2)) - hu(a1(5 * i)) + (n * e * e * Sin(b) * Sin(b) * da / a + m * (2 - e * e * Sin(b) * Sin(b)) * Sin(b) * Cos(b) * df / (1 - f)) / m
p(2 * i) = 1
Next
Dim u(), nni(), xx()
ReDim u(1 To 7), nni(1 To 28), xx(1 To 7)
nni() = 法系数(b1(), p())
u() = 常数项(b1(), p(), l1())
nni() = ni(nni(), 2 * t)
xx() = bx(nni(), u())
Dim a3()
ReDim a3(1 To 3 * n2)
For i = 1 To n2
l = hu(a2(3 * i - 1))
b = hu(a2(3 * i))
n = bn(a, e, f)
m = bm(a, e, b)
a3(2 * i - 1) = (-Sin(l) * xx(1) / Cos(b) + Cos(l) * xx(2) / Cos(b)) / n + Tan(b) * Cos(l) * xx(4) + Tan(b) * Sin(l) * xx(5) - xx(6)
a3(2 * i) = (-Sin(b) * Cos(l) * xx(1) - Sin(b) * Sin(l) * xx(2) + Cos(b) * xx(3)) / m - Sin(l) * xx(4) + Cos(l) * xx(5) - n * e * e * Sin(b) * Cos(b) * xx(7) / m + (n * e * e * Sin(b) * Cos(b) * da / a + m * (2 - e * e * Sin(b) * Sin(b)) * Sin(b) * Cos(b) * df / (1 - f)) / (m + n)
a3(2 * i - 1) = l + a3(2 * i - 1)
a3(2 * i) = b + a3(2 * i)
a3(2 * i - 1) = du(a3(2 * i - 1))
a3(2 * i - 1) = Format(a3(2 * i - 1), "0.000000")
a3(2 * i) = du(a3(2 * i))
a3(2 * i) = Format(a3(2 * i), "0.000000")
Next
剩余7页未读,继续阅读
zhdel
- 粉丝: 1
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功
- 1
- 2
前往页