一次二阶矩法

已于 2022-08-08 14:12:06 修改
阅读量2.8k 收藏 13
点赞数 4
分类专栏： 不确定性量化分析 文章标签： 算法
于 2022-08-07 20:25:08 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44606558/article/details/126165354
版权

一次二阶矩法

非线性功能函数g(x)展开成泰勒级数保留一次项，按照可靠指标形成求解方程，可以得到求解可靠度的一次二阶矩法。此方法简单、常用、易于掌握。
一次二阶矩法分中心点法和设计验算点法。中心点法不考虑随机变量的概率密度分布（不足）。

设计验算点法

原理

设独立正态分布变量组合在这里插入图片描述为结构极限状态方程及极限状态面在这里插入图片描述上的一点，在该点按泰勒级数展开保留一次项，得
在这里插入图片描述
同时可以得到在这里插入图片描述的均值以及标准差
在这里插入图片描述
在这里插入图片描述
结构可靠性指标为
在这里插入图片描述
令
在这里插入图片描述
处于极限状态面上及值为0，且把Xi用标准化变量在这里插入图片描述替换，并除以在这里插入图片描述得到

在这里插入图片描述
令
在这里插入图片描述

由上面推导，上式可整理成：
在这里插入图片描述
该方程为Y空间内的法线式超平面方程如图（二维情形），可靠度指标β为该空间内坐标原点到该极限状态面的最短距离，p*点称为设计验算点
在这里插入图片描述
设计验算点在Y空间中的坐标为
在这里插入图片描述
转化为原始空间的坐标为
在这里插入图片描述
可以利用如下迭代方法得到可靠度指标β
在这里插入图片描述

当量正态法(JC法)

当处理的X变量为非正态分布时，我们可以把它处理成相应的当量正态化变量X‘。JC法中进行该处理需满足在验算点处在这里插入图片描述和在这里插入图片描述的累积分布函数和概率密度函数分别对应相等，即
在这里插入图片描述
则当量正态化变量的均值和标准差为：
在这里插入图片描述
完成替换以后，我们就可以利用前面所述方法求得相应的可靠度指标。

映射变换法

该方法的思想是利用累积分布函数值相等的映射，将非正态分布随机变量变换为正态分布随机变量。
对于随机变量在这里插入图片描述有概率密度函数在这里插入图片描述和累计分布函数在这里插入图片描述，我们可以作一定变换（等概率变换）将它映射为标准正态变量Y：

在这里插入图片描述
由此可得
在这里插入图片描述

其中Y是标准正态向量，且可以推导出：
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
其中
在这里插入图片描述
同前面，我们可以用该方法，把其他分布的变量转化为标准正态变量。然后运用第一节方法得到结构可靠度指标。基本流程如下：
在这里插入图片描述

