中南大学matlab课后答案-第六章
发布时间：2020-12-16
实验指导 二 1，
（1）>> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];B=[,,]';x=A\B （2）>> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];B=[,,]';x=A\B （3）>> cond(A 2，
（1）建立函数文件，命令如下; Function fx=funx(x fx=x^41+x^3+1; $
调用fzero函数求根，命令如下; >> z=fzero(@funx,-1 (2 建立函数文件，命令如下; Function fx=sin(x fx=x-(sin(x/x;
调用fzero函数求根，命令如下; >> z=fzero(@sin,
(3 建立函数文件，命令如下; function q=myfun(p
x=p(1;
、
y=p(2; z=p(3;
q(1=sin(x+y^2+log(z-7; q(2=3*x+2^y-z^3+1; q(3=x+y+z-5;
调用fsolve函数求根，命令如下;
>> options=optimset('Display','off';x=fsolve(@myfun,[1,1,1]',options 3,
(1 建立函数文件,命令如下; function yp=funt(t,y
>
yp=-y*+sin(10*t;
求微分方程，程序如下：
>> t0=0;tf=5;y0=1;[t,y]=ode23(@funt,[t0,tf],y0 (2 建立函数文件,命令如下; function yp=funr(t,y
yp=cos(t-(y/(1+t^2;
求微分方程，程序如下：
>> t0=0;tf=5;y0=1;[t,y]=ode23(@funr,[t0,tf],y0 4,
建立函数文件命令如下:
{
function fx=mymax(x fx=-1*(1+x^2/(1+x^4;
求最大值，程序如下:
>> [x,y]=fminbnd(@mymax,0,2 5,
编写目标函数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,0,0;,,1,0;,,,1]; b=[400,440,484,]; 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=A\b ！
矩阵分解：>> A=[2,3,5;3,7,4;1,-7,1];b=[10,3,5]';[L,U]=lu(A;x=U\(L\b
(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=A\b
矩阵分解：>> 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\(L\b 2,
(1建立函数文件，命令如下：
function fx=fun1(x fx=3*x+sin(x-exp(x;
调用fzero函数求根，命令如下; >> y=fzero(@fun1,
、
(2 建立函数文件，命令如下： function fx=fun2(x
fx=1-(1/x+5;
调用fzero函数求根，命令如下; >> y=fzero(@fun2,1
(3 建立函数文件，命令如下;
function q=fun3(p x=p(1; y=p(2; q(1=x^2+y^2-9;
'
q(2=x+y-1;
调用fsolve函数求根，命令如下;
>> options=optimset('Display','off';x=fsolve(@fun3,[1,1]',options 3,
(1建立函数文件，命令如下：
function ydot=fun5(t,y
ydot(1=(2-3*y(2-2*t*y(1/(1+t^2; ydot(2=y(1; ydot=ydot';
求解微分方程，命令如下： '
>> t0=0;tf=5;x0=[0,1];[t,y]=ode45(@fun5,[t0,tf],x0;[t,y] (2 (1建立函数文件，命令如下：
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,y] 4, 】
建立函数文件命令如下:
function fx=max(x fx=-1*(sin(x+cos(x^2;
求最大值，程序如下:
>> [x,y]=fminbnd(@mymax,0,pi 5,
编写目标函数
