帮我用matlab编程，平面框架结构，层高3.3m，共八层，每榀横向宽度为5m，共两榀周期为1s，增加隔震支座后框架结构周期为3s，求解在el-centro波作用下的前200阶振型，并计算各振型对应周期

好的，我可以为您提供一些编程思路和代码。

首先，我们需要定义框架结构的几何参数：

h = 3.3; % 层高
L = 5; % 横向宽度
N = 8; % 楼层数
a = 2*L; % 周期
T0 = 1; % 原周期
T = 3; % 隔震后周期

然后，我们需要定义材料参数和荷载：

E = 2.1e11; % 弹性模量
I = L^4/12; % 惯性矩
rho = 7850; % 密度
q = 1; % 单位长度荷载

接下来，我们可以使用MATLAB中的pdepe函数求解前200阶振型：

Nt = 200; % 求解的振型数
x = linspace(0, a, 101); % 离散化的x坐标
t = linspace(0, 10, 201); % 离散化的时间坐标

% 定义偏微分方程
function [c,f,s] = pdefun(x,t,u,DuDx)
    c = E*I/rho/h^2;
    f = DuDx;
    s = -q/h;
end

% 定义边界条件
function [pl,ql,pr,qr] = bcfun(xl,ul,xr,ur,t)
    pl = ul(1);
    ql = 0;
    pr = ur(1);
    qr = 0;
end

% 求解偏微分方程
sol = pdepe(0, @pdefun, @(x) [0;0], @bcfun, x, t);

% 取得前Nt个振型
modes = sol(:,1:Nt);

% 绘制前Nt个振型
plot(x, modes(end,:));

最后，我们可以计算各振型对应的周期：

periods = 2*pi./sqrt(T0./sol(:,1:Nt));