但是Pg的有关表达式是这样的，dPg/dt+Pg=-100*(1-f)，f是表达式1

非常抱歉，我在之前的回答中写错了表达式。在这里提供一个更新的YALMIP编程程序，用于求解最优的`t2`值：

% 定义变量
t = sdpvar(1); % 时间变量
w = sdpvar(1); % 角速度变量
Pg = sdpvar(1); % 电网有功功率变量

% 定义常量
A = 0.1; % 常数
W = 0.2; % 常数
M = 0.3; % 常数
Thet = 0.4; % 常数
C = 0.5; % 常数

% 定义表达式
f = M*exp(t*A)*cos(t*W+Thet)+C; % 表达式1
t1 = argmin(f, t); % 表达式2

dw_dt = (0.248*w-1.16*w-50*derivative(f, t)-20*(1-f))/(2.1*0.875); % 表达式3
Pwe = 97.81*w^3; % 表达式4
dPg_dt = derivative(Pg, t); % 表达式5
Pun = 30-Pg-Pwe; % 表达式6
fsecmax = 1-Pun*exp(A)-(f-Pun*M); % 表达式7

% 定义约束条件
constr = [t2>t1, w>1.47, w<2.52]; % 约束条件

% 定义目标函数
obj = fsecmax; % 目标函数

% 定义限制条件
cons = [dPg_dt + Pg == -100*(1-f)]; % 表达式5

% 求解最优解
optimize([constr, cons], obj);

% 输出结果
disp(value(t2));

在这个程序中，我们首先定义三个变量`t`、`w`和`Pg`，分别表示时间、角速度和电网有功功率。然后我们定义了一些常量，包括`A`、`W`、`M`、`Thet`和`C`。接着，我们定义了一些表达式，即题目中给出的公式，如表达式1、2、3、4、5、6和7。

接下来，我们定义了一个限制条件，即表达式5。我们可以通过将它放到一个限制条件集合中来实现，即`cons = [dPg_dt + Pg == -100*(1-f)];`。

最后，我们通过调用YALMIP的`optimize`函数来求解最优解。在这个程序中，我们使用了`[constr, cons]`来表示约束条件和限制条件的合集。最终，我们输出最优解的值，即`disp(value(t2));`。