数学实验强课件--复习2.pptVIP

课件 * 课件 第二章复习 习题解答 结束 课件 ?29? 结束 第二章习题 1、用图形放大法求解方程 x sin x = 1. 并观察 该方程有多少个根。 程序代码如下： clear x=-10:0.05:10; y=0; plot(x,y) hold on ezplot(x*sin(x)-1,[-10,10]) 课件 显示： ?30? 结束 由图观察，有无穷多个实根。 课件 ?31? 结束 2、将方程x5 +5x3– 2x + 1 = 0 改写成各种等价 的形式进行迭代，观察迭代是否收敛，并给出 解释。 clear x=[-1:.01:1]; y=x.^5+5*x.^3-2*x+1; plot(x,y) grid on f=solve(x^5+5*x^3-2*x+1=0) f1=double(f) 求精确解的程序代码如下： 课件 ?32? 结束 f = [-.76845316155665081805627362000827] [ -.16186872229712980224734834175246e-1 -2.3181992385209164662001197641760*i] [ -.16186872229712980224734834175246e-1 +2.3181992385209164662001197641760*i] [ .40041345300803838925287164417938 -.28601637969992331846741572632061*i] [ .40041345300803838925287164417938 +.28601637969992331846741572632061*i] f1 = -0.7685 -0.0162 - 2.3182i -0.0162 + 2.3182i 0.4004 - 0.2860i 0.4004 + 0.2860i 显示： 课件 ?33? 结束 由图观察，实根在 -0.8— -0.7 之间。 课件 ?34? 结束 clear x(1)=1;y(1)=1;z(1)=1; for k=1:100 x(k+1)=(x(k)^5+5*x (k) ^3+1)/2; y (k+1) =((2*y (k) -y (k) ^5-1)/5)^(1/3); z (k+1) =1/(2-z (k) ^4-5*z (k) ^2); end x, y,z 三种迭代格式求近似解 的程序代码如下： x(100) = Inf y(100) = 0.4004 + 0.2860i z (100)= -0.0667 结果： 第一、三种迭代格式不收 敛，仅第二种迭代格式收 敛于一个虚根。 课件 ?35? 结束 改进第一个迭代公式如下： 牛顿迭代公式如下： 这两个迭代公式完全相同！ 课件 ?35? 结束 牛顿迭代公式求近似解的程序代码如下： clear, x(1)=1; for k=1:100 x(k+1)=(4*x(k)^5+10*x (k) ^3-1)/ (5*x(k)^4+15*x (k) ^2-2); end x 显示： x = Columns 1 through 7 1.0000 0.7222 0.4946 0.1667 0.6035 0.3679 -3.8969 Columns 8 through 14 -3.0369 -2.3402 -1.7814 -1.3474 -1.0363 -0.8508 -0.7791 Columns 15 through 21 -0.7687 -0.7685 -0.7685 -0.7685 -0.7685 -0.7685 -0.7685 收敛于一个实根-0.7685，且速度很快。 课件 ?36? 结束 3、求下列方程的根 1）xsin(x)+cos(x)-sin(x)-2x=0, 在[-1，1]上的近似 解。 2）判定方程 x7+x5+1=0 有几个实根? 3) 找出函数 f(x)=x3-6x2-2x+12 的所有零点。 课件 ?37? 结束 clear ezplot(x*sin(x)+cos(x)-sin(x)-2*x ,[-1,1]) s=solve(x*sin(x)+cos(x)-sin(x)-2*x) vpa(s,7) y=numeric