| 查看: 3085 | 回复: 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 |
回复此楼» 猜你喜欢
投稿咨询
已经有7人回复
-
今天放榜没戏了吧
已经有11人回复
-
我面上完蛋了
已经有5人回复
-
让我中一个面上吧！
已经有12人回复
-
基金啊基金
已经有4人回复
-
放榜前的不淡定
已经有13人回复
-
时间戳今天，20号变了
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有24人回复
-
估计是周四
已经有3人回复
-
科研孤儿太难了
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
hytao2012
铁杆木虫 (正式写手)
木头虫子
- 应助: 53 (初中生)
- 金币: 6609.8
- 散金: 115
- 红花: 11
- 帖子: 540
- 在线: 208.4小时
- 虫号: 2050091
- 注册: 2012-10-08
- 性别: GG
- 专业: 固体力学
【答案】应助回帖
★
jjdg: 金币+1, 辛苦了 2014-08-04 15:04:21
jjdg: 金币+1, 辛苦了 2014-08-04 15:04:21
|
我也写了一个小程序，感觉好像是因为你的方程太复杂，解不出来啊，注释的是我写了个简单的方程，就可以了。 clear; clc; syms N g D Ph Q R T To f Y K f=1-exp(-N*g*(D/Ph)^1.5*(exp(-Q/(R*T))*R*T^2/Q*(1-exp(Q*(To-T)/R/T^2)))^1.5); % f=D^2+Q^2; K=sum((f-Y)^2); dD=diff(K, D); dQ=diff(K, Q); Result=solve(dD==0, dQ==0, D, Q) % Result.D % Result.Q |
9楼2014-08-04 09:04:21
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
40