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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
141/2返回列表12下一页
查看: 2454  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wenzimu

新虫 (初入文坛)

[求助] 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)连续变化的前提下求出对应S的值，并拟合出相应的曲线；
            2.  P=50时，T在(0,30)连续变化的前提下求出对应S的值，并拟合出相应的曲线；
     请问这种情况在1stopt或matlab可以实现吗，请高手帮忙，谢谢！（类似以下的图形）

拟合图
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-30 22:20:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★
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发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

csgt0

荣誉版主 (著名写手)

彩色挂图
不是啊
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发消息 送TA红花 TA的回帖
普通回帖

wenzimu

新虫 (初入文坛)
谢谢  ~~~我看看哈！
赞一下 回复此楼
2楼2012-05-31 08:30:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenzimu

新虫 (初入文坛)
谢谢   我看看哈
回复此楼
3楼2012-05-31 08:31:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

感谢参与，应助指数 +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发消息 送TA红花 TA的回帖

wenzimu

新虫 (初入文坛)
引用回帖:
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发消息 送TA红花 TA的回帖

wenzimu

新虫 (初入文坛)
引用回帖:
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 ...

可以请教一下，您是通过什么软件计算的吗？还有，可以在P=20时，T在[0,30]的前提下计算出所有的S值吗？谢谢！
赞一下 回复此楼
6楼2012-05-31 12:18:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenzimu

新虫 (初入文坛)
引用回帖:
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发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图
引用回帖:
8楼: Originally posted by wenzimu at 2012-05-31 19:12:44
真的不好意思，还要麻烦您一下~由于第一次接触matlab，所以很多看不大懂，比如代码中的y1和y2的范围您是如何确定的呢？...

你这是个三次方程，有三个解，y1取得是S1的第一行，也就是第一个解，而S1是通过将T代入g1得到的，g1是通过解方程得到的S关于T的函数。
赞一下 回复此楼
showmethemoney
9楼2012-06-01 09:22:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenzimu

新虫 (初入文坛)
引用回帖:
9楼: Originally posted by csgt0 at 2012-06-01 09:22:31
你这是个三次方程，有三个解，y1取得是S1的第一行，也就是第一个解，而S1是通过将T代入g1得到的，g1是通过解方程得到的S关于T的函数。...

谢谢您的回复，那就是说y1和y2的语句写成：
y1=S1(1,subs(g1));
y2=S2(1,sub(g2));   
是这样吗？麻烦您再看一下，谢谢！
赞一下 回复此楼
10楼2012-06-01 10:49:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wenzimu 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 +3 鹿沫笙 2026-03-02 3/150 2026-03-02 18:10 by barlinike
[考研] 中国科学技术大学材料与化工281求调剂，有科研和获奖经历 +3 wsxw 2026-03-02 3/150 2026-03-02 18:08 by lature00
[考研] 289求调剂 +6 yang婷 2026-03-02 6/300 2026-03-02 18:05 by 无际的草原
[考研] 299求调剂 +4 kkcoco25 2026-03-02 4/200 2026-03-02 18:04 by barlinike
[考研] 268求调剂 +5 简单点0 2026-03-02 6/300 2026-03-02 18:03 by lature00
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 14/700 2026-03-02 17:44 by 好好好1233
[考研] 接收调剂 +5 津萌津萌 2026-03-02 12/600 2026-03-02 17:04 by 津萌津萌
[考研] 材料284求调剂，一志愿郑州大学英一数二专硕 +12 想上岸的土拨鼠 2026-02-28 12/600 2026-03-02 16:18 by youmomaoyan
[考研] 化工专硕348，一志愿985求调剂 +6 弗格个 2026-02-28 9/450 2026-03-02 14:09 by liyongv
[考博] 诚招农业博士 +3 心欣向荣 2026-02-28 3/150 2026-03-02 13:33 by 时间不狗
[考研] 0805总分292，求调剂 +8 幻想之殇 2026-03-01 8/400 2026-03-02 12:51 by 无际的草原
[考研] 272求调剂 +7 材紫有化 2026-02-28 7/350 2026-03-02 12:48 by 无际的草原
[考研] 281求调剂 +5 2026计算机_诚心 2026-03-01 8/400 2026-03-02 11:05 by 汪！？！
[考研] 284求调剂 +10 天下熯 2026-02-28 11/550 2026-03-02 11:03 by 无际的草原
[考研] 材料类求调剂 +11 wana_kiko 2026-02-28 14/700 2026-03-02 08:46 by 聪明的大松鼠
[考研] 材料复试调剂 +4 学材料的点 2026-03-01 5/250 2026-03-02 08:26 by houyaoxu
[基金申请] 成果系统访问量大，请一小时后再尝试。---NSFC啥时候好哦，已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考研] 307求调剂 +5 wyyyqx 2026-03-01 5/250 2026-03-01 15:21 by Fff-1
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
信息提示
关闭
请填处理意见
关闭