24小时热门版块排行榜    

查看: 2678  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cupidaxun

木虫 (著名写手)

[求助] 求解一元三次方程

已知(a-x)(b-x)(c-x)=A,求x值,请教用什么方法或工具能快速解出x?非常感谢!
回复此楼

» 猜你喜欢

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

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

cupidaxun

木虫 (著名写手)

引用回帖:
2楼: Originally posted by 月只蓝 at 2013-06-10 10:07:13
1stopt、MATLAB、Maple,EXCEL都可以做。

1stopt软件小巧,程序也非常简单:
代码:
Parameters x;
Function (1-x)*(2-x)*(3-x)=4;

迭代数: 16
计算用时(时:分:秒:毫秒): 00:00:00:61
计算中止原因: 达 ...

非常感谢提供的工具和方法,只是这些工具都没有用过,stopt正在研究中,请问用Excel计算能否实现?(比如,我列表给出a,b,c,A数据,对应的x,就可以在另一列计算出来)
3楼2013-06-10 12:33:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
cupidaxun: 金币+5, ★★★★★最佳答案, 非常感谢! 2013-06-10 12:30:17
1stopt、MATLAB、Maple,EXCEL都可以做。

1stopt软件小巧,程序也非常简单:
代码:
Parameters x;
Function (1-x)*(2-x)*(3-x)=4;

迭代数: 16
计算用时(时:分:秒:毫秒): 00:00:00:61
计算中止原因: 达到收敛判定标准
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
函数表达式 (1-x)*(2-x)*(3-x)-(4)
目标函数值: 0
x: 0.203678096741042

====== 计算结束 ======果:



对于a,b,c,A是已知参数的情况下,MATLAB程序如下,附图1以a=1,b=2,c=3,A=4为例子的输入和输出结果
function solvenole
clear all;clc
format long
global a b c A

a=input('请输入已知常数a=');
b=input('请输入已知常数b=');
c=input('请输入已知常数c=');
A=input('请输入已知常数A=');


x1=fsolve(@fun299792,[a])


%--------------------------------
function y=fun299792(x)
global a b c A
y=(a-x)*(b-x)*(c-x)-A;

Maple也可以做,见附图2
求解一元三次方程
附图1.jpg


求解一元三次方程-1
附图2.jpg

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2013-06-10 10:07:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

★ ★
csgt0: 金币+2, 谢谢 2013-06-14 10:42:06
引用回帖:
3楼: Originally posted by cupidaxun at 2013-06-10 12:33:15
非常感谢提供的工具和方法,只是这些工具都没有用过,stopt正在研究中,请问用Excel计算能否实现?(比如,我列表给出a,b,c,A数据,对应的x,就可以在另一列计算出来)...

可以,列个公式就可以
x一实根,2虚根
x=solve('(a-x)*(b-x)*(c-x)=A','x')

x =

                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           a/3 + b/3 + c/3 - ((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)/((((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)^3 + (A/2 + ((a + b + c)*(a*b + a*c + b*c))/6 - (a + b + c)^3/27 - (a*b*c)/2)^2)^(1/2) - A/2 - ((a + b + c)*(a*b + a*c + b*c))/6 + (a + b + c)^3/27 + (a*b*c)/2)^(1/3) + ((((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)^3 + (A/2 + ((a + b + c)*(a*b + a*c + b*c))/6 - (a + b + c)^3/27 - (a*b*c)/2)^2)^(1/2) - A/2 - ((a + b + c)*(a*b + a*c + b*c))/6 + (a + b + c)^3/27 + (a*b*c)/2)^(1/3)

a/3 + b/3 + c/3 + ((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)/(2*((((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)^3 + (A/2 + ((a + b + c)*(a*b + a*c + b*c))/6 - (a + b + c)^3/27 - (a*b*c)/2)^2)^(1/2) - A/2 - ((a + b + c)*(a*b + a*c + b*c))/6 + (a + b + c)^3/27 + (a*b*c)/2)^(1/3)) - ((((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)^3 + (A/2 + ((a + b + c)*(a*b + a*c + b*c))/6 - (a + b + c)^3/27 - (a*b*c)/2)^2)^(1/2) - A/2 - ((a + b + c)*(a*b + a*c + b*c))/6 + (a + b + c)^3/27 + (a*b*c)/2)^(1/3)/2 - (3^(1/2)*(((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)/((((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)^3 + (A/2 + ((a + b + c)*(a*b + a*c + b*c))/6 - (a + b + c)^3/27 - (a*b*c)/2)^2)^(1/2) - A/2 - ((a + b + c)*(a*b + a*c + b*c))/6 + (a + b + c)^3/27 + (a*b*c)/2)^(1/3) + ((((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)^3 + (A/2 + ((a + b + c)*(a*b + a*c + b*c))/6 - (a + b + c)^3/27 - (a*b*c)/2)^2)^(1/2) - A/2 - ((a + b + c)*(a*b + a*c + b*c))/6 + (a + b + c)^3/27 + (a*b*c)/2)^(1/3))*i)/2

a/3 + b/3 + c/3 + ((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)/(2*((((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)^3 + (A/2 + ((a + b + c)*(a*b + a*c + b*c))/6 - (a + b + c)^3/27 - (a*b*c)/2)^2)^(1/2) - A/2 - ((a + b + c)*(a*b + a*c + b*c))/6 + (a + b + c)^3/27 + (a*b*c)/2)^(1/3)) - ((((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)^3 + (A/2 + ((a + b + c)*(a*b + a*c + b*c))/6 - (a + b + c)^3/27 - (a*b*c)/2)^2)^(1/2) - A/2 - ((a + b + c)*(a*b + a*c + b*c))/6 + (a + b + c)^3/27 + (a*b*c)/2)^(1/3)/2 + (3^(1/2)*(((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)/((((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)^3 + (A/2 + ((a + b + c)*(a*b + a*c + b*c))/6 - (a + b + c)^3/27 - (a*b*c)/2)^2)^(1/2) - A/2 - ((a + b + c)*(a*b + a*c + b*c))/6 + (a + b + c)^3/27 + (a*b*c)/2)^(1/3) + ((((a*b)/3 + (a*c)/3 + (b*c)/3 - (a + b + c)^2/9)^3 + (A/2 + ((a + b + c)*(a*b + a*c + b*c))/6 - (a + b + c)^3/27 - (a*b*c)/2)^2)^(1/2) - A/2 - ((a + b + c)*(a*b + a*c + b*c))/6 + (a + b + c)^3/27 + (a*b*c)/2)^(1/3))*i)/2
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
4楼2013-06-10 18:02:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询


csgt0: 金币+1, 欢迎交流 2013-06-14 10:42:19
引用回帖:
4楼: Originally posted by dbb627 at 2013-06-10 18:02:47
可以,列个公式就可以
x一实根,2虚根
x=solve('(a-x)*(b-x)*(c-x)=A','x')

x =

                                                                                                             ...

对“一实根,2虚根”说法有异议。这完全取决于a,b,c,A。如果a,b,c,A都是实数,必然有一个实数解,其他两个则不一定。a,b,c复数, A=0,则是三个复数解。A=0,a,b,c是实数,就有三个实数解。再者,这样冗长的解析解在实用中不知道有什么意义。
Excel计算容易实现求实数根,见附件。复数运算我还不会。如可以实现复数运算,也可以用newton法解

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : CubicEqaution.xls
  • 2013-06-11 09:23:36, 14.5 K

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

5楼2013-06-11 09:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见