24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 479  |  回复: 2

dl355538623

木虫 (正式写手)

扬帆远航

[求助] mathematica 解方程 已有1人参与

给位虫友好:
       我是非计算机专业。现在想解决方程组,找到合适的因子。苦于没有方法,希望给位能给予大力帮助。谢谢。

已知y1=130 x^3.61; y2=30 x^2.62; y3=80 x^3.19; 且符合方程 .
已知 x= 0.5, 1, 1.5, 2.0, 2.5, 3.0.

能否通过mathematica写出程序,给出A和B的值或则可能的范围?或则符合什么样的分布?谢谢各位了。
回复此楼

» 猜你喜欢

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

抓住时间、勿自寻烦恼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
拟合呗:
CODE:
Clear[y1, y2, y3]
x = {0.5, 1., 1.5, 2.0, 2.5, 3.0};

data = Transpose@{130 x^3.61, 30 x^2.62, 80 x^3.19, ConstantArray[0., Length@x]};

coef = FindFit[data, -y3^a + b y1^a + (1 - b) y2^a, {a, b}, {y1, y2, y3}]
(* {a -> 5.33194*10^-16, b -> 0.658877} *)

具体效果怎么样搞不清楚。(确切地说是懒得验。)可能会面临非线性拟合时的初值敏感问题。不过残差好像还是满小的:
CODE:
{y1, y2, y3} = Most@Transpose@data;
-y3^a + b y1^a + (1 - b) y2^a /. coef

(*
{1.11022*10^-16, -1.11022*10^-16, -1.11022*10^-16, 1.11022*10^-16, -1.11022*10^-16, 0.}
*)

小木虫Mathematica版块已毁(当然原本也不咋的),建议大家前往百度贴吧或Stackexchange。
2楼2015-03-19 20:02:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
dl355538623: 金币+10, ★★★很有帮助, 谢谢你的帮助。非常感谢, 2015-03-28 15:48:02
把精度提了一下,这个结果大概更可靠些。不过a已经小到了这份上直觉上应该无所谓了:
CODE:
Clear[y1, y2, y3]

x = {0.5, 1., 1.5, 2.0, 2.5, 3.0};

data = Transpose@{130 x^3.61, 30 x^2.62, 80 x^3.19, ConstantArray[0., Length@x]};

coef = FindFit[Rationalize[data, 0], -y3^a + b y1^a + (1 - b) y2^a, {a, b}, {y1, y2, y3}, WorkingPrecision -> 40]

(*
{a -> 1.920883336444913919341922550921764678017*10^-34,
b -> 0.6588766973338608221513612419841153771431}
*)

又画了一下你的6组数据所代表aOb平面内的线,感觉上面的结果应该是可信的:
CODE:
{y1, y2, y3} = Most@Transpose@data;
eqs = -y3^a + b y1^a + (1 - b) y2^a == 0 // Thread
ContourPlot[Evaluate@Rationalize[eqs, 0], {a, -10^-34, 3 10^-34}, {b, 0.5, 0.8}, WorkingPrecision -> 40, PlotPoints -> 40]

mathematica 解方程
小木虫Mathematica版块已毁(当然原本也不咋的),建议大家前往百度贴吧或Stackexchange。
3楼2015-03-19 20:19:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dl355538623 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本科211,293分请求调剂 +6 莲菜就是藕吧 2026-04-03 7/350 2026-04-04 11:12 by 莲菜就是藕吧
[考研] 材料专业383求调剂 +8 郭阳阳阳成 2026-04-03 8/400 2026-04-04 10:29 by Rednal.
[考研] 数一英一274机械调剂 +4 星陨流霞 2026-04-04 5/250 2026-04-04 09:35 by liucky
[考研] 化工调剂303分,过四级 +28 栖梧待风 2026-04-02 28/1400 2026-04-03 21:40 by qlm5820
[考研] 303求调剂 +10 DLkz1314. 2026-03-30 10/500 2026-04-03 18:03 by Jimmyandyou
[考研] 一志愿北京交通大学材料工程总分358 +4 cs0106 2026-04-03 4/200 2026-04-03 13:41 by 百灵童888
[考研] 324求调剂 +12 想上学求调 2026-04-02 12/600 2026-04-03 11:15 by linyelide
[考研] 325分化学调剂 +5 15771691647 2026-04-02 5/250 2026-04-03 09:58 by ChemPharm
[考研] 330分求调剂 +15 qzenlc 2026-03-29 15/750 2026-04-02 21:41 by 百灵童888
[考研] 土木304求调剂 +4 兔突突突, 2026-04-02 5/250 2026-04-02 21:16 by 兔突突突,
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
[考研] 275学硕081000服从调剂到其他专业,保不住本专业了 +7 一只小小水牛 2026-04-02 8/400 2026-04-02 14:23 by alice-2022
[考研] 学硕化学工程与技术,一志愿中国海洋大学320+求调剂 +8 披星河 2026-04-02 8/400 2026-04-02 14:12 by oooqiao
[考研] 0710生物学求调剂 +9 manman511 2026-04-01 9/450 2026-04-02 10:00 by zxl830724
[考研] 一志愿安徽大学计算机科学与技术学硕,331分求调剂 +5 蒋昌鹏qtj 2026-04-01 5/250 2026-04-02 08:10 by fxue1114
[考研] 求调剂,一志愿南京师范大学计算机专硕,初试373,六级通过, +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 材料与化工(0856)304求B区调剂 +8 邱gl 2026-03-30 16/800 2026-04-01 17:58 by 邱gl
[考研] 求调剂:一志愿:南京大学 专业:0705 总分320 ,本科985,四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
[考研] 材料工程专硕求调剂 +10 hyl3153942 2026-03-29 10/500 2026-03-31 16:31 by hypershenger
信息提示
请填处理意见