重庆大学数学模型数学实验作业二.docVIP

开课学院、实验室：数统学院 实验时间 ： 课程 名称 数学实验 实验项目 名 称 方程求解 实验项目类型 验证 演示 综合 设计 其他 指导 教师 肖剑 成 绩 实验目的 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。 实验内容 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。 基础实验 一、问题重述 1．用图形放大法求解方程x sin(x) = 1. 并观察该方程有多少个根。 2．将方程x5 +5x3- 2x + 1 = 0 改写成各种等价的形式进行迭代，观察迭代是否收敛，并给出解释。 3．求解下列方程组 直接使用MATLAB命令：solve()和fsolve()对方程组求解。 4．编写用二分法求方程根的函数M文件。 5. 设非线性方程组为 其中已知，随机产生数据后，用fsolve解这个方程组。 6. 使用fsolve计算方程组的解时，为验证初值是否对解有影响，采用随机产生的100组随机数作为初始值，依次进行求解。 二、实验过程 1.作图得y=x sin(x)-1得其在区间[-10,10]的图像，可知原方程在[-10,10]上有8个根。如图： 下图是运用放大法求其在[0,2]上的解，可知解大概为x=1.114.如图： 编程Untitled1（见附件）。 2.画图y=x5 +5x3- 2x + 1,可得零点大概位于区间[-2,0]。即存在收敛的迭代函数，将y=0变形得x= (1/2)x5 +(5/2)x3 + 1/2,这是y1,以及x=2/x^3-5/x-1/x^4,这是y2，利用加速迭代使y1变形后得到x=(-4x5-10x3 + 1)/(2-5x4-15x2),这是y3.用MATLAB编程Untitled2（见附件），输出为（仅罗列等间距的10） Untitled2= c =2.0000 -2.5000 -0.8333 c =-0.6109 -0.8050 -0.7685 c =-0.7685 -0.7685 -0.7685 c =-0.7685 -0.7685 -0.7685 c =-0.7685 -0.7685 -0.7685 c =-0.7685 -0.7685 -0.7685 c =-0.7685 -0.7685 -0.7685 c =-0.7685 -0.7685 -0.7685 c =-0.7685 -0.7685 -0.7685 c =-0.7685 -0.7685 -0.7685 故迭代求得解为x= -0.7685,且y1、y2、y3均收敛，但y3收敛最快。 3.（1）直接利用solve()编程Untitled3（见附件）得 x1=0.5671 x2=0.5671 (2) 建立方程组的M-函数文件san.m（见附件）,然后直接利用MATLAB运行如下程序 y=fsolve(‘san’,[1,1,1],1) 得到如下结果 y = 0.0000 1.5492 0.0000 4.编辑函数m文件mysolve.m（见附件）。当然，此法也有其局限性，待定函数必须具有确定的显式表达式，同时得保证函数在给定区间有零点，且在端点函数值以及区间中点函数值不为0的情况下区间（a,b）只能有奇数个零点，且只能求出一个零点；若端点（区间中点）函数值为0，则输出端点（区间中点）值；且这个m文件只能求出一个零点。所以在命令窗口输入指令时最好先用MATLAB作图大致确定函数零点情况。验证：在命令窗口输入X=mysolve(@(x)x^2-3,0,3)，运行得到The equation root is X =1.7320。 显然X^2=3,是x^2-3=0的一个零点。 5．编辑函数-M文件h.m(见附件)，在MATLAB命令行窗口键入Untitled5.m（见附件），得到解。由于方程系数随机产生，故每次运行结果都不同。 6. 编辑函数-M文件liu.m(见附件)，在MATLAB命令行窗口键入Untitled6.m（见附件），得到解如下（仅列出部分）： a =5.6357 7.98