matlab讲义要点
前期知识：高等数学，机械原理，液压传动，控制原理
学习目标：了解基本原理、掌握基本应用、熟练使用matlab 、精通机械工程的建模和仿真 本课计划：课堂学习、课后练习、
任务：要看课堂笔记
引例：斐波纳西数列的递推公式为)2()1()(--+=n n n F F F ，通项表达式是什么？ 高等数学上册第48页。斐波纳西数列在优化设计和股票分析中有用处。
第1章 了解Matlab
本章要求：了解Matlab 的功能、组成
第1节 Matlab 的界面
1.命令窗口：输入命令和显示运行结果和寻求帮助的窗口。 第1个例子求不定积分
问题:在编写代码时不能输入或者不能显示汉字
解决办法：将use custom font 换成use desktop font
先定义一个符号变量x syms x 设2
11
)(x
x f +=
求不定积分
?dx x f )(
matlab 求解：int(1/(1+x .^2)) 详见l1_bdjf.m
最重要的问题一：工作路径 查看当前工作路径的命令是 pwd
matlab 默认的路径为安装好的目录下work ，为了需要我们需要改换路径。 改换路径的方法有：（1）采用DOS 命令 mkdir('根目录名称','新目录名称')
例：mkdir('d:\','mywork') 如果d 盘下没有mywork 即创建，如有就会给出警告。 进入新建文件夹 cd d:\mywork
（2）采用matlab 命令 editpath ，pathtool （3）通过matlab 界面 [file]菜单->set path
难点：我想将打开MATLAB 时的默认工作路径改为F:\Program\MATLAB\WorkSpace\ ,只需要在原来的默认路径(bin)下创建一个名为startup.m 的文件,内容为相对路径 cd ..\..\WorkSpace\ 或绝对路径 cd F:\Program\MATLAB\WorkSpace\
即可。再次打开MA TLAB 时便会自动执行startup.m 文件,将工作路径转至WorkSpace 下。
最重要的问题二：工作路径设置不能设置在有汉字的目录下或汉字文件夹，不支持汉字运算。
汉字用的是Unicode 编码一个字符占两个字节，字母用的是ASC Ⅱ编码，一个字母占一个字节。到2008版才能处理汉字。
初学者容易出现的错误就是把别人的程序拷在带汉字的文件夹下，运行出现错误。 如果我们已知某个文件名，但忘了在哪个文件夹下，可以用which 命令如which FUN what 命令：M-files in the current directory 思考题：what 和dir 的区别？
2.工作空间：显示数据的变量信息，包括变量名、字节大小、变量类型等。 输入 load wind 和load cities 加载了后缀为mat 的wind 和cities 数据文件 在命令窗口输入who 就可以列出空间的变量
在命令窗口输入whos 可以列出名称、大小和类型 whos -file 文件名.mat 可以查看加载前的数据信息。 3.历史记录：
显示所有在命令窗口输入的执行过的命令，清除历史的方法有两种
4.帮助
（1）help 函数名或命令名，可以显示对应的帮助信息 例如 help sin 可以得知sin 实现正弦函数的计算功能 例如对自编代码的解释，help l1_bdjf （2）联机帮助
（3）演示帮助
在命令窗口输入demo 或demos
第2节 工作环境
安装时必带虚拟java 机，在和其它程序如vc 、adams 等混合或联合编程时需要进行必要的配置。matlab 在工作时会调用它内部所带的函数、动态链接库等。这些函数有的是m 有的是c 文件、有的是mex 、有的是p 文件等。matlab 为了减轻人们的编程负担，编制了很多应用函数和API 函数。当然还有很多工具箱。 例如：
figure;peaks;查看当前时间 clock 查看当前日期 date xpbombs
funtool 函数运算工具
课后练习：
通过菜单preference 选项进行多种设置例如字体、背景颜色等。
第3节 Matlab 的历史
20世纪70年代，美国新墨西哥大学计算机科学系主任Cleve Moler 教授为方便学生在学习《数值计算》课程时矩阵运算而开发的FORTRAN 程序库和接口程序，取名MA TALAB 。1983年Cleve Moler 和John Little 等用C 语言开发了第二代MA TALAB 专业版，具有数值计算及数据图形化功能。
1984年成立MathWorks 公司，推出了商业版的MA TALAB4.0版。随后功能不断改进和完善。
第4节 Matlab 基本功能
1.数值计算 （1）普通计算
算术运算：+（plus ）、-（minus ）、*(times)、/(右除—即右为除数左为被除数)、\(左除—左数为除数)、幂方^(power) 注意：matlab 是以数组为基础运算的，也是以复数为基础运算的。 【例】用MATLAB 计算38- 能得到–2吗？8开立方有几个根？ 复数开方)sin (cos )]atan sin()atan [cos(22θθi r a
b
i a b b a bi a z +=++=
+=
r(cosθ+isinθ)的n 次方根
=r^(1/n) * (cosθ+isinθ)^ (1/n)=r^(1/n) * e^(i*(2k*pi+θ)/n)=r^(1/n) * (cos((2k*pi+θ)/n)+isin((2k*pi+θ)/n)) 其中的k=0,1,2,...(n-1)，所以:有n 个复数根
%%%%%%%%%%%%%%%%%%%%%%%%%%% % 一般程序
%学习angle( )函数 a=-8; r=a^(1/3) m=[0,1,2];
R=abs(a)^(1/3);
Theta=(angle(a)+2*pi*m)/3; rrr=R*exp(i*Theta)
t=0:pi/20:2*pi;x=R*sin(t);y=R*cos(t); plot(x,y,'b:'),grid hold on
plot(rrr([2,3]),'o','MarkerSize',15,'Color','b')
axis([-3,3,-3,3]),axis square
hold off
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
那如何利用matlab进行-8开立方呢？
syms x;f=x^3+8;y=solve(f);y1=y(1)
知识点：y1是什么类型的数据？ischar(y1);isnumeric(y1);iscell(y1);如何知道类型：class(y1)
练习：8开立方
在这里要复习一下复数开n次方的意义。
（2）符号运算
定义分别用sym和syms 思考：sym与syms的不同
clc;clear;syms a b;z=a+b
如何定义一个复数
clc;clear;a=sym('a ', 'real');b=sym('b', 'unreal');z=a+b; real(y)
clc;clear;a=sym('a ', 'real');b=sym('b', 'real');z=a+b*i; real(y)
clc;clear;syms a b real;y = a+b*i; real(y)
知识点：将数组转化为符号矩阵A=[3 1.1 2;2 4 1.5;pi 0.5 5];B=sym(A)
符号矩阵转化为数值矩阵C=double(B)
思考：非数值型的符号矩阵能不能转？n=['a' 'b' 'c' 'd'; 'b' 'c' 'd' 'a'; 'c' 'd' 'a' 'b'; 'd' 'a' 'b' 'c'];D=sym(n) 2.数据处理
l2_clown.m
Matlab的命令：
1.dos命令。已经学过建立文件夹，设置文件夹，
dir;delete 文件名;edit;exit; 查看主机ip: dos('ipconfig') 查看主机名：dos('hostname')
关机dos('shutdown /s /f -t 3600');取消dos('shutdown /a')
2.帮助命令；
3.操作命令(部分汇总)
一．常用的窗口命令
help 启动联机帮助文件显示
what 列出当前目录下的有关文件
lookfor 对help信息中的关键词查找
which 找出函数与文件所在的目录名
demo 运行MATLAB的演示程序
path 设置或查询MATLAB的路径
二．有关文件及其操作的语句
cd 改变当前的工作目录
dir 列出当前目录的内容
getenv 获得环境参数
fopen 打开文件
fclose 关闭文件
三．启动与退出的命令
quit,exit 退出MATLAB环境
startup MATLAB自启动文件
matlabrc 启动主程序
四．管理变量工作空间的命令
who 简要列出工作空间变量名
whos 详细列出工作空间变量名
load 从文件中读入变量
save 列出工作空间中变量存盘
clear 删除内存中的变量与函数
五．对命令窗口控制的常用命令
cedit 设置命令行编辑与回调的参数
clc 清除命令窗口中的显示
format 设置输出格式
echo 显示文件中的MATLAB命令
more 控制命令窗口的输出页面
Matlab语言的特点：
1.高级性
比通用编程语言c等更高级，定义不是很严格。例如c没有圆周率，1/0会出错，没有inf eps NaN等。
知识点：NaN的用途
bdw=0.01;tpw=0.15;pos=[1/2+bdw,2/3+bdw,1/2-2*bdw,1/3-bdw-tpw];
figure('NumberTitle','off','Units','normalized','Position',pos,'Color',[1,1,1]); set(gcf,'Name','试验窗');
P=ones(16,16)*NaN;
P(1,5)=1;P(2:14,4)=1;P(15,5)=1;P(16,6:11)=1;P(15,12)=1;P(9:14,13)=1;
P(2:9,6)=1;P(9,7)=1;P(8,8)=1;P(9,9)=1;P(8,10)=1;P(9,11)=1;P(8,12)=1;
P(2:14,5)=2;P(10:15,6:11)=2;P(9,8)=2;P(9,10)=2;P(9:14,12)=2;
set(gcf,'Pointer','Custom','PointerShapeCData',P,'PointerShapeHotSpot',[2,5])
2.通用性
应用领域广泛，和其他专业软件有接口
3.专用工具箱；
4.GUI
执行程序有丰富的界面，可以生成.exe文件，例zhizhenbiao.exe
Matlab的编程规则：
1.文件格式
.m .fig .mdl .mat .p .asv .mex .c .h .lib .exe .dll
例如hello.c
知识点：查询当前工作目录下的文件种类例如i=exist('hello') 结果i=3. 或者用i=exist('hello', '种类名')
'种类名'= 'var','builtin','dir', or 'file'
type('hello')
.p为预编译文件，可以保护系统中的文件，也可以隐藏代码。生成p代码文件的方法是将现有的<文件> .m ，在命令行输入pcode <文件> .m 有没有后缀都可以，见例l3_pcode
2.文件名
1) 文件名命名要用英文字符，第一个字符不能是数字和特殊字符，不可以是汉字或希腊字母
2) 文件名不要取为matlab的一个固有函数,m文件名的命名尽量不要是简单的英文单词，最好是由大小写英文/数字/下划线等组成。原因是简单的单词命名容易与matlab内部函数名同名，结果会出现一些“莫名其妙”的错误。
3) 文件存储路径一定为英文
4) m文件起名不能为两个单词，如three phase，应该写成three_phase或者ThreePhase
3.变量与常量
Matlab变量及命名规则
Matlab中所有的变量都是用矩阵形式来表示的，即所有的变量都表示一个矩阵或者一个向量。其命名规则如下：?变量名区分大小写
?变量名长度最多63 个字符
?变量名必须以英文字母开头，后可字母、数字和下划线
?变量名中不许出现标点符号
?保留字不能用作变量名
如果对它们重新赋值，则失去原来意义。
4.保留字与关键字：for end if while function return elseif case switch break等。abs bone box close conv conj disp echo find floor global grid hold inv jet keyboard log magic notebook（起动word）on ones pwd（显示当前工作目录）quiver rand rank randperm simulink（起动simulink）sim（运行mdl模型）text title version。
第5节Matlab的数值类型与数据结构
1.数值类型
基本类型为数组（Array），在存储时按列存储，与FORTRAN方式相同，与C语言有较大差别（数组下标从1开始，在编程时容易出错的地方
．．．．．．．．．．．）。
（1）数量（Scalar）、矢量（Vector）、矩阵（Matrix）、字符串（String）、单元阵列（Cell array）和结构（Structure）
（2）整型数据相关函数
ceil:朝正无穷大方向取整
class：返回对象的数据类型
fix：朝零方向取整
floor：朝负无穷大方向取整
isa 检查是否给定类的对象isa('a','char')
isinteger
isnumeric
round：四舍五入取整
mod 模数求余
rem 求余数
sign 符号函数
（3）逻辑数据：True —用1表示，False —用0表示，常用作程序结构的判断条件
1. 逻辑运算符~ & | Xor（异或，如两输入相同则输出为假，如输入一真一假输出为真）
2. 关系运算符= = ~ = > >= < <=
3. 常见逻辑函数
iskeyword(cell)
动态规划-图论
§1动态规划模型 如图所示，给定一个线路网络，两点之间连线上的数字表示 两点间距离，试求一条从A到E的路线，使总距离为最短。Mattlab求解： 首先利用Excel建立两个工作表edge和n分别存储图的上三 角阵和顶点数量。其中edge= 99999 5 2 99999 99999 99999 99999 99999 99999 99999 99999 99999 3 7 99999 99999 99999 99999 99999 99999 99999 99999 6 3 99999 99999 99999 99999 99999 99999 99999 99999 99999 6 99999 99999 99999 99999 99999 99999 99999 99999 3 8 99999 99999 99999 99999 99999 99999 99999 99999 1 99999 99999 99999 99999 99999 99999 99999 99999 99999 3 99999 99999 99999 99999 99999 99999 99999 99999 7 99999 99999 99999 99999 99999 99999 99999 99999 99999 n=9,然后在Matlab调入以上数据。同时将自编的动态规划 软件“dynamic.m”调入当前目录之中，在Matlab命令窗口
输入dynamic,回车后则在窗口显示出路径Path 和距离distance §2 最小生成树 例1 某工厂要架设局域网联通工厂各个部门。已知工厂有7个部门，各个部门间铺设网线的距离如上图所示，计算出铺设网线的最短距离。 Matlab 的算法： 首先，将上图的邻接矩阵存储为G ,顶点数存储为N ；即：G= 99999 50 60 99999 99999 99999 99999 50 99999 99999 65 40 99999 99999 60 99999 99999 52 99999 99999 45 99999 65 52 99999 50 30 42 99999 40 99999 50 99999 70 99999 99999 99999 99999 30 70 99999 99999 99999 99999 45 42 99999 99999 99999 2 5 3 1 4 7 6 50 60 45 65 52 40 50 70 30 42
MATLAB复习知识点汇总
MATLAB 复习 1、实现符号函数 运算功能的函数m 文件为： function y=sgn(x) if x<0 y1=-1; elseif x==0 y1=0; else y1=1; end y=y1; 2、求满足1+2+3+…n<100的最大正整数n 的MATLAB 程序为： sum=0;n=0; while sum<100 n=n+1; sum=sum+n; end 1,0,sgn()0,0, 1,0x y x x x >?? ===??-
sum=sum-n; n=n-1; n,sum 3、m文件有两种形式，一种称为命令文件（Script File），另一种称为函数文件（Function File），两种文件的扩展名都是m。 4、反馈控制系统品质要求：稳定性、准确性、快速性 5、按不同系统的特征方程式，可将自动控制系统分为线性系统和非线性系统。 6、MATLAB控制相关的工具箱 （1）控制系统工具箱 （2）系统辨识工具箱 （3）模型预测控制工具箱 （4）鲁棒控制工具箱 （5）神经网络工具箱 7、MATLAB系统的构成 （1）开发环境（2）数学函数库（3）MATLAB语言（4）图形处理系统（5）应用程序接口 8、常用工具箱 （1）控制类工具箱（2）应用数学类工具箱（3）信号处理类工具箱（4）其他常用工具箱 9、MATLAB语句形式：>>变量=表达式 10、MATLAB常用命令 quit 关闭MATLAB
exit 关闭MATLAB clc 清除MATLAB 命令窗口中的所有显示内容 clear 清除工作空间中保存的所有变量 11、MATLAB 基本数据类型：双精度数组、字符串数组、元胞数组、结构数组 12、矩阵的创建 （1）直接输入法（2）通过数据文件创建矩阵（3）通过m 文件创建矩阵（4）通过函数创建矩阵（5）冒号法 13、冒号法 ［1］冒号法构造向量 冒号表达式的一般格式为：向量名＝初值：步长：终值 ［2］冒号法构造矩阵 一般格式为： A(:,j ):表示矩阵A 的第j 列； A(i ,:):表示矩阵A 的第i 行。 A(i,j)表示取矩阵A 的第i 行第j 列交叉位置的元素 14、矩阵的运算 /(矩阵的右除)或\(矩阵的左除) A -1=inv(A) 矩阵的逆 15、求下面方程组的根 1231231 2323532255316 x x x x x x x x x +-=?? -+=??--=?
matlab的RBF-BP神经网络讲义
matlab的RBF BP神经网络讲义 一、RBF神经网络 1985年，Powell提出了多变量插值的径向基函数(Radical Basis Function，RBF)方法， 1988年，Moody和Darken提出了一种神经网络结构，即RBF神经网络。 RBF网络是一种三层前向网络，其基本思想是：（1）用RBF作为隐单元的“基”构成隐含层空间，将输入矢量直接(即不需要通过权连接)映射到隐空间（2）当RBF的中心点确定后，映射关系也就确定（3）隐含层空间到输出空间的映射是线性的。 newrb()函数 功能 建立一个径向基神经网络 格式 net = newrb(P，T，GOAL，SPREAD，MN，DF) 说明 P为输入向量，T为目标向量，GOAL为圴方误差，默认为0，SPREAD为径向基函数的分布密度，默认为1，MN为神经元的最大数目，DF为两次显示之间所添加的神经元神经元数目。 例子： 设[P,T]是训练样本，[X,Y]是测试样本； net=newrb(P,T,err_goal,spread); %建立网络 q=sim(net,p); e=q-T; plot(p,q); %画训练误差曲线 q=sim(net,X); e=q-Y; plot(X,q); %画测试误差曲线 二、BP神经网络 训练前馈网络的第一步是建立网络对象。函数newff()建立一个可训练的前馈网络。这需要4个输入参数。 第一个参数是一个Rx2的矩阵以定义R个输入向量的最小值和最大值。 第二个参数是一个设定每层神经元个数的数组。 第三个参数是包含每层用到的传递函数名称的细胞数组。 最后一个参数是用到的训练函数的名称。 举个例子，下面命令将创建一个二层网络。它的输入是两个元素的向量，第一层有三个神经元(3)，第二层有一个神经元(1)。 第一层的传递函数是tan-sigmoid，输出层的传递函数是linear。 输入向量的第一个元素的范围是-1到2[-1 2]，输入向量的第二个元素的范围是0到5[0 5]，训练函数是traingd。 net=newff([-1 2; 0 5],[3,1],{'tansig','purelin'},'traingd'); 这个命令建立了网络对象并且初始化了网络权重和偏置，因此网络就可以进行训练了。 我们可能要多次重新初始化权重或者进行自定义的初始化。 下面就是初始化的详细步骤。 在训练前馈网络之前，权重和偏置必须被初始化。初始化权重和偏置的工作用命令init来实
MatLab讲义
2011年数学中国国赛培训讲座 Matlab的基础及数学建模中的应用 周吕文:zhou.lv.wen@https://www.wendangku.net/doc/a8135190.html, 大连大学数学建模工作室&中国科学院力学研究所 2011年7月
第一部分 MatLab基础 1 简单介绍 MATLAB是Matrix Laboratory“矩阵实验室”的缩写。MatLab语言是由美国的Clever Moler博士于1980年开发的，初衷是为解决“线性代数”课程的矩阵运算问题。1984年由美国 MathWorks公司推向市场，历经十多年的发展与竞争，现已成为国际公认的最优秀的工程应用开发环境。MATLAB功能强大、简单易学、编程效率高，深受广大科技工作者的欢迎。 在数学建模竞赛中，由于只有短短的三到四天，而论文的评判不仅注重计算的结果更注重模型的创造性等很多方面，因此比赛中把大量的时间花费在编写和调试程序上只会喧宾夺主，是很不值得的。使用MATLAB 可以很大程度上的方便计算、节省时间，使我们将精力更多的放在模型的完善上，所以是较为理想的。 这里快速的介绍一下MATLAB与数学建模相关的基础知识，并列举一些简单的例子，很多例子都是源于国内外的数学建模赛题。希望能帮助同学们在短时间内方便、快捷的使用MATLAB 解决数学建模中的问题。当然要想学好MatLab更多的依赖自主学习，一个很好的学习MatLab的方法是查看MatLab的帮助文档： z如果你知道一个函数名，想了解它的用法，你可以用'help'命令得到它的帮助文档：>>help functionname z如果你了解含某个关键词的函数，你可以用'lookfor'命令得到相关的函数： 2 基本命令与函数 基本运算 z变量的赋值 实数赋值>> x=5; 复数赋值>> x=5+10j; (或>>x=5+10i) z向量的一般值方法 行向量赋值：>>x=[1 2 3]; (或x=[1, 2 ,3]) 列向量赋值：>>y=[1;2;3]; 矩阵的赋值：>>x=[1 2 3; 4 5 6; 7 8 9]; z常用矩阵（zeros ones eye） n行m列0矩阵：>>x=zeros(n,m); n行m列1矩阵：>>x=ones(n,m); n 阶的单位阵：>>y=eye(n); z矩阵行列操作 >> A=[1 2 3;4 5 6;7 8 9] A= 1 2 3 4 5 6 7 8 9 >>x=A(1,3) %取第一行的第三列元素 x= 3
天津大学matlab讲义-应用基础第一章
MATLAB应用基础 赵国瑞 天津大学电子信息工程学院 计算机基础教学部 2000.3 制作
概述 MATLAB是世界流行的优秀科技应用软件之一。具有功能强大(数值计算、符号计算、图形生成、文本处理及多种专业工具箱)、界面友好，可二次开发等特点。 自1984年由美国MathWorks公司推向市场以来，先后发布了多个版本，1993年发布4.0版，1996年发布5.0版，1999年发布5.3版。目前发布的为6.5版。 MATLAB有专业和学生版之分。二者功能相同，但计算规模和计算难度有差别。 在国内外，已有许多高等院校把MATLAB列为本科生、研究生必须掌握的基本技能。我校自1999年列为研究生选修课程。而且有很多教师、研究生把它作为进行科研的重要工具。 国内关于MATLAB的书籍很多，如： 《精通MATLAB 5.3》张志涌等编著北京航空航天大学出版社，2000.8 《科学计算语言MATLAB简明教程》杜藏等编著南开大学出版社，1999.6 《精通MATLAB 5》张宜华编写清华大学出版社，1999.6 《精通MATLAB--综合辅导与指南》 Duane Hanselman、Bruce Littlefield编著李人厚等译较西安交通大学出版社，1998.1 等等 本课程主要介绍MATLAB 5.3的基本功能和基础知识。至于其包含的多种工具箱，如仿真工具箱、解非线性方程（组）工具箱、优化工具箱等，应通过本学习后，结合各专业自己进一步学习和使用。 第1章MATLAB基础 1.1 源文件（M－文件） 分为两类：函数文件和非函数文件。 都用扩展名.M 1.1.1函数文件 格式1（无返回值函数） function函数名(输入表) %称为函数头 函数体 例如： function box(opt_box); %BOX Axis box. % BOX ON adds a box to the current axes. % BOX OFF takes if off. % BOX, by itself, toggles the box state. % % BOX sets the Box property of the current axes. % % See also GRID, AXES. % Copyright (c) 1984-98 by The MathWorks, Inc. % $Revision: 1.5 $ $Date: 1997/11/21 23:32:59 $
动态计划求解方法的Matlab实现及应用[]
动态规划求解方法的Matlab实现及应用[1].txt我自横刀向天笑，笑完我就去睡觉。你的手机比话费还便宜。路漫漫其修远兮，不如我们打的吧。第 %卷第 ,期信息工程大学学报 S>:+% <>+, !""’年 >月 T>8D3F: >C 53C>DEFB2>3 G3?23@@D23? 032H@DA2BI 6@N+!""’ !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! 动态规划求解方法的 !"#$"%实现及应用 于斌，刘姝丽，韩中庚 <信息工程大学信息工程学院，河南郑州 #’"""!） 摘要：文章对动态规划问题的求解方法进行了分析研究，根据问题的特点、难点和关键点做了 针对性的处理，然后用 !"#$"%做了实现尝试，从而实现了“最佳组队”和“最短路线”等问题的 求解。实践证明所采用方法和程序都是有效的。 关键词：动态规划；基本方程；!"#$"%实现；最佳组队 中图分类号：* !!&+,文献标识码：-文章编号：&%.& $ "%.,’
$ "# !"#$"% &’"$(>"#(*+ *, #-’ ./+"0(1 23*43"00(+4 5663*"1-"+7 8#9 566$(1"#(*+ /0 123，450 6789:2，。-< =7>3?9?@3? <53AB2B8B@ >C 53C>DEFB2>3 G3?23@@D23?，53C>DEFB2>3 G3?23@@D23? 032H@DA2BI，=7@3?J7>8 #’"""!，K723F） 5%9#3"1#：1I F3F:IJ23? F3L 23H@AB2?FB23? B7@ LI3FE2M ND>?DFEE23? FNND>FM7，F3 @CC@MB2H@ L2AN>AF: 7FA O@@3 L>3@
MATLAB讲义
第一章基础准备及入门 什么是MATLAB？ MATLAB是MathWorks公司于1984年推出的数学软件，是一种用于科学工程计算的高效率的高级语言。MATLAB最初作为矩阵实验室（Matrix Laboratory），主要向用户提供一套非常完善的矩阵运算命令。随着数值运算的演变，它逐渐发展成为各种系统仿真、数字信号处理、科学可是化的通用标准语言。 在科学研究和工程应用的过程中，往往需要大量的数学计算，传统的纸笔和计算机已经不能从根本上满足海量计算的要求，一些技术人员尝试使用Basic,Fortran,C\C++等语言编写程序来减轻工作量。但编程不仅需要掌握所用语言的语法，还需要对相关算法进行深入分析，这对大多数科学工作者而言有一定的难度。与这些语言相比， MATLAB的语法更简单，更贴近人的思维方式。用MATLAB编写程序，犹如在一张演算纸上排列公式和求解问题一样高效率，因此被称为“科学便笺式”的科学工程计算语言。 MATLAB由主包和功能各异的工具箱组成，其基本数据结构是矩阵。正如其名“矩阵实验室”，MATLAB起初主要是用来进行矩阵运算。经过MathWorks 公司的不断完善，时至今日，MATLAB已经发展成为适合多学科、多工作平台的功能强大的大型软件。 本章有两个目的：一是讲述MATLAB正常运行所必须具备的基础条件；二是简明系统地介绍高度集成的Desktop操作桌面的功能和使用方法。 本章的前两节分别讲述：MATLAB的正确安装方法和MATLAB 环境的启动。因为指令窗是MATLAB最重要的操作界面，所以本章用第 1.3、1.4 两节以最简单通俗的叙述、算例讲述指令窗的基本操作方法和规则。这部分内容几乎对MATLAB各种版本都适用。 MATLAB6.x 不同于其前版本的最突出之处是：向用户提供前所未有的、成系列的交互式工作界面。了解、熟悉和掌握这些交互界面的基本功能和操作方法，将使新老用户能事半功倍地利用MATLAB去完成各种学习和研究。为此，本章特设几节用于专门介绍最常用的交互界面：历史指令窗、当前目录浏览器、工作空间浏览器、内存数组编辑器、交互界面分类目录窗、M文件编辑/调试器、及帮助导航/浏览器。 本章是根据MATLAB6.5版编写的，但大部分内容也适用于其他6.x版。 1.1M ATLAB的安装和内容选择
动态规划 销售人员分配问题(matlab编程)
数学规划课程设计 题目：销售人员费配问题 姓名： 学号： 成绩： 2011年6月
销售人员费配问题 摘要：动态规划程序设计是对解最优化问题的一种途径、一种方法，而不是一种特殊算法，本论文通过对动态规划的基本概念和基本思路，并利用Matlab对动态规划中的销售人员分配问题进行了分析，然后利用Matlab语言进行了程序设计和计算，是复杂问题简单化，避免了繁琐的计算，从而使问题能跟方便地得到解决。 关键词：动态规划销售人员分配问题Matlab语言
一、问题重述 某企业甲、乙、丙三个销售市场，其市场的利润与销售人员的分配有关，现有6个销售人员， 二、问题分析 首先我们对设备的分配规定一个顺序，即先考虑分配给甲市场，其次乙市场，最后丙市场，但分配时必须保证企业的总收益最大。 将问题按分配过程分为三个阶段，根据动态规划逆序算法，可设： 1、阶段数k=1，2，3（即甲、乙、丙三个市场的编号分别为1，2，3）； 2、状态变量x k 表示分配给第k 个市场至第3个市场的人员数（即第k 阶段初尚未分配的人员数）； 3、决策变量u k 表示分配给第k 市场的人员数； 4、状态转移方程：x k+1=x k -u k ; 5、g k (u k )表示u k 个销售人员分配到第k 个市场所得的收益值，它由下表可查得； 6、f k (x k )表示将x k 个销售人员分配到第k 个市场所得到的最大收益值，因而可得出递推方程： f k (x k )= 6 ,...,1,0max =k u [ g k (u k )+ f k+1(x k -u k )]，k=1，2，3 f 4(x 4)=0 三、问题求解 1）k=3时，市场丙的分配方案和总收益. 最大收益：f 3(x 3)=6 ,...,1,0max 3=u [g 3(x 3)]
Matlab基础知识点汇总
MATLAB讲义 第一章 MATLAB系统概述 1.1 MATLAB系统概述 MATLAB（MATrix LABoratory）矩阵实验室的缩写，全部用C语言编写。 特点： （1）以复数矩阵作为基本编程单元，矩阵运算如同其它高级语言中的语言变量操作一样方便，而且矩阵无需定义即可采用。 （2）语句书写简单。 （3）语句功能强大。 （4）有丰富的图形功能。如plot,plot3语句等。 （5）提供了许多面向应用问题求解的工具箱函数。目前，有20多个工具箱函数，如信号处理、图像处理、控制系统、系统识别、最优化、神经网络的模糊系统等。 （6）易扩充。 1.2 MATLAB系统组成 （1）MATLAB语言 MATLAB语言是高级的矩阵、矢量语言，具有控制流向语句、函数、数据结构、输入输出等功能。同时MATLAB又具有面向对象编程特色。MATLAB语言包括运算符和特殊字符、编程语言结构、字符串、文件输入/输出、时间和日期、数据类型和结构等部分。 （2）开发环境 MATLAB开发环境有一系列的工具和功能体，其部分具有图形用户界面，包括MATLAB桌面、命令窗口、命令历史窗口、帮助游览器、工作空间、文件和搜索路径等。 （3）图形处理 图形处理包括二维、三维数据可视化，图像处理、模拟、图形表示等图形命令。还包括低级的图形命令，供用户自由制作、控制图形特性之用。 （4）数学函数库 有求和、正弦、余弦等基本函数到矩阵求逆、求矩阵特征值和特征矢量等。 MATLAB数学函数库可分为基本矩阵和操作、基本数学函数、特殊化数学函数、线性矩阵函数、数学分析和付里叶变换、多项式和二重函数等。 （5）MATLAB应用程序接口（API） MATLAB程序可以和C/C++语言及FORTRAN程序结合起来，可将以前编写的C/C++、FORTRAN语言程序移植到MATLAB中。 1.3 MATLAB的应用围包括： MATLAB的典型应用包括： ●数学计算 ●算法开发 ●建模、仿真和演算 ●数据分析和可视化 ●科学与工程绘图 ●应用开发（包括建立图形用户界面） 以矩阵为基本对象 第二章 Matlab基础
最优化方法的Matlab实现(公式(完整版))
第九章最优化方法的MatIab实现 在生活和工作中，人们对于同一个问题往往会提出多个解决方案，并通过各方面的论证从中提取最佳方案。最优化方法就是专门研究如何从多个方案中科学合理地提取出最佳方案的科学。由于优化问题无所不在，目前最优化方法的应用和研究已经深入到了生产和科研的各个领域，如土木工程、机械工程、化学工程、运输调度、生产控制、经济规划、经济管理等，并取得了显著的经济效益和社会效益。 用最优化方法解决最优化问题的技术称为最优化技术，它包含两个方面的内容： 1）建立数学模型即用数学语言来描述最优化问题。模型中的数学关系式反映了最优化问题所要达到的目标和各种约束条件。 2）数学求解数学模型建好以后，选择合理的最优化方法进行求解。 最优化方法的发展很快，现在已经包含有多个分支，如线性规划、整数规划、非线性规划、动态规划、多目标规划等。 9.1 概述 利用Matlab的优化工具箱，可以求解线性规划、非线性规划和多目标规划问题。 具体而言，包括线性、非线性最小化，最大最小化，二次规划，半无限问题，线性、非线性方程（组）的求解，线性、非线性的最小二乘问题。另外，该工具箱还提供了线性、非线性最小化，方程求解，曲线拟合，二次规划等问题中大型课题的求解方法，为优化方法在工程中的实际应用提供了更方便快捷的途径。 9.1.1优化工具箱中的函数 优化工具箱中的函数包括下面几类： 1 ?最小化函数
2.方程求解函数 3.最小—乘（曲线拟合）函数
4?实用函数 5 ?大型方法的演示函数 6.中型方法的演示函数 9.1.3参数设置 利用OPtimSet函数，可以创建和编辑参数结构；利用OPtimget函数，可以获得o PtiOns优化参数。 ? OPtimget 函数 功能：获得OPtiOns优化参数。 语法：
Matlab入门基础知识整理
MATLAB基础知识 MATLAB的主要功能 1．数值计算和符号计算功能 例如，求解线性方程组 在MATLAB命令窗口输入命令： a=[2,3,-1;3,-5,3;6,3,-8]; b=[7;8;9]; x=inv(a)*b 也可以通过符号计算来解此方程 syms x1 x2 x3 [x1,x2,x3]=solve(2*x1+3*x2-x3-7,3*x1-5*x2+3*x3-8,6*x1+3*x2-8*x3-9) 2．绘图功能 例如，分别绘制函数y=300sinx/x和y=x2的曲线 x=-20:0.1:20; plot(x,300*sin(x)./x,':',x,x.^2); 3．程序设计语言功能 MATLAB是解释性语言，程序执行速度较慢，而且不能脱离MATLAB环境而独立运行。MathWorks公司提供了将MATLAB源程序编译为独立于MATLAB集成环境运行的EXE文件以及将MATLAB程序转化为C语言程序的编译器。4．扩展功能 MATLAB包含两部分内容：基本部分和各种可选的工具箱。 基本部分构成了MATLAB的核心内容，也是使用和构造工具箱的基础。 工具箱扩展了MATLAB的功能。 功能性工具箱 学科性工具箱 MATLAB的集成开发环境 命令窗口（Command Window）用于输入命令并显示除图形以外的所有执行结果。一般来说，一个命令行输入一条命令，命令行以回车或分号结束 但一个命令行也可以输入若干条命令，各命令之间以逗号或分号隔开 如果一个命令行很长，一个物理行之内写不下，可以在第1个物理行之后加上3个小黑点并按下回车键，然后接着下一个物理行继续写命令的其他部分。3个小黑点称为续行符，即把下面的物理行看作该行的逻辑继续。例如： z=1+1/(1*2)+1/(1*2*3)+1/(1*2*3*4)+ … 1/(1*2*3*4*5) 工作空间（Workspace）是MATLAB用于存储各种变量的内存空间。 当前目录（Current Directory）是指MATLAB运行时的工作目录，只有在当前目录下的文件、函数才可以被运行或调用。 如果没有特殊指明，数据文件也将存放在当前目录下。 命令历史窗口 自动保留自安装起所有用过的命令的历史记录，并且还标明了使用时间，从而方便用户查询。而且，通过双击命令可进行历史命令的再运行。 MATLAB的帮助功能 进入MATLAB帮助界面可以通过以下方法。 ●单击MATLAB主窗口工具栏中的Help按钮。
基于Matlab的动态规划程序实现
动态规划方法的Matlab 实现与应用 动态规划(Dynamic Programming)是求解决策过程最优化的有效数学方法，它是根据“最优决策的任何截断仍是最优的”这最优性原理，通过将多阶段决策过程转化为一系列单段决策问题，然后从最后一段状态开始逆向递推到初始状态为止的一套最优化求解方法。 1．动态规划基本组成 (1) 阶段 整个问题的解决可分为若干个阶段依次进行，描述阶段的变量称为阶段变量，记为k (2) 状态 状态表示每个阶段开始所处的自然状况或客观条件，它描述了研究问题过程的状况。各阶段状态通常用状态变量描述，用k x 表示第k 阶段状态变量，n 个阶段决策过程有n+ 1个状态。 (3) 决策 从一确定的状态作出各种选择从而演变到下一阶段某一状态，这种选择手段称为决策。描述决策的变量称为决策变量，决策变量限制的取值范围称为允许决策集合。用()k k u x 表示第k 阶段处于状态k x 时的决策变量，它是k x 的函数。用()k k D x Dk(xk)表示k x 的允许决策的集合。 (4) 策略 每个阶段的决策按顺序组成的集合称为策略。由第k 阶段的状态k x 开始到终止状态的后部子过程的策略记为{}11(),(),,()k k k k n n u x u x u x ++ 。可供选择的策略的范围称为允许策略集合，允许策略集合中达到最优效果的策略称为最优策略。从初始状态* 11()x x =出发，过程按照最优策略和状态转移方程演变所经历的状态序列{ } **** 121,,,,n n x x x x + 称为最优轨线。 (5) 状态转移方程 如果第k 个阶段状态变量为k x ，作出的决策为k u ，那么第k+ 1阶段的状态变量1k x +也被完全确定。用状态转移方程表示这种演变规律，记为1(,)k k k x T x u +=。 (6) 指标函数 指标函数是系统执行某一策略所产生结果的数量表示，是衡量策略优劣的数量指标，它定义在全过程和所有后部子过程上，用()k k f x 表示。过程在某阶段j 的阶段指标函数是衡量该阶段决策优劣数量指标，取决于状态j x 和决策j u ，用(,)j j j v x u 表示。 2．动态规划基本方程 (){} 11()min ,,(),()k k k k k k k k k k f x g v x u f x u D x ++=∈???? Matlab 实现 (dynprog.m 文件) function [p_opt,fval]=dynprog (x,DecisFun,SubObjFun,TransFun,ObjFun) % x 是状态变量，一列代表一个阶段的所有状态； % M-函数DecisFun(k,x) 由阶段k 的状态变量x 求出相应的允许决策变量; % M-函数SubObjFun(k,x,u) 是阶段指标函数, % M-函数ObjFun(v,f) 是第k 阶段至最后阶段的总指标函数 % M-函数TransFun(k,x,u) 是状态转移函数, 其中x 是阶段k 的某状态变量, u 是相应的决策变量; %输出 p_opt 由4列构成，p_opt=[序号组;最优策略组;最优轨线组;指标函数值组]； %输出 fval 是一个列向量，各元素分别表示p_opt 各最优策略组对应始端状态x 的最优函数值。
MATLAB软件基础知识讲义(doc 74页)
MATLAB软件基础 §1MATLAB 概述 MATLAB 是MATrix LABoratory（“矩阵实验室”）的缩写，是由美国MathWorks 公司开发的集数值计算、符号计算和图形可视化三大基本功能于一体的，功能强大、操作简单的语言。是国际公认的优秀数学应用软件之一。 20世纪80年代初期，Cleve Moler与John Little等利用C语言开发了新一代的MATLAB语言，此时的MATLAB语言已同时具备了数值计算功能和简单的图形处理功能。1984年，Cleve Moler与John Little等正式成立了Mathworks公司，把MA TLAB 语言推向市场，并开始了对MATLAB工具箱等的开发设计。1993年，Mathworks公司推出了基于个人计算机的MATLAB 4.0版本，到了1997年又推出了MATLAB 5.X版本（Release 11），并在2000年又推出了最新的MATLAB 6版本（Release 12）。 现在，MATLAB已经发展成为适合多学科的大型软件，在世界各高校，MATLAB已经成为线性代数、数值分析、数理统计、优化方法、自动控制、数字信号处理、动态系统仿真等高级课程的基本教学工具。特别是最近几年，MATLAB在我国大学生数学
建模竞赛中的应用，为参赛者在有限的时间内准确、有效的解决问题提供了有力的保证。 概括地讲，整个MATLAB系统由两部分组成，即MATLAB 内核及辅助工具箱，两者的调用构成了MATLAB的强大功能。MATLAB语言以数组为基本数据单位，包括控制流语句、函数、数据结构、输入输出及面向对象等特点的高级语言，它具有以下主要特点： 1）运算符和库函数极其丰富，语言简洁，编程效率高，MATLAB 除了提供和C语言一样的运算符号外，还提供广泛的矩阵和向量运算符。利用其运算符号和库函数可使其程序相当简短，两三行语句就可实现几十行甚至几百行C或FORTRAN的程序功能。 2）既具有结构化的控制语句（如for循环、while循环、break 语句、if语句和switch语句），又有面向对象的编程特性。 3）图形功能强大。它既包括对二维和三维数据可视化、图像处理、动画制作等高层次的绘图命令，也包括可以修改图形及编制完整图形界面的、低层次的绘图命令。 4）功能强大的工具箱。工具箱可分为两类：功能性工具箱和学科性工具箱。功能性工具箱主要用来扩充其符号计算功能、图
Matlab仿真常用知识及参考程序
实验一 三、 5.设计M 文件计算: x=0:0.1:10 当sum>1000时停止运算，并显示求和结果及计算次数。 参考程序： x=0:0.1:10; i=1; sum=x(i)*x(i)-2*x(i); while(sum<=1000) i=i+1; sum=sum+x(i)*x(i)-2*x(i); end if(sum>1000) sum=sum-x(i)*x(i)+2*x(i); i=i-1; end sum i 实验二 二、6、编写验证魔方矩阵的函数文件，输出要求如下： （1） 如果输入矩阵的维数小于3，输出显示’error’ （2） 如果输入矩阵的不是方阵，输出显示’the size of matrix X must be N-by-N matrix’ （3） 显示行、列和及其对角线求和后的值，并判断其和是否相同。若不同，显示‘No ’,相同显示‘Yes ’。 参考程序： function mofang_test(data) [Row,Com]=size(data); dimen=min(Row,Com); if(dimen<3) msgbox('error','Warning'); elseif(Row~=Com) msgbox('the size of matrix X must be N-by-N matrix','Warning'); else )x 2x (sum i 1000i 2 i -=∑=
j=Com:-1:1; SumRC=0; SumCR=0; for i=1:Row SumRC=SumRC+data(i,i); SumCR=SumCR+data(j(i),i); end if( SumRC~=SumCR) msgbox('No-1','Sorry'); else SumRow=zeros(Row); for i=1:Row for j=1:Com SumRow(i)=SumRow(i)+data(i,j); end end SumCom=zeros(Com); for j=1:Com for i=1:Row SumCom(j)=SumCom(j)+data(i,j); end end Value=SumRow(1); flag=0; for i=1:Row if(Value~=SumRow(i)) msgbox('No-2','Sorry'); flag=1; break; end end if(flag~=1) for i=1:Com if(Value~=SumCom(i)) msgbox('No-2','Sorry'); flag=1; break; end end end if(flag==0) msgbox('Yes','Prosit'); end end
图论算法及matlab程序的三个案例
图论实验三个案例 单源最短路径问题 Dijkstra 算法 Dijkstra 算法是解单源最短路径问题的一个贪心算法。其基本思想是，设置一个顶点集合S 并不断地作贪心选择来扩充这个集合。一个顶点属于集合S 当且仅当从源到该顶点的最短路径长度已知。设v 是图中的一个顶点，记()l v 为顶点 v 到源点v 1的最短距离， ,i j v v V ?∈，若 (,)i j v v E ?，记i v 到j v 的权ij w =∞。 Dijkstra 算法： ① 1{}S v =,1()0l v =；1{}v V v ??-,()l v =∞,1i =,1{}S V v =-； ② S φ=，停止，否则转③； ③ ()min{(),(,)} j l v l v d v v =， j v S ∈，v S ?∈； ④ 存在 1 i v +，使 1()min{()} i l v l v +=，v S ∈； ⑤ 1{} i S S v +=， 1{} i S S v +=-，1i i =+，转②； 实际上，Dijkstra 算法也是最优化原理的应用：如果12 1n n v v v v -是从1v 到 n v 的最短路径，则 12 1 n v v v -也必然是从1v 到 1 n v -的最优路径。 在下面的MATLAB 实现代码中，我们用到了距离矩阵，矩阵第i 行第j 行元 素表示顶点i v 到j v 的权ij w ，若i v 到j v 无边，则realmax ij w =，其中realmax 是 MATLAB 常量，表示最大的实数+308)。 function re=Dijkstra(ma)
实验二 MATLAB基础知识(二)
Experiment 1. Fundamental Knowledge of Matlab (II) 【Experimental Purposes】 1、熟悉并掌握MATLAB的工作环境。 2、运行简单命令，实现数组及矩阵的输入输出，了解在MATLAB下如何绘图。【Experimental Principle】 1. Vectors Let's start off by creating something simple, like a vector. Enter each element of the vector (separated by a space) between brackets, and set it equal to a variable. For example, to create the vector a, enter into the MATLAB command window (you can "copy" and "paste" from your browser into MATLAB to make it easy): a = [1 2 3 4 5 6 9 8 7] MATLAB should return: a = 1 2 3 4 5 6 9 8 7 To generate a series that does not use the default of incrementing by 1, specify an additional value with the colon operator (first:step:last). In between the starting and ending value is a step value that tells MATLAB how much to increment (or decrement, if step is negative) between each number it generates. To generate a vector with elements between 0 and 20, incrementing by 2(this method is frequently used to create a time vector), use t = 0:2:20 t = 0 2 4 6 8 10 12 14 16 18 20 Manipulating vectors is almost as easy as creating them. First, suppose you would like to add 2 to each of the elements in vector 'a'. The equation for that looks like: b = a + 2 b = 3 4 5 6 7 8 11 10 9 Now suppose, you would like to add two vectors together. If the two vectors are the same length, it is easy. Simply add the two as shown below: c = a + b c = 4 6 8 10 12 14 20 18 16 Subtraction of vectors of the same length works exactly the same way. MATLAB sometimes stores such a list in a matrix with just one row, and other times in a matrix with just one column. In the first instance, such a 1-row matrix is called a row-vector; in the
动态规划_销售人员分配问题(matlab编程)
一、问题重述 某企业甲、乙、丙三个销售市场，其市场的利润与销售人员的分配有关，现有6个销售人员，分配到各市场所获利润如下表示，试问应如何分配销售人员才能使总利润最大？ 二、问题分析 首先我们对设备的分配规定一个顺序，即先考虑分配给甲市场，其次乙市场，最后丙市场，但分配时必须保证企业的总收益最大。 将问题按分配过程分为三个阶段，根据动态规划逆序算法，可设： 1、阶段数k=1，2，3（即甲、乙、丙三个市场的编号分别为1，2，3）； 2、状态变量x k 表示分配给第k 个市场至第3个市场的人员数（即第k 阶段初尚未分配的人员数）； 3、决策变量u k 表示分配给第k 市场的人员数； 4、状态转移方程：x k+1=x k -u k ; 5、g k (u k )表示u k 个销售人员分配到第k 个市场所得的收益值，它由下表可查得； 6、f k (x k )表示将x k 个销售人员分配到第k 个市场所得到的最大收益值，因而可得出递推方程： f k (x k )= 6 ,...,1,0max =k u [ g k (u k )+ f k+1(x k -u k )]，k=1，2，3 f 4(x 4)=0 三、问题求解 1）k=3时，市场丙的分配方案和总收益. 最大收益：f 3(x 3)=6 ,...,1,0max 3=u [g 3(x 3)]
最大收益：f 2(x 2)=2 max u [g 2(u 2)+ f 3(x 3)]= 2 max u [g 2(u 2)+ f 3(x 2- u 2 )] 最大收益：f 1(x 1)=1 max u [g 1(u 1)+ f 2(x 1- u 1)]= max[g 1(u 1)+ f 2(4- u 1)] 为此，我们可以用Matlab 语言编程使问题能跟方便地得到解决，其算法设计如下图：
(整理)matlab 动态规划讲义.
第四章动态规划 §1 引言 1.1 动态规划的发展及研究内容 动态规划（dynamic programming）是运筹学的一个分支，是求解多阶段决策问题的最优化方法。20世纪50年代初R. E. Bellman等人在研究多阶段决策过程(multistep decision process)的优化问题时，提出了著名的最优性原理（principle of optimality），把多阶段过程转化为一系列单阶段问题，逐个求解，创立了解决这类过程优化问题的新方法—动态规划。1957年出版了他的名著《Dynamic Programming》，这是该领域的第一本著作。 动态规划问世以来，在经济管理、生产调度、工程技术和最优控制等方面得到了广泛的应用。例如最短路线、库存管理、资源分配、设备更新、排序、装载等问题，用动态规划方法比用其它方法求解更为方便。 虽然动态规划主要用于求解以时间划分阶段的动态过程的优化问题，但是一些与时间无关的静态规划（如线性规划、非线性规划），只要人为地引进时间因素，把它视为多阶段决策过程，也可以用动态规划方法方便地求解。 应指出，动态规划是求解某类问题的一种方法，是考察问题的一种途径，而不是一种特殊算法（如线性规划是一种算法）。因而，它不象线性规划那样有一个标准的数学表达式和明确定义的
一组规则，而必须对具体问题进行具体分析处理。因此，在学习时，除了要对基本概念和方法正确理解外，应以丰富的想象力去建立模型，用创造性的技巧去求解。 例1 最短路线问题 下面是一个线路网，连线上的数字表示两点之间的距离（或费用）。试寻求一条由A到G距离最短（或费用最省）的路线。 例2 生产计划问题 工厂生产某种产品，每单位（千件）的成本为1（千元），每次开工的固定成本为3（千元），工厂每季度的最大生产能力为6（千件）。经调查，市场对该产品的需求量第一、二、三、四季度分别为2，3，2，4（千件）。如果工厂在第一、二季度将全年的需求都生产出来，自然可以降低成本（少付固定成本费），但是对于第三、四季度才能上市的产品需付存储费，每季每千件的存储费为0.5（千元）。还规定年初和年末这种产品均无库存。试制定一个生产计划，即安排每个季度的产量，使一年的总费用（生产成本和存储费）最少。 1.2 决策过程的分类