| 查看: 4853 | 回复: 8 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
origin 非线性自定义拟合问题~~方程组编写及运行~~
|
|||||
|
小弟专业是有机化学，现在课题设计到一个用origin进行非线性拟合的问题~~挣扎了很久也不会~~数学程序底子薄~~急求帮助~~目前只知道应该在origin里进行自定义编辑公式~~可是这个公式就是编不了，请高手帮忙~~里面的β是我要求的量，里面的T（z）、Z是我的实验测得的数据量，其他的都为已知的量，其中L eff和L是近似相等的，可以当做一个量。还有公式里面有积分，是不是还要把要将公式进行一番推导变换呢？帮我看看~~还有这个公式怎么才能编辑出来~~用origin C 语言怎么表示？如能帮小弟大忙，感激不尽~~如果来北京~~小弟请客~~ 公式见附图中公式（1）、（2）、（3）~~  计算截图.JPG |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有158人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有122人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求高手相助：origin非线性拟合 新建公式怎么编辑
已经有8人回复
- 毕不了业了，发帖求助Origin8.0多元非线性拟合的问题
已经有9人回复
- origin非线性拟合
已经有5人回复
- 急：求救Origin8.5自定义非线性模拟的问题
已经有7人回复
- origin8.5自定义拟合公式并拟合
已经有7人回复
- 请问有人会用origin拟合数据用自定义的公式么？急急！
已经有8人回复
- origin的非线性拟合的重复性问题，已经R-square 的标准
已经有5人回复
- origin进行自定义非线性拟合
已经有15人回复
- 求助origin非线性拟合实验数据
已经有12人回复
- 如何用origin进行非线性拟合？
已经有9人回复
- origin进行非线性拟合如何选择函数或者自定义合适的函数
已经有19人回复
- matlab非线性参数拟合问题
已经有7人回复
- 求助，origin7.5 中的非线性拟合参数为什么一直在变，而且误差很大！！！
已经有4人回复
- 求助：origin8.0 自定义非线性拟合，参数拟合值与设定初始值一样
已经有18人回复
- origin非线性拟合
已经有4人回复
- 关于origin里面的非线性拟合的问题
已经有6人回复
- 求助非线性拟合问题
已经有5人回复
- 求助origin高手 非线性拟合问题
已经有6人回复
- origin中非线性拟合中logistic模型的疑问？
已经有6人回复
- 请教关于origin中自定义函数非线性回归拟合问题！
已经有10人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68178.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
jiqiang: 金币+50, ★有帮助, 还是很谢谢你~~~ 2013-05-31 23:11:14
csgt0: 金币+2, 谢谢 2013-06-03 14:39:36
jiqiang: 金币+50, ★有帮助, 还是很谢谢你~~~ 2013-05-31 23:11:14
csgt0: 金币+2, 谢谢 2013-06-03 14:39:36
|
对于第一组数据，按z0= 0.0015117处理 参数 Beta = 0.0000013251423784 各常数数值都很大，初值不好给出，MATLAB程序列如下，仅供参考，高版本的1stopt跑出来效果应该会好很多。 %-------------------------------------------------------------------------------------------------- function feixianxingnihe314 %clear all;clc format long zspan=[-40 -37.9375 -35.875 -33.8125 -31.75 -29.6875 -27.625 -25.5625 -23.5 -21.4375 -19.375 -17.3125 -15.25 -13.1875 -11.125 -9.0625 -7 -4.9375 -2.875 -0.8125 1.25 3.3125 5.375 7.4375 9.5 11.5625 13.625 15.6875 17.75 19.8125 21.875 23.9375 26 28.0625 30.125 32.1875 34.25 36.3125 38.375 40.4375 ]; Texp=[ 0.99345 0.99 0.98806 0.98477 0.98216 0.98184 0.97636 0.97683 0.96763 0.96138 0.9594 0.9472 0.92474 0.89896 0.8736 0.85005 0.80012 0.7539 0.74086 0.678 0.60586 0.68747 0.7964 0.85913 0.89399 0.90719 0.93361 0.95553 0.96915 0.97888 0.98453 0.98644 0.98746 0.98978 0.99194 0.99773 0.99861 1.00266 1.00016 1.00033 ]; %E的数据，在此输入 k0=1e-6; lb=0; ub=inf; %------------------------------------------------------------------------- % 使用函数lsqnonlin()进行参数估计 OPTIONS=optimset('MaxFunEvals',1000); [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,OPTIONS,Texp); %ci = nlparci(k,residual,jacobian); %residual; fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n') fprintf('\n\t参数 D = %.16f',k(1)) y=KineticsEqs(k); R2=1-sum((Texp-y').^2)./sum((Texp'-mean(y)).^2); fprintf('\n\t相关系数之平方R^2 = %.16f',R2); figure plot(zspan,KineticsEqs(k),'b',zspan,Texp,'or'),legend('计算值','实验值','Location','Best') %------------------------------------------------------------------------- function f = ObjFunc(k,Texp) f=KineticsEqs(k)'-Texp; %------------------------------------------------------------------------ function xt = KineticsEqs(k) zspan=[-40 -37.9375 -35.875 -33.8125 -31.75 -29.6875 -27.625 -25.5625 -23.5 -21.4375 -19.375 -17.3125 -15.25 -13.1875 -11.125 -9.0625 -7 -4.9375 -2.875 -0.8125 1.25 3.3125 5.375 7.4375 9.5 11.5625 13.625 15.6875 17.75 19.8125 21.875 23.9375 26 28.0625 30.125 32.1875 34.25 36.3125 38.375 40.4375 ]; I0=1.492e12; Leff=1; z0= 0.0015117; %z0=4.2764e-16; %q0=k(1)*I0*Leff./(1+zspan.^2./(z0^2)); n=length(zspan); for i=1:n; q0=k(1)*I0*Leff./(1+zspan(i).^2./(z0^2)); fun=@(x) log(1+q0.*exp(-x.^2)); integ(i)=quad(fun,-1e5,1e5); xt(i)=1./sqrt(pi)./q0.* integ(i); end  附图1.jpg |
7楼2013-05-31 15:25:44
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68178.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
jiqiang: 金币+20, 请求帮助~~最好能是origin拟合~~小弟是学化学的，常用是的origin软件~~其他的软件更不会了 2013-05-30 23:22:47
csgt0: 金币+1, 谢谢 2013-05-31 11:41:20
感谢参与，应助指数 +1
jiqiang: 金币+20, 请求帮助~~最好能是origin拟合~~小弟是学化学的，常用是的origin软件~~其他的软件更不会了 2013-05-30 23:22:47
csgt0: 金币+1, 谢谢 2013-05-31 11:41:20
|
你把需要拟合的数据也发上来吧。 常用的最强大的三款曲线拟合软件：1stopt，MATLAB，Origin应该都能解决你的问题。 1stopt就拟合效果而言，几乎是最强大的，可以直接拟合带积分项的公式。 MATLAB自己写个小程序，就能把积分项包含进拟合公式中，或者考虑用MATLAB自带的cftool工具箱。 |
2楼2013-05-30 08:52:41
3楼2013-05-30 23:17:52
| 这是要用的数据文件~~ |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 计算数据信息.doc
- 附件 2 : 数据1.xls
- 附件 3 : 数据2.xls
- 附件 4 : 数据3.xls
2013-05-30 23:20:38, 70 K
2013-05-30 23:20:39, 9.5 K
2013-05-30 23:20:41, 10.5 K
2013-05-30 23:20:42, 9.5 K
4楼2013-05-30 23:20:52
