MATLAB考核期末(含答案).docVIP

一、微分方程求解。（15分） 设，，，用ode45方法求微分方程的解。要求: 1）仿真时间t=30秒 2）结果绘制在一张图中，包括曲线，一阶曲线，二阶曲线，三阶曲线 3）用图例命令分别说明四条曲线为“”，“”，“” ，“” 4）定义横坐标为“时间”，纵坐标为“输出”，图形标题名称为“微分方程的解” 解：系统方程为 , 这是一个单变量三阶常微分方程。将上式写成一个一阶方程组的形式，这是函数ode45调用规定的格式。 令： 函数文件程序： function ydot=myfun1(t,y) ydot=[y(2);y(3);1-8*y(1)-2*y(3)-4*y(2)]; 主文件程序： t=[0 30]; y0=[0;1;0]; [tt,yy]=ode45(@myfun1,t,y0); y=(1-yy(:,3)-2*yy(:,2)-4*yy(:,1))/8; plot(tt,y,r,tt,yy(:,1),k,tt,yy(:,2),-g,tt,yy(:,3),-.b); legend(y-t,yˊ-t,yˊˊ-t,yˊˊˊ-t) title(微分方程的解) xlabel(时间) ylabel(输出) 二、数值计算，编程完成以下各题（共20分，每小题5分） 1、画出函数在区间[3, 5]的图形，求出该函数在区间[3, 5]中的最小值点和函数的最小值. 解：程序如下 x=3:0.05:5; y=(sin(5*x).^2).*exp(0.05*x.^2)-5*(x.^5).*cos(1.5*x)+1.5*abs(x+5.5)+x.^2.5; plot(x,y); fn=(sin(5*x).^2).*exp(0.05*x.^2)-5*(x.^5).*cos(1.5*x)+1.5*abs(x+5.5)+x.^2.5; [xmin,fmin]=fminbnd(fn,3,5) hold on plot(xmin,fmin,r*); Xlabel(x) Ylabel(f(x)) 经过运行后得到的图像截图如下： 运行后的最小值点=4.6051，= -8.3383e+003。 2、画出函数在[1，3]区间的图形，并用编程求解该非线性方程的一个根，设初始点为. 解：在MATLAB中提供了一个fzero函数，可以用来求单变量非线性方程的根，其调用格式为：z=fzero(fname,x0,tol,trace),其中fname是待求根的函数文件名，x0为搜索的起点。一个函数可能有多个根，但是fzero 函数只给出离x0最近的那个根。Tol控制结果的精度，缺省时取tol=esp，trace指定迭代信息是否在运算中显示，为1时显示，为0时不显示，缺省时取trace=0.下面就用fzero函数来求解。 函数文件： function y=fun2_2(x) y=((cos(x)).^2).*exp(-0.3*x)-2.5*abs(x); 主程序： clc xdomain=[1 3]; fplot(@fun2_2,xdomain); x=fsolve(@fun2_2,2,optimset(Display,off)) 运行程序后,函数在[1，3]区间的图形如下： 运行后求得该方程的一个根为z=0.3256。 3、求矩阵，使其满足方程，设初始解向量为. 解：设，则， 又因为，所以写成方程组的形式如下： 显然，这是一个非线性方程组，下面用fsolve函数求解器编程求解。 %在新建中写出函数文件fun2_3.m并保存。 function y=fun2_3(x) y=[x(1).^3+2*x(1).*x(2).*x(3)+x(2).*x(3).*x(4)-1,(x(1).^2).*x(2)+(x(2).^2).*x(3)+ x(1).*x(2).*x(4)+x(2).*x(4).^2-2,... (x(1).^2).*x(3)+x(1).*x(3).*x(4)+x(2).*x(3).^2+x(3).*x(4).^2-3,x(1).*x(2).*x(3)+2*x(2).*x(3).*x(4)+x(4).^3-4]; %主程序矩阵求解juzhen.m x0=[1 1 1 1]; fsolve(@fun2_3,x0) %给定初值，求解器求解。 运行后的结果为： ans =-0.1291 0.8602 1.2903 1.1612 即是a=-0.1291 ,b=0.8602,c=1.2903,d=1.1612. 所以所求的矩阵为 。 4、已知非线性方程组如下，编程求方程组的解，设初始点为[1 0.5 -1]. 解：%在新建中建立函数文件fun2_4