| 查看: 937 | 回复: 3 | |||
zw073916木虫 (著名写手)
|
[求助]
求助Matlab编程，万分感谢 已有1人参与
|
|
有2组数： x: 0 10 40 110 155 215 310 410 505 605 700 800 900 y: 0 0.00468 0.00704 0.00996 0.01132 0.01308 0.01568 0.01868 0.02148 0.02468 0.02828 0.03292 0.03888 用方程 y’=(p(1)*sinh((p(2)*130*(1-p(5)*(1-exp(-(p(4)*y/(p(5)*130))))))/(((1+p(6)*x)^(-1/3))*(1-p(3)*y)))) 进行拟合 求解参数P(1)-p(6)？ y'为y的导数，上面方程可以表示为y'=f(x,y)由于这是个隐函数，用origin拟合不了，据说可以用matlab用最小二乘法进行非线性拟合，可是我不会。哪位大牛能帮忙把程序写一下？ |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有52人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助各位大神，帮忙用MATLAB编程数出一个立方体内所有颗粒的个数
已经有0人回复
- 求matlab wolfe线搜索下DY共轭梯度法程序，新人求助，万分感谢
已经有1人回复
- 求助matlab牛人随机生成序列相邻两数不相等
已经有2人回复
- matlab动力学模型编程求助
已经有6人回复
- 求助matlab编程问题~
已经有1人回复
- 求助高手帮忙运用matlab分析数据，得出回归系数
已经有2人回复
- 求助代码
已经有0人回复
- 【求助】求助 matlab编程：可否实现清除文本文档数据
已经有5人回复
- matcom混合编程求助
已经有1人回复
- 【求助】向高手求助matlab编码
已经有2人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
zw073916: 金币+100, ★★★★★最佳答案 2016-05-16 08:41:48
感谢参与，应助指数 +1
zw073916: 金币+100, ★★★★★最佳答案 2016-05-16 08:41:48
2楼2016-05-14 23:15:42
zw073916
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1475.9
- 散金: 8397
- 红花: 1
- 帖子: 1236
- 在线: 144.2小时
- 虫号: 3433901
- 注册: 2014-09-22
- 性别: GG
- 专业: 金属材料的力学行为
|
我按照你给我的代码修改后，怎么计算不了 clear all;clc format long tspan=[20 50 120 165 225 320 420 515 615 710 810 910]-10; yexp=[0.00468 0.00704 0.00996 0.01132 0.01308 0.01568 0.01868 0.02148 0.02468 0.02828 0.03292 0.03888]'; k0=[1 1 1 1 1 1]; %%%请注意这里，初值的选取 y0=0; lb=-[1 1 1 1 1 1]*1e3; ub=[1e12 1e9 1e6 1000 10 1]; yy=[y0 yexp']; [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp); fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n') fprintf('\t待拟合参数 k1 = %.6f\n',k(1)) fprintf('\t待拟合参数 k2 = %.6f\n',k(2)) fprintf('\t待拟合参数 k3 = %.6f\n',k(3)) fprintf('\t待拟合参数 k4 = %.6f\n',k(4)) fprintf('\t待拟合参数 k5 = %.6f\n',k(5)) fprintf('\t待拟合参数 k6 = %.6f\n',k(6)) fprintf(' \t残差平方和= %.6f\n\n',resnorm) ts=0:1:max(tspan); [ts ys]=ode45(@KineticsEqs,ts,y0,[],k); [ttt XXsim] = ode45(@KineticsEqs,tspan,y0,[],k); y=XXsim(2:end); xexp=yexp; R2=1-sum((xexp-y).^2)./sum((xexp-mean(y)).^2); fprintf('\n\t决定系数R-Square = %.6f',R2); figure(1) plot(ts,ys,'b',tspan,yy,'or'),legend('计算值','实验值','Location','best'); yr=y-yexp; figure(2) plot(tspan(2:end),yr,'r*',[-100 1000],[0 0]),axis([-100 1000 -50 50]); figure(3) plot(yexp,y,'ro',[0 0.05],[0 0.05],'b-'); %--------------------------------------------------------- function f = ObjFunc(k,tspan,y0,yexp) [t Xsim] = ode45(@KineticsEqs,tspan,y0,[],k) ; ysim = Xsim(2:end); size(ysim); size(yexp); f=ysim-yexp; %---------------------------------------------------------- function dydt = KineticsEqs(t,y,k) beta(1)=k(1); beta(2)=k(2); beta(3)=k(3); beta(4)=k(4); beta(5)=k(5); beta(6)=k(6); dydt = (beta(1)*sinh((beta(2)*130*(1-beta(5)*(1-exp(-(beta(4)*y/(beta(5)*130))))))/(((1+beta(6)*t)^(-1/3))*(1-beta(3)*y)))); |
3楼2016-05-16 18:12:58
zw073916
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1475.9
- 散金: 8397
- 红花: 1
- 帖子: 1236
- 在线: 144.2小时
- 虫号: 3433901
- 注册: 2014-09-22
- 性别: GG
- 专业: 金属材料的力学行为
4楼2016-05-16 18:14:03
30