codePorter-faceCV
关注
  • 4
    点赞
  • 踩
  • 13
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
专栏目录
数据预处理：数据的正态化和标准变换
weixin_43899514的博客
03-25 2万+
一、正态数据 正态数据的分布是一个钟型分布，数据点的分布是中间高两边低的，例如一个学校里学生的身高。特别高和特别低的都少，高度居中的人较多。如下图（其中x轴为数据的取值，y轴为对应取值的频率） 二、数据需要正态化的原因 很多统计学的理论，如线性回归、假设检验和参数估计等，很多地方需要数据的分布为正态。而现实中的数据不一定服从正态或者是标准正态，因此需要对数据进行正态化，才可以使用已有的一些模型。有时候，模型的表现不好，可能跟数据的预处理有很大的关系。 三、两个指标判断是否是正态 数据一般有两个指标，样本均
可靠性方法：改进一次二阶矩法+Rackwitz-Fiessler方法的Matlab源代码
02-26
可靠性算法，改进一次二阶矩法的Matlab源代码，结合Rackwitz-Fiessler方法，能够考虑任意分布的随机变量，里面包含部分测试例子，可直接在Matlab软件中调用执行，文件中包含详细的注释。
可靠性分析算法：一次二阶矩法的Matlab源代码
02-26
可靠性算法，一次二阶矩法的Matlab源代码，里面包含部分测试例子，可直接在Matlab软件中调用执行，文件中包含详细的注释。
什么是一阶矩和二阶矩？
热门推荐
kyle1314608的博客
07-24 2万+
期望的公式扩展 一阶矩就是期望值，换句话说就是平均数(离散随机变量很好理解，连续的可以类比一下)。举例：xy坐标系中，x取大于零的整数，y1, y2, ...,yn 对应x=1, 2,..., n的值，现在我要对y求期望，就是所有y累加除以n，也就是y的均值。 此时y的均值我可以在坐标系中画一条线，我会发现所有的点都在这条线的两边。如果是中心矩我就会用每个值减去均值z=yn-y均作为...
yjej.rar_一次二阶_一次二阶 可靠_一次二阶矩_二阶矩_可靠指标法
07-14
编制改进一次二阶矩法计算可靠指标的程序 编制改进一次二阶矩法计算可靠指标的程序
用Matlab改进一次二阶矩法程序.doc
07-05
用Matlab改进一次二阶矩法程序
一阶二次矩.zip_一次二阶矩_一阶矩_可靠性_可靠性算法_算法 可靠性
07-14
一次二阶矩法，可靠性算法中的一种，适合初学者学习
11 二阶矩方法和Lovasz局部引理
hanhan的博客
06-12 933
文章目录11 二阶矩方法和Lovasz局部引理11.1 The Second-Moment Method——二阶矩方法11.1.1 二阶矩方法定理11.1.2 二阶矩方法的应用——随机图阈值11.2 Lovasz Local Lemma——Lovasz局部引理11.2.1 LLL定理11.2.2 LLL定理证明11.3 Asymmetric LLL 11 二阶矩方法和Lovasz局部引理 11.1 The Second-Moment Method——二阶矩方法 11.1.1 二阶矩方法定理 首先写出二阶矩方
Monte Carlo Method(蒙特·卡罗方法)
weixin_30951389的博客
04-17 3049
0-故事： 　　蒙特卡罗方法是计算模拟的基础，其名字来源于世界著名的赌城——摩纳哥的蒙特卡罗。 　　蒙特卡罗一词来源于意大利语，是为了纪念王子摩纳哥查理三世。蒙特卡罗（MonteCarlo）虽然是个赌城，但很小，估计跟北京的一条街差不多大。 　　其思想来源于著名的蒲丰投针问题(提出用投针实验的方法求圆周率π)。而后20世纪40年代美国在第二次世界大战中研制原子弹的“曼哈顿计划”计划的...
matlab一次二阶矩法正态分布,一次二阶矩法.ppt
weixin_39867142的博客
03-19 1511
设计验算点法求可靠指标评述 验算点法对极限状态方程中服从正态分布的随机变量计算结果尚可，而非正态分布误差较大。 JC法求可靠指标 问题的提出：验算点法对极限状态方程中服从正态分布的随机变量计算结果尚可，而非正态分布误差较大。 在工程结构可靠度分析中，永久荷载一般服从正态分布，风压、雪压、楼面活载服从其它类型(如极值Ⅰ型等)，截面抗力R服从对数正态分布。因此在极限状态方程中，常包含非正态变量分布的基...
一次二阶矩法程序设计
10-15
一次二阶矩法就是一种在随机变量的分布尚不清楚的情况下，采用只有均值和标准差的数学模型去求解结构可靠度的方法。由于该法将功能函数Z=g(x1, x2,……, xn)在某点用泰勒级数展开，使之线性化，然后求解结构的可靠度，因此称为一次二阶矩法。 C#是一种安全的、稳定的、简单的，由C和C++衍生出来的面向对象的编程语言。它在继承C和C++强大功能的同时去掉了一些它们的复杂特性。C#综合了VB简单的可视化操作和C++的高运行效率，以其强大的操作能力、优雅的语法风格、创新的语言特性和便捷的面向组件编程的支持成为.NET开发的首选语言。 传统的手工计算方法已不再适应今天高速发展的社会，必须采用一定的工具才会使得工程技术或科研人员的工作效率大大提高。利用程序语言来计算结构可靠度，对工程计算具有很重要的作用。本文中将以均值一次二阶矩法为理论基础，利用C#语言开发一个计算结构可靠指标的程序。
可靠性方法：改进一次二阶矩法的Matlab源代码
02-26
可靠性算法，改进一次二阶矩法的Matlab源代码，里面包含部分测试例子，可直接在Matlab软件中调用执行，文件中包含详细的注释。
高斯分布归一化、期望、二阶矩、方差推导证明
qq_37549266的博客
07-15 1万+
目录 写在前面的唠叨： 归一化推导证明： 期望（一阶矩）推导证明： 二阶矩推导证明： 方差推导证明： 写在前面的唠叨： 最近这段时间一直在研究深度学习之类的东西，虽然如今对几种常见的神经网络都有了很好的了解，用起来也比较顺手，但是越学也越觉得瓶颈越来越明显了，最大的问题觉得还是数学基础不行，学习那些常见的模型已经把线性代数的知识捡的差不多了，而到了想自己设计模型的时候，才忽然发...
jc法matlab程序,科学网-利用matlab编程实现jc法来求解可靠性-杜宝雷的博文
weixin_35733790的博客
03-17 182
这个应该算是课堂作业了，利用当量正态法来求解以下问题，求管道的可靠性指标，管材承载能力R服从对数正态分布，环向应力服从正态分布，地面载荷服从极值I型分布。ContentsJC METHODAUTHOR:[email protected]:2014-04-02functionexer()%EXER Summary of this function goes here% Detailed...
基于二阶矩阵的最优化问题（二）（附matlab代码）
qq_42035274的博客
05-31 1815
基于二阶矩阵的优化问题（二）非精确搜索确定步长沃尔夫（wolfe）条件Armijo条件曲率下降条件强沃尔夫（wolfe）条件Goldstein条件回溯法（backtracking） 非精确搜索确定步长 在之前的文章中，我们已经讨论了如何从XkX_kXk​更新到Xk+1X_{k+1}Xk+1​，这篇文章我们主要不精确搜索中确定步长的方法 沃尔夫（wolfe）条件 wolfe条件分为Armijo条件和曲率下降条件，只要函数f有下界，所求的α就满足wolfe和强wolfe条件。 Armijo条件 Armijo条件
概率-什么是一阶矩，二阶矩？
Allenlzcoder的博客
11-14 1万+
什么是一阶矩，二阶矩？       根据S.M.罗斯的概率论教程，一阶矩指E[X]，即数列X的均值称为一阶矩。       以此类推，E[Xn] ，n≥1，称为X的 n阶矩，也就是二阶矩、三阶矩...       参考：       1、图灵数学·统计学丛书08-概率论基础教程(第7版)-美S·M·罗斯-郑忠国(译)-人民邮电出版.
二阶矩过程
HeinSven的博客
01-03 1万+
随机变量序列： 依概率收敛于X：，记为（a.s.） 依概率1收敛于X：，记为（p） 依分布收敛于X：分布函数列弱收敛于，随机变量X的分布函数F(x)，记为   二阶矩空间： H空间上的范数： 随机变量序列的均方极限： 定义：设，若 称随机变量序列{}均方收敛于随机变量X，称X为的均方极限，记为 ps：当极限运算直接接触随机变量时，用符号l.i.m 柯西序列： 若...
蓝桥杯算法心得——矿石样本分析（贪心+两数之和）
最新发布
晴天学长
05-19 320
链表的时间复杂度要知道查找（o（1））（o（n））一般都用ArrayList，除非是修改。
一次二阶矩法matlab计算程序
07-22
### 回答1： 二阶矩法是一种常用的数值积分方法，用于计算定积分的近似值。在MATLAB中，可以通过编写程序来实现二阶矩法的计算。 以下是一次二阶矩法的MATLAB计算程序： ```matlab function result = second_order_quadrature(a, b, n, f) h = (b - a) / n; % 计算步长 x = a:h:b; % 构造等距节点 y = f(x); % 计算节点上的函数值 sum1 = 0; % 存储第一部分求和结果 sum2 = 0; % 存储第二部分求和结果 for i = 2:n sum1 = sum1 + y(i); % 第一部分求和 sum2 = sum2 + f((x(i - 1) + x(i)) / 2); % 第二部分求和 end result = h / 6 * (y(1) + 2 * sum1 + 4 * sum2 + y(end)); % 计算最终结果 end ``` 在上面的程序中，需要定义函数 `f(x)` 来表示被积函数，并传入积分区间的上下限 `a` 和 `b` 以及划分的子区间数 `n`。程序会自动计算出步长 `h`，并构造等距节点数组 `x`。然后，根据二阶矩法的公式，对每个子区间进行求和计算，最后得到近似的积分结果。 例如，如果要计算函数 `f(x) = x^2` 在区间 `[0, 2]` 上的积分近似值，并将区间划分成 10 个子区间，可以调用下面的代码： ```matlab result = second_order_quadrature(0, 2, 10, @(x) x^2); disp(result); ``` 程序会输出近似的积分结果。 ### 回答2： 二阶矩法是一种数值积分方法，用于近似计算函数的积分值。在Matlab中，可以通过编写计算程序来实现二阶矩法的计算。 具体而言，我们可按照以下步骤编写程序： 首先，要定义被积函数f(x)。可以在程序中使用匿名函数的形式，例如 f = @(x) x^2。 然后，我们需要定义积分区间[a, b]和积分剖分节点数n。 接下来，计算每个子区间的宽度h，即 h = (b - a) / n。 初始化积分值sum为0。 使用for循环遍历每个子区间： - 计算每个子区间的中点xm，即 xm = a + (i - 0.5) * h。 - 计算每个子区间的二阶矩值M，即 M = (f(xm + h/2) - f(xm - h/2)) / h。 - 将M乘以子区间宽度h，并累加到sum上，即 sum = sum + h * M。 循环结束后，输出sum作为最终的数值积分结果。 下面是一个示例的二阶矩法计算程序： ```matlab f = @(x) x^2; % 定义被积函数 a = 0; % 积分区间左端点 b = 1; % 积分区间右端点 n = 100; % 积分剖分节点数 h = (b - a) / n; % 子区间宽度 sum = 0; % 积分值初始化 for i = 1:n xm = a + (i - 0.5) * h; % 子区间中点 M = (f(xm + h/2) - f(xm - h/2)) / h; % 二阶矩值 sum = sum + h * M; % 累加到积分值 end disp(sum); % 输出积分结果 ``` 以上程序计算了函数f(x) = x^2在区间[0, 1]上的数值积分结果。您也可以根据实际需要修改函数和积分区间，以及设定不同的剖分节点数来计算不同函数在不同区间上的数值积分值。 ### 回答3： 二阶矩法是一种常用的数值积分方法，用于求解积分问题。它的基本思想是利用多个子区间上的二阶多项式来逼近积分函数，并通过求解线性方程组得到积分结果。 下面是一个用MATLAB编写的一次二阶矩法计算程序的示例： ``` function result = trapezoidal(f, a, b, n) % f：被积函数（在MATLAB中定义为函数句柄） % a：积分下限 % b：积分上限 % n：子区间数 h = (b - a) / n; % 子区间宽度 x = a:h:b; % 子区间节点 y = f(x); % 子区间上的函数值 % 计算权重 w = ones(size(y)); w(2:end-1) = 2; % 计算积分 result = (h/2) * sum(w .* y); end ``` 该程序定义了一个名为`trapezoidal`的函数，它接受被积函数`f`、积分下限`a`、积分上限`b`和子区间数`n`作为输入参数。函数首先计算出子区间宽度`h`，然后在积分上限和下限之间生成等间距的子区间节点。然后，计算每个子区间上的函数值，并根据二阶矩法的权重规则计算出权重。最后，利用求和函数`sum()`和点乘计算得到积分结果。 这个程序可以用于求解一维函数的定积分问题，并且通过调整子区间数`n`的大小可以改变数值积分方法的精度。在使用时，需要根据具体的被积函数和积分区间进行相应的修改。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值