《数学建模与数学实验》上机报告

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

《数学建模与数学实验》上机报告(第 1 次)

一、上机训练目的、题目或内容(简述综述)等

题目一:数学软件(MathType5.2、MATLAB 、Maple、Mathematica4.0、LINGO8.0)安装调试;基本命令使用(变量赋值、定义函数、过程控制、绘图命令、拟合、线性规划、非线性规划);高等数学实验(绘图,极限,求导,积分,解微分方程);线性代数实验(矩阵基本运算,线性方程组求解,解超定方程组,优化命令)。调试运行给定的两个程序:

题目二:

1、以两种方式打开MATLAB 工作窗口,进入MATLAB 6.0 的工作环境,并尝试用不同的方式退出。(这个在报告里面说明方法就可以)

2、尝试、熟悉MATLAB 6.0 的各栏菜单以及各个工具栏的功能。(自己掌握,报告里面就不写了)

3、绘制函数y=cos(5x+2)/sin(3x+1) 的图像,并求解当x=2 时的函数值。

4、练习并熟练掌握MATLAB 的帮助命令,学会利用MATLAB 的帮助信息。

5、求矩阵A=的行列式、逆的特征根;B=,解方程BX=

6、两个矩阵A=B=将矩阵改为3行3列的矩阵,作加、减、乘和除(左

除,右除)运算,同事运用数组运算法则进行运算,比较二者计算结果有何异同。

二、数学模型或求解分析或算法描述程序命令图形等

题目一:

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)

题目二:

3. x=2;

y=cos(5*x+2)./sin(3*x+1)

x=[-10:0.01:10];

y=cos(5*x+2)./sin(3*x+1);

plot(x,y)

5. A=[8 1 -5 1;9 -3 0 -6;-5 2 -1 2;0 4 -7 6]

det(A)

inv(A)

[V,D]=eig(A)

6.B=[1 3 5;2 4 6; 7 8 9;1 0 1]

b=[1;2;3;1]

X=B\b

三、结果或结论或截图

四、结果分析讨论或评价、推广、小结等

运用数学软件求解相关数学问题更加简便,省去了繁琐的计算过程,所求的结果也准确可靠。可以广泛推广到各个领域应用。

《数学建模与数学实验》上机报告(第 2 次)

一、上机训练目的、题目或内容(简述综述)等

题目三:

1、创建一个表达式Ζ= ,并求当x=1,y=2是的z值。

2、按照的步长间隔绘制函数在0时的曲线。

3、用曲面画图命令surf表现函数的图像。

4、作函数的图标,并用所有的修饰命令进行修饰。

5、计算积分:

(1)(2)

6、求函数的导数

(1)(2)

题目四:

1、用while循环写一个1~200之间的整数和。

2、编写一个M文件,画出下列分段函数所表示的曲面。

=

二、数学模型或求解分析或算法描述程序命令图形等

题目三:

1、 x=1;

y=2;

z=(sqrt(4*x^2+1)+0.5457*exp(-0.75*x^2-3.75*y^2-1.5*x))/(2*sin(3*y)-1)

2、 x=0:0.1:1;y=x.*exp(-x);plot(x,y);

3、[x,y]=meshgrid(-1:0.02:1);

z=x.^2+y.^2;

surf(x,y,z)

axis equal

4、y=tan(sin(x))-sin(tan(x));

plot(x,y)

xlabel('x')

ylabel('y')

title('这是一个图','fontsize',20);

5、(1)syms x

y=int((exp(x)-1)^4)*exp(x);

yy=int(y,x,0,1)

(2)syms x

y=(2*x)/(x^2+1);

yy=int(y,x,0,inf)

6、(1)syms x

f=sym(2^(x/log10(x)));

diff(f)

(2)syms x a

f=(1/2*a)*(log10(sqrt(a^2+x^2)/a+x)-a/(a+x));

diff(f)

题目四:

start=1

en=200

sum=0

while start<=en

sum = start + sum;

start = start + 1;

end

sum

三、结果或结论或截图

题目三:

1、

2、

3、

相关文档
最新文档