求解方程111

clear
clc
close all

%%%********************************************************
U = 2;
H1 =0 ; % -0.511927637861460
H2 = -0.356812993614989;
H3 = -0.179577655772018;
H4 = 0.718894052877446;

A1 = 0.0501965897343734;
A2 = -0.0127991489739105; % 
A3 = 0.0316187993091640; %
A4 = 0.00273551358771455;

K1h = 3.14863993587337;
K1a = 7.41980044165051;

w1h = 7.67960914900707;
w1a = 18.0970741381280;

B = 0.82;             % 模型宽度
I = 0.666;            % 质量惯性矩
m = 14.92;            % 每延米质量
wh = 1.247*pi*2;      % 竖弯频率  其实是角速度
wa = 2.918*pi*2;      % 扭转频率  其实是角速度
ksh = 0.000;          % 竖弯阻尼比
ksa = 0.000; &nbs