万能的ncf吧求能解超越方程的大神！！

yukai881208 · 发表于 2012-9-4 10:26:51

设计时不小心遇到超越方程，，说简单点就是知道弧长，弦长，求弧高{:09:} 于是解到最后得到这货
sin(8HL/(4H^2+x^2))=4Hx/(4H^2+x^2) h 和l 为已知量求x ，，，，据说是能通过刘顿二分法通过程序
求得近似值，，期待能用c++或vb解决的大神的粗线！！{:19:}

wangjunxiong · 发表于 2012-9-4 11:58:00

本帖最后由 wangjunxiong 于 2012-9-4 15:46 编辑

验算你的公式，公式是有问题的！例如，我随意画个圆，截取一小部分，测量得：H=27.95 ，L=25.32 ,那么 X =5.07 。代入公式的左边 sin(8HL/(4H^2+x^2))=0.031359 ，右边 4Hx/(4H^2+x^2) = 0.179915，这明显不相等。

artrustee · 发表于 2012-9-4 15:16:43

化成f(x)=(d+x^2)/(c*x)*sin(a/b+x^2)-1
其中，a=8HL,b=4H^2,c=4H,d=4H^2
x的定义域自己算一下，然后就是求f（x）=0的可行解
加油吧~

wangjunxiong · 发表于 2012-9-4 15:44:24

用别人的公式解决了。

登录/注册后可看大图

wangjunxiong · 发表于 2012-9-4 16:27:22

用rhinoscript 也做了一个，运行速度还是可以的！修改一些参数就可以为自己所用了。

登录/注册后可看大图

yukai881208 · 发表于 2012-9-4 16:32:57

wangjunxiong 发表于 2012-9-4 15:44

登录/注册后可看大图

用别人的公式解决了。

咦我当时验算60度角时貌似是正确的，，难道验算错了囧，，，，可否看下前辈的原文件这2张图看的还是不甚理解，，，另外好像问过2次问题都是前辈解答的- -真是不好意思了，，，

yukai881208 · 发表于 2012-9-4 16:34:12

artrustee 发表于 2012-9-4 15:16

登录/注册后可看大图

化成f(x)=(d+x^2)/(c*x)*sin(a/b+x^2)-1
其中，a=8HL,b=4H^2,c=4H,d=4H^2
x的定义域自己算一下，然后就是 ...

嗯感谢大神的回答~：）~

wangjunxiong · 发表于 2012-9-4 16:53:40

本帖最后由 wangjunxiong 于 2012-9-4 17:00 编辑

yukai881208 发表于 2012-9-4 16:32

登录/注册后可看大图

咦我当时验算60度角时貌似是正确的，，难道验算错了囧，，，，可否看下前辈的原文件这2张图看的还是不甚 ...

点击图片就可以放大看了！蛮清晰的！
原谅我！我没有保存。

那个rhinoscript的代码很低级，效率非常低。记得读书的时候有个叫二分法，这个办法应该会提高效率。我是懒人，就不深入搞了。

yukai881208 · 发表于 2012-9-4 16:59:57

wangjunxiong 发表于 2012-9-4 16:53

登录/注册后可看大图

点击图片就可以放大看了！蛮清晰的！
原谅我！我没有保存。

嗯我也确实要自己多理解一下！不能当纯伸手党！还是感谢前辈了！！！

artrustee · 发表于 2012-9-4 17:09:01

yukai881208 发表于 2012-9-4 16:34

登录/注册后可看大图

嗯感谢大神的回答~：）~

= =我不是什么大神。。。

artrustee · 发表于 2012-9-4 20:41:22

wangjunxiong 发表于 2012-9-4 16:53

登录/注册后可看大图

点击图片就可以放大看了！蛮清晰的！
原谅我！我没有保存。

其实方程出来了，有很多数值解法。。。

maniux · 发表于 2012-9-4 21:57:17

这种果断mathematica啊，费力写什么script……

panhao1 · 发表于 2012-9-7 10:44:32

求 f（θ）=θ*L/2/c 和 f（θ）=sin（θ/2) 两条曲线的交点

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

[VB & C#] 万能的ncf吧 求能解超越方程的大神！！

[VB & C#] 万能的ncf吧求能解超越方程的大神！！