NCF参数化建筑论坛

标题: 万能的ncf吧求能解超越方程的大神！！ [打印本页]

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

作者: wangjunxiong 时间: 2012-9-4 11:58
本帖最后由 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
化成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
用别人的公式解决了。

[attach]34749[/attach]

[attach]34750[/attach]

作者: wangjunxiong 时间: 2012-9-4 16:27
用rhinoscript 也做了一个，运行速度还是可以的！修改一些参数就可以为自己所用了。

[attach]34751[/attach]

[attach]34752[/attach]

作者: yukai881208 时间: 2012-9-4 16:32

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

登录/注册后可看大图

用别人的公式解决了。

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

作者: yukai881208 时间: 2012-9-4 16:34

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
本帖最后由 wangjunxiong 于 2012-9-4 17:00 编辑

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

登录/注册后可看大图

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

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

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

作者: yukai881208 时间: 2012-9-4 16:59

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

登录/注册后可看大图

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

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

作者: artrustee 时间: 2012-9-4 17:09

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

登录/注册后可看大图

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

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

作者: artrustee 时间: 2012-9-4 20:41

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

登录/注册后可看大图

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

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

作者: maniux 时间: 2012-9-4 21:57
这种果断mathematica啊，费力写什么script……

作者: panhao1 时间: 2012-9-7 10:44
求 f（θ）=θ*L/2/c 和 f（θ）=sin（θ/2) 两条曲线的交点

欢迎光临 NCF参数化建筑论坛 (http://ncf-china.com/)