| 查看: 3190 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
King_A金虫 (小有名气)
|
[求助]
求matlab编程代码，求帮助，谢谢！ 已有3人参与
|
||
|
有上述公式22，其中除f,D和Q外，其它均为已知量，其中T（t）为一维数组，此外还已知另一数组Y(t)。 现想进行一下计算：定义K=∑（f-Y(t)）^2，分别对K求D和Q的一阶偏导数，并另其一阶偏导数等于0，解方程组，最终得到D和Q 的值， 最后将计算所得的D和Q值返回到公式22中，画出(T,f)曲线。 自行研究多日，却不得其所，实在无奈，只得到此求助大神帮帮，在此先表谢意。  22.png |
回复此楼» 猜你喜欢
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有12人回复
-
极限拉扯，如何抉择
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
12楼2014-12-07 10:42:17
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20659.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.7小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2014-07-30 22:27:35
King_A
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 711
- 散金: 201
- 帖子: 214
- 在线: 97.7小时
- 虫号: 1130184
- 注册: 2010-10-23
- 专业: 金属材料的合金相、相变及
3楼2014-07-31 08:49:24
King_A
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 711
- 散金: 201
- 帖子: 214
- 在线: 97.7小时
- 虫号: 1130184
- 注册: 2010-10-23
- 专业: 金属材料的合金相、相变及
|
下面是我自行编辑的，请大神帮忙看下哪有错误。 clear all syms a b c d r N g x y f=1-exp(-N*g*(a/b)^1.5*(exp(-c/(r*x))*(r*(x^2)/c)*(1-exp(c*(d-x)/(r*(x^2)))))^(3/2)) k=sum(f-y)^2 z=diff(k,'a') h=diff(k,'b') [a,b]=solve('-(3*N*g*exp(-N*g*(a/b)^(3/2)*(-(r*x^2*exp(-c/(r*x))*(exp((c*(d - x))/(r*x^2)) - 1))/c)^(3/2))*(a/b)^(1/2)*(y + exp(-N*g*(a/b)^(3/2)*(-(r*x^2*exp(-c/(r*x))*(exp((c*(d - x))/(r*x^2)) - 1))/c)^(3/2)) - 1)*(-(r*x^2*exp(-c/(r*x))*(exp((c*(d - x))/(r*x^2)) - 1))/c)^(3/2))/b=0','(3*N*a*g*exp(-N*g*(a/b)^(3/2)*(-(r*x^2*exp(-c/(r*x))*(exp((c*(d - x))/(r*x^2)) - 1))/c)^(3/2))*(a/b)^(1/2)*(y + exp(-N*g*(a/b)^(3/2)*(-(r*x^2*exp(-c/(r*x))*(exp((c*(d - x))/(r*x^2)) - 1))/c)^(3/2)) - 1)*(-(r*x^2*exp(-c/(r*x))*(exp((c*(d - x))/(r*x^2)) - 1))/c)^(3/2))/b^2=0') 下面是结果： Warning: The solutions are parametrized by the symbols: z = C_ z11 = C_ z4 = ({0, d} union Dom::ImageSet(-(c*(((c + 8*PI*d*k*r*I)/c)^(1/2) - 1)*I)/(4*PI*k*r), k, Z_ minus {0}) union Dom::ImageSet((c*(((c + 8*PI*d*k*r*I)/c)^(1/2) + 1)*I)/(4*PI*k*r), k, Z_ minus {0})) minus {0} z5 = ({0, d} union Dom::ImageSet(-(c*(((c + 8*PI*d*k*r*I)/c)^(1/2) - 1)*I)/(4*PI*k*r), k, Z_ minus {0}) union Dom::ImageSet((c*(((c + 8*PI*d*k*r*I)/c)^(1/2) + 1)*I)/(4*PI*k*r), k, Z_ minus {0})) minus {0} > In solve at 190 a = z z4 z5 b = 1 - exp(-N*g*(a/b)^(3/2)*((r*z^2*exp(-c/(r*z)))/c - (r*z^2*exp((c*d)/(r*z^2) - c/(r*z))*exp(-c/(r*z)))/c)^(3/2)) z11 1 - exp(-N*g*(a/b)^(3/2)*((r*z5^2*exp(-c/(r*z5)))/c - (r*z5^2*exp((c*d)/(r*z5^2) - c/(r*z5))*exp(-c/(r*z5)))/c)^(3/2)) |
4楼2014-07-31 10:16:45
