求积分公式或VBA代码

YANG6815475 · 发表于 2011-9-15 21:21

求积分公式或VBA代码！

fjmxwrs · 发表于 2011-9-15 22:24

回复 YANG6815475 的帖子

什么问题？

YANG6815475 · 发表于 2011-9-16 09:24

回复 fjmxwrs 的帖子

计算公路匝道点位坐标的复化辛普森公式是一个积分公式计算的，在Excel函数中无此函数，想编制此函数或相应的VBA太难了！！！这个问题不解决，无法编计算程序。

wayy · 发表于 2011-9-16 09:25

没有附件也没有问题说明，如何给你解答？

fjmxwrs · 发表于 2011-9-16 12:38

回复 YANG6815475 的帖子

你得讲清楚你的计算是怎么个过程，或上附件写清结果样子。

YANG6815475 · 发表于 2011-9-17 22:49

ZD-2程序(匝道中边桩坐标计算程序) 积分法.rar (30.86 KB, 下载次数: 8)

2011-9-17 22:48 上传

点击文件名下载附件
匝道中边桩坐标计算程序 (二十二) ZD-2程序(匝道中边桩坐标计算程序) 积分法
⑴有关资料
采用复化辛普森公式：积分法。
⒈ 匝道是组成高等级公路立交的基本单元，其形式千变万化，就线形而言，也是由直线段、回旋曲线段、圆曲线段组成。但是，组成立交的匝道涉及线形的曲率变化特点，利用复化辛普森公式导证了计算公路匝道点位坐标的通用公式。
　⒉、公路匝道坐标计算公式推导
　　任一点切线方位角
αi＝αA＋(Рi＋РA）／2 ×（DK i －DKA）（弧度）    （1）
式中：Рi─曲线上任一点的曲率，
рi＝рA＋（РB－РA）／（DKB－DKA）×（DK i －KDA）（2）

РA、РB─曲线上始、终点的曲率，等于1 / R，R为半径（米），
DKA 、DKB─曲线上始、终点的里程〈m〉，
DK i ─曲线上任一点的里程〈m〉，
αA─曲线上始点的方位角（弧度），
αi─曲线上任一点的方位角，
设P=（РB－РA）／（DKB－DKA），Li=DK i －KDA，
（1）式化为
αi＝αA＋(РA+P* Li＋РA）／2× Li
αi＝αA＋РA* Li +P* Li2  / 2 （弧度）          ⑶
为方便计算，令变量X= Li，则⑶式化为
αi＝αA＋РA* X +P* X 2  / 2                      ⑷
  ⑷根据复化辛普森公式，任一点的坐标为
                     ⑸
　　                                           （6）
  由⑸、⑷式得
X=XA+∫DKADki（Cos（P/2*X2+PA*X+аA）dl
y=XA+∫DKADki（Sin（P/2*X2+PA*X+аA）dl       ⑹

在CASIO4800P计算器中积分，占用K、L、M、N，分别存贮a、b、n、积分结果，因此积分时不能使用这些变量。
⑵ ZD-2程序设计
AC  MODE  5  1  Z  D  ALPHA — 2  EXE  1
A"KA"∶G"XA"∶C"YA"∶D"QA"∶                         始点桩号，坐标，切线方位角
R"PA"∶E "KB "∶F"PB"∶B∶ O∶                         始、终点曲率，                                                                                                    Lb1 1∶｛ S ｝∶ S"KP"∶ S＜A  => S＝A∶≠=>H＝S-A∶ P＝（F-R）÷（E-A）∶
Q=π÷180∶J =D+（RH+PH^2÷2）÷Q ∶
J“QB”=360 Frac((J+360) ÷ 360 )  ▲  Rad∶
I"X" = G+∫（cos（PX^2÷2+RX+D Q , 0 , H , 4 ▲
U"Y" = C+∫（sin（PX^2÷2+RX+D Q, 0 , H , 4 ▲ Deg∶
I"X1" =I-B cos（J+O）  ▲ U"Y1" =U-B sin（J+O） ▲
I"X2" =I+2B cos（J+O）▲ U"Y2" =U+2B sin（J+O）▲ Goto  1 EXE
  字母代表意义
S─所求点桩号(m)，
A"KA"、G"XA"、C"YA"、─ 始点桩号、坐标,
D"QA"─始点切线方位角，
R"PA"、F"PB"─始、终点曲率，左转输入负值，右转输入正值，为曲线半径的倒数(1/R)，直线输入零（或1/E50），
B─公路半宽度(m)，
O─横断面与中桩点切线的左夹角(度)，
J─所求点切线方位角(度)，
X、Y ─中桩坐标(m)，
X1、Y1 ─左边桩坐标(m)，
X2、Y2 ─右边桩坐(m)，
X=XA+∫DKADki（Cos（P/2*X2+PA*X+аA）dl
y=XA+∫DKADki（Sin（P/2*X2+PA*X+аA）dl       ⑹

YANG6815475 · 发表于 2011-9-17 23:03

以下列积分公式做成函数公式或VBA：
X=XA+∫DKA Dki（Cos（P/2*X2+PA*X+аA）dl
y=XA+∫DKA Dki（Sin（P/2*X2+PA*X+аA）dl       ⑹
化为数学公式：
y1=XA+∫DKA Dki（Cos（P/2*X2+PA*X+аA）dl
y2=XA+∫DKA Dki（Sin（P/2*X2+PA*X+аA）dl       ⑹
式中：XA、P、PA、аA为常数，
DKA、 Dki为定积分上、下区间，等于a、b值。
X2  是x的平方。详见文件。
这是用卡西欧4800P计算器的积分功能编制的程序。

YANG6815475 · 发表于 2011-9-18 09:50

实际上就是请高手编制以下公式的函数：
∫DKA Dki（Cos（P/2*X2+PA*X+аA）dl
∫DKA Dki（Sin（P/2*X2+PA*X+аA）dl ⑹
式中：XA、P、PA、аA为常数，
DKA、 Dki为定积分上、下区间，等于a、b值。
X2 是x的平方。详见文件。

高 · 发表于 2011-9-18 11:36

让大家帮你做公式并不难，但是你自己的公式太复杂了，我都不知道能不能看懂。
最好的方法有两个：
1、把你的问题分解成多个问题，你自己根据大家的答案再把若干公式合在一起，成为你想要的公式
2、自己多学习、多积累，自己以后做公式
3、很有耐心或对你的问题比较专业的人帮你

YANG6815475 · 发表于 2011-9-18 20:38

这是一个比较难的问题，我查阅了许多资料，在Excel函数中没有积分函数公式。这个积分公式见下：
实际上就是请高手编制以下公式的函数：
∫DKA Dki（Cos（P/2*X2+PA*X+аA）dl
∫DKA Dki（Sin（P/2*X2+PA*X+аA）dl ⑹

式中：XA、P、PA、аA为常数，
DKA、 Dki为定积分上、下区间，等于a、b值。
X2 是x的平方。详见文件。
积分式子是一个一元2次式，对我这样的水平无法解决的，对Excel函数很熟的可能解决。电脑软件不知有没有计算定积分的。不知怎的，计算器都有解积分的功能，而电脑却没有？！

		自动登录	找回密码
密码			注册

求积分公式或VBA代码

浏览过的版块