主要內(nèi)容：matlab參數(shù)識(shí)別應(yīng)用，主要適用于微分方程、微分方程組參數(shù)識(shí)別、simulink模型參數(shù)識(shí)別，領(lǐng)域不限。1 使用matlab識(shí)別微分方程參數(shù)以及微分方程組(多個(gè)微分方程)參數(shù)

2 使用matlab調(diào)用simulink并識(shí)別simulink模型的參數(shù)(m函數(shù)與simulink交互)

內(nèi)容為本人在學(xué)習(xí)過程中總結(jié)的知識(shí)，拿出來與大家分享，希望大家多多討論。

下邊貼出一部分源碼，其它完整內(nèi)容在附件的rar壓縮包中，大家可以自己拿去聯(lián)系或者修改

%首先 在參數(shù)設(shè)置的data import /export下設(shè)置input和output，這樣直接可以用simulink 的in1和out1端口

%輸入和輸出數(shù)據(jù)，與workspace進(jìn)行數(shù)據(jù)交互。當(dāng)然也可以用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]';%將時(shí)間和數(shù)據(jù)同時(shí)寫入，要求如此。

options=[]

t=x'; %與simulink configuration paramter下的import 名稱一致。

u=t;

[t,x1,ydata1]=sim('pm',tspan,options,xdata);

% save('xdata.mat','xdata');

ydata=ydata1;%將第一次計(jì)算真實(shí)值賦給ydata；

%%

%注意點(diǎn)1，ydata1也就是實(shí)時(shí)仿真值需要更新；2，參數(shù)k要能進(jìn)行傳遞。

%調(diào)試心得：只要仿真時(shí)能夠?qū)崿F(xiàn)上述兩點(diǎn)，即可得到正確結(jié)果。

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, 下載次數(shù): 10)

2014-6-23 23:14 上傳

待識(shí)別模型,k1,k2參數(shù)

2014-6-23 23:25 上傳

點(diǎn)擊文件名下載附件

8.18 KB, 下載次數(shù): 807

matlab識(shí)別simulink參數(shù)

2014-6-23 23:25 上傳

點(diǎn)擊文件名下載附件

45.06 KB, 下載次數(shù): 803

微分方程擬合

2014-6-23 23:25 上傳

點(diǎn)擊文件名下載附件

29.53 KB, 下載次數(shù): 2071

微分方程組擬合

總結(jié)

以上是生活随笔為你收集整理的matlab参数识别,教程：matlab参数识别（包括微分方程和simulink）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。