【答案】应助回帖

» 猜你喜欢

北伊利诺伊大学(Northern Illinois University)机械系招收全奖博士生已经有0人回复
北伊利诺伊大学(Northern Illinois University)机械系招收全奖博士生已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有269人回复
第五届粉末衍射晶体结构解析及Rietveld精修高级研讨会已经有44人回复
FDTD软件里的时域图像画图已经有2人回复
复杂量子系统的新物态和新机理讲习班已经有0人回复
导电原子力显微镜cAFM如何测试mcAFM(即磁导AFM)性能？已经有0人回复
感谢小木虫，已完成目标已经有0人回复
零点能已经有2人回复
科研求助已经有0人回复
【25申博】博导推荐求助已经有4人回复

一入机械深似海

1楼2018-07-04 17:45:53

独孤神宇

版主 (知名作家)

应助: 489 (硕士)
贵宾: 0.008
金币: 29618.8
散金: 802
红花: 121
沙发: 1
帖子: 5357
在线: 850.4小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

为何不试试直接拟合圆弧曲线而选择分段

发自小木虫Android客户端

数值计算

2楼2018-07-04 18:46:45

646452271

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 24.5
帖子: 22
在线: 7.3小时
虫号: 2174212
注册: 2012-12-08
性别: GG
专业: 传动机械学

引用回帖:

2楼: Originally posted by 独孤神宇 at 2018-07-04 18:46:45
为何不试试直接拟合圆弧曲线而选择分段

因为已知条件好像没法做一个连续函数，我是用几何知识确定的第二个圆弧和第三个圆弧圆心，但是不能直接算出来，因为算出来就会近似取值，交点就不相切，所以只能选择这个方式来求曲线。

赞一下(1人)

一入机械深似海

3楼2018-07-04 22:57:56

求学者hyt

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 1165.9
散金: 5
红花: 2
帖子: 230
在线: 12.2小时
虫号: 4116943
注册: 2015-10-03
专业: 结构工程

感觉你把分段函数单独写成一个m文件，直接调用就好了，不存在什么复杂的。拟合是有误差的，你还要进行误差分析。

发自小木虫Android客户端

4楼2018-07-05 07:04:30

646452271

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 24.5
帖子: 22
在线: 7.3小时
虫号: 2174212
注册: 2012-12-08
性别: GG
专业: 传动机械学

引用回帖:

4楼: Originally posted by 求学者hyt at 2018-07-05 07:04:30
感觉你把分段函数单独写成一个m文件，直接调用就好了，不存在什么复杂的。拟合是有误差的，你还要进行误差分析。

额，还能这样，那应该如何求导呢？感觉有点不太清楚

发自小木虫Android客户端

一入机械深似海

5楼2018-07-05 12:04:04

独孤神宇

版主 (知名作家)

应助: 489 (硕士)
贵宾: 0.008
金币: 29618.8
散金: 802
红花: 121
沙发: 1
帖子: 5357
在线: 850.4小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

引用回帖:

3楼: Originally posted by 646452271 at 2018-07-04 22:57:56
因为已知条件好像没法做一个连续函数，我是用几何知识确定的第二个圆弧和第三个圆弧圆心，但是不能直接算出来，因为算出来就会近似取值，交点就不相切，所以只能选择这个方式来求曲线。...

你图1中已经绘制出了三段圆弧，你可以将整个圆弧的数据整理一下附上来，说不定可以用一个公式得到较好的拟合效果呢。

数值计算

6楼2018-07-05 14:42:02

求学者hyt

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 1165.9
散金: 5
红花: 2
帖子: 230
在线: 12.2小时
虫号: 4116943
注册: 2015-10-03
专业: 结构工程

引用回帖:

5楼: Originally posted by 646452271 at 2018-07-05 12:04:04
额，还能这样，那应该如何求导呢？感觉有点不太清楚
...

如果函数连续，切可导，那就用最基本的导数定义啊
比如你有f(x)的m文件
导数f’(x)=(f(x+h)-f(x))/h
h取很小就可以了

发自小木虫Android客户端

7楼2018-07-06 16:09:48

646452271

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 24.5
帖子: 22
在线: 7.3小时
虫号: 2174212
注册: 2012-12-08
性别: GG
专业: 传动机械学

引用回帖:

6楼: Originally posted by 独孤神宇 at 2018-07-05 14:42:02
你图1中已经绘制出了三段圆弧，你可以将整个圆弧的数据整理一下附上来，说不定可以用一个公式得到较好的拟合效果呢。...

不好意思，前两天出差了，才回来，请见谅。圆弧的数据其实在我给的第二个图中就有了，第一个圆弧半径是三百，圆心位置我选取的是(0,-285)，然后在距离Y轴10mm的位置和这段圆弧存在一交点，过这一点做圆弧1的相切圆弧，这段圆弧的半径为80mm，之后，在距离Y轴25.35mm点线和圆弧二存在交点，过该点做圆弧3，半径为13.

赞一下(1人)

一入机械深似海

8楼2018-07-09 15:25:43

646452271

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 24.5
帖子: 22
在线: 7.3小时
虫号: 2174212
注册: 2012-12-08
性别: GG
专业: 传动机械学

引用回帖:

7楼: Originally posted by 求学者hyt at 2018-07-06 16:09:48
如果函数连续，切可导，那就用最基本的导数定义啊
比如你有f(x)的m文件
导数f’(x)=(f(x+h)-f(x))/h
h取很小就可以了
...

对的，但是我现在做的是分段函数，在数据上这一点应该是相切的，但是我不知道如果做成连续函数是否相切，因为我在画这个plot函数的时候切点是我个人定义的，不知道如果做成.m文件是否能够同样适用

一入机械深似海

9楼2018-07-09 15:31:52

求学者hyt

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 1165.9
散金: 5
红花: 2
帖子: 230
在线: 12.2小时
虫号: 4116943
注册: 2015-10-03
专业: 结构工程

★
独孤神宇: 金币+1 2018-07-10 16:39:54

引用回帖:

9楼: Originally posted by 646452271 at 2018-07-09 15:31:52
对的，但是我现在做的是分段函数，在数据上这一点应该是相切的，但是我不知道如果做成连续函数是否相切，因为我在画这个plot函数的时候切点是我个人定义的，不知道如果做成.m文件是否能够同样适用...

间隔足够小是可以的求出切线斜率的，或者还有另一个方法，就是你既然有分段函数了，可以分段求导，然后再进行下一步运算。

另外我说的所谓的m文件意思是定义一个
function f=f(x)
这种，之后你用时调用即可。而且数值求导也比较方便。

发自小木虫Android客户端

10楼2018-07-10 12:17:16