大家好,又见面了,我是你们的朋友全栈君。
主要内容:matlab参数识别应用,主要适用于微分方程、微分方程组参数识别、simulink模型参数识别,领域不限。1 使用matlab识别微分方程参数以及微分方程组(多个微分方程)参数
2 使用matlab调用simulink并识别simulink模型的参数(m函数与simulink交互)
内容为本人在学习过程中总结的知识,拿出来与大家分享,希望大家多多讨论。
下边贴出一部分源码,其它完整内容在附件的rar压缩包中,大家可以自己拿去联系或者修改
%首先 在参数设置的data import /export下设置input和output,这样直接可以用simulink 的in1和out1端口
%输入和输出数据,与workspace进行数据交互。当然也可以用xdata1=sim(‘pm’,tspan,options,ut);
clc;
clear;
global k1
global k2
global xdata;
k1=3;
k2=-4;
kc=[k1,k2]
tspan=[0,3];
load_system(‘pm’)
opt=simget(‘pm’);
x=min(tspan):opt.FixedStep:max(tspan);
xdata=[x;x]’;%将时间和数据同时写入,要求如此。
options=[]
t=x’; %与simulink configuration paramter下的import 名称一致。
u=t;
[t,x1,ydata1]=sim(‘pm’,tspan,options,xdata);
% save(‘xdata.mat’,’xdata’);
ydata=ydata1;%将第一次计算真实值赋给ydata;
%%
%注意点1,ydata1也就是实时仿真值需要更新;2,参数k要能进行传递。
%调试心得:只要仿真时能够实现上述两点,即可得到正确结果。
k0=[-1,-1];
F=@(k,tspan)mysf(k,tspan);
yd=F(k0,tspan)
plot(x’,yd,’or’,x’,ydata,’b*’)
lb=[-10,-10];
ub=[10,10];
options=optimset(‘TolFun’,1e-8,’TolX’,1e-8,’MaxFunEvals’,100, ‘Algorithm’,’trust-region-reflective’, ‘display’, ‘iter’);
[kp,resnorm]=lsqcurvefit(F,k0,tspan’,ydata,lb,ub,options)
yest=F(kp,tspan);
plot(x’,yest,’or’,x’,ydata,’b*’)
pm.jpg
(13.7 KB, 下载次数: 10)
2014-6-23 23:14 上传
待识别模型,k1,k2参数
2014-6-23 23:25 上传
点击文件名下载附件
8.18 KB, 下载次数: 807
matlab识别simulink参数
2014-6-23 23:25 上传
点击文件名下载附件
45.06 KB, 下载次数: 803
微分方程拟合
2014-6-23 23:25 上传
点击文件名下载附件
29.53 KB, 下载次数: 2071
微分方程组拟合
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。
发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/196481.html原文链接:https://javaforall.cn