King_A

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 711
散金: 201
帖子: 214
在線: 97.7小時(shí)
蟲號(hào): 1130184
注冊(cè): 2010-10-23
專業(yè): 金屬材料的合金相、相變及

[求助] 求matlab編程代碼，求幫助，謝謝！已有3人參與

有上述公式22，其中除f,D和Q外，其它均為已知量，其中T（t）為一維數(shù)組，此外還已知另一數(shù)組Y(t)。
現(xiàn)想進(jìn)行一下計(jì)算：定義K=∑（f-Y(t)）^2，分別對(duì)K求D和Q的一階偏導(dǎo)數(shù)，并另其一階偏導(dǎo)數(shù)等于0，解方程組，最終得到D和Q 的值，
最后將計(jì)算所得的D和Q值返回到公式22中，畫出(T,f)曲線。

自行研究多日，卻不得其所，實(shí)在無(wú)奈，只得到此求助大神幫幫，在此先表謝意。

22.png

回復(fù)此樓

» 猜你喜歡

材料，紡織，生物（0856、0710），化學(xué)招生啦已經(jīng)有8人回復(fù)
0703化學(xué)調(diào)劑，六級(jí)已過(guò)，有科研經(jīng)歷已經(jīng)有9人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有6人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有5人回復(fù)
生物學(xué)071000 329分求調(diào)劑已經(jīng)有3人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有12人回復(fù)
288求調(diào)劑，一志愿華南理工大學(xué)071005 已經(jīng)有3人回復(fù)
307求調(diào)劑已經(jīng)有3人回復(fù)
312求調(diào)劑已經(jīng)有8人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

matlab 編程小問(wèn)題，求幫助，謝謝。已經(jīng)有13人回復(fù)
求高手給做一道題，用MATLAB或1stOpt，請(qǐng)把結(jié)果直接給出來(lái)，我的軟件就不行，謝謝已經(jīng)有5人回復(fù)
matlab動(dòng)力學(xué)模型編程求助已經(jīng)有6人回復(fù)
求一個(gè)matlab詳細(xì)代碼已經(jīng)有5人回復(fù)
求回歸方程參數(shù)的最優(yōu)化算法matlab代碼已經(jīng)有17人回復(fù)
matlab編程齒輪已經(jīng)有9人回復(fù)
請(qǐng)高人指點(diǎn)一個(gè)matlab編程問(wèn)題，謝謝已經(jīng)有5人回復(fù)
新手matlab求助，謝謝各位已經(jīng)有10人回復(fù)
求助用matlab或者（1st0pt)編程解方程組已經(jīng)有3人回復(fù)
MATLAB高手幫我看看這個(gè)數(shù)據(jù)存的格式，謝謝已經(jīng)有3人回復(fù)
求基于DCT域的數(shù)字圖像水印Matlab代碼，謝謝各位大神！已經(jīng)有3人回復(fù)
請(qǐng)matlab高手看一下這個(gè)問(wèn)題怎樣編程才好，附程序代碼，謝謝！已經(jīng)有9人回復(fù)
matlab編程大家看這兩個(gè)式子該如何編寫呢？已經(jīng)有3人回復(fù)
請(qǐng)問(wèn)如何用Matlab編程做這個(gè)圖？謝謝大家！已經(jīng)有5人回復(fù)
如何用matlab軟件或者用C及C++自己手工編程尋找轉(zhuǎn)折點(diǎn) 已經(jīng)有12人回復(fù)
求助matlab編程已經(jīng)有13人回復(fù)
【求助】用matlab編程，需要將數(shù)據(jù)自檢分類（高手進(jìn)來(lái)指點(diǎn)下）已經(jīng)有5人回復(fù)
【求助】請(qǐng)教一個(gè)matlab程序，謝謝幫忙！！已經(jīng)有15人回復(fù)
【資源】Matlab編程基礎(chǔ)及應(yīng)用(視頻教程) 已經(jīng)有50人回復(fù)

1樓 2014-07-30 08:39:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

hytao2012

鐵桿木蟲 (正式寫手)

木頭蟲子

應(yīng)助: 53 (初中生)
金幣: 6341.8
散金: 115
紅花: 11
帖子: 479
在線: 206.6小時(shí)
蟲號(hào): 2050091
注冊(cè): 2012-10-08
性別: GG
專業(yè): 固體力學(xué)

【答案】應(yīng)助回帖

★
jjdg: 金幣+1, 辛苦了 2014-08-04 15:04:21

我也寫了一個(gè)小程序，感覺(jué)好像是因?yàn)槟愕姆匠烫珡?fù)雜，解不出來(lái)啊，注釋的是我寫了個(gè)簡(jiǎn)單的方程，就可以了。

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人)

回復(fù)此樓

9樓2014-08-04 09:04:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

jerkwin

專家顧問(wèn) (正式寫手)

專家經(jīng)驗(yàn): +14
應(yīng)助: 454 (碩士)
金幣: 20699.1
散金: 148
紅花: 81
帖子: 813
在線: 2648.3小時(shí)
蟲號(hào): 1023452
注冊(cè): 2010-05-19
專業(yè): 理論和計(jì)算化學(xué)
管轄: 分子模擬

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

做最小二乘擬合吧

贊一下

回復(fù)此樓

2樓2014-07-30 22:27:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

