wenzimu

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學

[求助] 1stopt或matlab如何求解以下的非線性方程并擬合出相應曲線？

有一非線性方程為：
            -5022.5+1.3535*S-8.985*10E(-5)*S^2+0.1257*T-5.574*10E(-5)*S*T+5.106*S^2*T+4P/S+0.0028P-23.079*P*S=0

   想令1.  P=20時，T在(0,30)連續(xù)變化的前提下求出對應S的值，并擬合出相應的曲線；
         2.  P=50時，T在(0,30)連續(xù)變化的前提下求出對應S的值，并擬合出相應的曲線；
   請問這種情況在1stopt或matlab可以實現(xiàn)嗎，請高手幫忙，謝謝�。愃埔韵碌膱D形）

擬合圖

回復此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有166人回復
哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有5人回復
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復
求調(diào)劑已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復
光學工程學碩調(diào)劑信息已經(jīng)有25人回復
歡迎加入課題組已經(jīng)有0人回復
散金幣，求好運，祝面上順利！已經(jīng)有29人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較已經(jīng)有4人回復
擬合非線性方程系數(shù)的精度問題已經(jīng)有5人回復
如何使用matlab求解非線性方程組的所有整數(shù)解？已經(jīng)有9人回復
matlab求解非線性方程組已經(jīng)有16人回復
非線性方程組matlab擬合已經(jīng)有9人回復
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！已經(jīng)有18人回復
用1stopt進行曲線擬合每次給出結果不一樣已經(jīng)有6人回復
【原創(chuàng)】多項式和非多項式曲線擬合的MATLAB實現(xiàn)視頻教程已經(jīng)有10人回復
【求助】如何使得matlab擬合的曲線強制經(jīng)過一個點已經(jīng)有8人回復
【求助】求教matlab解非線性方程組已經(jīng)有9人回復

1樓 2012-05-30 22:20:08

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

回帖置頂 ( 共有1個 )

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

★ ★ ★
wenzimu: 金幣+1 2012-05-31 14:39:33
dbb627: 金幣+2, 感謝應助 2012-05-31 22:51:05
wenzimu: 回帖置頂 2012-06-01 12:30:24

