matlab求解-a(发展系数)的方法。
2015-10-09 10:22
309 查看
<灰色系统理论及应用>刘思峰第7版,P165~166,求解-a(发展系数)的方法。b:灰色作用量。
clear all;clc;close all
a=-0.1;
X0=[];
for i=0:5
X0= [X0,exp(-a*i)];
end
X1=cumsum(X0);
Z1=[0];
for i=2:6
Z1=[Z1,(X1(i)+X1(i-1))/2];
end
a=(sum(Z1(2:6))*sum(X0(2:6))-(6-1)*sum(Z1(2:6).*X0(2:6)))/((6-1)*sum(Z1(2:6).^2)-sum(Z1(2:6))^2);
b=( sum(X0(2:6))*sum(Z1(2:6).^2) - sum(Z1(2:6))* sum(Z1(2:6).*X0(2:6))) /((6-1)*sum(Z1(2:6).^2)-sum(Z1(2:6))^2)
求解出来:a是个正值 ,前面加个负号,变成-a(发展系数)。
clear all;clc;close all
a=-0.1;
X0=[];
for i=0:5
X0= [X0,exp(-a*i)];
end
X1=cumsum(X0);
Z1=[0];
for i=2:6
Z1=[Z1,(X1(i)+X1(i-1))/2];
end
a=(sum(Z1(2:6))*sum(X0(2:6))-(6-1)*sum(Z1(2:6).*X0(2:6)))/((6-1)*sum(Z1(2:6).^2)-sum(Z1(2:6))^2);
b=( sum(X0(2:6))*sum(Z1(2:6).^2) - sum(Z1(2:6))* sum(Z1(2:6).*X0(2:6))) /((6-1)*sum(Z1(2:6).^2)-sum(Z1(2:6))^2)
求解出来:a是个正值 ,前面加个负号,变成-a(发展系数)。
相关文章推荐
- matlab_BP
- python,matlab等字符串处理
- Matlab 的fspecial函数用法
- 目标检测光流法(四):matlab下光流法工具箱
- matlab里.p文件的一些说明
- matlab计时函数
- Matlab取整函数(fix, floor, ceil, round)的使用
- matlab 逻辑数组及其应用
- matlab 嵌套循环的学习
- matlab impixel的用法
- 将OpenCV里SIFT的keypoint坐标以及描述子向量输出给MATLAB
- MATLAB并行编程
- matalb 画plot图
- matlab中的max函数
- matlab中的元胞数组保存问题
- MATLAB并行计算工具箱使用
- 如何在Zemax与Matlab间通信
- Matlab 日期频次统计
- matlab中的sparse和full函数(稀疏矩阵和非稀疏矩阵之间的转换)
- Writing Fast Matlab code 3-5