数学实验报告6.docVIP

数学实验报告 实验序号： 实验16 日期： 年 月 日 班级 姓名 学号 实验 名称 方程和方程组的求解 问题 背景 描述 在MATLAB中，由函数solve( )、null( )、fsolve( )，fzero等来解决线性方程（组）和非线性方程（组）的求解问题，其具体格式如下： X=solve(‘eqn1’,’eqn2’,…,’eqnN’,’var1’,’var2’,…,’varN’) ???? X=fsolve(fun,x0,options) 函数solve用来解符号方程、方程组，以及超越方程，如三角函数方程等非线性方程。参数’eqnN’为方程组中的第N个方程，’varN’则是第N个变量。 函数null(A)则用来解线性方程组AX=O的基础解系，实际是求系数矩阵A的零空间，在null函数中可加入参数’r’，表示有理基。通过求系数矩阵的秩和增广矩阵的秩，可以判定方程组是否有解，以及是否需要求基础解系。 实验 目的 熟悉MATLAB软件中关于求解方程和方程组的各种命令，掌握利用MATLAB软件进行线性方程组、非线性方程、非线性方程组的求解。 实验所用软件及版本：MATLAB7.0.1 实验 过程 记录 1、建立M文件sy161.m，并在文件内输入以下命令： clear format rat A=sym([1,1,3,-1;0,2,-1,1;1,1,2,2;1,-1,1,-1]); b=sym([-2;1;4;0]); B=[A,b]; n=length(b) RA=rank(eval(A)) RB=rank(eval(B)) if(RA==RBRA==n) ? X=eval(A\b)??? elseif(RA==RBRAn) ??? C=eval(A\b) ??? D=null(eval(A),r)?? ??? syms k1 k2 ??? X=k1*D(:,1)+k2*D(:,2)+C ??? else ??? fprintf(No Solution for the Equations) end 并在命令窗口运行文件sy161.m 2、在命令窗口输入 X=solve(x-exp(-x)=0,x) 并按回车执行命令。 3、在命令窗口输入 fplot(5*x^2*sin(x)-exp(-x),[0,10]) 按回车执行命令，画出图形发现方程的共有四个根，分别在 0、3、6、9附近。 在命令窗口输入 x1=fsolve(5*x^2*sin(x)-exp(-x),0) 并回车； 在命令窗口输入 x2=fsolve(5*x^2*sin(x)-exp(-x),3) 并回车； 在命令窗口输入 x3=fsolve(5*x^2*sin(x)-exp(-x),6) 并回车； 在命令窗口输入 x4=fsolve(5*x^2*sin(x)-exp(-x),9) 并回车。 4、先建立一个存放函数的m文件(sy164.m)： function y=sy164(x) y(1)=x(1)-0.7*sin(x(1))-0.2*cos(x(2)) y(2)=x(2)-0.7*cos(x(1))-0.2*sin(x(2)) 接着，建立在命令窗口输入: format short x0=[0.1,0.1] fsolve(@sy164,x0) 实验 结果 报告 与 实验 总结 1、n = 4 RA = 4 RB = 4 X = 13/8 -5/8 -3/8 15/8 可知该方程组有唯一解，唯一解为X= （13/8，-5/8，-3/8，15/8）T 2、方程的解为： X =lambertw(1) 3、函数 y=5*x^2*sin(x)-exp(-x)在[0,10]内的图像为： 方程组的四个根为： x1 =0.5018 x2=3.1407 x3 =6.2832 x4=9.4248 4、方程组的解为x=0.4442 y=0.7715。 思考 与 深入 教师 评语