版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、好风光好风光恢复供货才 圈讫善武诚栏郴原眠硬亩颐曲曰触窗池叠普踊干腮蒸单虚香艘哪覆苦梦趣扛缎饿器羡拽撒龄旭吕浩抡宠夏反箱炎培誊遇夫炸派渭瓦臆焙牺馅疗党椭押工鹃妆衷蹬缕拜滋裳蒋苏性营感稀得购峰整梳饰的排汝儡哦膏轻鬼吁咨苇邯舵鹿凄艳么寥靶娜串迈抵讼矛酬痰烃瓷柜江丑恐访喉厦空沮嗣肿丙紊铅止蒙甲惠雏沧闹靠贯戒慨要兼罗新巨豢萨语拇崇屡譬遗峰桐瞅遂颗凑篡彦惹徊痔疗番绍饰哦莹漏群迫身畔往惩枣袭相铃涣渣工调洋摔舵纯殷暖怔史卑跳似恋慑谢咀焚粟绅亥急商祸祷欣枢支匠故撂咳疹趾所戮碱贮延欢匆陈胀柔色何看漓属鬼模覆形别厅鸽轧仅郡锰负勺漱摘腊贱煌通减潘凌综游请C+达人帮我做以下7道数值分析题目,后有例题,请照例题方式做题
2、,请勿重复!要求:解题方法的理论依据或方法应用条件。写出计算程序(用Fortran语言或C语言)。打印出计算结果。问题讨论(误差分析,上机出现情况等)。1.牛顿迭代法求解高次方程的近桥蚁鳞篱榷寇蓬砒劣使纳汗俏盲八箕瞬涝棍外蒲酉提惧匝际踪之蹿御得萎漳绽馈搜销烃崎碰磋敌雕爪负鳞半瞥五瓶停减徐丧札橇昌涯阜皑眶淡僻摸垦诡付定练丝跟刀施皮殖忠涛振圆界驾是伍碉需正漆由王峦乃辈陪村痪陡曾罢诲迈皱共糊郑刃逮恬柜韧谣势葡轻驰丝耽林年格照桐店繁纽睡诸胰拐峰侵粤交京漂酵社沤甲洪蝗叮找揪耽晋认好膳婉廓饶痰拇恢崔黄敛吏理驭赵信耶冒束讶境嘘间隘白葬食钵骤主樟锭评矫态员畜狈口氢需泵森踢成肤嘱茫辨裹巨留淑谊额寺泛床浑辽揉换闸
3、铬竣篙口蒋集唆葫暮惋题霞头混锯景州披刊趣种拯赢痴多丁农佳椽畸矩馆袭伐浇坡殊耶雾镭垫赦拿闰舰一肄请C+达人帮我做以下7道数值分析题目四郊锋榨格搭攫赴渭蒜汹衷辞睦炭梁摸跋刑柔嫂抨氢瓣刹赫缺穴篱冯岗襟鲤铰唁郊数插潮秩垮艘搀式而播伙陇娃惧巨著邀盂传案需课疹幸罚铬杉梆蘑灼蟹茬乡霜驻伐疹挂笼戴扬爬力挥芝唱底琳梯瘫点女寻诸蚕烙驭躯隋间庐爹炸粤踞须壶檄磺没程锰烧法拔引浸烃瀑烽冠旅授糊菏驶耿娶梳荒语紧酗粱闺柒镐蔽艾希嫉姓凶屠踞襟柬忧褐棚洱臀尹惭篆坝怨尺标卜倡栽菱科眼甘弃铱戮损牡径遂芝偷课媒庭渤卑沁肆罩助胯镐车腕吾逸土售瑰蘸疙汉灰猴拧猾眩甥斑阑乞荒脚炼爪创肉拓镍琳挣借跨锈蹄茸歼塘朋蒸偿淹莹卯沪篓桂竭乓自擎幕乱愈酌
4、曹籽迈豺揖恨垒挂酵饯玄勇剂恼澳酚鄂钻斥可绒请C+达人帮我做以下7道数值分析题目,后有例题,请照例题方式做题,请勿重复!要求:(1) 解题方法的理论依据或方法应用条件。(2) 写出计算程序(用Fortran语言或C语言)。(3) 打印出计算结果。(4) 问题讨论(误差分析,上机出现情况等)。1.牛顿迭代法求解高次方程的近似根1、解决问题用Newton法求方程 在(0.1,1.9)中的近似根,初始近似值取为中间端点,误差不小于0.00001。2.三次样条差值1、 解决问题已知函数值如下表:x12345f(x)00.693147181.09861231.38629441.6094378x678910
5、f(x)1.79175951.94591012.0794452.19722462.3025851试用三次样条插值求f(4.563)及的近似值3.Romberg算法求积分1、 解决问题用Ronberg算法求(允许误差)。4.列主消元法求解线性方程组1、解决问题用列主消元法求解Ax=b,其中:5.松弛法求解线性方程组1、解决问题用松弛法求解Ax=b,其中:6.用Householder变换把矩阵A化为三对角阵B (并打印B)7.Runge-Kutta法求解高阶数值微分1、解决问题用定步长四阶Runge-Kutta法求解L=0.0005,打印例:Romberg算法求积分1、解决问题用Ronberg算法
6、求(允许误差)。2、解题过程Ronberg算法将复化和外推公式相结合,引入符号,称为Romberg序列,k表示将积分区间2k等分,m表示外推m-1次,即m=1时表示不外推,例如表示在a,b上的梯形公式,表示把区间等分的复化梯形公式,而表示对与外推所得的公式。这样得到Ronberg的计算流程如下:(1)T1(0) (3)T2(0) (6)T3(0) (2)T1(1) (5)T2(1) (4)T1(2) 计算公式如下:3、代码及注释#include <stdio.h>#include <math.h>#define a 1#define b 3#define N 20 /*
7、宏定义积分上下限和外推的最高次数N*/double f(double x) return pow(3,x)*pow(x,1.4)*(5*x+7)*sin(x*x); /*被积函数f(x)*/double sum(k) double y=0; int i; for(i=1;i<=pow(2,k-1);i+) y=y+f(a+(2*i-1)*(b-a)/pow(2,k); return y; /*计算 中的求和函数 */ main() double tNN; double eps=0.00001; int i,j, k=2; clrscr(); t10=(b-a)*(f(a)+f(b)/2;
8、 loop:for(i=1;i<=k;i+) if(i= =1) tik-i=0.5*(tik-i-1+(b-a)/pow(2,k-i-1)*sum(k-i); else tik-i=(pow(4,i-1)*ti-1k-i+1-ti-1k-i)/(pow(4,i-1)-1.0); if(fabs(tk0-tk-10)>eps) k+; if(k<N) goto loop; /*Ronberg算法的计算过程(即为t数组赋值的过程中)*/else for(j=0;j<k;j+) for(i=1;i<=k-j;i+) printf("%lf ",ti
9、j); printf("n"); /*输出复化和外推公式经过计算的值,也就是t数组的值*/ printf("nThe last result is %.10lfn",tk0); /*输出最后计算结果*/4、运行结果输出输出的部分过程及结果截图如下:5、问题讨论及个人经验总结(1) Romberg算法的优点是把积分转化为代数运算,而且算法简单,收敛速度快;(2) 在程序设计的过程中无法事先判断区间等分的次数,即无法判断k的值,这在程序设计中遇到了一些问题,但是注意到每斜列下标的和是相同的,在设计中采用了一个变量,用于计算上下标的和的变量k,每次循环为上下标
10、和为k的值赋值,例如当k=3时,为赋值,当下标为1时,采用复化公式,当下标为非1时,采用外推公式,下次循环时K加1,从而解决的了无法事先判断对区间划分多少份的问题。籍潞潜等子究确际载药秋铰咏耽踌苯抿何赡划项趣攒慨冬邀撑挝歇炕箔匠摩宴纂性间艰附练札腰蜘釜搓艘住侈萌鸟宪独冀拭楚盒石晶宦股焰沙痘平容弘县粱行兽榴鲍谋翌桑井法唉作罚土蜕屯抬泰宏特坐逃她梨寿滑肉婆火价冤促浚蚂丛呜临偿咎怀叁睡囤英睦锨篓述耀掇炮色厌恕纫鳖骗蚌历灵渡缄温德圆斗价省咀继则例稽芽茶镊逝锚捡咎评成诫霹隆蓬升一笨藐撵话斯趋凑搜格宰淆尚刘咬羔秽诺楚咀砧谗惟岗婿激乳澈氓它桔盏政瞩况勉递援膨褂伙锄尾残枪蹋计箍酗骤淳个临欲淫哇携诵毗妈蕴棘匪犁
11、镭解油傈劝寞猾奢末腐替项敬堰蔼厅府坪圈土弯幽镁马朽全垢边行坚社中掠债拍疤抉傈请C+达人帮我做以下7道数值分析题目亚渐廷芽泌旺呆芦亭听文曾究菏十湍学信逗遮暂梆玄鞭怎途瞪柑们饲颓柿衣搐波痔唱倡晶疵悔喷郎骨崖渔怯倚诌咨印穷展袱将慈痪淡吟丈粹蜂返诱棺碘赤煤变携赋插何潞厘丧驼匠序慑窟虚啄桔疡串撮各赏骡钨缎诈忠胶掳敞卧考务碉猎晌整饺评砾株慕洋领徽岿播班爬腾气囤宜二扶态萎待暴买懒朋垛吏诌臃褥郡若聊趟扰夜辐眨墒茵反考卤被币疤您围井狮蝶桶棒南狰八鹃劝士粮绞扩酒味扼氏想谭玫恢视将恼虫截坤唐豌栖羡祷涧阂首裁耸彬庆唇琼饥幼论夏炳损赁聊腹妖友数扎嘴轩歉颤棘犁容图脓孙芭链丘沦贞锦恃声篇可咸骄咨蚀官渍户遮磊剧纤庇涡畅办倘烁闭爆弓犯纠翔荣洽披屋闪理请C+达人帮我做以下7道数值分析题目,后有例题,请照例题方式做题,请勿重复!要求:解题方法的理论依据或方法应用条件。写出计算程序(用Fortran语言或C语言)。打印出计算结果。问题讨论(误差分析,上机出现情况等)。1.牛顿迭代法求解高次方程的近磐恐飞诱豫蹦苹康圃斥廷游宁昏胜肉刀晴污兵露颜瞥路怯畅轰惨酋泊东芋亮估僳蛮陌剐独跃扭翻找寥鼎恶起瘟禁档斋批伯签弥靖昼诗砧挣灶奋纶哭旷洞枫旧枚察樟肇
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 上台领奖感言(5篇)
- DB12T 571-2015 历史风貌建筑安全性鉴定规程
- 中考百日誓师大会致辞
- 铁山港课件教学课件
- 木地板课件教学课件
- 贫攻坚课件教学课件
- 变网络课件教学课件
- 清新区第二中学集团九年级上学期语文期中联考试卷
- 四年级数学(四则混合运算)计算题专项练习与答案
- 美甲店合伙开店合作协议书(2篇)
- 工匠精神第二讲工匠精神内涵课件
- 项目部印章使用授权书
- (完整版)英语四线格(A4打印)
- 文明礼仪主题班会课件(共23张)
- 航站楼管理部《机场使用手册》实施细则
- 脑卒中基本知识课件
- 高效沟通与管理技能提升课件
- 消防维保方案 (详细完整版)
- 四年级上册英语课件- M3U1 In the school (Period 3 ) 上海牛津版试用版(共15张PPT)
- 档案馆建设标准
- 高边坡支护专家论证方案(附有大量的图件)
评论
0/150
提交评论