[发明专利]一种仿真数学模型参数对的估量方法在审

专利信息
申请号: 201810212438.4 申请日: 2018-03-15
公开(公告)号: CN108491353A 公开(公告)日: 2018-09-04
发明(设计)人: 张留伟 申请(专利权)人: 信阳师范学院
主分类号: G06F17/11 分类号: G06F17/11;G06F17/18
代理公司: 西安铭泽知识产权代理事务所(普通合伙) 61223 代理人: 韩晓娟
地址: 464000 河*** 国省代码: 河南;41
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 仿真数学模型 平方和 仿真建模技术 最小二乘原理 解微分方程 参数估计 参数结果 残余误差 初始时刻 仿真系统 精度参数 模型参数 目标函数 实验数据 数学模型 误差理论 优化仿真 优化目标 状态确定 数值解 最优化 最优解
【权利要求书】:

1.一种仿真数学模型参数对的估量方法,其特征在于:该估量方法的具体步骤如下:

S1:假设一个仿真系统中有A、B两个参数,在t时刻,参数A的数目为x(t),参数B的数目为y(t),其中x(t)、y(t)均大于0,x(t)、y(t)之间满足以下变化规律

初始条件为:x(t0)=α5,y(t0)=α6,其中αk(1≤k≤6)为模型的待定参数,则方程的通解为

α1lny(t)-α2y(t)+α3lnx(t)+α4x(t)+C

其中,C为积分常数,由初始条件决定;

S2:在确定模型参数之前,先确定参数对最优解的优化目标,根据误差理论中的最小二乘原理,将残余误差平方和ε作为目标函数,建立最优化模型,即

其中,Xtj为微分方程的n值解在t=t0,t1,L,tm,n刻的值。

S3:根据初始时刻t0的状态确定参数α5=x0,α6=y0,之后将三组实验数据代入到公式α1lny(t)-α2y(t)+α3lnx(t)+α4x(t)+C中,可初步设定α2α3的值,求出参数αk(k=1,…,4)和积分常数C的值;

S4:根据参数值用数值的方法解微分方程,得到t0,t1,…,tm时刻对应的数值解Xtj,Ytj(j=0,1,…,m),将数值解Xtj,Ytj(j=0,1,…,m)与实验数据xj,yj(j=0,1,…,m)的参与误差平方和ε作为目标函数,不断循环调整αk(k=1,…,4),最终得到使ε最小的αk(k=1,…,4)值,并建立最优仿真数学模型。

2.根据权利要求1所述的一种仿真数学模型参数对的估量方法,其特征在于:所述步骤S4中,在进行具体的计算时:

S41:由初始条件直接确定αk(k=5,6)的值;

S42:在估计值附近对αk(k=2,3)循环搜索,每次搜索时,针对特定的αk(k=2,3)值,由三组实验数据初步确定参数αk(k=1,4)和积分常数C;

S43:用四阶龙格-库塔公式提供开始值,用Adams四阶预测-校正系统求解微分方程组的数值解,计算数值解Xtj,Ytj(j=0,1,…,m)与实验数据xj,yj(j=0,1,…,m)的均方差ε;

S44:重复步骤S41~S43,直到找到本次搜索的最小误差平方和ε及其对应的αk(k=1,…,4)值;

S45:缩小搜索范围,减小搜索步长,重复步骤S41~S44,直至满足ε的精度要求。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于信阳师范学院,未经信阳师范学院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201810212438.4/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top