基本语法
c=c(:) %将矩阵c转化为一维
zeros(10,25); %创建5行10列
ones(i);%产生全为1的矩阵
eye(2) %产生一个两行两列的单位矩阵 对角线全为1
x=linesapce(0,2,5) %产生0到2,5个步长的向量
Inf 无穷大
flipud:上下旋转 fliplr:左右旋转 rot90:逆时针旋转90
max(x) %计算出x向量中最大值
max(z, 2) %将z中列向量的值与2进行比较
max(z, [], 2) %将z中行向量的值进行比较
max(z, 2, 2) %将z中行向量的值与2进行比较
sum %求和 默认列向量相加
sum(x, 2)指定行向量相加
sum(x(:))则将所有数相加
abs %取绝对值;
sqrt %开根号
mod() %求余
axis([1,1,1,1,1,1]) %固定坐标轴显示范围
axis equal; %将坐标轴单位长度调整一致
一、第一章
- linprog(f,A,b,Aeq,beq,lb,ub)
f:价值
向量(目标函数的系数)
A,b:线性不等式的约束
Aeq,beq:线性等式约束
lb,ub:上下线性约束 - matlab中没有求最大值,应加负号转为求最小值问题
- 投资收益和风险问题
模型一:固定风险水平,优化收益
模型二:固定盈利水平,极小化风险
二、第二章
- 整数规划问题:将线性规划问题的变量限制为整数(尽管整数规划因为限制变量而增加困难,但是整数解是有限个,为枚举法提供了方便)
- 在一定计算量下,使用蒙特卡洛法可以得到满意解(选取足够多的随机数)
- 指派工作问题可以选择匈牙利算法、拍卖算法等
- 如指派问题等0-1整数规划问题,可以利用bintprog求解
三、一致性检验的步骤
- 计算一致性指标CI=(λmax -n)/(n-1)
- 查找对应的平均随机一致性指标RI
- 计算一致性比例CR
CR=CI/RI
【如果CR < 0.1,则可以判断矩阵的一致性可以接受,否则应该进行1修正】
四、特征值法求权重
- 检验判断矩阵一致性是否可以接受
CR<0.1 - 求出矩阵A的最大特征值以及对应的特征向量
λmax 、特征向量A
[V, D] = eig(A) %A为矩阵 V为特征向量 D为特征值
Max_eig = max(max(D)) %求出特征值中最大值,即求最大特征值
- 对求出的特征值进行归一化即可得到权重
A[1]/A - A[n]/A
在这里插入代码片
五、层次分析法AHP解法
- 找出评价目标:为小明同学选择最佳的旅游景点
- 为达到目标选择方案:三种,分别是去苏杭、北戴河、桂林
- 评价标准或者指标:景色、花费、居住、饮食、交通
六、构造一个简单的计算评分方法
(x - max) / (max - min)
七、判断多项式拟合阶数
for i=1:5 %最高进行5阶判断
y2=polyfit(x,y,i); %构造多项式
Y=polyval(y2,x);%计算拟合函数在x处的值。
if sum((Y-y).^2)<0.1
c=i %得到最合适阶数
break;
end
end