24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 1998  |  回复: 9
【悬赏金币】回答本帖问题,作者646452271将赠送您 5 个金币

646452271

铁虫 (初入文坛)

[求助] 初学Matlab,函数图像用plot画出来了,但是不知道拟合,希望可以得到帮助已有1人参与

各位大大好,我最近在做一个铁路上的项目,需要用到Matlab进行计算,因为是刚开始学习这个软件,有些地方不是很清楚,希望能够得到各位大佬的帮助。

       我做的这个项目第一步是需要建立钢轨曲线,曲线是三段相切圆弧得到。 因为只是仿真,不用完整铁路模型,因此第一段圆弧圆心的位置我自己确定。 第一、二、三段圆弧半径,以及圆弧两个交点横坐标都给出了,就是图上红圈标注的六个尺寸。这条曲线我通过之前的学习利用分段函数已经做出来了。如下图一表示。

但是现在我需要做的任务是已知点l1的横坐标,同时知道红色阴影面积s1=s2=s3=....=a。连续求 l2,l3 ,l4...的坐标。因为 a 会取很小,所以l2,l3,l4...应该是需要用到循环来做。

      而我面临的问题就是两个:
      1、我现在做的函数是分段函数,如果带入循环计算考虑区间肯定会相当复杂,这条线我认为是需要连续求导,因此需要拟合。所以就需要先导出数据点再导入拟合,这一部分我还是学的有点云里雾里的,不太会用。
      2、第二个最难的就是这个循环函数了,因为确实项目时间有点短,所以希望求助大家。当然我也不是伸手党,如果大神觉得一步一步教我太慢或者觉得我这样伸手要结果不应该,也希望能为我指条路,比如有什么类似的帖子或者教程,我去学也行,只是希望能尽快做出这个函数。

     再次感谢各位能够阅读我的帖子,也希望能够得到你们的帮助,谢谢。

初学Matlab,函数图像用plot画出来了,但是不知道拟合,希望可以得到帮助
QQ截图20180704172032.png


初学Matlab,函数图像用plot画出来了,但是不知道拟合,希望可以得到帮助-1
QQ截图20180704172936.png


初学Matlab,函数图像用plot画出来了,但是不知道拟合,希望可以得到帮助-2
QQ图片20180704173435.jpg
回复此楼

» 猜你喜欢

一入机械深似海
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

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

发自小木虫Android客户端
数值计算
2楼2018-07-04 18:46:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

646452271

铁虫 (初入文坛)

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

因为已知条件好像没法做一个连续函数,我是用几何知识确定的第二个圆弧和第三个圆弧圆心,但是不能直接算出来,因为算出来就会近似取值,交点就不相切,所以只能选择这个方式来求曲线。
一入机械深似海
3楼2018-07-04 22:57:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

求学者hyt

银虫 (小有名气)

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

发自小木虫Android客户端
4楼2018-07-05 07:04:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

646452271

铁虫 (初入文坛)

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

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

发自小木虫Android客户端
一入机械深似海
5楼2018-07-05 12:04:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

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

你图1中已经绘制出了三段圆弧,你可以将整个圆弧的数据整理一下附上来,说不定可以用一个公式得到较好的拟合效果呢。
数值计算
6楼2018-07-05 14:42:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

求学者hyt

银虫 (小有名气)

引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

646452271

铁虫 (初入文坛)

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

不好意思,前两天出差了,才回来,请见谅。圆弧的数据其实在我给的第二个图中就有了,第一个圆弧半径是三百,圆心位置我选取的是(0,-285),然后在距离Y轴10mm的位置和这段圆弧存在一交点,过这一点做圆弧1的相切圆弧,这段圆弧的半径为80mm,之后,在距离Y轴25.35mm点线和圆弧二存在交点,过该点做圆弧3,半径为13.
一入机械深似海
8楼2018-07-09 15:25:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

646452271

铁虫 (初入文坛)

引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

求学者hyt

银虫 (小有名气)


独孤神宇: 金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 646452271 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] AC审稿 +3 马winner 2024-04-26 7/350 2024-04-27 23:33 by topedit
[硕博家园] 聊天 +13 暮色恋伊人 2024-04-22 14/700 2024-04-27 23:26 by @tanzelin
[硕博家园] 有没有想聊天的 +5 兔跑儿2024 2024-04-27 6/300 2024-04-27 22:49 by 简风23
[考博] 博士招生-211农业院校 +5 NYC917 2024-04-26 10/500 2024-04-27 22:31 by NYC917
[论文投稿] 求论文投稿期刊推荐 15+3 mihudawang 2024-04-26 5/250 2024-04-27 22:30 by lizhengke06
[找工作] 普通院校药学硕士,做合成的,感觉找不到工作 +12 pom戴墨镜 2024-04-24 20/1000 2024-04-27 22:08 by saenpe
[硕博家园] 博士白读了 +49 Da_Meng_Zi 2024-04-21 54/2700 2024-04-27 20:12 by 慕容家族
[基金申请] 优青已经开始了么? +3 majunge000 2024-04-27 3/150 2024-04-27 18:51 by spring_1987
[高分子] 请问UV灯是365nm的,那么选光引发剂的波长选多少的?要完全一致吗? +5 engledd2004 2024-04-21 5/250 2024-04-27 18:32 by quirinus
[找工作] 江苏理工学院要慎选 +3 jjchenshui 2024-04-27 3/150 2024-04-27 17:45 by 阿司匹林888
[考博] 材料方向24博士申请/一作SCI三篇 +3 白天不碰 2024-04-24 5/250 2024-04-27 17:41 by 安塔瓦拉多
[考博] 25年博士申请 +7 Changzixuan 2024-04-25 13/650 2024-04-27 14:11 by 我属驴核动力驴
[有机交流] 环肽的合成 +3 徐来不惊 2024-04-25 6/300 2024-04-27 13:54 by xxzhangci
[基金申请] "颜宁:基础研究应顶天立地"能做到基础研究同时顶天立地的才是牛人 +6 zju2000 2024-04-24 6/300 2024-04-27 10:17 by kcmn1000
[论文投稿] 投稿RSC旗下杂志突然间看不到投稿状态,里面投稿的文章信息不见了,有遇到过的吗? 50+3 dlying 2024-04-26 3/150 2024-04-27 09:58 by bobvan
[考博] 申博求助 +4 dskabdh 2024-04-24 11/550 2024-04-26 15:54 by dskabdh
[基金申请] 国社科项目,你们学校都限额申报吗? +7 屡战屡败 2024-04-21 10/500 2024-04-25 12:10 by 屡战屡败
[电化学] 耗材发问 +4 Happy C 2024-04-22 4/200 2024-04-25 11:03 by 普通小虫
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +4 我属驴核动力驴 2024-04-22 5/250 2024-04-24 08:47 by 晓目崇
[考博] 研二光催化6月底4篇2区 +7 wjtab 2024-04-22 11/550 2024-04-23 06:59 by byron2012
信息提示
请填处理意见