24小时热门版块排行榜    

查看: 2571  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 6/300 2026-09-02 23:59 by rM1TE0WVDIIY
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +14 Kittylucky 2026-08-27 15/750 2026-09-02 19:35 by liumei0601
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
信息提示
请填处理意见