《数学模型及数学软件》上机训9个题.docVIP

《数学模型及数学软件》上机报告（第1次）参考模板 专业班级姓名学号2007年9月日星期结果分析评价专业班级姓名学号2010年9月日星期结果分析评价MathType5.2、MATLAB 、Maple、Mathematica4.0、LINGO8.0）安装调试；基本命令使用（变量赋值、定义函数、过程控制、绘图命令、拟合、线性规划、非线性规划）；高等数学实验（绘图，极限，求导，积分，解微分方程）；线性代数实验（矩阵基本运算，线性方程组求解，解超定方程组，优化命令）。调试运行给定的两个程序： 1) c=[6,3,4]; A=[0,1,0]; b=[50]; Aeq=[1,1,1]; beq=[120]; vlb=[30,0,20]; vub=[]; [x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub) 2) function f=fun3(x); f=-x(1)-2*x(2)+(1/2)*x(1)^2+(1/2)*x(2)^2 x0=[1;1]; A=[2 3 ;1 4]; b=[6;5]; Aeq=[];beq=[]; VLB=[0;0]; VUB=[]; [x,fval]=fmincon(fun3,x0,A,b,Aeq,beq,VLB,VUB) 题目二 1、以两种方式打开 MATLAB 工作窗口，进入MATLAB 6.0 的工作环境，并尝试用不同的方式退出。 （这个在报告里面说明方法就可以） 2、尝试、熟悉 MATLAB 6.0 的各栏菜单以及各个工具栏的功能。（自己掌握，报告里面就不写了） 3、绘制函数y=cos(5x+2)/sin(3x+1) 的图像，并求解当 x=2 时的函数值。 4、练习并熟练掌握 MATLAB 的帮助命令，学会利用 MATLAB的帮助信息。 5、求矩阵A=的行列式、逆和特征根 6、两个矩阵 将矩阵A改为3行3列的矩阵，作加、减、乘和除（左除，右除）运算，同时运用数组运算法则进行运算，比较二者的计算结果有何异同。 题目三 1、创建一个表达式 ，并求当 x=1, y=2 时的z 值。 2、按照的步长间隔 绘制函数 在0≤x≤1时的曲线。 3、用曲面画图命令surf表现函数z=的图像 4、作函数y=tan(sin(x))-sin(tan(x))的图形，并用所有的修饰命令进行图像修饰。 5、计算积分： (1) (2) 6、求函数的导数 (1) (2) 题目四 1、用 while 循环求 1～200 之间的整数之和 2、编写一个 M 文件，画出下列分段函数所表示的曲面。 题目五 利用Matlab求解线性规划问题% min fx % s.t .(约束条件)： Ax=b % (等式约束条件)： Aeqx=beq % lb=x=ub linprog函数的调用格式如下： x=linprog(f,A,b) x=linprog(f,A,b,Aeq,beq) x=linprog(f,A,b,Aeq,beq,lb,ub) x=linprog(f,A,b,Aeq,beq,lb,ub,x0) x=linprog(f,A,b,Aeq,beq,lb,ub,x0,options) [x,fval]=linprog(…) [x, fval, exitflag]=linprog(…) [x, fval, exitflag, output]=linprog(…) [x, fval, exitflag, output, lambda]=linprog(…) 其中： x=linprog(f,A,b)返回值x为最优解向量。 x=linprog(f,A,b,Aeq,beq) 作有等式约束的问题。若没有不等式约束，则令A=[ ]、b=[ ] 。 x=linprog(f,A,b,Aeq,beq,lb,ub,x0,options) 中lb ,ub为变量x的下界和上界，x0为初值点，options为指定优化参数进行最小化。 Options的参数描述： Display显示水平。 选择’off’ 不显示输出；选择’Iter’显示每一 步迭代过程的输出；选择’final’ 显示最终结果。 MaxFunEvals 函数评价的最大允许次数 Maxiter 最大允许迭代次数 TolX x处的终止容限 [x,fval]=linprog(…) 左端 fval 返回解x处的目标函数值。 [x,fval,exitflag,outp