| 查看: 1791 | 回复: 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值并算相关系数。求帮助。。急求。。快被整疯了  |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
现代”学阀”该如何界定
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
国社科系统bug了，是不是要放榜了？
已经有9人回复
-
申博发邮件
已经有9人回复
-
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有9人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有8人回复
-
上海工程技术大学激光智能制造课题组招收博士研究生
已经有8人回复
津津爱睡觉
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1100.8
- 散金: 33
- 红花: 1
- 帖子: 145
- 在线: 60小时
- 虫号: 2029058
- 注册: 2012-09-25
- 性别: GG
- 专业: 信号理论与信号处理
3楼2015-07-20 14:18:59
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68175.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.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
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68175.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-07-20 14:25:09