matlab微分方程组_matlab求微分方程特解

2022-10-04 15:36:56 浏览数 (1)

大家好,又见面了,我是你们的朋友全栈君。

主要内容: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

0 人点赞