24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3491  |  回复: 7

青沙山

木虫 (小有名气)

[求助] Mathematica中含参数的方程求解时如何自动舍去虚数解或不符合条件的解 已有1人参与

本人用的是9.0版,求解一个含参数的方程时若添加解的条件时,总是解不出来,求高手指点?

我想添加个解的不等式条件或者解必须是实数的条件,见下图:
Mathematica中含参数的方程求解时如何自动舍去虚数解或不符合条件的解
s.png

代码如下:
Solve[{-6 lx^2 (x - Subscript[lx, 1])^2 +
    12 lx (x - Subscript[lx, 1])^2 Subscript[lx, 1] -
    4 x^3 Subscript[lx, 2] +
    12 lx (x - Subscript[lx, 1])^2 Subscript[lx, 2] + 4
\!\(\*SubsuperscriptBox[\(lx\), \(1\), \(3\)]\) Subscript[lx, 2] +
    2 lx^2
\!\(\*SubsuperscriptBox[\(lx\), \(2\), \(2\)]\) -
    4 lx Subscript[lx, 1]
\!\(\*SubsuperscriptBox[\(lx\), \(2\), \(2\)]\) + 6
\!\(\*SubsuperscriptBox[\(lx\), \(1\), \(2\)]\)
\!\(\*SubsuperscriptBox[\(lx\), \(2\), \(2\)]\) - 4 lx
\!\(\*SubsuperscriptBox[\(lx\), \(2\), \(3\)]\) + 4 Subscript[lx, 1]
\!\(\*SubsuperscriptBox[\(lx\), \(2\), \(3\)]\) +
\!\(\*SubsuperscriptBox[\(lx\), \(2\), \(4\)]\) == 0,
  0 < Subscript[lx, 1] < x < Subscript[lx, 1] + Subscript[lx, 2] <
   lx}, x, Reals]
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
青沙山: 金币+5, ★★★很有帮助 2016-01-08 19:55:03
……下标害的。对现有代码而言最简单的改法是:
CODE:
<< Notation`

Symbolize[ParsedBoxWrapper[\(lx \_ _\)]](*这行用相应的面板输入会比较好看,具体看看帮助*)

Solve[{-6*lx^2*(x - Subscript[lx, 1])^2 +
    12*lx*(x - Subscript[lx, 1])^2*Subscript[lx, 1] - 4*x^3*Subscript[lx, 2] +
         12*lx*(x - Subscript[lx, 1])^2*Subscript[lx, 2] +
    4*Subscript[lx, 1]^3*Subscript[lx, 2] + 2*lx^2*Subscript[lx, 2]^2 -
         4*lx*Subscript[lx, 1]*Subscript[lx, 2]^2 +
    6*Subscript[lx, 1]^2*Subscript[lx, 2]^2 - 4*lx*Subscript[lx, 2]^3 +
         4*Subscript[lx, 1]*Subscript[lx, 2]^3 + Subscript[lx, 2]^4 == 0,
  0 < Subscript[lx, 1] < x < Subscript[lx, 1] + Subscript[lx, 2] < lx}, x, Reals]

至于为什么下标不好,解释起来字太多这里就用中文不重复了,有空看看这:http://mathematica.stackexchange.com/a/18395/1871
楼主实在理解不了就记住:没事别用下标。

» 本帖已获得的红花(最新10朵)

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

青沙山

木虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by xzczd at 2016-01-08 00:25:38
……下标害的。对现有代码而言最简单的改法是:

<< Notation`

Symbolize(*这行用相应的面板输入会比较好看,具体看看帮助*)

Solve


至于为什么下标不好,解释起来字太多这里就用中文不重复了,有 ...

谢谢,小白表示仍然不懂呢!

如下图:
Mathematica中含参数的方程求解时如何自动舍去虚数解或不符合条件的解-1

第一行看不懂,能给个图吗?

若不用下标,请问还有什么方式表示一个参数序列呢?
尤其是我想对这个参数序列求和,如下例中定义的函数:
Mathematica中含参数的方程求解时如何自动舍去虚数解或不符合条件的解-2

fx[x_, i_] := -qx*f[x] + \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(u = 1\), \(i\)]\((
\*SubscriptBox[\(rx\), \(u\)]*
\*SuperscriptBox[\(f[x -
\*UnderoverscriptBox[\(\[Sum]\), \(v = 1\), \(u\)]
\*SubscriptBox[\(lx\), \(v\)]]\), \(0\)])\)\)
3楼2016-01-08 10:42:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)

引用回帖:
3楼: Originally posted by 青沙山 at 2016-01-08 10:42:56
谢谢,小白表示仍然不懂呢!

