中南大学matlab课后答案-第六章

实验指导二1，（1） A=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6;B=0.95,0.67,0.52;x=AB（2） A=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6;B=0.95,0.67,0.53;x=AB（3） cond(A)2，（1）建立函数文件funx.m，命令如下; Function fx=funx(x) fx=x41+x3+1;调用fzero函数求根，命令如下; z=fzero(funx,-1)(2) 建立函数文件sin.m，命令如下; Function fx=sin(x) fx=x-(sin(x)/x;调用fzero函数求根，命令如下; z=fzero(sin,0.5)(3) 建立函数文件myfun.m，命令如下; function q=myfun(p)x=p(1);y=p(2);z=p(3);q(1)=sin(x)+y2+log(z)-7;q(2)=3*x+2y-z3+1;q(3)=x+y+z-5;调用fsolve函数求根，命令如下; options=optimset(Display,off);x=fsolve(myfun,1,1,1,options)3,(1) 建立函数文件funt.m,命令如下; function yp=funt(t,y)yp=-y*(1.2+sin(10*t); 求微分方程，程序如下： t0=0;tf=5;y0=1;t,y=ode23(funt,t0,tf,y0)(2) 建立函数文件funr.m,命令如下; function yp=funr(t,y)yp=cos(t)-(y)/(1+t2);求微分方程，程序如下： t0=0;tf=5;y0=1;t,y=ode23(funr,t0,tf,y0)4,建立函数文件mymax.m命令如下:function fx=mymax(x)fx=-1*(1+x2)/(1+x4);求最大值，程序如下: x,y=fminbnd(mymax,0,2)5,编写目标函数M文件fop.m，命令如下:function f=fop(x)f=-1*(x(1)(1/2)+x(2)(1/2)+x(3)(1/2)+x(4)(1/2);设定约束条件，并调用fmincon函数求解此约束最优化问题，程序如下：x0=200,200,200,200;A=1,0,0,0;1.21,1,0,0;1.331,1.21,1,0;1.4641,1.331,1.21,1;b=400,440,484,532.4;Ib=0,0,0,0;options=optimset(Display,off);x,y=fmincon(fop,x0,A,b,Ib,options)思考练习1,(1)矩阵求逆： A=2,3,5;3,7,4;1,-7,1;b=10,3,5;x=inv(A)*b矩阵除法： A=2,3,5;3,7,4;1,-7,1;b=10,3,5;x=Ab矩阵分解： A=2,3,5;3,7,4;1,-7,1;b=10,3,5;L,U=lu(A);x=U(Lb)(2)矩阵求逆： A=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2;b=-4,13,1,11;x=inv(A)*b 矩阵除法： A=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2;b=-4,13,1,11;x=Ab 矩阵分解： A=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2;b=-4,13,1,11; L,U=lu(A);x=U(Lb)2,(1)建立函数文件fun1.m，命令如下：function fx=fun1(x)fx=3*x+sin(x)-exp(x);调用fzero函数求根，命令如下; y=fzero(fun1,1.5)(2) 建立函数文件fun2.m，命令如下： function fx=fun2(x)fx=1-(1/x)+5;调用fzero函数求根，命令如下; y=fzero(fun2,1)(3) 建立函数文件myfun3.m，命令如下;function q=fun3(p)x=p(1);y=p(2);q(1)=x2+y2-9;q(2)=x+y-1;调用fsolve函数求根，命令如下; options=optimset(Display,off);x=fsolve(fun3,1,1,options)3,(1)建立函数文件fun5.m，命令如下：function ydot=fun5(t,y)ydot(1)=(2-3*y(2)-2*t*y(1)/(1+t2);ydot(2)=y(1);ydot=ydot;求解微分方程，命令如下： t0=0;tf=5;x0=0,1;t,y=ode45(fun5,t0,tf,x0);t,y(2) (1)建立函数文件fun5.m，命令如下：function ydot=fun6(t,y)ydot(1)=cos(t)+(5*y(1)*cos(2*t)/(t+1)2-y(2)-y(3)/(3+sin(t);ydot(2)=y(1);ydot(3)=y(2);ydot=ydot;求解微分方程，命令如下： t0=0;tf=5;x0=0,1;t,y=ode45(fun5,t0,tf,x0);t,y4,建立函数文件max.m命令如下:function fx=max(x)fx=-1*(sin(x)+cos(x2);求最大值，程序如下: x,y=fminbnd(mymax,0,pi)5,编写目标函数M文件fop.m，命令如下:function f=topm(x)f=-1*x*(3-2*