| 查看: 1620 | 回复: 3 | ||||
津津爱睡觉金虫 (小有名气)
|
[求助]
matlab解含参隐函数方程 已有1人参与
|
|
我的方程是x=1/(x(2).*49.005^2).*[log((y./x(1)).*((x(1)-49.005)./(y-49.005))-49.005.*(1./y-1./x(1))]..x1、x2是要拟合的参数 实验数据：xdata=[0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.5,3,3.5,4,4.5,5,5.5,6,6.5] ydata=[83.94986152,69.18056238,63.78006779,62.16487246,60.65377495,59.89822619,58.46520206,57.7096533,56.95410455,55.769572,54.95441885,53.89832959,53.14278084,52.14041949,51.64595481,51.34373531,51.11791018,50.85514711]..我需要拟合出了方程中的参数x1,x2然后根据xdata求y值并算相关系数。求帮助。。急求。。快被整疯了  |
回复此楼» 猜你喜欢
265求调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有21人回复
-
295材料求调剂，一志愿武汉理工085601专硕
已经有4人回复
-
0856调剂，是学校就去
已经有9人回复
-
298-一志愿中国农业大学-求调剂
已经有9人回复
-
求调剂
已经有4人回复
-
学校已经提交到NSFC，还能修改吗？
已经有7人回复
-
材料学硕297已过四六级求调剂推荐
已经有6人回复
-
295复试调剂
已经有7人回复
-
本人考085602 化学工程 专硕
已经有20人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
感谢参与，应助指数 +1
|
先拟合出： x2 = 0.000124354608284474 x1 = 71.0611809557734 相关系数(R): 0.99453642822085 相关系数之平方(R^2): 0.989102707058286 决定系数(DC): 0.989102707058287 代入原方程，再用MATLAB解方程： function solve_eqs clear all;clc format long global x y0=50; x_data=[0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.25,2.5,3,3.5,4,4.5,5,5.5,6,6.5]; for i=1:length(x_data) x=x_data(i); [y(i) fval(i)]=fsolve(@eqs,y0); end [x_data ;y; fval]' function f=eqs(y) global x f=1/(0.000124354608284474*49.005^2).*(log((y./71.0611809557734).*((71.0611809557734-49.005)./(y-49.005)))-49.005.*(1./y-1./71.0611809557734))-x; 计算结果： x y 方程残差 0 71.061023088663632 0.000011398222637 0.250000000000000 67.988292938180209 0.000000000096047 0.500000000000000 65.534689215008385 0.000000331981567 0.750000000000000 63.529133478234002 0.000000000137998 1.000000000000000 61.859615079896656 0.000000036989439 1.250000000000000 60.449405535313012 0.000003469065679 1.500000000000000 59.244034566066560 0.000000108193081 1.750000000000000 58.203506004930126 0.000000000309269 2.000000000000000 57.297823894942496 0.000000002283940 2.250000000000000 56.503949371087856 0.000000000099167 2.500000000000000 55.803867449155170 0.000000000001618 3.000000000000000 54.630610121986017 0.000000002101571 3.500000000000000 53.693164548810202 0.000000017349405 4.000000000000000 52.934038863155749 0.000000762365175 4.500000000000000 52.312777591278476 0.000000000526119 5.000000000000000 51.800007824280115 0.000000000001291 5.500000000000000 51.373858765788839 0.000000363302177 6.000000000000000 51.017693708877445 0.000000025319006 6.500000000000000 50.718627743671938 0.000000000780623 |
2楼2015-07-20 13:04:09
津津爱睡觉
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1100.8
- 散金: 33
- 红花: 1
- 帖子: 145
- 在线: 60小时
- 虫号: 2029058
- 注册: 2012-09-25
- 性别: GG
- 专业: 信号理论与信号处理
3楼2015-07-20 14:18:59
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-07-20 14:25:09
