wenzimu

新蟲 (初入文壇)

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

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

有一非線性方程為：
            -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ù)變化的前提下求出對應(yīng)S的值，并擬合出相應(yīng)的曲線；
         2.  P=50時，T在(0,30)連續(xù)變化的前提下求出對應(yīng)S的值，并擬合出相應(yīng)的曲線；
   請問這種情況在1stopt或matlab可以實現(xiàn)嗎，請高手幫忙，謝謝�。愃埔韵碌膱D形）

擬合圖

回復(fù)此樓

» 猜你喜歡

歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有27人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有291人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有12人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名已經(jīng)有0人回復(fù)
山西大同大學(xué)物理學(xué)碩士研究生招收調(diào)劑生已經(jīng)有0人回復(fù)
[調(diào)劑信息]211智能人工感知方向國家青年特聘專家課題組招收調(diào)劑研究生已經(jīng)有0人回復(fù)
［調(diào)劑信息］211智能人工感知方向國家青年特聘專家課題組招收調(diào)劑研究生已經(jīng)有0人回復(fù)
[調(diào)劑信息]211智能人工感知方向國家青年特聘專家課題組招收調(diào)劑研究生已經(jīng)有0人回復(fù)

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

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

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

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

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

★ ★ ★
wenzimu: 金幣+1 2012-05-31 14:39:33
dbb627: 金幣+2, 感謝應(yīng)助 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人)

回復(fù)此樓

showmethemoney

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

不是啊

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,:))

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

贊一下(1人)

回復(fù)此樓

showmethemoney

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

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

普通回帖

wenzimu

新蟲 (初入文壇)

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

謝謝 ~~~我看看哈！

贊一下

回復(fù)此樓

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

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

wenzimu

新蟲 (初入文壇)

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

謝謝我看看哈

回復(fù)此樓

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

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

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

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

贊一下

回復(fù)此樓

showmethemoney

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

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

wenzimu

新蟲 (初入文壇)

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

引用回帖:

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 ...

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

贊一下

回復(fù)此樓

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

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

wenzimu

新蟲 (初入文壇)

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

引用回帖:

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

贊一下

回復(fù)此樓

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

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

wenzimu

新蟲 (初入文壇)

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

引用回帖:

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的范圍您是如何確定的呢？

贊一下

回復(fù)此樓

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

引用回帖:

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

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

贊一下

回復(fù)此樓

showmethemoney

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

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

wenzimu

新蟲 (初入文壇)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

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

[考研] 321求調(diào)劑 +10	y-yh 2026-04-01	12/600	2026-04-02 18:15 by y-yh
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +15	GENJIOW 2026-03-31	16/800	2026-04-02 17:56 by cyh—315
[考研] 349求調(diào)劑 +8	zwjjjjjj 2026-03-31	8/400	2026-04-02 17:54 by 無懈可擊111
[論文投稿] chinese chemical letters英文版投稿求助 120+4	Yishengeryi 2026-03-30	6/300	2026-04-02 17:19 by Yishengeryi
[考研] 求調(diào)劑 +7	Aniyaio 2026-04-02	7/350	2026-04-02 16:42 by zzsw+
[考研] 0703化學(xué) +4	goldtt 2026-04-02	6/300	2026-04-02 16:05 by goldtt
[考研] 學(xué)碩化學(xué)工程與技術(shù)，一志愿中國海洋大學(xué)320+求調(diào)劑 +8	披星河 2026-04-02	8/400	2026-04-02 14:12 by oooqiao
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +9	雙馬尾痞老板2 2026-04-01	9/450	2026-04-02 12:14 by oooqiao
[考研] 0856材料與化工調(diào)劑，339 +13	10213207 2026-03-31	13/650	2026-04-02 11:35 by michael2011
[考研] 348環(huán)境工程調(diào)劑 +3	吳彥祖24k 2026-04-01	3/150	2026-04-02 09:14 by nanaliuyun
[考研] 一志愿346上海大學(xué)生物學(xué) +3	上海大學(xué)346調(diào)劑 2026-04-01	3/150	2026-04-02 08:36 by w蟲蟲123
[考研] 材料調(diào)劑 +11	一樣YWY 2026-03-31	11/550	2026-04-01 22:25 by zhouyuwinner
[考研] 315求調(diào)劑 +6	akie... 2026-03-28	7/350	2026-03-31 16:48 by asdfzly
[考研] 277跪求調(diào)劑 +8	1915668 2026-03-27	13/650	2026-03-31 14:58 by 王亮_大連醫(yī)科大
[考研] 286求調(diào)劑 +6	Faune 2026-03-30	6/300	2026-03-31 14:37 by jp9609
[考研] 一志愿食品科學(xué)與工程083200求調(diào)劑 +4	XQTJZ 2026-03-30	4/200	2026-03-31 04:10 by fmesaito
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 085602 化工專碩 338分求調(diào)劑 +12	路癡小琪 2026-03-27	12/600	2026-03-28 15:41 by L135790
[考研] 340求調(diào)劑 +5	jhx777 2026-03-27	5/250	2026-03-28 04:18 by fmesaito