找回密码
 欢迎注册
楼主: 数学星空

[讨论] 两参数曲面的相交曲线计算问题

[复制链接]
发表于 2024-11-10 17:00:15 | 显示全部楼层
另外如果牛顿迭代法,为什么引入反三角函数呢,直接用\(f(t)=A \cos(t)\sin(t)+B\cos(t)+C\sin(t)+D\)不好吗?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2024-11-10 17:14:56 | 显示全部楼层
mathe 发表于 2024-11-10 17:00
另外如果牛顿迭代法,为什么引入反三角函数呢,直接用\(f(t)=A \cos(t)\sin(t)+B\cos(t)+C\sin(t)+D\)不好 ...

怎么直接利用f(x)构造合适的迭代公式呢?

点评

t_0任意选择,使用迭代$t_{n+1}=t_n - \frac{f(t_n)}{f'(t_n)}$  发表于 2024-11-12 09:55
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2024-11-11 09:14:01 | 显示全部楼层
数学星空 发表于 2024-11-10 13:35
关于两个一般椭圆台曲面相交的计算,现遇到两个很麻烦的计算问题需要探讨一下:

对第一个椭圆台方程:

第一问:解的是一个二元二次方程
  1. Solve[{s^2+c^2==1,1/2Sin[2 Subscript[t, 1]](Subscript[m, 0]+c Subscript[m, 1])+Cos[Subscript[t, 1]] (Subscript[n, 0]+c s Subscript[n, 1]+c Subscript[n, 2]+s Subscript[n, 3])+Sin[Subscript[t, 1]] (Subscript[p, 0]+c Subscript[p, 1])+Subscript[q, 0]+c s Subscript[q, 1]+c Subscript[q, 2]+s Subscript[q, 3]==0},{c,s}]
复制代码

点评

和mathe的万能公式代换是一样的 ,消元后是一个四次代数方程~若不依赖数学软件计算还是不方便  发表于 2024-11-11 09:44
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2024-11-12 09:57:59 | 显示全部楼层
牛顿迭代法用单变量还是多变量没什么区别,单变量我们会计算f(x)/f'(x), 而换成多变量,那么1/f'(x)需要替换为雅克比矩阵(一阶偏导)求逆。
主要问题是牛顿迭代法无法保证收敛,这个可以通过多次随机选择合适的起始点,直到找到可以收敛的初始值
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2024-11-13 19:45:59 | 显示全部楼层
mathe 发表于 2024-11-12 09:57
牛顿迭代法用单变量还是多变量没什么区别,单变量我们会计算f(x)/f'(x), 而换成多变量,那么1/f'(x)需要替 ...

对于下面的方程,如何构造一个合适的迭代公式,使计算更快速的收敛到零点?(在区间0~2pi内对于任意给定的t2,求t1?)
((1500*cos(t2)*sqrt(3) - 150000)*sin(t1) -| 150*cos(t2)*sqrt(3)*sin(t2) + 20000*cos(t2)*sqrt(3) + 475000*sin(t2) - 2000000)*cos(t1) + (-13500*cos(t2)*sqrt(3) - 37500000*sqrt(3) + 1575000*cos(t2) + 1350000)*sin(t1) + 1500*cos(t2)*sqrt(3)*sin(t2) + 37500000*sin(t2)*sqrt(3) + 50000*cos(t2)*sin(t2) - 12500000*sqrt(3) - 3375000*cos(t2) - 1275000*sin(t2)=0

点评

是否可以用渐进表达第一个初始值?  发表于 7 天前
牛顿法初一般情况始值好像没有特别好的选择  发表于 7 天前
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-21 17:50 , Processed in 0.022111 second(s), 17 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表