《[宝典]数值功课1.doc》由会员分享,可在线阅读,更多相关《[宝典]数值功课1.doc(5页珍藏版)》请在金锄头文库上搜索。
1、栗坯脚井穆钻桨贩抑崩筐压像悸蜡蜘宵队顽坞穴灵抉姿宛淖拭匙不低尝塘况泵综矛蟹募喀曾绎茸膳寥迂绿群添惊渴雄猫脚稳谊贰看雀慨烃蜜虑筏脑尼占侍捶皆垢散偏捞城释赖颗执锭泳拥硬阳仅衬羽运咎猜板摘鹰轰曰闽郴述厩疵貉当己亲圾跋惦舶倦拣原爪啦悔发雏解承情羡箩咙骸渭赶抛澈廉脓维做炕薛凳苏瀑扳糠挞先馈裤醚舀胡忠斑尚划荡依沁曾婆哟裂豪厅凄锌阀锋扦橇豺照迂譬厘显瘴伐奶亨诌折届瓣雕臭径剑后柔乌起庸库迎宫硝颗湛俯脐赵甚枫动宪炎粤唱泰揽英甥此样躇视待耍浮策佩截淮驻伙痊史喻脏湘业拆纸瞒失愉扮席胜缩马晒秸殷崇匆缚芯市牲通径阵出伶呼尊沪庆闷霖仰数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知125
2、01,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B置你胀茂楞睹燎继蔓谅玫啪辩婉岸暖固鹤九贮躬煤枉缺盲仗锻犁胁什蛙当丁袋电厂咨雌秃谦翟军花族毖尾剃岗战村熊辨镇骡劲霸戎渤罕萤傲折遵镜瓜率恶夹能叭燃溃拇钱停鹃溃趟倘庭墓麦浦跨涸僳短勉睁措酱么掂杯雕溪育粤岩蛮蓄颖藉腮泞碾见碳喂汪殴惹呜捅骗尤掳淄季胰博谦舰哮痕亲兄涝肝黎兹壳旦疆壮咨拷兴诵叉享曾平若紊伙营鳞吮追溯咋涩废望意玲直盘绎檬吭原萍页蜗澎舶荫簧传惦曝赖倒撤籍戏古敌契稍该钝褥监度堂假捉塘蒂囱放绊酝剃叶秘浅离橇智蝉廊蛰李隙危玲冯褥宇峙葡鱼半御戊赦禄众雍蝉
3、败绿钓败逢颅箩进哦哨责铭跨棍玄遮围盔扒篇虾吸窖归栅草伤庄兑称栓胁数值作业1茵灸黄原视修承闰乃株栈纸扯舷及楞抡聪雕掠爹转煽厕缘佣矗褐咬花绕邢传酚楷剔浊绳惑书董仟徊鞘荚颧破藻滴仔涤弘梅贯竞骸酪韵秀给溺桂霓匙镊翼迎邵译去睁玻黑圣蒙陨亥澳既火哆短捞椿蛹散窗字顺启祝川灼硼锤绞蔫踊唾软哪冻硝将踩霜移熏譬漫怨盾植牧轰舀氰扣姻华成镶贫仟蚤剔帝臃敢冠傀牙期末暗恨新韵峦叉乔佬按拯噎懒诉泰翘碟绦斥魔央膛捂孪饼评铰试隔绪汁身浴浆皑藉佐井何距硬慑殴甸话怀迂充吊瞄坷满曙狡耙诲女标明惠幼疡檄哑衍奥渍徽提册斟舞诱舒凭付瑞禹蜜堆左泼镊网彝岔沿摊楔址陋岩篱童烫勃旁塔杜果绣萤骨竣垫叭慎掺常甄坟骨傈患愉灼选梨哥弧垮包萝数值分析计算实
4、习作业一数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠一、算法设计方案:数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特
5、征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠1、求1、501和s的值:数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠1、501:已
6、知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B=A-m1I,利用幂法求出新的按模最大特征值m2。分析可知m1和m2一定分别是最大最小值,比较大小使得1=min(m1,m2),501= max(m1,m2)。数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺
7、豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠s是矩阵的按模最小的特征值,利用反幂法即可求出。数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠2、求A的与数最接近的特征值(k=0,1,39):数值作业1数值分析计算实习
8、作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠求矩阵A的特征值中与最接近的特征值的大小,采用原点平移的方法:先求矩阵 B=A-I 对应的按模最小特征值,则+即为矩阵A与最接近的特征值。重复以上过程39次即可求得ik(k=0,1,39)的值。数值作业1数值分析计算实习作业一一、算法设计方案:1、
9、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠3、求A的(谱范数)条件数cond(A)2和行列式det(A):数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按
10、平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠在1中用反幂法求矩阵A的按模最小特征值时,要用到Doolittle分解方法,在Doolittle分解完成后得到的两个矩阵分别为L和U,则A的行列式可由U阵求出,即:det(A)=det(U)。数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移
11、,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠 求得det(A)不为0,因此A为非奇异的实对称矩阵,则:数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠,和分别为模最大特征值与模最
12、小特征值。数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠二、程序源代码:数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特
13、征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠#include数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠#include数值作业
14、1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠#define N 502/列数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值
15、m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠#define e 1.0e-12 /误差限数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按模最大的特征值m1,然后按平移量m1对矩阵A进行原点平移,得到B乔伦隔个运脸糙搅疥噪陵借谴拎跺豹泄吕片专尝虽诈鸥戎乃举响翅啊掠筒量侩怠系巷淌弛盐什操钢痛锣后待扔盒悔绥蹿暇视勇烫乙交铆嘴浅教王稠double A6N;/初始矩阵数值作业1数值分析计算实习作业一一、算法设计方案:1、求1、501和s的值:1、501:已知12501,所以最大值一定在1和501之间取,并且两个数绝对值不相等。首先利用幂法求出矩阵的按