24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2789  |  回复: 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)曲线。

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

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hytao2012

铁杆木虫 (正式写手)

木头虫子


jjdg: 金币+1, 好建议 2014-08-02 21:19:53
你对变量的命名和文字叙述的不一样,求导的时候都搞错了吧。为什么不定义一致了,那样多清晰。
5楼2014-07-31 11:11:16
已阅   回复此楼   关注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 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 电子信息279求调剂,有书读就行 +4 wwwooden 2026-04-08 7/350 2026-04-10 11:34 by wwwooden
[考研] 269求调剂 +10 啊啊我我 2026-04-07 10/500 2026-04-10 10:38 by 高维春
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 296求调剂 +5 汪!?! 2026-04-09 5/250 2026-04-09 17:47 by 柠檬不酸zy
[考研] 085501机械英二77总分294求调剂,接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
[考研] 材料与化工专硕306分找合适调剂 +27 沧海轻舟e 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[考研] 281求调剂 +10 椰子蘑菇 2026-04-06 10/500 2026-04-08 11:43 by zzucheup
[考研] 计算机408|在校多次国家级竞赛获奖|申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 304求调剂(085602,过四级,一志愿985) +25 化工人999 2026-04-04 26/1300 2026-04-07 22:06 by hemengdong
[考研] 277求调剂 数一104分 +9 瓶子PZ 2026-04-05 14/700 2026-04-07 17:52 by 蓝云思雨
[考研] 0854求调剂 +9 亨氏番茄沙司 2026-04-06 10/500 2026-04-07 14:37 by shdgaomin
[考研] 304求调剂 +4 luoye0105 2026-04-05 4/200 2026-04-06 21:05 by 木子君1218
[考研] 一志愿南航,数一英一学硕317求调剂!! +6 Acaciad 2026-04-04 6/300 2026-04-06 12:13 by 考研学校招点人
[考研] 296求调剂 +3 汪!?! 2026-04-05 5/250 2026-04-05 17:38 by 蓝云思雨
[考研] 295求调剂 +4 A你好研究生 2026-04-04 5/250 2026-04-04 22:46 by yu221
[考研] 316求调剂 +9 墨辰_Orion926 2026-04-04 9/450 2026-04-04 21:35 by lbsjt
信息提示
请填处理意见