clc
clear
syms P T S
f=-5022.5+1.3535*S-8.985*10^(-5)*S.^2+0.1257*T-5.574*10^(-5)*S*T+5.106*S.^2*T+4*P./S+0.0028*P-23.079*P.*S;
P=20;
g1=solve(subs(f),S)
P=50;
g2=solve(subs(f),S);
T=1:.1:30;
S1=subs(g1);
S2=subs(g2);
y1=S1(1,

;
y2=S2(1,

;
plot(T,y1,T,y2)

有三個解，好像第一個比較符合你的

贊一下(1人)

回復此樓

showmethemoney

7樓2012-05-31 13:37:49

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

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

不是啊

CODE:

y1=S1(1,:);

y2=S2(1,:);

plot(T,y1,T,y2)

我只不過把它突出一下，用y1,y2表示出來，不用也可以，就是取S1和S2中的3行的第1行

CODE:

plot(T,S1(1,:),T,S2(1,:))

當然你要的可能不是第一行，也可能是第2、3行，把1改成2或者3就行了

贊一下(1人)

回復此樓

showmethemoney

11樓2012-06-01 10:59:12

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

普通回帖

wenzimu

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學

謝謝 ~~~我看看哈！

贊一下

回復此樓

2樓2012-05-31 08:30:26

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

wenzimu

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學

謝謝我看看哈

回復此樓

3樓2012-05-31 08:31:36

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

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

我怎么算出來S沒那么大值，比如P=20，T=2時，有
S=
54.147511926504082426692461168811
0.015906120687527972318210073677971 - 5.1427877848475078473632223503473e-39*i
- 9.0957857707756256775368611954232 + 5.1427877848475078473632223503473e-39*i

贊一下

回復此樓

showmethemoney

4樓2012-05-31 09:48:44

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

wenzimu

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學

引用回帖:

4樓: Originally posted by csgt0 at 2012-05-31 09:48:44
我怎么算出來S沒那么大值，比如P=20，T=2時，有
S=
54.147511926504082426692461168811
0.015906120687527972318210073677971 - 5.1427877848475078473632223503473e-39*i
- 9.09578577077562567753686 ...

上面的圖只是相似的，并不是該方程的值，我想通過這個方程得到跟那個圖一樣的結果，而不是上面圖中的值是這個方程的值，給您造成誤解，真不好意思！

贊一下

回復此樓

5樓2012-05-31 12:15:43

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

wenzimu

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學

引用回帖:

可以請教一下，您是通過什么軟件計算的嗎？還有，可以在P=20時，T在[0,30]的前提下計算出所有的S值嗎？謝謝！

贊一下

回復此樓

6樓2012-05-31 12:18:33

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

wenzimu

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學

引用回帖:

7樓: Originally posted by csgt0 at 2012-05-31 13:37:49
clc
clear
syms P T S
f=-5022.5+1.3535*S-8.985*10^(-5)*S.^2+0.1257*T-5.574*10^(-5)*S*T+5.106*S.^2*T+4*P./S+0.0028*P-23.079*P.*S;
P=20;
g1=solve(subs(f),S)
P=50;
g2=solve(subs(f),S);
T=1:.1:30; ...

真的不好意思，還要麻煩您一下~由于第一次接觸matlab，所以很多看不大懂，比如代碼中的y1和y2的范圍您是如何確定的呢？

贊一下

回復此樓

8樓2012-05-31 19:12:44

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

引用回帖:

8樓: Originally posted by wenzimu at 2012-05-31 19:12:44
真的不好意思，還要麻煩您一下~由于第一次接觸matlab，所以很多看不大懂，比如代碼中的y1和y2的范圍您是如何確定的呢？...

你這是個三次方程，有三個解，y1取得是S1的第一行，也就是第一個解，而S1是通過將T代入g1得到的，g1是通過解方程得到的S關于T的函數(shù)。

贊一下

回復此樓

showmethemoney

9樓2012-06-01 09:22:31

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

wenzimu

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學

引用回帖:

9樓: Originally posted by csgt0 at 2012-06-01 09:22:31
你這是個三次方程，有三個解，y1取得是S1的第一行，也就是第一個解，而S1是通過將T代入g1得到的，g1是通過解方程得到的S關于T的函數(shù)。...

謝謝您的回復，那就是說y1和y2的語句寫成：
y1=S1(1,subs(g1));
y2=S2(1,sub(g2));
是這樣嗎？麻煩您再看一下，謝謝！

贊一下

回復此樓

10樓2012-06-01 10:49:23

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

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

[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 268求調(diào)劑 +8	一定有學上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 304求調(diào)劑 +8	小熊joy 2026-03-14	8/400	2026-03-17 17:29 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 085600材料與化工求調(diào)劑 +5	緒幸與子 2026-03-17	5/250	2026-03-17 16:40 by laoshidan
[考研] 271求調(diào)劑 +12	生如夏花… 2026-03-11	14/700	2026-03-17 10:56 by lovewei0727
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 機械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 308 085701 四六級已過求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7	mysdl 2026-03-11	9/450	2026-03-13 22:43 by JourneyLucky
[考研] 四川大學085601材料工程專碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 工科，求調(diào)劑 +3	我887 2026-03-11	3/150	2026-03-13 21:39 by JourneyLucky
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 求調(diào)劑 +5	一定有學上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 314求調(diào)劑 +7	無懈可擊的巨人 2026-03-12	7/350	2026-03-13 15:40 by JourneyLucky
[考研] 工科調(diào)劑 +4	Jiang191123！ 2026-03-11	4/200	2026-03-13 15:15 by Miko19
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321