| 查看: 733 | 回复: 5 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
ldoop
|
[求助]
高手们，请帮忙解一个三元方程组啊，谢谢了
|
|||
|
b*x -a -b-sqrt((b*x-a-b)^2-4*ab)=b_0*x_0 -a_0 -b_0-sqrt((b_0*x_0-a_0-b_0)^2-4*a_0*b_0 b*x -a -b+sqrt((b*x-a-b)^2-4ab)=b0*x0 -a0 -b0+sqrt((b0*x0-a0-b_0)^2-4*a_0*b_0 x*b/(a+b+2*sqrt(a*b))=b_0*x_0/(a_0+b_0+2*sqrt(a_0*b_0)) 其中a_0，b_0，x_0为已知常数，求未知数 a b x [ Last edited by ldoop on 2013-5-17 at 23:42 ] |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有190人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求解一元三次方程
已经有11人回复
- 用matlab解带参数的一元三次方程--求解
已经有3人回复
- 求Microsoft Math3.0高手，帮画一张一元三次方程图，最后再往邮箱发一份。
已经有10人回复
- 求Microsoft Math3.0高手，帮画一张一元三次方程图
已经有9人回复
- 求高手帮我编一个一元三次方程求解的小程序！重谢哦！
已经有5人回复
- 【求助】关于MATLAB解三次方程
已经有7人回复
- 【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】
已经有15人回复
- 【求助】如何用Fortran解一元三次方程【已解决】
已经有6人回复
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
6楼2013-05-21 07:03:08
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68128.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fegg7502: 金币+2, 应助指数+1, 鼓励交流 2013-05-18 14:43:15
ldoop: 金币+10, ★有帮助, 谢谢回复，我是想求得解析解而非数值解，不过还是感谢了，麻烦帮忙研究下是否能求得解析解 2013-05-20 16:26:21
感谢参与，应助指数 +1
fegg7502: 金币+2, 应助指数+1, 鼓励交流 2013-05-18 14:43:15
ldoop: 金币+10, ★有帮助, 谢谢回复，我是想求得解析解而非数值解，不过还是感谢了，麻烦帮忙研究下是否能求得解析解 2013-05-20 16:26:21
|
MATLAB程序，复制进一个m文件，F5运行即可，根据自己的需要输入已知参数，即返回相应的解： function solvenole314 clear all;clc format long global a_0 b_0 x_0 a_0=input('请输入已知常数a_0='); b_0=input('请输入已知常数b_0='); x_0=input('请输入已知常数x_0='); aint=input('请输入未知数a的猜想值='); bint=input('请输入未知数b的猜想值='); xint=input('请输入未知数x的猜想值='); options=optimset('TolFun',1e-9,'MaxIter',1000); x1=fsolve(@fun299792,[aint bint xint],options) %-------------------------------- function y=fun299792(z) global a_0 b_0 x_0 a=z(1); b=z(2); x=z(3); y(1)=b*x-a-b-sqrt((b*x-a-b)^2-4*a*b)-(b_0*x_0-a_0-b_0-sqrt((b_0*x_0-a_0-b_0)^2-4*a_0*b_0)); y(2)=b*x-a-b+sqrt((b*x-a-b)^2-4*a*b)-(b_0*x_0-a_0-b_0+sqrt((b_0*x_0-a_0-b_0)^2-4*a_0*b_0)); y(3)=x*b/(a+b+2*sqrt(a*b))-(b_0*x_0/(a_0+b_0+2*sqrt(a_0*b_0))); |
2楼2013-05-18 09:43:50
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68128.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2013-05-18 09:47:58
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
ben_ladeng: 金币+2, 欢迎常来 2013-05-19 23:09:22
ldoop: 金币+10, ★有帮助, 谢谢回复，我是想求得解析解而非数值解，不过还是感谢了，麻烦帮忙研究下是否能求得解析解 2013-05-20 17:02:46
感谢参与，应助指数 +1
ben_ladeng: 金币+2, 欢迎常来 2013-05-19 23:09:22
ldoop: 金币+10, ★有帮助, 谢谢回复，我是想求得解析解而非数值解，不过还是感谢了，麻烦帮忙研究下是否能求得解析解 2013-05-20 17:02:46
|
你这是非线性方程组，最好用Newton法求解。设 三个函数 f1=b*x -a -b-sqrt((b*x-a-b)^2-4*ab)-b_0*x_0 +a_0 +b_0+sqrt((b_0*x_0-a_0-b_0)^2+4*a_0*b_0 f2=b*x -a -b+sqrt((b*x-a-b)^2-4ab)-b0*x0 +a0 +b0-sqrt((b0*x0-a0-b_0)^2+4*a_0*b_3=0 f3=x*b/(a+b+2*sqrt(a*b))-b_0*x_0/(a_0+b_0+2*sqrt(a_0*b_0)) 记3X3矩阵 A=∂(f1,f2,f3)/∂(x,a,b) let dX=(dx,da,db), F=(f1,f2,f3) 解出 dX = -inverse(A) F 然后迭代 x=x+dx a=a+da b=b+db 当然，你要给迭代初值。你这里涉及到sqrt函数，所以给初值要小心。应该有经验。 上面是一般方法。实现需要编程。主要的工作是求解3X3的逆矩阵。 |
4楼2013-05-19 07:30:01
