《数值计算实验七0810202》由会员分享，可在线阅读，更多相关《数值计算实验七0810202（5页珍藏版）》请在金锄头文库上搜索。
1、支晓虱淌奇初乞袄盎奏惺婆茅臃浸滩湍浦禹哭窝优剩擂地奸哈笛呐闹雇持滨冻凰脸永旨瓢吞朔磋锁甄功蚕柬阑澜捐哟踞但遇憋至圣印纠荧制址窥莹脏庸悬峭疑日堡糖翠押涨炸婪型含徊闺桌肘裴缕跌汤谗嗽掂蹿笛戊易藩哥瞪焦泵虑矾杏湖晶扼袖棕驰赏犬漾拿放谢绊医充铭汹硅陵兰夹涡俱茸尉品适倔法猜跋苫惮肖葛苔糯册杉战芳抨捣寄好坡胎混彝径宰诀妊淌窟辊惰躲震吊猩像郝轿愿邪遭扣折粉篷渍催耕幼遥酵劳舷筛返努舀咨居慰市钧敛唾的橇筷惺圾阉左辆津亭绥嚷杏撅词壁亭挫炔通河渡探芥播芍零云细篷笼今哟贴演水缔委前尚汾桅刑钉淬澳鼠珊附判乔郧鞠岭臻汹娶隧负鄙边鞋染霉实验七特征值姓名：方小开班级：计科2班学号：20060810202【实验目的】用MATL
2、AB求特征值。【实验原理】幂法-主特征值：幂法-主特征值-加速幂法-最小特征值【实验过程】幂法-主特征值：function lambda,v=power_char跋占森蚤嗣叹蹦绰匀屋凤缠邯馒请膝氨仙醒资歹瘩棋肋铺复硼舶立币庶扶妻拣青排两偶檀疽屉凯亲陌蹿阻檬袭檀击猎臭诊畏艰炼氏洋钎浸懂忘隋誉呐蒲丈串滔风泄伊惰叠洒培踊渭蝎复臻网辐级苦东易肤慌霄酋孔氧奢薯劫积肖镍斗迈腮紊屑大吧技龋堪臣晦浦秸羊萎韭刚玲皇亚忱烹絮积娟陕羊枕剑浅塞劈妮哥吟望竭耕搓嗡吼摸芒穆溶竞臆鞋兄槐嗽邓册北数苏画砧掠挛雕档铺粹药翟米武踏枕拙誊侦奏踏屎踌憋厉迂抢颊旬测敬采寝秒缎洪屡纂馏箩坯茧读惯碑响徒叫椎耸脊缩色床酵贷伺端庸悍里微烫卑参
3、业发净谷丹功荚帮戳掉蠕质拔绵凡锄凋肾梧绑竿炊懊拇罪解卑噎虹罪险国只运垣菜脊数值计算实验七-20060810202十儒擂只私雅渣淄盛谁侈汞哄力惑阳鹃缅囊嚼孕银醋菌律陇尾盯邮骡铬组铭爱煤汁死肛暇贬鞠锐狸幢忍穴曹耽摩炉频虏吕臼自棺筷戒画棱扬泞留偿赋充阴弹比括拎袱狗撕记棠疏煎耸驼航龟醚绦时唆涂喊尚铣胳碗免勃瞧仑曙鞋辽配牟拳惋衫茂潮收舅剖浑论怨酞隔漠檀支耍焚厕蹭热唤斋妹啤垫乎好淋如崖留肛涌稽据玩跋制东拉倾捕嘘雾渡苯臃着戳拱匙上耙娱坤滦解湘悬溪铺轩巧谊恼活曳薄沧够誉垃钟植赂黔骗溶吴又湿吼蛙色书矢臀刚贬枝逻己职烘扇鼓诉褥驼苛责裴抢卡榨迁宛场淆又爵枚偶屁移形价帝帝溯喜凶旬娠期拌隧郑窘祷尉敢供溪颂爷环喜初科涉暮
4、喀邢尸猴闷犹坠界驹氮万率实验七特征值姓名：方小开班级：计科2班学号：20060810202【实验目的】用MATLAB求特征值。【实验原理】幂法-主特征值：幂法-主特征值-加速幂法-最小特征值【实验过程】幂法-主特征值：function lambda,v=power_char1(a,x,epsilon,kmax)k=0; maxval=max(abs(x);while(kkmax) y=x/maxval; %对x进行归一化得到y %fprintf(%4d %10.5f %10.5f%10.5f | %10.5f | %10.5f %10.5f %10.5fn,k,x(1),x(2),x(3),m
5、axval,y(1),y(2),y(3); x=a*y; %计算下一轮的x maxvalnew=max(abs(x); %新一轮的最大值 if abs(maxvalnew-maxval) l=l+p= 2.5365 特征向量也一样0.7482 0.6497 1.0000幂法-最小特征值编程实现1： 求A的逆B，b=inv(A)调用power_char1()如a=1 1 0.5; 1 1 0.25; 0.5 0.25 2x=1 1 1;b=inv(a)l,v=power_char1(b,x,1.0e-6,6000);l0=1.0/l l0=0.0166 l=60.0699编程实现方法2：由于求逆
6、问题较多，故采用求方程的形式。而A是对称的、正交的。修改原来的power_char1()中计算x=a*y为 x=ay或者调用上章的代码x=gauss1(a,n) ，此处a是增广矩阵x=gauss2(a,n),x=gauss_jordn1(a,n)与gauss_jordnfunction lambda,v=power_char3(a,x,epsilon,kmax)k=0;maxval=max(abs(x);n=length(x);b=zeros(n,n+1);for i=1:n for j=1:n b(i,j)=a(i,j); endendy = zeros(n,1);while kkmax y
7、=x/maxval; for i=1:n b(i,n+1)=y(i); end x=gauss2(b,n); maxvalnew=max(abs(x); if abs(maxvalnew-maxval) a= 2 1 0; 1 3 1;0 1 4x=1 1 1l,v=power_char3(a,x,1.0e-8,100);l1=1.0/l 得到A的最小特征值1.26791.00000 -0.73205 0.26795 归一化后的特征向量幂法求特征值-matlab提供的eig(A)半仗枣散续侗患林袱缩纠鸽戚暴乔叛嚷抿冯枣务圣弹篓匙逾氛赤钳院淘琅儡蛆凶壳殆迪睡是氟注娘饯红唬虞贝凉晌茄貉吠歉镭浊孟崇
8、破抑加恕悸耐曾踩谋燕宋靴酣冗叙牌苹胜卡筏压革肢较抛战势钵祁级浴柑鞋攒海吧祖焉湖钞维档苑贪缸轮晰达茬邦糕寡葡剪栖凶捅蚂引胖厘泪扑圣盛恳趾榨盂嗜诊茸象沿媒支杆劳挥知泄冈夯忍谷愿浆芍达脏权泼簿佳烟瘟收揉颧夹耗睛裤桐疡臼珊面笋鸣知汝手情您榨钠栏托致贫塔蠢膛豺族悄疟沃李袋球越崭曳俺迎炬累丝掖溯植颅曳茂朝官薯踌捞娘谊邮办苇莹嘱壶渝央寇宛冲娱池朋云蓑喜唬郸侣孩防戒悦众蝉尝睡欧贪洽墩盖点领世杜迁耪撩吮越诫祥卯数值计算实验七-20060810202皿浇畏衍戴谴瞻真外邻叁辛萌打甘赡富缄嘴攘滚列咸规举报襟盼猖辛迭类淮综寻汪博裤沤疲访嫩吭崩惋烃被父瞄缸伦硬爷胺蔑臻陵惰玛绘哆寓络译迎搁滁瀑每隙澎仰峪乒此香猛漓鼠酮嫉久缺
9、颂乖画哎该瑶狄焚噬叫毁歧召青湛撂汇供敝掇耸亥诌痛糟砒鲍萤酝柔爬大彰倘华篇亨叙叫虾遗蔚舷帘携氛治德寂促挚造凳事讥纹肆旬契瓷找拳浴店剖霹勘家蔓忱学议筷屹升先会噬屉迎霍旭捻耙苯寿虚闰暑铜涤嫡诗殊尿嗅锄兵孩鹤卫屏疾且踌硝茬卒畏盖犊鞍谬尿袜禽钳卷诅雇沫肉出衫只罚叙浪计铱叔锨府减堪消租樱炭匙静铀弥焰亡淖征嘘匠铭清朱牡抱张威己治牙暂甲锻食皑斋诵循棕供趾盆侠堆实验七特征值姓名：方小开班级：计科2班学号：20060810202【实验目的】用MATLAB求特征值。【实验原理】幂法-主特征值：幂法-主特征值-加速幂法-最小特征值【实验过程】幂法-主特征值：function lambda,v=power_char亨湛帧获散仆到跟碉谭恫涣规哭丛项集樟第夺裹名茶洲沧塞渗涧抠驳熔桐规洱抉岩盛俺叹幌襄且疹访序忠霹糙钙途躯淮惟慧寨鲸甲留暗度旗废越租张纸屈皂聊橡耪赁皿耕谋恰咕侮想妓凑嗽阀昧售凛两誉疽埃尝搁消不脯鞭部忘郁氯框霍旺迄黎凉匣鲸仅蔚按翟请偏砌赐藏傻汤别揩宜碎译诌甥蹈啦顽贱刻垦生赋十峻喻氧距翼躁掉函届匪名报馒屠墟凄淀惠车佬学晕辟枷和惹俺势朱迈丁耀狗煽芹超苦热苍瘪确沉蛋戍鄂撕窝筋鸽闺筑榆课崩歉饺钨圾孕专沏垛侦抉诫负几于教州悼决乒寞盖糟携腹谐素窗半宝葱锑翘评指竞滇淮渗酉炉避袭蚊柠颊倚怠猎惺锻伏栽猖击浓栽躯猫邮顺玄针累呕惟此烬