King_A

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 711
散金: 201
帖子: 214
在線: 97.7小時(shí)
蟲號(hào): 1130184
注冊(cè): 2010-10-23
專業(yè): 金屬材料的合金相、相變及

引用回帖:

2樓: Originally posted by jerkwin at 2014-07-30 22:27:35
做最小二乘擬合吧

最小二乘擬合達(dá)不到我的要求啊，我還需要具體的D和Q的值。最小二乘擬合只能得到方程，得不到具體D和Q

贊一下

回復(fù)此樓

3樓2014-07-31 08:49:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

King_A

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 711
散金: 201
帖子: 214
在線: 97.7小時(shí)
蟲號(hào): 1130184
注冊(cè): 2010-10-23
專業(yè): 金屬材料的合金相、相變及

下面是我自行編輯的，請(qǐng)大神幫忙看下哪有錯(cuò)誤。
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')

下面是結(jié)果：
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))

贊一下

回復(fù)此樓

4樓2014-07-31 10:16:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hytao2012

鐵桿木蟲 (正式寫手)

木頭蟲子

應(yīng)助: 53 (初中生)
金幣: 6341.8
散金: 115
紅花: 11
帖子: 479
在線: 206.6小時(shí)
蟲號(hào): 2050091
注冊(cè): 2012-10-08
性別: GG
專業(yè): 固體力學(xué)

★
jjdg: 金幣+1, 好建議 2014-08-02 21:19:53

你對(duì)變量的命名和文字?jǐn)⑹龅牟灰粯樱髮?dǎo)的時(shí)候都搞錯(cuò)了吧。為什么不定義一致了，那樣多清晰。

贊一下

回復(fù)此樓

5樓2014-07-31 11:11:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

King_A

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 711
散金: 201
帖子: 214
在線: 97.7小時(shí)
蟲號(hào): 1130184
注冊(cè): 2010-10-23
專業(yè): 金屬材料的合金相、相變及

引用回帖:

5樓: Originally posted by hytao2012 at 2014-07-31 11:11:16
你對(duì)變量的命名和文字?jǐn)⑹龅牟灰粯�，求�?dǎo)的時(shí)候都搞錯(cuò)了吧。為什么不定義一致了，那樣多清晰。

只是臨時(shí)代替一下，因?yàn)檫€在嘗試

贊一下

回復(fù)此樓

6樓2014-08-02 19:30:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hc315

鐵蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 30.5
散金: 50
紅花: 1
帖子: 116
在線: 60.6小時(shí)
蟲號(hào): 1341899
注冊(cè): 2011-07-10
性別: GG
專業(yè): 運(yùn)籌與管理

引用回帖:

3樓: Originally posted by King_A at 2014-07-31 08:49:24
最小二乘擬合達(dá)不到我的要求啊，我還需要具體的D和Q的值。最小二乘擬合只能得到方程，得不到具體D和Q...

對(duì)D和Q做最小二乘擬合，就可以得到D和Q的值啊

贊一下

回復(fù)此樓

7樓2014-08-03 13:20:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

King_A

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 711
散金: 201
帖子: 214
在線: 97.7小時(shí)
蟲號(hào): 1130184
注冊(cè): 2010-10-23
專業(yè): 金屬材料的合金相、相變及

引用回帖:

7樓: Originally posted by hc315 at 2014-08-03 13:20:53
對(duì)D和Q做最小二乘擬合，就可以得到D和Q的值啊...

計(jì)算要比擬合更準(zhǔn)確點(diǎn)吧

贊一下

回復(fù)此樓

8樓2014-08-04 08:19:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hc315

鐵蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 30.5
散金: 50
紅花: 1
帖子: 116
在線: 60.6小時(shí)
蟲號(hào): 1341899
注冊(cè): 2011-07-10
性別: GG
專業(yè): 運(yùn)籌與管理

★
jjdg: 金幣+1, 辛苦了 2014-08-04 15:04:32

引用回帖:

8樓: Originally posted by King_A at 2014-08-04 08:19:07
計(jì)算要比擬合更準(zhǔn)確點(diǎn)吧...

根據(jù)最小二乘的計(jì)算原理，如果你的式子嚴(yán)格滿足輸入變量和輸出函數(shù)值的映射關(guān)系，那么擬合出來(lái)的系數(shù)值應(yīng)該等同于你計(jì)算出來(lái)的值

贊一下

回復(fù)此樓

10樓2014-08-04 14:46:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 312求調(diào)劑 +7	陌宸希 2026-03-16	8/400	2026-03-18 09:54 by macy2011
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 293求調(diào)劑 +6	世界首富 2026-03-11	6/300	2026-03-17 17:04 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 【0856】化學(xué)工程（085602）313 分，本科學(xué)科評(píng)估A類院校化學(xué)工程與工藝，誠(chéng)求調(diào)劑 +7	小劉快快上岸 2026-03-11	8/400	2026-03-17 16:57 by ruiyingmiao
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475 +5	Becho359 2026-03-11	5/250	2026-03-14 11:35 by 哦哦123
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 310求調(diào)劑 +3	【上上簽】 2026-03-11	3/150	2026-03-13 16:16 by JourneyLucky
[考研] 314求調(diào)劑 +7	無(wú)懈可擊的巨人 2026-03-12	7/350	2026-03-13 15:40 by JourneyLucky
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655