应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 求matlab编程代码，求帮助，谢谢！
51/1返回列表
查看: 2581  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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)曲线。

自行研究多日，却不得其所，实在无奈，只得到此求助大神帮帮，在此先表谢意。

求matlab编程代码，求帮助，谢谢！
22.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-07-30 08:39:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hytao2012

铁杆木虫 (正式写手)

木头虫子

【答案】应助回帖

★
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
赞一下(1人) 回复此楼
9楼2014-08-04 09:04:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

jerkwin

专家顾问 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
做最小二乘拟合吧
赞一下 回复此楼
2楼2014-07-30 22:27:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

King_A

金虫 (小有名气)
引用回帖:
2楼: Originally posted by jerkwin at 2014-07-30 22:27:35
做最小二乘拟合吧

最小二乘拟合达不到我的要求啊，我还需要具体的D和Q的值。最小二乘拟合只能得到方程，得不到具体D和Q
赞一下 回复此楼
3楼2014-07-31 08:49:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

King_A

金虫 (小有名气)
下面是我自行编辑的，请大神帮忙看下哪有错误。
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