数学实验实验二答案

实验05 数学规划模型㈡（2学时）

（第4章 数学规划模型）

1.（求解）汽车厂生产计划（LP，整数规划IP）p101~102

(1) (LP)在模型窗口中输入以下线性规划模型

max z = 2x1 + 3x2 + 4x3 s.t. 1.5x1 + 3x2 + 5x3 ≤ 600

280x1 + 250x2 + 400x3 ≤ 60000

x1, x2, x3 ≥ 0

并求解模型。

★(1) 给出输入模型和求解结果（见[101]）：

model: TITLE汽车厂生产计划（LP）; !文件名：p101.lg4; max=2*x1+3*x2+4*x3; 1.5*x1+3*x2+5*x3<600; 280*x1+250*x2+400*x3<60000; end (2) (IP)在模型窗口中输入以下整数规划模型

max z = 2x1 + 3x2 + 4x3 s.t. 1.5x1 + 3x2 + 5x3 ≤ 600

280x1 + 250x2 + 400x3 ≤ 60000

x1, x2, x3均为非负整数

1

并求解模型。

LINGO函数@gin见提示。

★(2) 给出输入模型和求解结果（见[102]模型、结果）：

model: TITLE汽车厂生产计划（IP）; !文件名：p102.lg4; max=2*x1+3*x2+4*x3; 1.5*x1+3*x2+5*x3

实验05 数学规划模型㈡（2学时）

（第4章 数学规划模型）

1.（求解）汽车厂生产计划（LP，整数规划IP）p101~102

(1) (LP)在模型窗口中输入以下线性规划模型

max z = 2x1 + 3x2 + 4x3 s.t. 1.5x1 + 3x2 + 5x3 ≤ 600

280x1 + 250x2 + 400x3 ≤ 60000

x1, x2, x3 ≥ 0

并求解模型。

★(1) 给出输入模型和求解结果（见[101]）：

model: TITLE汽车厂生产计划（LP）; !文件名：p101.lg4; max=2*x1+3*x2+4*x3; 1.5*x1+3*x2+5*x3<600; 280*x1+250*x2+400*x3<60000; end (2) (IP)在模型窗口中输入以下整数规划模型

max z = 2x1 + 3x2 + 4x3 s.t. 1.5x1 + 3x2 + 5x3 ≤ 600

280x1 + 250x2 + 400x3 ≤ 60000

x1, x2, x3均为非负整数 并求解模型。

LINGO函数@gin见提示。

★(2) 给出输入模型和求解结果（见[102]模型、结果）：

model: TITLE汽车厂生产计划（IP）; 1

!文件名：p102.lg4; max=2*x1+3*x2+4*x3; 1.5*x1+3*x

359242463.docx数据描述（基本数据类型及运算符）答案

实验二数据描述（基本数据类型及运算符）答案

编程及调试实例2-1改正错误后的程序 #include void main( ) {

int celsius, fahr; fahr = 100;

celsius = 5 * (fahr - 32) / 9;

printf(\ %d, celsius = %d\\n\ }

编程及调试练习2-1：源程序 （1）

#include void main( ) {

int celsius, fahr;

fahr = 150; /*只改变fahr的值*/

celsius = 5 * (fahr - 32) / 9; /*celsius = 5 * (fahr - 32) / 9与上面的实例完全一样 */ printf(\ }

运行程序后，输出： fahr = 150, celsius = 65

（2）

#include void main( ) {

int celsius, fahr;

fahr = 150;

实 验 六

实验项目名称：优化

实验时间： 2010-5-26、2010-6-2、2010-6-9 实验地点：理学实验楼525 实验目的：

1、掌握Matlab优化工具箱的基本用法，利用优化工具包解线性规划和非线性规划的问题；对不同算法作初步分析、比较； 2、练习实际问题的非线性最小二乘拟合。

3、初步了解实际问题的线性规划和非线性规划的模型建立。 实验内容： 1、

教材205页习题2;求解min ex1(4x12+2x22+4x1x2+2x2+1)，初值（-1，1），对不同算法的结果分析，比较。

编写M-文件 fun.m: function f=fun(x)

f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1); 输入M文件wliti.m如下： x0=[-1,1];

x=fminunc('fun',x0); y=fun(x) 2、

教材205页习题4;有一组数据（ti,yi），i=1,…,33,其中ti=10(i-1),yi由表给出。现要用这组数据拟合函数f(x,t)=x1+x2e-x4t+x3e-x5t中的参数x,初值可选为(0,5,1.5,-1,0.01,0.02)，用GN和LM

实验二 卷积实验

一、实验目的

1．熟悉并验证卷积的性质

2．利用卷积生成新的波形，建立波形间的联系 3．验证卷积定理

