数学建模 实验报告(免积分).docVIP

数学建模实验报告 1．用图形放大法求解方程 x sin(x) = 1。并观察该方程有多少个根。 解：方程变换为y=x.*sin(x)-1。代码：x=-100:0.01:100; y=x.*sin(x)-1; plot(x,y); 结果： 注：从图中可以看出有无数个根。再取一定范围： x=-4:0.01:4; y=x.*sin(x)-1; plot(x,y); grid on; %画坐标图 再放大， x=-1.2:0.01:1.2; y=x.*sin(x)-1; plot(x,y); grid on; %画坐标图 2．将方程x5 +5x3- 2x + 1 = 0 改写成各种等价的形式进行迭代，观察迭代是否收敛，并给出解释。用其他方法求解方程得到方程的全部根。 解：方程可以变为三种不同形式迭代， （1）（a）x=(2*x-5*x^3-1)^(1/5) （b）y=(((-y)^5+2*y-1)/5)^(1/5) （c） z=(z^5+5*z^3+1)/2 程序：x=1;y=1;z=1; for k=1:100 x=(2*x-5*x^3-1)^(1/5); y=(((-y)^5+2*y-1)/5)^(1/5); z=(z^5+5*z^3+1)/2 x,y,z; end 最后结果为：x = 2.0162 - 0.8223i y = 0.5303 + 0.1491i z = Inf 函数x,y,z分别对应方程（a）（b ）（c），从结果可以看出方程（ c）不收敛，结果趋于无穷大，方程(a)(b)收敛。 （2）先图形放大法： 程序；x=-1:0.01:1; y=x.^5 +5*x.^3-2*x+1; plot(x,y) grid on, 可以看出取值在（-0.8，-0.7），取初值 n=0,1,2,3……… 程序：x=-077;y=-0.77;z=-0.77; for k=1:30 x=(-4*x^5-10*x^3+1)/(2-5*x^4-15*x^2); y=(2*y^6+4*y^2-3*y)/(5*y^3+3*y^5+2*y-2); z=(8*z^2-2*z)/(z^5+5*z^3+6*z-1); x,y,z; end 迭代最后结果为： x = -0.7685 y = -0.7685 即方程根为-0.7685. （3）