i am trying to solve bouc-wen hysteresis model anybody help me to solve this problem

Question

SABARIANAND D V 2018-4-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/395144-i-am-trying-to-solve-bouc-wen-hysteresis-model-anybody-help-me-to-solve-this-problem

评论： Alireza Ansari 2021-11-8

采纳的回答： Hongtao Li

i am trying to solve the bouc-wen hysteresis model using ordinary differential equations help me solve this problem

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Hongtao Li 2018-5-24

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/395144-i-am-trying-to-solve-bouc-wen-hysteresis-model-anybody-help-me-to-solve-this-problem#answer_321854

编辑：Hongtao Li 2018-5-24

在 MATLAB Online 中打开

%% Reproduce Figs. 7-8 in SMYTH, Andrew; WU, Meiliang. Multi-rate Kalman filtering for the data fusion of %%displacement and acceleration response measurements in dynamic system monitoring. Mechanical %%Systems and Signal Processing, 2007, 21.2: 706-723.

%% Data source from %% http://strongmotioncenter.org/vdc/scripts/earthquakes.plx

C = textscan(fID, '%f','HeaderLines',51);
xgt = cell2mat(C);
Fsa = 200;
gt = ( 0:1:length(xgt)-1 ) / Fsa;
%%%Solve ordinary differential equation
%%%https://se.mathworks.com/help/matlab/ref/ode45.html
[t, y] = ode45(@(t,y)funBoucWen(t,y, gt, xgt), gt, [1 -2 0]);
%%%Build a new matlab script file and save as funBoucWen.m
function dydt = funBoucWen(t,y, gt, xgt)
% dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];
m = 1;
c = 0.3;
k = 9;
beta = 2;
gamma = 1;
n = 2;
xgt = interp1(gt, xgt, t);
dydt = [y(2); -c / m * y(2) - k/m * y(3) - xgt;...
    y(2) - beta * abs(y(2) ) * ( abs(y(3)) )^(n-1) * y(3) - gamma * y(2)...
    * ( abs(y(3)) )^n ];
end