代码:
c_0a=21;c_0b=3.5;k_0=46.9;c_1a=283;c_1b=2.95;k_1=5;x_0=14.3;
alpha_a=140;alpha_b=695;gamma=363;beta=363;A=301;n=2;enta=190;
alpha=alpha_a+alpha_bv;
c_1=c_1a+c_1bv;
c_0=c_0a+c_0b*v;
c_0=100;
output=0.0;
switch Flag
case 'getK' %鑾峰彇鍒氬害
case 'getF' %鑾峰彇鎶楀姏
K1=1/(c_0+c_1)*(alpha*zi+c_0*x_v+k_0*(x_d-yi));
L1=-gamma*abs(x_v-K1)*zi*(abs(zi))^(n-1)-beta*(x_v-K1)*(abs(zi))^n+A*(x_v-K1);
K2=1/(c_0+c_1)*(alpha*(zi+dt/2*L1)+c_0*x_v+k_0*(x_d-(yi+dt/2*K1)));
L2=-gamma*abs(x_v-K2)*(zi+dt/2*L1)*(abs(zi+dt/2*L1))^(n-1)-beta*(x_v-K2)*(abs(zi+dt/2*L1))^n+A*(x_v-K2);
结果:
错误使用 *
内部矩阵维度必须一致。
出错 Bouc_Wen_Func (line 37)
L2=-gammaabs(x_v-K2)(zi+dt/2L1)(abs(zi+dt/2L1))^(n-1)-beta(x_v-K2)(abs(zi+dt/2L1))^n+A*(x_v-K2);
出错 Solution_of_SDOF_under_Earthquake (line 79)
Fd_j=Bouc_Wen_Func('getF',0.0,dt,i*dt,disp_new,velo_new); %Bouc_Wen模型阻尼�?