24小时热门版块排行榜    

查看: 2567  |  回复: 0

fengxing24

新虫 (初入文坛)

[求助] matlab曲线拟合过端点求拟合函数,新手求 指导

本人新手,做论文急需拟合下面两组数据的函数,偶函数,函数为要求过tdata的两个端点,即A [23.8188 0]和B [0 -15.654]两点,且在A点一阶导数为正无穷,即垂直于x轴;B点一阶导数为0,即垂直于y轴;
matlab曲线拟合过端点求拟合函数,新手求 指导

下面我是做的非线性拟合,不知道怎么做才过端点且满足端点的导数??恳求帮助指导
函数文件curvefun1.m
unction f=curvefun1(x,tdata)
% f=x(1)+x(2)*exp(-0.02*x(3)*tdata).^2+x(4)*tdata.^2
f=x(1)+x(2)*exp(x(3)*tdata.^2)+x(4)*exp(x(5)*tdata.^2)+x(6)*tdata.^2%f=A+BEXP(C*X^2)+D*EXP(E*X^2)+F*x^2

下面是计算数据[tdata y]

%%按照直角坐标拟合
tdata=[23.8188
23.8342
23.7987
23.7929
23.8527
23.7062
23.6544
23.5103
23.2279
23.0524
23.0524
22.4726
22.3073
21.9304
21.621
21.1967
20.9019
20.6114
20.2232
19.7369
19.3375
19.1869
18.6443
17.8836
17.2772
17.1087
16.5737
16.0344
15.7027
15.0864
14.7988
14.5615
14.0628
13.4977
13.3248
12.8461
12.393
12.1607
11.9807
11.5081
11.2076
10.9033
10.822
10.2442
9.6234
9.5353
8.9276
8.2954
8.155
7.4794
7.0256
6.6483
6.3191
5.6674
4.9366
4.9366
4.7496
4.1897
3.7401
3.3996
3.1933
2.9862
2.4351
2.1442
1.3893
1.374
0.6974
0.1026
0


];
y=[0
-0.3221
-1.1519
-2.0705
-3.0347
-3.2926
-3.9812
-4.8063
-6.0083
-6.6225
-6.6225
-8.043
-8.3916
-8.9705
-9.5804
-10.1079
-10.538
-10.8858
-11.3962
-11.8158
-12.3589
-12.4146
-13.0249
-13.685
-14.1091
-14.155
-14.4614
-14.8013
-14.9063
-15.2473
-15.2502
-15.4975
-15.7695
-15.9686
-15.9417
-16.0339
-16.1741
-16.1872
-16.2183
-16.2189
-16.275
-16.2563
-16.2056
-16.2188
-16.1824
-16.2529
-16.2532
-16.2354
-16.1712
-16.2036
-16.2777
-16.2136
-16.1762
-16.1061
-16.0701
-16.0701
-16.1101
-16.079
-15.9136
-15.9923
-15.9638
-15.9691
-15.9098
-15.8816
-15.8233
-15.9487
-15.8259
-15.7887
-15.654
];
x0=[0.2,0,0,0.0005,0.005,-1];
x=lsqcurvefit(\'curvefun1\',x0,tdata,y)
f=curvefun1(x,tdata);
m=f-y0
plot(tdata,y,\'o\',tdata,f,\'r\');
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 fengxing24 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +4 myrtle 2026-08-26 4/200 2026-08-31 08:43 by 逸域隐帆
[基金申请] 中青基了要发朋友圈吗? +6 349506619 2026-08-28 6/300 2026-08-31 08:43 by kingkocxr
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 08:28 by xiangy672
[基金申请] 面上意见出来了 +10 黄鸟于飞Chao 2026-08-29 19/950 2026-08-31 08:01 by blueearth171
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见