二、实验内容`

信号的卷积是针对时域信号处理的一种分析方法。信号的卷积一般用于求取信号通过某系统后的响应。在信号与系统中，我们通常求取某系统的单位冲激响应，所求得的h[n]可作为系统的时域表征。任意系统的系统响应可用卷积的方法求得：

y[n]?x[n]?h[n]

根据实验原理提示编写以下程序：

（1）MATLAB提供了一个内部函数conv()来计算两个有限长序列的卷积。conv()函数假定两个序列都从n?0开始。

给出序列x=[3, 11, 7, 0, -1, 4, 2]和h=[2, 3, 0, -5, 2, 1]，求两者的卷积y。

ans =

6 31 47 6 -51 -5 41 18 -22 -3 8 2

将函数conv()稍加扩展为函数conv_m()，它可以对n从任意取值开始的序列求卷积。

格式如下：

function [y, ny]=conv_m(x, nx, h, nh) % 信号处理的改进卷积程序 % [y, ny]=conv_m(x, nx, h, n

教 师 教 案

( 2008 —2009 学年第 二 学期 )

课 程 名 称：数学实验 授 课 学 时：32 授 课 班 级：2007级 任 课 教 师：钟尔杰 教 师 职 称：副教授 教师所在学院：应用数学学院

电子科技大学

数学实验教案 第 1 页 共 11 页

课程名称 课程编号 必修课 课程类型 选修课 授课方式 是否采用多媒体 学时分配 教材 数学实验 21100520 授课专业班级 年级 2007级 修课人数 二年级 140人 学科基础课 (? ) ； 专业基础课 ( ) ； 专业课 ( ) 专 选 课 ( ) ； 公选课 ( ) 考核方式 是否采用双语 考 试 (? ) 考 查 ( ) 否 理论课 (? ) ；实践课 (? ) 是 课堂讲授 20 学时； 实践课 12 学时 名称 数学实验讲义 1、数学实验 2、数学实验 3、线性代数与空间解析几何 4、微积分（上、下） 5、数学实验简明教程 6、实用数值计算方法 作者 数学学院 李尚志 等 谢云荪 等 黄廷祝 等 电子科大数学系 电子科大数学系 李丽，王振领 张志涌等 出版社及出版时间 讲义2009年 高等教育出版社

一、选择题

1、空间曲线绘图命令是：（ ）。

(A) plot2 (B) plot3 (C) surf (D) plot 2、在matlab中，命令hold off 表示是：（ ）。

(A) 在图中消隐分隔线。 (B) 使系统处于可放大状态。 (C) 保留当前图形和它的轴，使其后图形放在当前图形上。 (D) 此后图形指令运作将抹掉当前窗口中的旧图形，然后画上新图形。

3、生成三维表面图的命令是：（ ）。

(A) surf (B) meshc (C) mesh (D) meshz 4、下面哪个命令的功能是是系统处于可放大状态（ ）。 (A) grid on (B) hold on (C)zoom on (D) hidden on 5、在MATLAB的若干通用操作指令中，清除内存中的变量的是（ ） （A）clear （B） cla （C） clf （D）clc 6、下列表示可以作为MATLAB的变量的是( ) （A） abcd-2 （B）xyz_2# （C）@h （D）X_1_a

一、选择题

1、空间曲线绘图命令是：（ ）。

(A) plot2 (B) plot3 (C) surf (D) plot 2、在matlab中，命令hold off 表示是：（ ）。

(A) 在图中消隐分隔线。 (B) 使系统处于可放大状态。 (C) 保留当前图形和它的轴，使其后图形放在当前图形上。 (D) 此后图形指令运作将抹掉当前窗口中的旧图形，然后画上新图形。

3、生成三维表面图的命令是：（ ）。

(A) surf (B) meshc (C) mesh (D) meshz 4、下面哪个命令的功能是是系统处于可放大状态（ ）。 (A) grid on (B) hold on (C)zoom on (D) hidden on 5、在MATLAB的若干通用操作指令中，清除内存中的变量的是（ ） （A）clear （B） cla （C） clf （D）clc 6、下列表示可以作为MATLAB的变量的是( ) （A） abcd-2 （B）xyz_2# （C）@h （D）X_1_a

浙 江 海 洋 学 院

数 学 实 验 报 告 二

学生姓名： 学 号： 指导教师： 实验时间: 报告评分：

一、实验室名称：实验室

二、实验项目名称：抛射曲线实验 三、实验原理：

炮弹飞行的简单数学模型为抛射曲线，其参数方程为

?x?v0cos??t??12 y?v0sin??t?gt?2?其中v0为初始速度，? 为发射角，g 是重力加速度（约为 9.8米/秒2）。

四、实验目的：

熟悉MATLAB的集成环境，掌握MATLAB的基本运算和数学表达式计算方法及其画图。

五、实验内容：

对已知射程数据，为使炮弹击中目标需计算发射角；反之对已知的发射角数据，也需要计算射程。设v0 = 515米/秒，对不同目标计算射程、发射角和炮弹飞行时间。 六、实验步骤及操作：

弹道数据包括射程、发射角及炮弹飞行时间。令参数方程中y 的表达式为零可确定炮弹飞行时间 t(? )，代入 x 的表达式可导出射程计算公式。

2v0sin? g22v0sin?v0（2） 射程计算公式：x1?v0cos???sin2?

《数学实验》实验报告

2010-2011学年第1学期

学生姓名： 学 号： 院 部： 数理学院 专 业： 班 级： 任课教师：

实验报告1

实验目的：

熟悉Mathematica软件包的使用。

实验内容：

1、用两种方式编写如下自定义函数，并求其导数f’(x)在x=-2.0,x=1.0,x=5.0处的值

2、分别用

Plot3D, ParametricPlot3D

函数画出x2?y2?z2?1（0?x?1,0?y?1）的图像。

3、用Mathematica实现一个四人追逐问题，给出结果并划出追逐路线（如下图）。

实验要求：

1、撰写实验报告；

2、写出试验过程中所使用的Mathematica程序或语句和计算结果；

实验步骤：

1、 用两种方式编写如下自定义函数，并求其导数f’(x)及f’(x)在

x=-2.0,x=1.0,x=5.0处的值

方法1

f[x_]:=E^x Sin[x]/;x<=0;f[x_]:=Cos[x]/;0E;

g[x_]=D[f