| 查看: 2462 | 回复: 6 | ||
| 【悬赏金币】回答本帖问题，作者menthalj将赠送您 10 个金币 | ||
[求助]
求教：matlab怎么自定义积分方程来拟合参数。已有2人参与
|
||
想用下面的方程来拟合一组数据，不知道怎么编程。 可以在curve fitting里直接输吗？麻烦大神能帮忙给出这个公式的程序。。。跪谢 QQ截图20180813152225.jpg |
回复此楼» 本帖已获得的红花（最新10朵）
menthalj» 猜你喜欢
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有92人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
超导三元凸包图
已经有0人回复
-
最近遇到这样一个问题
已经有10人回复
-
论文接收，散金币了！
已经有171人回复
-
vasp算居里温度
已经有0人回复
-
郑州大学光电器件方向招收一名博士
已经有6人回复
2楼2018-08-13 15:27:49
独孤神宇
版主 (知名作家)
- 应助: 488 (硕士)
- 贵宾: 0.008
- 金币: 29512.8
- 散金: 802
- 红花: 120
- 沙发: 1
- 帖子: 5333
- 在线: 850.2小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
如果是为了学习，可以参考下面例子 ： http://www.ilovematlab.cn/forum. ... &tid=514429 若是为了计算出参数，可以将数据一并给出。 |
3楼2018-08-14 21:19:20
hzlhm
至尊木虫 (著名写手)
- 应助: 372 (硕士)
- 金币: 16128.8
- 红花: 51
- 帖子: 2560
- 在线: 571.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+1, 鼓励交流 2018-08-17 07:47:33
独孤神宇: 金币+1 2018-08-17 10:01:12
menthalj: 金币+10 2018-08-23 15:40:51
感谢参与，应助指数 +1
月只蓝: 金币+1, 鼓励交流 2018-08-17 07:47:33
独孤神宇: 金币+1 2018-08-17 10:01:12
menthalj: 金币+10 2018-08-23 15:40:51
|
用拟合工具箱（curve fitting）目前还不能解决带有积分方程形式拟合问题。但可以用matlab自带的拟合函数（如nlinfit）来解决。解决该问题，首先自定义积分方程函数，其主要代码为 syms t y=c(1)*int(exp(c(2)/t^2),0,x) 然后，在主程序中调用其函数 x=[。。。] y=[。。。] c0=[c10 c20] c= nlinfit(x,y,func,c0) |
4楼2018-08-15 11:42:21
5楼2018-08-23 15:43:28
6楼2018-08-23 15:45:51
hzlhm
至尊木虫 (著名写手)
- 应助: 372 (硕士)
- 金币: 16128.8
- 红花: 51
- 帖子: 2560
- 在线: 571.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
【答案】应助回帖
|
分析你给出的x,y数据，是无法拟合出其自定义积分函数的系数。所以需对你的数据作一个变换，即x1=ln(x)，y1=ln(y)。然后再按给出的数学模型进行拟合计算。 x=[。。。]; y=[。。。]; x1=log(x); y1=log(y); x0=[-25 -32] [c,resnorm,residual,exitflag]=lsqcurvefit(@func,x0,x1,y1); C1=c(1);C2=c(2); str = ['C1=',num2str(C1),'，C2=',num2str(C2)]; fprintf('%s\n',str); x2=x1; y2=func(c,x2); B=[y1' y2']; disp(' 原始数据 拟合数据') disp(B) err=max(double(y1-y2)); str = ['最大误差=',num2str(err)]; fprintf('%s\n',str); R2=R2_func(y1,y2); x3=exp(x2); y3=exp(y2); plot(x,y,'rp',x3,y3,'k-'),grid on xlabel('x'),ylabel('y')  lhm1.jpg |
7楼2018-08-26 16:01:00
50