标题:[讨论][挑战]解方程程序
只看楼主
乌鸦丘比特
Rank: 1
等 级:新手上路
威 望:2
帖 子:625
专家分:0
注 册:2004-7-19
 问题点数:0 回复次数:24 
[讨论][挑战]解方程程序
题目描述
f(x)=a[n]*x^n+a[n-1]*x^(n-1)+......a[0]*x^0
是一个一元N次的多项式,问题是编写一个程序,能够解出方程f(x)=0的所有解,精确到0,01;
程序越完善越好哦.
输入:
一个整数N
接下来输入N个实数,分别对应a[n],a[n-1].....a[0],系数为0也照常输入

请大家说说自己的思路.
如果可能试下写个程序.

该问题简化版本:N=3.如果任意N不会可以先做这个..
搜索更多相关主题的帖子: 解方程 挑战 
2006-07-19 14:40
穆扬
Rank: 1
等 级:禁止发言
帖 子:1910
专家分:0
注 册:2006-6-1
得分:0 
提示: 作者被禁止或删除 内容自动屏蔽

2006-07-19 23:52
乌鸦丘比特
Rank: 1
等 级:新手上路
威 望:2
帖 子:625
专家分:0
注 册:2004-7-19
得分:0 

是的,但是这第一个也不太容易解


我喜欢创造,一只扑腾着翅膀向天空飞翔的乌鸦
2006-07-20 10:00
ysol
Rank: 1
等 级:新手上路
帖 子:107
专家分:0
注 册:2006-6-11
得分:0 
第一解让计算机试根吧,比如随便找2个数a,b
f(a)>0,f(b)<0,a,b间就必有一解,然后用插值法逐步逼近。
关于插值方法最简单的就是2分法了,就是速度不敢恭维了,不过对于计算机来讲只要N不超过1000就绰绰有余了

个人意见仅供参考!!!
2006-07-25 23:37
穆扬
Rank: 1
等 级:禁止发言
帖 子:1910
专家分:0
注 册:2006-6-1
得分:0 
提示: 作者被禁止或删除 内容自动屏蔽

2006-07-25 23:48
ysol
Rank: 1
等 级:新手上路
帖 子:107
专家分:0
注 册:2006-6-11
得分:0 
无解还用试吗????
2006-07-25 23:50
乌鸦丘比特
Rank: 1
等 级:新手上路
威 望:2
帖 子:625
专家分:0
注 册:2004-7-19
得分:0 
无解需要判断无解

我喜欢创造,一只扑腾着翅膀向天空飞翔的乌鸦
2006-07-26 08:05
穆扬
Rank: 1
等 级:禁止发言
帖 子:1910
专家分:0
注 册:2006-6-1
得分:0 
提示: 作者被禁止或删除 内容自动屏蔽

2006-07-26 08:56
穆扬
Rank: 1
等 级:禁止发言
帖 子:1910
专家分:0
注 册:2006-6-1
得分:0 
提示: 作者被禁止或删除 内容自动屏蔽

2006-07-26 08:56
ysol
Rank: 1
等 级:新手上路
帖 子:107
专家分:0
注 册:2006-6-11
得分:0 
以下是引用穆扬在2006-7-26 8:56:10的发言:

题目要求所有解

无解算不算所有解??????你的那点数学常识不要在我面前显摆,起码在你连我的基本意思都没有搞清楚之前不要说一些不经过大脑考虑的话。

无解的判断实际上可以放在寻找ab的过程中,设定经过一定数量级的2分插值还不能满足f(a)>0,f(b)<0基本上可以判断在
区间[a,b]上无实根,方程的求解过程关键在于预测该方程的解区间也就是a,b的值的确定。这只是简单的解方程思想。

如果要精确的求解恐怕要用到拉格朗日多项式插值公式和矩阵,特别是Vandermonde矩阵或者Toeplitz矩阵,可惜Vandermonde问题本身是一个严重病态的,解搂主的方程个人感觉最好用Toeplitz矩阵,不过问题是相当的........................哈哈!楼主要是有兴趣可以共同探讨!!

[此贴子已经被作者于2006-7-26 13:35:19编辑过]

2006-07-26 13:32



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-79169-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.163931 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved