精华内容
下载资源
问答
  • 研究的结果是对粒子群优化算法定义了一个与迭代步有关的压缩因子,随着迭代步不断增大压缩因子逐渐减小,使得在算法初期,压缩因子较大,提高算法的全局搜索能力,在算法后期,压缩因子较小,提高算法的局部搜索能力,另外,...
  • 根据ISO12213-2:2009中Fortran程序编写,迭代过程有改动。该标准中混合物能量参数U公式B.6第一项缺少系数2,国标17747.2-1999中是正确的。
  • 使用c#窗体应用实现天然气压缩因子的计算 第2部分:用摩尔组成进行计算,支持修改温压和各成分占比,若使用的单片机内存足够,可进行移植
  • 国标 天然气压缩因子的计算 GB17747.1-2011
  • 天然气压缩因子的计算
  • 天然气压缩因子计算的软件非常非常的好用哦!
  • 新算法基于压缩因子粒子群算法,利用双重压缩因子;第一个压缩因子用来调节全局和局部搜索模型;第二压缩因子利用时变的加速因子,进一步平衡全局和局部最优值对粒子种群升级的影响;通过对基本粒子群算法,压缩因子...
  • SRK方程混合物压缩因子和焓值计算的MATLAB程序
  • 考虑DCT阈中量化效应,采用8X8的DCT模板,选择量化步长,结合缩放因子进行图像压缩
  • 设置压缩因子改变储存大小,并设置分辨率和DPI
  • 计算天然气压缩因子的公式,核实无误,放心使用
  • mat9_rk_压缩因子_源码

    2021-09-29 08:24:38
    用matlab编程计算混合气体的压缩因子,这个算法只列举了三种混合气体的例子,可继续增加
  • 粒子群算法改进——压缩因子

    千次阅读 2020-11-09 15:31:31
    在研究完粒子群算法中有关惯性权重的优化之后,我们把目光转向速度更新公式的后两项,根据之前所学可知:个体学习因子c1和社会学习因子c2决定了粒子本身经验信息和其他粒子的经验信息对粒子运行轨迹的影响,其反映了...

    前言概述

    粒子速度更新公式如下:

    vid = wvid-1 + c1r1(pbestid-xid)+ c2r2(gbestd-xid)

    在研究完粒子群算法中有关惯性权重的优化之后,我们把目光转向速度更新公式的后两项,根据之前所学可知:个体学习因子c1和社会学习因子c2决定了粒子本身经验信息和其他粒子的经验信息对粒子运行轨迹的影响,其反映了粒子群之间的信息交流。
    较大的c1值,会使粒子过多地在自身的局部范围内搜索
    较大的c2值,则又会促使粒子过早收敛到局部最优值。

    为了有效地控制粒子的飞行速度,使算法达到全局搜索与局部搜索两者间的有效平衡,Clerc构造了引入收缩因子的PSO模型,采用了压缩因子,这种调整方法通过合适选取参数,可确保PSO算法的收敛性,并可取消对速度的边界限制

    参数设置

    • 个体学习因子c1 = 2.05 (应用较多)
    • 社会学习因子 c2 = 2.05(应用较多)
    • 惯性权重 w = 0.9
    • C = c1+c2 = 4.1
    • 收缩因子在这里插入图片描述

    综上,速度更新公式如下
    **vid = ==K== ·(wvid-1 + c1r1(pbestid-xid)+ c2r2(gbestd-xid) )**

    示例代码

    求解函数y = 7cos(5x) + 4*sin(x)在[-5,5]内的最大值

    …… ……
    %% 初始化参数
    c1 = 2.05;  % 每个粒子的个体学习因子
    c2 = 2.05;  % 每个粒子的社会学习因子
    C = c1 + c2;
    w = 0.9;  % 惯性权重
    fai = 2/abs((2-C-sqrt(C^2-4*C)));% 收缩因子
    
    %% 初始化粒子的位置和速度
    …… ……
    
    %% 计算适应度
    …… ……
    
    %% 迭代K次来更新速度与位置
    fitnessbest = ones(K,1);  % 初始化每次迭代得到的最佳的适应度
    for d = 1:K  % 开始迭代,一共迭代K次
        for i = 1:n   % 依次更新第i个粒子的速度与位置
           v(i,:) = fai * (w*v(i,:) + c1*rand(1)*(pbest(i,:) - x(i,:)) + c2*rand(1)*(gbest - x(i,:)));  % 更新第i个粒子的速度
            x(i,:) = x(i,:) + v(i,:); % 更新第i个粒子的位置
            % 如果粒子的位置超出了定义域,就对其进行调整
         		…… ……
          …… ……
    end
    
    …… ……
    
    function y = Obj_fun1(x)
        y = 7*cos(5*x) + 4*sin(x);
    end
    

    运行结果如下

    在这里插入图片描述在这里插入图片描述对比未经优化的粒子群算法
    在这里插入图片描述在这里插入图片描述可知,求得的最优解更接近实际最优解

    参考文献:Eberhart R C . Comparing inertia weights and
    constriction factors in optimization[C]//
    Proceedings of the 2000 IEEE Congress on
    Evolutionary Computation, La Jolla, CA. IEEE, 2000.

    展开全文
  • 基于21中摩尔成分进行操作状态下天然气压缩因子的计算。
  • 输气管道压缩因子计算,石油大学,输气管道设计与管理
  • Matlab与LabVIEW混合编程的天然气压缩因子软件设计.pdf
  • 目的:采用压缩因子粒子群算法实现双指数拟合。 function x_opt = PSO_ExpFit2(t, Et) %{ 函数功能:压缩因子粒子群算法实现指数拟合:y = a1*exp(-x/b1) + a2*exp(-x/b2) + c 输入: t:自变量; Et:因变量; ...

    目的:采用压缩因子粒子群算法实现双指数拟合。

    function x_opt = PSO_ExpFit2(t, Et)
    %{
    函数功能:压缩因子粒子群算法实现指数拟合:y = a1*exp(-x/b1) + a2*exp(-x/b2) + c
    输入:
      t:自变量;
      Et:因变量;
    输出:
      x_opt:最优解;
    调用格式:
    clear;clc;
    t = 0.2*(1:3000)';
    data = 400*exp(-t/5) + 600*exp(-t/200) + 10*rand(size(t));
    tic;
    p = PSO_ExpFit2(t, data)
    toc;
    fit = p(1)*exp(-t/p(2)) + p(3)*exp(-t/p(4)) + p(5);
    plot(t, data, t, fit, 'LineWidth', 2)
    legend('model', 'PSO');
    %} 
    % 初始值
    D = 5;                      % 粒子维度,未知数个数;
    IterMax = 500;      % 最大迭代次数;
    Vmin = -1;               % 速度最小值;
    Vmax = 1;               % 速度最大值;
    c1 = 2.05;               % 学习因子1,认知部分;
    c2 = 2.05;               % 学习因子2,社会部分;
    N = 100;               % 种群个数;
    Xmin = [0, -2, 0, -2, 0];    % 下限
    max_sig = max(Et);
    Xmax = [1.5*max_sig, 4, 1.5*max_sig, 4, 0.5*max_sig];  % 上限;
    % 计算
    phy = c1 + c2;
    lamda = 2 / abs(2 - phy - sqrt(phy^2 - 4*phy));   % 压缩因子
    % 初始化种群个体(限定位置和速度) 
    x = rand(N, D).*repmat(Xmax - Xmin, N, 1) + repmat(Xmin, N, 1);
    v = rand(N, D)*(Vmax - Vmin) + Vmin;
    % 初始化个体最优位置和最优值 
    p = x;
    pbest = ones(1, N);
    for i = 1 : N
       pbest(i) = Obj_Fit(x(i, :), t ,Et);
    end
    % 初始化全局最优位置和最优值 
    g = ones(1, D);
    gbest = inf;
    for i = 1 : N
       if pbest(i) < gbest
          g = p(i, :);
          gbest = pbest(i);
       end
    end
    % 迭代直到满足精度或者迭代次数 
    fx = zeros(1, N);
    for i = 1 : IterMax
       for j = 1 : N
          % 更新个体最优位置和最优值 
          fx(j) = Obj_Fit(x(j, :), t ,Et);
          if fx(j) < pbest(j)
             p(j, :) = x(j, :);
             pbest(j) = fx(j);
          end
          % 更新全局最优位置和最优值 
          if pbest(j) < gbest
             g = p(j, :);
             gbest = pbest(j);
          end
          % 更新位置和速度值  
          v(j, :) = lamda*v(j, :) + c1*rand*(p(j, :) - x(j, :)) + c2*rand*(g - x(j, :));
          x(j, :) = x(j, :) + v(j, :);
          % 边界条件处理 
          for ii = 1 : D
             if (v(j, ii) > Vmax) || (v(j, ii) < Vmin)
                v(j, ii) = rand*(Vmax - Vmin) + Vmin;
             end
             if (x(j, ii) > Xmax(ii)) || (x(j, ii) < Xmin(ii))
                x(j, ii) = rand*(Xmax(ii) - Xmin(ii)) + Xmin(ii);
             end
          end
       end
    end
    x_opt = g;
    x_opt(2) = 10^(x_opt(2));
    x_opt(4) = 10^(x_opt(4));
    end
    
    % 目标函数
    function fitError = Obj_Fit(p0, t ,Et)
    A1 = p0(1);
    B1 = p0(2);
    A2 = p0(3);
    B2 = p0(4);
    A0 = p0(5);
    f = A1*exp(-t/10^B1) + A2*exp(-t/10^B2) + A0;  
    fitError = norm(Et - f);
    end

     

    展开全文
  • 天然气压缩因子计算

    热门讨论 2011-09-14 13:33:29
    基于AGA8 92DC的天然气压缩因子计算软件
  • 一个带压缩因子粒子群算法的matlab程序,用于求解多维无约束优化问题
  • 摘要:介绍了求解 BWRS状态方程中压缩因子的数值方法 ,在求解方面重点介绍了弦截法。对计算方法和求解方法进行对比 ,找到了求解 BWRS状态方程中压缩因子的简便快捷途径 ,给出了详细的求解过程;并结合计算实例对结果...
  • 天然气压缩因子
  • 主成分分析和因子分析的介绍、区别和联系

    万次阅读 多人点赞 2019-07-14 00:46:00
    本文介绍一下主成分分析、因子分析,进而介绍它们之间的区别和联系。 两个方法的推导我也还有一些没有完全理解,因此中间有些理解可能有误,请大家批评指正 主成分分析 主成分分析:将多个有一定相关性的指标进行...

    本文介绍一下主成分分析、因子分析,进而介绍它们之间的区别和联系。
    两个方法的推导我也还有一些没有完全理解,因此中间有些理解可能有误,请大家批评指正

    主成分分析

    主成分分析:将多个有一定相关性的指标进行线性组合,以最少的维度解释原数据中尽可能多的信息为目标进行降维,降维后的各变量间彼此线性无关。
    从其介绍中我们可以知道,主成分分析实际上是将相关性很强的变量通过一定的线性组合合并成一个综合变量,达到降维的效果,因此对于主成分分析只对那些变量间相关性强的变量才能有效降维。

    • 推导角度一:
      原来样本有p个变量x1···xp,现首先期望通过一个线性组合在这里插入图片描述使得相应的方差最大,接下来,类似的再通过一个线性组合找到第二个新变量在这里插入图片描述,使其在原样本信息中除被第一维度y1解释掉的部分,剩余部分信息的方差最大,依次类推下去,直到能够解释足够多的原样本空间信息(一般累计方差占80%以上),舍弃剩余的一些信息,以达到降维的效果。
      详细的推导见下面两张图片(划得有点乱了,实在抱歉啦)
      在这里插入图片描述
      图1

    在这里插入图片描述
    图2

    • 推导角度二:
      对于正交空间的样本点,寻找一个超平面,使得
      • 最近重构性:样本点到这个超平面的距离足够近;
      • 最大可分性:样本点在这个超平面上的投影能尽可能分开;
        周志华老师的书中表明这两种方法得到的是公式是等价的。
        在这里插入图片描述
        图3
        在这里插入图片描述
        图4
        我觉得周老师书中的最近重构性也可以理解为丢失掉的信息尽可能少,因为得到新的正交坐标系的基向量{w1,w2,w3````xd},然后丢弃掉一些坐标维度后,将源空间的样本点投影到新空间得到新空间样本点Zi,然后基于Zi重构xi,最后计算两者的距离,使其距离和最小,实际上就是通过距离来度量损失掉的信息,令丢失掉的信息最少,即距离最大,得到新的坐标维度。

    因子分析

    ————————————————————————————————————————
    首先,先介绍下因子分析的整个思路:
    因子分析实际上是提取出反应原变量的一些共性因子,同时希望尽可能少的损失样本信息。
    其推导的思路是这样的:1、建立一个模型:X=AF+e;2、希望求出A,于是对模型取方差,得在这里插入图片描述,若没有特殊因子项,然后使用样本的协方差阵去估计X的协方差,进而就可以求出A了;3、但是实际上我们希望得到少数几个能反应共性变异的因子,因此舍弃特征根小的那几个,即尽可能的少损失一些信息,并以此达到降维的效果,实际上损失的这一部分也就是特殊因子了。
    —————————————————————————————————————————
    接下来详细介绍这个过程:

    1. 因子分解

    因子分析概括来说就是要得到一些公共因子,利用这些因子来解释原来的样本空间数据,即将X进行分解,X=AF,即在这里插入图片描述
    .目标也是尽可能保留x中的信息。

    1. 因子旋转

    在完成分解,得到公共因子后,如果各个公共因子的实际意义不明显的话,一般还需要再进行坐标旋转,坐标旋转就是不改变因子所携带的各变量的信息量,但是改变各个因子所携带的信息量,使得各个因子的实际意义明显。
    比如最常见的正交旋转(Varimax),最大方差正交旋转法,使得因子载荷的方差尽可能大,即信息分布尽可能不均匀,也就是信息要集中分布于几个不同的因子上。

    1. 因子得分计算
      得到意义明确的公共因子之后,由于我们之前实际上是得到了因子载荷矩阵A,即得到了下述方程组在这里插入图片描述,但是我们会发现,方程数量p多余我们要求的值F的数量,这样一来我们需要用一些手段来得到近似F,比如回归法等等,来计算出因子得分,也就是各个样本各个因子的实际值。
      详细推导见下述,也画的有点多了,见谅
      因子分解的推导
      在这里插入图片描述
      图5
      在这里插入图片描述
      图6
      在这里插入图片描述
      图7
      因子旋转的具体过程
      在这里插入图片描述
      图8
      在这里插入图片描述
      图9
      因子得分的计算
      在这里插入图片描述
      图10
      ——————————————————————————————————————
      总结一下主成分分析的步骤和因子分析的步骤:
      主成分分析分析过程:

      1.将原始数据标准化,以消除变量间在数量级和量纲上的不同。
      2. 求标准化数据的相关矩阵。(标准化数据的协方差矩阵就是原数据的相关系数矩阵)
      3. 求相关矩阵的特征值和特征向量。
      4. 计算方差贡献率与累计方差贡献率:每个主成分的贡献率代表了原数据总信息量的百分比。
      5. 确定主成分:设C1,C2,C3…Cp为p个主成分,其中前m个主成分的的总信息量(方差贡献度)不低于80%时,课提取前m个主成分来反映原评价对象。
      6. 用原指标的线性组合来计算各主成分得分:以各主成分对原指标的相关系数为权,将各个主成分表示为原指标的线性组合,而主成分的经济意义则由权数较大的指标的综合意义来确定(一般主成分可以不做解释,这是因子分析擅长的)
      7. 综合得分:
      在这里插入图片描述
      8. 得分排序:利用总得分得到得分名次

    因子分分析分析过程:

    1. 将原始数据标准化,以消除变量间在数量级和量纲上的不同。
    2. 求标准化数据的相关矩阵。(标准化数据的协方差矩阵就是原数据的相关系数矩阵)
    3. 求相关矩阵的特征值和特征向量。
    4. 计算方差贡献率与累计方差贡献率。
    5. 确定因子:设F1,F2,F3…Fp为p个主成分,其中前m个因子的的总信息量(方差贡献度)不低于80%时,课提取前m个因子来反映原评价对象。
    6. 因子旋转:若所获得的m个因子无法确定或其实际意义不是很明显,这时需将因子进行旋转以获得较为明显的实际意义。
    7. 用原指标的线性组合来计算各因子得分:采用回归估计法,Bartlett估计法计算因子得分
    8. 综合得分:以各因子的方差贡献率为权,由各因子的线性组合得到综合评价指标函数
      在这里插入图片描述
    9. 得分排序:利用总得分得到得分名次
      ———————————————————————————————————————
      说一下我的理解:主成分分析和因子分析的区别和联系(可能理解有误,希望大家能指正)
      1、主成分的目的是将原来的几个相联系的变量通过线性组合压缩成少数的几个综合变量,加以考察,但是这几个综合变量往往不具有明显的现实意义,而因子分析则是期望将一些相互联系的变量压缩,也是通过线性组合得到一些公共因子,这些公共因子往往有明显的现实意义。一个在于将原来的变量综合(Y=AX),一个在于将变量分解,提取公共因素(X=AF)。
      2、主成分分析的目标是以尽可能少的综合变量去蕴含原来数据中尽可能多的信息,而因子分析的目标是尽可能包含原数据尽可能多的信息量,没有要求要以尽可能少的因子。也就是说主成分分析的关键在于“变异数”的问题,希望综合指标尽可能少,而因子分析在于共变异数的问题,也就是希望变量能够提取出原数据中共有信息(共有方差/变异)。
      对于其中说的共有方差/变异有些难理解,我是这样来理解的。
      在这里插入图片描述
      因子分解如上图所示,这个等式X=AF+e可以理解成一个回归式,因变量是X,自变量是F,e是误差(随机变量),所以实际上该因子分解时在做一个提取影响X的自变量(公共因子)的过程,也就是提取共有方差。
      3、当因子分析的特殊因子方差贡献率为0时,主成分分析和因子分析(使用主成分法求公共因子)是完全等价的。两者的系数矩阵之间只是相差一个系数,这个系数是各个特征根的二次平方根。当主成分分析取标准化系数时,两个系数矩阵就完全一致了。
      在这里插入图片描述
      ——————————————————————————————————
      验证性因子分析
      上面所说的因子分析也就是常用的探索性因子分析,而实际因子分析还有一种验证性因子分析,验证性因子分析正如其名字,是用来验证一种结构的,举一个例子:根据理论我们可以构建出一个问卷,问卷中有很多题目,其中1-3题共同反映因子A,4-6题共同反映因子B,然后便可以使用验证性因子分析(CFA),去验证这种关系是否存在,以及问卷设计的信效度是否很好。具体的介绍可以我的这篇博文

    ————————————————————————————————————
    参考文献:
    王斌会.《多元统计分析及R语言建模》
    周志华.《机器学习》
    张文彤,董伟.《SPSS统计分析高级教程》

    展开全文
  • 该算法将原始图像分成8×8的图像块,利用特定压缩因子对应的量化表生成图像块的DCT量化矩阵。对每个图像块,利用重要DCT量化系数和相应的水印标识生成待嵌入水印信息,并基于密钥随机嵌入在其它图像块对应的量化系数...
  • 时变压缩因子粒子群算法.pdf
  • 本文使用双Yukawa流体混合物模型来描述混合胶体溶液,使用一阶平均球近似理论和Monte Carlo模拟研究混合胶体溶液的温度、组成、粒径大小和作用能量参数对其渗透压缩因子的影响。计算结果表明,有一阶平均球近似计算...
  • 压缩因子综合信息粒子群算法.pdf
  • 混合胶体溶液的渗透压缩因子研究,靳琳,于养信,混合胶体溶液的热力学性质(如渗透压等)对于血液生理机理、制药以及光学晶体的制备等都有着重要的意义。本文使用双Yukawa流体混合
  • 带自适应压缩因子粒子群优化算法.pdf
  • Matlab编程天然气压缩因子计算模型.doc
  • 通过研究发现,AAC音频比例因子的值会随着压缩次数的增加而逐渐减小。基于此,提出了一种基于比例因子转移概率差值统计特性的AAC音频压缩历史检测算法。实验结果表明,该算法能对多次压缩的AAC音频进行准确分类,其...

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 33,050
精华内容 13,220
关键字:

压缩因子查询