基于MATLAB的数值分析.docVIP

基于MATLAB的数值分析 数值分析是一门计算量大、算法多、实践性强的应用性数学学科，其中每一种算法都是源于实际计算的需要而寻求数学方法的产物。利用MATLAB软件的数学计算功能和友好的图形界面，对数值分析中涉及到的基本问题进行辅助设计，不仅高效、准确，而且可以直观的看出结果，便于对算法进行整体的描述与分析。基于上述思想，我们可以在MATLAB环境下进行数值分析的研究与应用。 一、幂法的MATLAB程序 用幂法计算矩阵的主特征值和对应的特征向量的MATLAB主程序 function [k,lambda,Vk,Wc]=mifa(A,V0,jd,max1) lambda=0;k=1;Wc =1; ,jd=jd*0.1;state=1; V=V0; while((k=max1)(state==1)) Vk=A*V; [m j]=max(abs(Vk)); mk=m; tzw=abs(lambda-mk); Vk=(1/mk)*Vk; Txw=norm(V-Vk); Wc=max(Txw,tzw); V=Vk;lambda=mk;state=0; if(Wcjd) state=1; end k=k+1;Wc=Wc; end if(Wc=jd) disp(请注意：迭代次数k,主特征值的近似值lambda,主特征向量的近似向量Vk,相邻两次迭代的误差Wc如下：) else disp(请注意：迭代次数k已经达到最大迭代次数max1,主特征值的迭代值lambda,主特征向量的迭代向量Vk,相邻两次迭代的误差Wc如下：) end Vk=V;k=k-1;Wc; 用幂法计算下列矩阵的主特征值和对应的特征向量的近似向量，精度.并把（1）和（2）输出的结果与例5.1.1中的结果进行比较. (1)； (2)；(3)；（4）. 解 （1）输入MATLAB程序 A=[1 -1;2 4]; V0=[1,1]; [k,lambda,Vk,Wc]=mifa(A,V0,0.00001,100), [V,D] = eig (A), Dzd=max(diag(D)), wuD= abs(Dzd- lambda), wuV=V(:,2)./Vk, 运行后屏幕显示结果 请注意：迭代次数k,主特征值的近似值lambda,主特征向量的近似向量Vk,相邻两次迭代的误差Wc如下： k = 33 lambda =3.00000173836804 Wc = 8.691862856124999e-007 Vk =-0.49999942054432 -0.70710678118655 0.44721359549996 -0.89442822756294 V = 1.00000000000000 0.70710678118655 wuV =-0.89442719099992 -0.89442719099992 Dzd = 3 wuD = 1.738368038406435e-006 由输出结果可看出，迭代33次，相邻两次迭代的误差Wc 8.69 19e-007，矩阵的主特征值的近似值lambda3.000 00和对应的特征向量的近似向量Vk （-0.500 00，1.000 00, lambda与例5.1.1中的最大特征值近似相等，绝对误差约为1.738 37e-006，Vk与特征向量 的第1个分量的绝对误差约等于0，第2个分量的绝对值相同.由wuV可以看出，的特征向量V(:,2) 与Vk的对应分量的比值近似相等.因此，用程序mifa.m计算的结果达到预先给定的精度. 二、输入MATLAB程序 B=[1 2 3;2 1 3;3 3 6]; V0=[1,1,1]; [k,lambda,Vk,Wc]=mifa(B,V0,0.00001,100), [V,D] = eig (B), Dzd=max(diag(D)), wuD= abs(Dzd- lambda), wuV=V(:,3)./Vk, 运行后屏幕显示结果 请注意：迭代次数k,主特征值的近似值lambda,主特征向量的近似向量Vk,相邻两次迭代的误差Wc如下： k = 3 lambda = 9 Wc =0 Dzd = 9 wuD = 0