19次方系数对方程根的影响课程设计报告书.docVIP

西安郵電學院实验X 19次方系数对方程根的影响 一、问题的提出 考虑代数多项式,显然它的全部根为1,2,…,20;如果19次方的系数有扰动,那么对方程的解有何影响. 二、实验内容 直接利用Matlab中的roots()和poly()函数,选择不同的扰动做数值计算. 主要步骤:ve=zeros(1,21);ve(2)=ess;roots(poly(1:20)+ve). Matlab程序:见附件bingtaiwenti.m 三、实验结果 1.取扰动分别为所得结果如下: 27.0817+5.03812i 21.3025+1.56717i 19.952+0i 27.0817-5.03812i 21.3025-1.56717i 19.2293+0i 19.5337+9.1664i 18.5028+3.6004i 17.6573+0.692896i 19.5337-9.1664i 18.5028-3.6004i 17.6573-0.692896i 13.8235+7.77167i 15.1651+3.76125i 15.4524+0.875524i 13.8235-7.77167i 15.1651-3.76125i 15.4524-0.875524i 10.7211+5.4609i 12.4866+2.88278i 13.3527+0.486992i 10.7211-5.4609i 12.4866-2.88278i 13.3527-0.486992i 8.91282+3.47317i 10.5225+1.71958i 11.8578+0i 8.91282-3.47317i 10.5225-1.71958i 11.0427+0i 7.69268+1.89884i 9.04487+0.59455i 9.9916+0i 7.69268-1.89884i 9.04487-0.59455i 9.00201+0i 6.75761+0.654714i 7.94891+0i 7.99952+0i 6.75761-0.654714i 7.00247+0i 7.00009+0i 5.95208+0i 5.99995+0i 5.99999+0i 5.00061+0i 5+0i 5+0i 4+0i 4+0i 4+0i 3+0i 3+0i 3+0i 2+0i 2+0i 2+0i 1+0i 1+0i 1+0i 2.绘制相应的图形,即将以上结果可视化. 四、实验结果分析 观测现象: 由结果可以观察到误差是不可避免的,并且扰动的减小不能使得所有根除的误差减小,还可以观察到处的误差总是最大,而靠近处的误差较小. 误差分析: 原问题数学上描述为: 的求解问题,也可等价为: 那么要考虑的问题是对方程的某个解的影响,我们不妨将看成的函数,从而可以求得此问题的条件数为: 可以求得20个根处的条件数为: 19次方系数变化对根的影响问题条件数 X=1,…,5 X=6,…,10 X=11,…,15 X=16,…,20 .172633e-14 2038.69 .886669e9 .315986e11 .859841e-8 76273.1 .347381e10 .235250e11 .114368e-3 .156701e7 .978210e10 .116152e11 .114955 .195843e8 .199945e11 .341541e10 25.5252 .159475e9 .296669e11 .452548e9 可见而靠近处的条件数非常小,而处的条件数最大.这与我们的观测结果是一致的. 由于条件数太大所以此问题是病态问题. 实验Y 用Jacobi法求对称矩阵的特征值及特征向量 一、实验内容 已下列矩阵为例,求对称矩阵的全部特征值及特征向量: 1. 2. 3. 二、方法步骤 1.在的非主对角线元素中,找出按模最大的元素; 2.计算平面旋转矩阵,其中的及由计算; 3.计算,(的初始值取单位阵); 4.如果(其中为的元素),则停止计算,所求特征值为:,特征向量:即得第列;否则令,重复以上各步. 三、实验结果 讨论的矩阵为: A = 5 2 4 2 5 3 4 3 1 要求误差为:err=1e-005 迭代次数:7 第1个特征值为:9.8057 相应的特征向量为:0.65053 0.57817 0.49248 第2个特征值为:3.0604 相应的特征向量为:-