计算方法的作业——请网友点评一下
本帖最后由 bugatti100P 于 2017-11-5 19:54 编辑这是我自己练手的网站，代码也放这里，图片可能加载不出来，可以转来这里看
或者去我的github主页bugatti100Peagle
https://github.com/Bugatti100Peagle/-/blob/master/03.jpeg
已经不指望老师了，而我又这么爱学习，所以向广大网友求助。这是我们学校《计算方法》的作业，我把题目和我写的代码放上来，请各位朋友看看有什么问题。代码格式问题，或者计算精度问题什么的都可以。欢迎来找bug.
!/*************************************************************************
! > File Name: kt0301_2.f90
! > Author: 青羽100P
! > Mail: 3443988515@qq.com
! > Created Time: 2017年10月11日 星期三 18时03分27秒
! ************************************************************************/
program kt0301_2
!拉格朗日主程序代码
implicit none
integer::pa(17),pa_1(4),i,z,error
real::tem(17),tem_1(4),chads
!使用数组存放数值
open(10,file='sample_unformatted_sequential.dat',f&
&orm='unformatted',access='sequential',status='old')
i=1
do
read(10,iostat=error),tem(i),pa(i)
if(error/=0)exit
i=i+1
end do
close(10)
pa(1)=1000
pa_1=(/pa(1),pa(4),pa(7),pa(10)/)
tem_1=(/tem(1),tem(4),tem(7),tem(10)/)
!挑出1000，700，400，200处的数值作为已知点
!-------------------------
open(11,file="chazhishuju.txt")
do z=1,17
call lar(4,pa_1,tem_1,pa(z),chads)
write(11,"(i2,f10.4,i5,f10.4,f10.4)"),z,tem(z),pa(z),chads,chads-tem(z)
end do
end program
include 'lar.f90'
!/************************************************************************
!> File Name: lar.f90
!> Author: 青羽100P
!> Mail: 3443988515@qq.com
!> Created Time: 2017年10月10日 星期二 21时59分40秒
! *************拉格朗日子程序代码***********************************************************/
subroutine lar(n,x,y,cha,ds)
implicit none
integer::n,x(n),cha,i,j
!n表示数组长度，x横坐标，cha表示要差值x的位置
real::y(n),larlc,ds
ds=0
larlc=1
do j=1,n
larlc=1
!保证不受上一次运算结果影响
do i=1,n
if(i==j)then
larlc=larlc*1.0
else
larlc=larlc*((cha-x(i)*1.0)/(x(j)*1.0-x(i)*1.0))
end if
!此处if 块，实现连乘算式，当算到自己是等于1就可以了，也可以直接用cycle
!每个算式后 *1.0 保证实型运算，否则会被省略，结果全为零
end do
ds=ds+larlc*y(j)
end do
print*,ds
end subroutine 本帖最后由 bugatti100P 于 2017-11-5 19:39 编辑
https://github.com/Bugatti100Peagle/-/blob/master/03.jpeg 弟子不必不如师，师不必贤于弟子。
lz心气颇高，可是大胆的设想往往需要小心的求证诶！
首先，题目中说用N次拉格朗日差值法差值，莫非N可以取正无穷？龙格现象就华丽丽的忽略了？随着N变化，差值误差如何变化？
其次，拉格朗日差值法就没有改进的？知网搜索过相关文献了吗？多项式差值系数为啥不能预先计算好？如何算？两分法求根？伴随矩阵求逆？矩阵求逆的数值方法（直接法、迭代法）？计算方法课那么多内容为啥不能融会贯通呢？
最后，计算方法重在方法的推导与误差分析，切忌舍本逐末侧重代码本身的奇技淫巧 谢谢，受教了 后面的那些方法我还没学，我再去看看课本，希望老师下次继续来点评。 bugatti100P 发表于 2017-11-12 23:37
后面的那些方法我还没学，我再去看看课本，希望老师下次继续来点评。
1、看完课本记得做习题；
2、前面已经说了，计算方法课程主要是学习推导过程与分析，编程只是辅助工具。希望点评计算方法作业，请去知乎或者stack exchange
页:
[1]