| 查看: 1621 | 回复: 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值并算相关系数。求帮助。。急求。。快被整疯了  |
回复此楼» 猜你喜欢
287求调剂
已经有5人回复
-
22408 344分 求调剂 一志愿 华电计算机技术
已经有3人回复
-
289求调剂
已经有3人回复
-
317求调剂
已经有10人回复
-
广西大学家禽遗传育种课题组2026年硕士招生（接收计算机专业调剂）
已经有3人回复
-
学校已经提交到NSFC，还能修改吗？
已经有8人回复
-
08工学调剂
已经有5人回复
-
281求调剂（0805）
已经有25人回复
-
085600材料与化工
已经有6人回复
-
265求调剂
已经有8人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +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
