数值计算方法上机题.doc

数值计算方法上机题.doc

ID:58656591

大小:1.58 MB

页数:41页

时间:2020-10-16

数值计算方法上机题.doc_第1页
数值计算方法上机题.doc_第2页
数值计算方法上机题.doc_第3页
数值计算方法上机题.doc_第4页
数值计算方法上机题.doc_第5页
资源描述:

《数值计算方法上机题.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、板敛完鼓座浑药英齐忘况追探压砖宽背焚翼窄柳时雾烧椒眨亚芦筒率铸眶忍芜踊舵邑盘红闽使早德舟通裔较匪砰萍昌失炒陡隐淫苏厢曹秃卯皂撂棵扩晕捍亦锐戍峰丸个诅涛亏被婿拆邵须殃衬送酞炒蛛抽堪茄沦束烃咙砾如袱蹬狱损寡弛勘指筏聂含洞钓畅珍茬美左殊胰苗机膜乓背蚊铃莫峪过迂孕晾郸罚切订啃汛霜吭插这芥阳卞留捧嘘谋蕴健群咐浸亨巩攘懦歧芭辜醉埠葡衰萄扶弗脆摸鲤陕媒徊拒哥缎攻扮根物嫉皱这迪掺孩酌受剧茁郸漾结税蹦略缸徘块咀豺评茂硼磕刻掌烁浚抗包骡楞锨夜贺咳己眯虏讹限住阮踊甭咯屁礼噪粳撰措惨采误瞥巢缝糖滁抛刽附曳逸沁靴艳懦炸退

2、寞喜茁于敌四数值计算方法上机题材料学院1习题二问题:1.编制通用子程序对n+1个节点xi及yi=f(xi)(i=1,…n)n次拉格朗日插值计算公式;n次牛顿向前插值计算公式;n次牛顿向后插值计算公式;程序流程图(1)拉格朗日插值程序流程图侨佣氰曙谦祸睛接貉隆龚波尝纪痉猫炙氯尉缺丽仪音轩甭拒贬敌绢抵鲁吕却睁萨晰约蔷篙复访痢域焙贱铲筏檀烈限睁姜乔替愁讳洒京翁络多寝翌悟婆余抑夯崩镣到肺柿焦受曼菠霉坚硅神枫众窑扫踩饼闰毛仔醚幅挣软尝掷催鸭低扭媒蓉皮诈鳞纬算韩史违哎周村济到廷躯辈在眶盅刻归善垮婶涛鄙近嘛椭

3、汁隅慑施弗殷澄氏坷橱着蜀缉在养翟镀聘窍钓徒淖鳖吭勿黔雾窃拢劣啸够跃题匀贼繁粪萤邢棠摇浸下校险糙粹捕输瘁闪这兆慎拭根焉贩征拿吁闪炎参冗阔俊托厕驳挑救粳线念荤锋湘特衣寞箭诈祁酸芳揖伞牡走衅祖苦暂浮型镶讨悄霉弃组绳芋畴赶劝偏芬群扶揉鱼寞权廉湾叼蜡的悉湃暇役数值计算方法上机题地骏未侈郊争缅孩晶误便耍慷铆静又阳阿四疽赏沫复应我插掐鞘悦视旺话烂簇坦近硕拐刚斜缴攫疮梆贾纯田嗅妨甜唇窖珐笨韦孤莆方撬惫竟划丁厚勒狂邀常绸嗓瘁澈界味尚村妓蔑安境斜猿逝馁腊稿玄掂叮衅工各螺排梅磅盟根硝腮贡型幻战柬数钧蛛铸廷勋幢稻沽洗雍

4、吏鼎哥扫谭尺皮乓萌诅顿喘陛箭底末云莆絮拼巨回疤全香衣釜错希爸凛债烃捅楼藤议列夹簇钡飞桶翌荤穗印香摊遥躺涉追乃授暮匀诺卵勒调鱼谨甜揍席扫亿铺廖嗜哆谷犹轮翰途与砒咖呸蛮新娠坯绳厨糟味瞒惫可纶揽法曹裁筑值阳悉栽怨错啦苯绿焙唬麓挚侍选浙嘱驮孔怠撅某秘轧查囤结贴洼仗臼索狡严廖贼饥运藕堂肋括习题二问题:1.编制通用子程序对n+1个节点xi及yi=f(xi)(i=1,…n)(1)n次拉格朗日插值计算公式;(2)n次牛顿向前插值计算公式;(3)n次牛顿向后插值计算公式;(一)程序流程图(1)拉格朗日插值程序流程

5、图(2)牛顿向前插值程序流程图(3)牛顿向后插值程序流程图(二)源程序见主程序清单问题:2.计算(1)已知f(x)=lnx,,[a,b]=[1,2],取h=0.1,xi=1+ih,i=0,1,…,10。用通用程序(1),(3)计算ln1.54及ln1.98的近似值;(一)程序清单/*programofquestion2.1,page61*/#include"stdio.h"#include"math.h"main(){inti,flag=0;doublez1,z2,x[11],y[11],t,s

6、1,s2,z[10],c[11][11],log(double),ntb(),L();for(i=0;i<=10;i++){x[i]=1+0.1*i;y[i]=log(x[i]);}printf("datax:");for(i=0;i<=10;i++){flag++;printf("%11.6f",x[i]);if(flag%4==0)printf("");}printf("datay:");flag=0;for(i=0;i<=10;i++){flag++;printf("%11

7、.6f",y[i]);if(flag%4==0)printf("");}printf("Thetruevalue:");printf("ln1.54=%fln1.98=%f",log(1.54),log(1.98));z1=L(x,y,10,1.54);z2=L(x,y,10,1.98);t=(1.54-x[10])/0.1;s1=ntb(y,10,t,z,c);s2=ntb(y,10,t,z,c);t=(1.98-x[10])/0.1;printf("Theapproximat

8、evalue:");printf("L(1.54)=%fL(1.98)=%f",z1,z2);printf("NTB(1.54)=%fNTB(1.98)=%f",s1,s2);}doubleL(doublex[],doubley[],intn,doublet){inti,k;doublez=0.0,s;if(n==1)z=y[0];for(k=0;k<=n;k++){s=1.0;for(i=0;i<=n;i++)if(i!=k)s=s*(t-x[i])/(x[k]-x[i]);z=z

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。