24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3468  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 青沙山 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 280求调剂 +7 兮兮夜夜 2026-04-09 10/500 2026-04-12 00:33 by 蓝云思雨
[考研] 290求调剂 +13 luoziheng 2026-04-10 13/650 2026-04-11 23:15 by labixiaoqiao
[硕博家园] 有没有学校材料专业收跨调(一志愿085410) +6 momo(上岸版) 2026-04-06 9/450 2026-04-11 22:38 by wj165256
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 电子信息279求调剂,有书读就行 +8 wwwooden 2026-04-08 11/550 2026-04-11 20:22 by cq2548
[考研] 求调剂 +6 电气300求调剂不 2026-04-08 6/300 2026-04-11 20:14 by 逆水乘风
[考研] 复试调剂 +9 积极向上; 2026-04-10 11/550 2026-04-11 09:25 by 猪会飞
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
[考研] 071000生物学调剂求助 +17 zzzzwww 2026-04-09 20/1000 2026-04-10 15:55 by 求调剂zz
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +20 慕绝cc 2026-04-09 24/1200 2026-04-10 10:22 by xujun0624
[考研] 08工学 309分求调剂 +6 Yin DY 2026-04-08 6/300 2026-04-10 09:18 by Delta2012
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 求调剂 +3 猪肉墩粉条cc 2026-04-08 4/200 2026-04-09 10:05 by 猪肉墩粉条cc
[考研] 085404,334分,求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
[考研] 304求调剂 +10 素年祭语 2026-04-06 17/850 2026-04-08 09:05 by 蓝云思雨
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 专硕085403,291分,有两篇专利,一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
信息提示
请填处理意见