%求解Bouc-Wen非線性滯回恢復模型的主函數
function?[tY]?=?bouc_wen（tspan）
%?tspan?若取二元數組，則該數組的元素決定解算的時間區間
%????????若是長度大于2的單調增或減數組，則該數組決定解算時間點
A=1;w=2*pi;phi=0;
alpha=1;beta=1.5;gamma=-0.5;n=1;
%利用4、5階龍格-庫塔法求解微分方程
%?Y?微分方程的解，此處為滯回恢復力
[tY]=ode45（@dzdttspan0[]alphabetagammanAwphi）;?
[x~]=jili（tAwphi）;????%位移
plot（xY）
end

%Bouc-Wen模型
%?Bouc-Wen模型由兩部分組成：
%?1）?非記憶部分——g（xdx/dt）
%?2）?記憶部分——z（x）
%?其中z（x）用以下所示的微分方程表達:
%

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1083??2015-04-06?21:58??bouc_wen.m

-----------?---------??----------?-----??----

?????????????????1083????????????????????1