| 查看: 2582 | 回复: 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 |
回复此楼» 猜你喜欢
A期刊撤稿
已经有5人回复
-
临港实验室与上科大联培博士招生1名
已经有8人回复
-
26申博自荐
已经有7人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有4人回复
-
带资进组求博导收留
已经有9人回复
-
求助大佬们，伤口沾上了乙腈
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 编程小问题，求帮助，谢谢。
已经有13人回复
- 求高手给做一道题，用MATLAB或1stOpt，请把结果直接给出来，我的软件就不行，谢谢
已经有5人回复
- matlab动力学模型编程求助
已经有6人回复
- 求一个matlab详细代码
已经有5人回复
- 求回归方程参数的最优化算法matlab代码
已经有17人回复
- matlab编程齿轮
已经有9人回复
- 请高人指点一个matlab编程问题，谢谢
已经有5人回复
- 新手matlab求助，谢谢各位
已经有10人回复
- 求助用matlab或者（1st0pt)编程解方程组
已经有3人回复
- MATLAB高手帮我看看这个数据存的格式，谢谢
已经有3人回复
- 求基于DCT域的数字图像水印Matlab代码，谢谢各位大神！
已经有3人回复
- 请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！
已经有9人回复
- matlab编程大家看这两个式子该如何编写呢？
已经有3人回复
- 请问如何用Matlab编程做这个图？谢谢大家！
已经有5人回复
- 如何用matlab软件或者用C及C++自己手工编程寻找转折点
已经有12人回复
- 求助matlab编程
已经有13人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【资源】Matlab编程基础及应用(视频教程)
已经有50人回复
7楼2014-08-03 13:20:53
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 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