24小时热门版块排行榜    

查看: 1122  |  回复: 4

passionfly

金虫 (小有名气)

[交流] 【求助】用最小二乘法求解线性方程组的Fortran代码 已有4人参与

For a standard linear system, y(n*1)=A(n*n)x(n*1)

If y is exact and A is well-conditioned, it is easy to calculate x.

However, if y has some disturbances or some errors, at the same time, A matrix is ill-conditioned. It is difficult to accurately obtain vector x. Alternatively, we can measure more components of y (becomes m*1, where m>n) and use the least square method to calculate the approximate x vector.

It is a common problem so I guess there should be some ready-to-use Fortran codes. Anybody knows where can I find this algorithms or codes? Any advice or suggestion is welcomed!

Thank you very much and with best wishes,
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

貳零壹壹

铜虫 (小有名气)

这个对我有难度,呵呵
2楼2010-10-30 20:04:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ptolomaeus

铁杆木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
cqsmath(金币+1):谢谢分享! 2010-11-08 21:40:44
LAPACK里有成熟的代码
http://www.netlib.org/lapack/
LAPACK是当今计算数学界最权威的软件包,实现了几乎所有成熟的数值代数算法(F和C)。
3楼2010-11-02 23:00:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

passionfly

金虫 (小有名气)

呵呵,非常感谢!
4楼2010-11-07 11:22:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer



小木虫: 金币+0.5, 给个红包,谢谢回帖
现在用的Fortran是哪个版本的?老夫学的时候还是Fortran 77呢,已经是将近30年前的事了。
5楼2013-11-01 16:11:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 passionfly 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见