如下图:


第一行看不懂,能给个图吗?

若不用下标,请问还有什么方式表示一个参数序列呢?
尤其是我想对这个参数序列求和,如下例中定义的函数:


fx := -q ...

……好吧我忘了这茬了。总之你执行我给的代码前,先按Ctrl+Shift+N把代码转成StandardForm。

至于后一个,直接用中括号啊:
CODE:
fx[x_, i_] := (-qx)*f[x] + Sum[ru[u]*f[x - Sum[lv[v], {v, 1, u}]]^0,
       {u, 1, i}]

» 本帖已获得的红花(最新10朵)

小木虫Mathematica版块已毁(当然原本也不咋的),建议大家前往百度贴吧或Stackexchange。
4楼2016-01-08 12:27:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青沙山

木虫 (小有名气)

送红花一朵
引用回帖:
4楼: Originally posted by xzczd at 2016-01-08 12:27:06
……好吧我忘了这茬了。总之你执行我给的代码前,先按Ctrl+Shift+N把代码转成StandardForm。

至于后一个,直接用中括号啊:

fx := (-qx)*f + Sum[ru*f[x - Sum]^0,
       {u, 1, i}]...

谢谢你!

Solve解出来有3个解,Reduce解出1个解,解出来的表达式好长!
好像Solve出来的结果仍然没有排除不符合条件的解……
学了3周Mathematica,依然搞不定我要的结果……
总之,非常感谢!
5楼2016-01-08 20:05:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)

引用回帖:
5楼: Originally posted by 青沙山 at 2016-01-08 20:05:53
谢谢你!

Solve解出来有3个解,Reduce解出1个解,解出来的表达式好长!
好像Solve出来的结果仍然没有排除不符合条件的解……
学了3周Mathematica,依然搞不定我要的结果……
总之,非常感谢!...

我觉得你这个方程弄解析解意义很有限,你仔细看看Solve给出的解,它是以高阶多项式的根的形式给出的,并且阶数已经大于4了,这基本说明它根本没有通常意义上的解析解。

» 本帖已获得的红花(最新10朵)

小木虫Mathematica版块已毁(当然原本也不咋的),建议大家前往百度贴吧或Stackexchange。
6楼2016-01-08 21:42:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青沙山

木虫 (小有名气)

送红花一朵
引用回帖:
6楼: Originally posted by xzczd at 2016-01-08 21:42:50
我觉得你这个方程弄解析解意义很有限,你仔细看看Solve给出的解,它是以高阶多项式的根的形式给出的,并且阶数已经大于4了,这基本说明它根本没有通常意义上的解析解。...

最近一直在弄这个,果然还是算不出来
这个带参数lx1,lx2的方程解出来x的值还要代入别的表达式中,然后通过相关方程组求参数lx1,lx2的近似解,应该如何处理?不断循环赋值求近似解吗?
7楼2016-01-08 23:07:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)

引用回帖:
7楼: Originally posted by 青沙山 at 2016-01-08 23:07:09
最近一直在弄这个,果然还是算不出来
这个带参数lx1,lx2的方程解出来x的值还要代入别的表达式中,然后通过相关方程组求参数lx1,lx2的近似解,应该如何处理?不断循环赋值求近似解吗?...

……你这么抽象地说我也很难给什么具体建议。运气好的话直接把上面的解代入你说的表达式就行了。(姑且问一句你应该不会不知道怎么用ReplaceAll(/.)吧……)Mathematica是有求方程数值解的函数的。你看看NSolve还有FindFit的帮助。
小木虫Mathematica版块已毁(当然原本也不咋的),建议大家前往百度贴吧或Stackexchange。
8楼2016-01-08 23:49:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 青沙山 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考博] 申博/考博 +3 啃面包的小书虫 2026-04-17 4/200 2026-04-17 23:54 by 阳阳阳^_^
[考研] 304求调剂 +7 castLight 2026-04-16 7/350 2026-04-17 20:05 by 关一盏灯cd
[考研] 294求调剂 +7 淡然654321 2026-04-17 8/400 2026-04-17 16:36 by wutongshun
[考研] 0854求调剂 +21 门路摸摸 2026-04-15 25/1250 2026-04-17 15:45 by qzxyhcsy
[考研] 279求调剂 +13 张番茄不炒蛋 2026-04-11 13/650 2026-04-17 10:38 by cuisz
[考研] 294求调剂 +14 淡然654321 2026-04-15 14/700 2026-04-16 21:01 by lpl364211
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 085801电气专硕272求调剂 +19 电气李 2026-04-13 21/1050 2026-04-15 13:37 by 黑科技矿业
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
信息提示
请填处理意见